annotate BC/test/batchcorrection_tests.R @ 3:2e3a23dd6c24 draft default tip

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