Mercurial > repos > ethevenot > univariate
diff runit/univariate_tests.R @ 3:140290de7986 draft
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) }