|
a |
|
b/tests/limma/test_dups.py |
|
|
1 |
import unittest |
|
|
2 |
|
|
|
3 |
import numpy as np |
|
|
4 |
|
|
|
5 |
from inmoose.limma import uniquegenelist, unwrapdups |
|
|
6 |
|
|
|
7 |
|
|
|
8 |
class Test(unittest.TestCase): |
|
|
9 |
def test_unwrapdups(self): |
|
|
10 |
M = np.arange(12).reshape((4, 3)) |
|
|
11 |
self.assertTrue(np.array_equal(unwrapdups(M, ndups=1), M)) |
|
|
12 |
self.assertTrue(np.array_equal(unwrapdups(M, ndups=1, spacing=2), M)) |
|
|
13 |
self.assertTrue( |
|
|
14 |
np.array_equal(unwrapdups(M), [[0, 3, 1, 4, 2, 5], [6, 9, 7, 10, 8, 11]]) |
|
|
15 |
) |
|
|
16 |
self.assertTrue( |
|
|
17 |
np.array_equal( |
|
|
18 |
unwrapdups(M, ndups=4), [[0, 3, 6, 9, 1, 4, 7, 10, 2, 5, 8, 11]] |
|
|
19 |
) |
|
|
20 |
) |
|
|
21 |
self.assertTrue( |
|
|
22 |
np.array_equal( |
|
|
23 |
unwrapdups(M, ndups=2, spacing=2), |
|
|
24 |
[[0, 6, 1, 7, 2, 8], [3, 9, 4, 10, 5, 11]], |
|
|
25 |
) |
|
|
26 |
) |
|
|
27 |
|
|
|
28 |
def test_uniquegenelist(self): |
|
|
29 |
M = np.arange(12).reshape((4, 3)) |
|
|
30 |
self.assertTrue(np.array_equal(uniquegenelist(M, ndups=1), M)) |
|
|
31 |
self.assertTrue(np.array_equal(uniquegenelist(M, ndups=1, spacing=2), M)) |
|
|
32 |
self.assertTrue(np.array_equal(uniquegenelist(M), [[0, 1, 2], [6, 7, 8]])) |
|
|
33 |
self.assertTrue(np.array_equal(uniquegenelist(M, ndups=4), [[0, 1, 2]])) |
|
|
34 |
self.assertTrue( |
|
|
35 |
np.array_equal( |
|
|
36 |
uniquegenelist(M, ndups=2, spacing=2), [[0, 1, 2], [3, 4, 5]] |
|
|
37 |
) |
|
|
38 |
) |
|
|
39 |
|
|
|
40 |
M = np.arange(6) |
|
|
41 |
self.assertTrue(np.array_equal(uniquegenelist(M, ndups=1), M)) |
|
|
42 |
self.assertTrue(np.array_equal(uniquegenelist(M, ndups=1, spacing=2), M)) |
|
|
43 |
self.assertTrue(np.array_equal(uniquegenelist(M), [0, 2, 4])) |
|
|
44 |
self.assertTrue(np.array_equal(uniquegenelist(M, ndups=3), [0, 3])) |
|
|
45 |
self.assertTrue(np.array_equal(uniquegenelist(M, ndups=3, spacing=2), [0, 1])) |
|
|
46 |
self.assertTrue( |
|
|
47 |
np.array_equal(uniquegenelist(M, ndups=2, spacing=3), [0, 1, 2]) |
|
|
48 |
) |