annotate batchcorrection-57edfd3943ab/runit/batchcorrection_tests.R @ 3:73892ef177e3 draft

Uploaded
author melpetera
date Tue, 02 May 2017 09:47:22 -0400
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
3
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
1 test_input_allLoessPool <- function() {
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
2
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
3 testDirC <- "input"
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
4 determineL <- FALSE
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
5 allLoessL <- TRUE
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
6 argLs <- list(method = "all_loess_pool",
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
7 span = "1")
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
8
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
9 if(!allLoessL)
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
10 argLs[["analyse"]] <- ifelse(determineL, "determine_bc", "batch_correction")
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
11
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
12 argLs <- c(defaultArgF(testDirC, determineL = determineL), argLs)
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
13 outLs <- wrapperCallF(argLs, allLoessL = allLoessL)
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
14
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
15 checkEqualsNumeric(outLs[['datMN']][1, 1], 25803076, tolerance = 1e-3)
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
16
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
17 }
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
18
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
19 test_input_allLoessSample <- function() {
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
20
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
21 testDirC <- "input"
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
22 determineL <- FALSE
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
23 allLoessL <- TRUE
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
24 argLs <- list(method = "all_loess_sample",
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
25 span = "1")
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
26
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
27 if(!allLoessL)
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
28 argLs[["analyse"]] <- ifelse(determineL, "determine_bc", "batch_correction")
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
29
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
30 argLs <- c(defaultArgF(testDirC, determineL = determineL), argLs)
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
31 outLs <- wrapperCallF(argLs, allLoessL = allLoessL)
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
32
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
33 checkEqualsNumeric(outLs[['datMN']][1, 1], 23402048, tolerance = 1e-3)
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
34
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
35 }
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
36
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
37 test_example1_allLoessSample <- function() {
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
38
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
39 testDirC <- "example1"
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
40 determineL <- FALSE
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
41 allLoessL <- TRUE
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
42 argLs <- list(method = "all_loess_sample",
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
43 span = "1")
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
44
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
45 if(!allLoessL)
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
46 argLs[["analyse"]] <- ifelse(determineL, "determine_bc", "batch_correction")
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
47
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
48 argLs <- c(defaultArgF(testDirC, determineL = determineL), argLs)
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
49 outLs <- wrapperCallF(argLs, allLoessL = allLoessL)
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
50
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
51 checkEqualsNumeric(outLs[['datMN']][1, 1], 21732604, tolerance = 1e-3)
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
52
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
53 }
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
54
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
55 test_example1_allLoessSampleSpan06 <- function() {
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
56
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
57 testDirC <- "example1"
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
58 determineL <- FALSE
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
59 allLoessL <- TRUE
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
60 argLs <- list(method = "all_loess_sample",
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
61 span = "0.6")
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
62
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
63 if(!allLoessL)
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
64 argLs[["analyse"]] <- ifelse(determineL, "determine_bc", "batch_correction")
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
65
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
66 argLs <- c(defaultArgF(testDirC, determineL = determineL), argLs)
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
67 outLs <- wrapperCallF(argLs, allLoessL = allLoessL)
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
68
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
69 checkEqualsNumeric(outLs[['datMN']][1, 1], 134619170, tolerance = 1e-3)
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
70
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
71 }
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
72
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
73 test_sacurine_allLoessPool <- function() {
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
74
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
75 testDirC <- "sacurine"
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
76 determineL <- FALSE
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
77 allLoessL <- TRUE
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
78 argLs <- list(method = "all_loess_pool",
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
79 span = "1")
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
80
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
81 if(!allLoessL)
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
82 argLs[["analyse"]] <- ifelse(determineL, "determine_bc", "batch_correction")
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
83
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
84 argLs <- c(defaultArgF(testDirC, determineL = determineL), argLs)
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
85 outLs <- wrapperCallF(argLs, allLoessL = allLoessL)
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
86
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
87 checkEqualsNumeric(outLs[['datMN']]["HU_neg_017", "M53T345"], 7902.366, tolerance = 1e-3)
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
88
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
89
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
90 }
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
91
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
92 test_sacurine_determinebc <- function() {
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
93
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
94 testDirC <- "sacurine"
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
95 determineL <- TRUE
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
96 allLoessL <- FALSE
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
97 argLs <- list(ref_factor = "batch",
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
98 span = "none")
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
99
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
100 if(!allLoessL)
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
101 argLs[["analyse"]] <- ifelse(determineL, "determine_bc", "batch_correction")
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
102
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
103 argLs <- c(defaultArgF(testDirC, determineL = determineL), argLs)
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
104 outLs <- wrapperCallF(argLs, allLoessL = allLoessL)
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
105
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
106 checkEqualsNumeric(outLs[['sumDF']]["M59T62", "batch.2.linear"], 3, tolerance = 1e-3)
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
107
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
108
73892ef177e3 Uploaded
melpetera
parents:
diff changeset
109 }