3aeae6f8fe0596f90f9cceaec800233e7b0df937
[synclust.git] / R / tests / t.utils.R
1 #test matrix [de]standardization
2 test.de_standardize = function()
3 {
4 n = 100
5 m = 10
6 M = matrix(rnorm(n*m,mean=2.0,sd=2.0),nrow=n,ncol=m)
7
8 std = synclust:::standardize(M)
9 #result is centered:
10 checkEquals(rep(0.0, m), colMeans(std$M))
11 #result is standardized:
12 checkEquals(rep(1.0, m), sqrt( unlist( apply(std$M, 2, var) ) ))
13
14 #rebuilt M == M:
15 M_rec = synclust:::destandardize(std)
16 checkEquals(M, M_rec)
17 }
18
19 #test neighborhoods remapping into one smaller component
20 test.remap_neighbors = function()
21 {
22 #connex comp : 1-2-5-8-10
23 #to remap into 1-2-3-4-5
24 NI = list(
25 "ix" = list(
26 c(2,8,10), #V(1)
27 c(1,5,8,10), #V(2)
28 c(4,7,11), #V(3)
29 c(3,6,9,12), #V(4)
30 c(2,10), #V(5)
31 c(4,7), #V(6)
32 c(3,6,9,12), #V(7)
33 c(1,2,10), #V(8)
34 c(4,7,11), #V(9)
35 c(1,2,5,8), #V(10)
36 c(3,9), #V(11)
37 c(4,7)), #V(12)
38 "ds" = list(
39 c(1.0,2.0,3.0), #1
40 c(1.0,2.0,3.0,4.0), #2
41 c(1.0,1.0,1.0), #3
42 c(1.0,1.0,1.0,1.0), #4
43 c(2.0,2.0), #5
44 c(1.0,1.0), #6
45 c(1.0,1.0,1.0,1.0), #7
46 c(2.0,3.0,1.0), #8
47 c(1.0,1.0,1.0), #9
48 c(3.0,4.0,2.0,1.0), #10
49 c(1.0,1.0), #11
50 c(1.0,1.0))) #12
51
52 indices = c(1,2,5,8,10)
53 locNI = synclust:::remapNeighbors(NI, indices)
54 checkEquals(2, length(locNI))
55 checkEquals(length(indices), length(locNI$ix))
56 checkEquals(length(indices), length(locNI$ds))
57
58 #formerly index 1 (now 1)
59 checkEquals(c(2,4,5), locNI$ix[[1]])
60 checkEquals(NI$ds[[1]], locNI$ds[[1]],)
61 #formerly index 2 (now 2)
62 checkEquals(c(1,3,4,5), locNI$ix[[2]])
63 checkEquals(NI$ds[[2]], locNI$ds[[2]])
64 #formerly index 5 (now 3)
65 checkEquals(c(2,5), locNI$ix[[3]])
66 checkEquals(NI$ds[[5]], locNI$ds[[3]])
67 #formerly index 8 (now 4)
68 checkEquals(c(1,2,5), locNI$ix[[4]])
69 checkEquals(NI$ds[[8]], locNI$ds[[4]])
70 #formerly index 10 (now 5)
71 checkEquals(c(1,2,3,4), locNI$ix[[5]])
72 checkEquals(NI$ds[[10]], locNI$ds[[5]])
73 }