Mercurial > repos > ethevenot > univariate

diff runit/univariate_tests.R @ 3:140290de7986 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/workflow4metabolomics/univariate.git commit 27bc6157f43574f038b3fb1be1f46ce4786e24b1
author ethevenot
date Sun, 30 Oct 2016 14:17:09 -0400
parents 09799fc16bc6
children
line wrap: on
line diff
--- a/runit/univariate_tests.R	Sat Aug 06 12:42:42 2016 -0400
+++ b/runit/univariate_tests.R	Sun Oct 30 14:17:09 2016 -0400
@@ -1,7 +1,24 @@
+test_input_anova <- function() {
+
+    testDirC <- "input"
+    argLs <- list(facC = "qual",
+                  tesC = "anova",
+                  adjC = "fdr",
+                  thrN = "0.05")
+
+    argLs <- c(defaultArgF(testDirC), argLs)
+    outLs <- wrapperCallF(argLs)
+
+    checkEqualsNumeric(outLs[["varDF"]]["v6", "qual_anova_fdr"], 1.924156e-03, tolerance = 1e-6)
+
+    checkEqualsNumeric(outLs[["varDF"]]["v4", "qual_anova_D.C_fdr"], 0.01102016, tolerance = 1e-6)
+
+}
+
 test_input_kruskal <- function() {
 
     testDirC <- "input"
-    argLs <- list(facC = "ageGroup",
+    argLs <- list(facC = "qual",
                   tesC = "kruskal",
                   adjC = "fdr",
                   thrN = "0.05")
@@ -9,7 +26,9 @@
     argLs <- c(defaultArgF(testDirC), argLs)
     outLs <- wrapperCallF(argLs)
 
-    checkEqualsNumeric(outLs[["varDF"]]["HMDB01032", "ageGroup_kruskal_senior.experienced_pva"], 0.1231246, tolerance = 1e-6)
+    checkEqualsNumeric(outLs[["varDF"]]["v4", "qual_kruskal_fdr"], 0.0008194662, tolerance = 1e-7)
+
+    checkEqualsNumeric(outLs[["varDF"]]["v6", "qual_kruskal_D.A_fdr"], 0.002945952, tolerance = 1e-7)
 
 }