1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
|
test_that("anyNA() method for SVT_SparseArray objects", {
## input of type() "integer"
m1 <- matrix(c(0L, 0L, 155L,
0L, 8L, -1L), nrow=2, byrow=TRUE)
svt1 <- as(m1, "SVT_SparseArray")
test_summarize_op1(m1, svt1, "anyNA")
m1[1, 2] <- NA
svt1 <- as(m1, "SVT_SparseArray")
test_summarize_op1(m1, svt1, "anyNA")
## input of type() "logical"
m2 <- matrix(c(FALSE, FALSE, TRUE,
FALSE, TRUE, TRUE), nrow=2, byrow=TRUE)
svt2 <- as(m2, "SVT_SparseArray")
test_summarize_op1(m2, svt2, "anyNA")
m2[1, 2] <- NA
svt2 <- as(m2, "SVT_SparseArray")
test_summarize_op1(m2, svt2, "anyNA")
## input of type() "double"
m3 <- matrix(c(0, 0, pi,
0, 0.25, 1e3), nrow=2, byrow=TRUE)
svt3 <- as(m3, "SVT_SparseArray")
test_summarize_op1(m3, svt3, "anyNA")
m3[1, 2] <- svt3[1, 2] <- NaN
test_summarize_op1(m3, svt3, "anyNA")
m3[1, 2] <- svt3[1, 2] <- NA
test_summarize_op1(m3, svt3, "anyNA")
## input of type() "complex"
m4 <- matrix(c(0, 0, pi,
0, 2-5i, 1e3), nrow=2, byrow=TRUE)
svt4 <- as(m4, "SVT_SparseArray")
test_summarize_op1(m4, svt4, "anyNA")
m4[1, 2] <- NaN # 1st type of "complex" NaN
svt4 <- as(m4, "SVT_SparseArray")
test_summarize_op1(m4, svt4, "anyNA")
m4[1, 2] <- NaN * 1i # 2nd type of "complex" NaN
svt4 <- as(m4, "SVT_SparseArray")
test_summarize_op1(m4, svt4, "anyNA")
m4[1, 2] <- NA
svt4 <- as(m4, "SVT_SparseArray")
test_summarize_op1(m4, svt4, "anyNA")
## input of type() "character"
m5 <- matrix(c("", "", "Hello",
"", "dear", "world"), nrow=2, byrow=TRUE)
svt5 <- as(m5, "SVT_SparseArray")
test_summarize_op1(m5, svt5, "anyNA")
m5[1, 2] <- NA
svt5 <- as(m5, "SVT_SparseArray")
test_summarize_op1(m5, svt5, "anyNA")
})
test_that("other summarization methods for SVT_SparseArray objects", {
## input of type() "integer"
m1 <- matrix(c( 0L, 0L, NA, 0L, NA,
NA, 0L, -3L, 1L, NA,
0L, 0L, 0L, 0L, 0L,
15L, 0L, 0L, 0L, NA), nrow=4, byrow=TRUE)
svt1 <- as(m1, "SVT_SparseArray")
test_summarize_op2(m1, svt1, "any")
test_summarize_op2(m1, svt1, "all")
test_summarize_op2(m1, svt1, "min")
test_summarize_op2(m1, svt1, "max")
test_summarize_op2(m1, svt1, "range")
test_summarize_op2(m1, svt1, "sum")
test_summarize_op2(m1, svt1, "prod")
test_summarize_op2(m1, svt1, "mean")
test_summarize_op2(m1, svt1, "var")
test_summarize_op2(m1, svt1, "sd")
m0 <- m1[0, ]
svt0 <- svt1[0, ]
expect_warning(min(svt0), "NAs introduced")
expect_warning(max(svt0), "NAs introduced")
expect_warning(range(svt0), "NAs introduced")
expect_identical(suppressWarnings(min(svt0)), NA_integer_)
expect_identical(suppressWarnings(max(svt0)), NA_integer_)
expect_identical(suppressWarnings(range(svt0)), rep(NA_integer_,2))
## input of type() "logical"
m2 <- is.na(m1)
svt2 <- as(m2, "SVT_SparseArray")
test_summarize_op2(m2, svt2, "any")
test_summarize_op2(m2, svt2, "all")
test_summarize_op2(m2, svt2, "min")
test_summarize_op2(m2, svt2, "max")
test_summarize_op2(m2, svt2, "range")
test_summarize_op2(m2, svt2, "sum")
test_summarize_op2(m2, svt2, "prod")
test_summarize_op2(m2, svt2, "mean")
test_summarize_op2(m2, svt2, "var")
test_summarize_op2(m2, svt2, "sd")
m0 <- m2[0, ]
svt0 <- svt2[0, ]
expect_warning(min(svt0), "NAs introduced")
expect_warning(max(svt0), "NAs introduced")
expect_warning(range(svt0), "NAs introduced")
expect_identical(suppressWarnings(min(svt0)), NA_integer_)
expect_identical(suppressWarnings(max(svt0)), NA_integer_)
expect_identical(suppressWarnings(range(svt0)), rep(NA_integer_,2))
})
test_that("summarization methods for 3D SVT_SparseArray objects", {
## input of type() "double"
a <- array(0, 6:4)
a[1, , 2] <- c(1e12, -1234.55, -2.1, -1, -0.55)
a[3, , 2] <- c(-0.55, 0, 1e-10, 0.88, 1)
a[5, , 2] <- c(pi, 10.33, 3.4567895e8, 300, 2009.01)
svt3 <- as(a, "SVT_SparseArray")
test_summarize_op1(a, svt3, "anyNA")
a[6, 3, 2] <- svt3[6, 3, 2] <- NA
a[6, 4, 2] <- svt3[6, 4, 2] <- NaN
test_summarize_op1(a, svt3, "anyNA")
expect_error(any(svt3), "does not support")
expect_error(all(svt3), "does not support")
test_summarize_op2(a, svt3, "min")
test_summarize_op2(a, svt3, "max")
test_summarize_op2(a, svt3, "range")
test_summarize_op2(a, svt3, "sum")
test_summarize_op2(a, svt3, "prod")
test_summarize_op2(a, svt3, "mean")
test_summarize_op2(a, svt3, "var")
test_summarize_op2(a, svt3, "sd")
})
|
