annotate maaslin-4450aa4ecc84/maaslin-4450aa4ecc84/src/test-Maaslin/test-Maaslin.R @ 6:ca61989bc3b4

Uploaded
author george-weingart
date Sun, 08 Feb 2015 23:39:43 -0500
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
6
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
1 c_strDir <- file.path(getwd( ),"..")
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
2
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
3 source(file.path(c_strDir,"lib","Constants.R"))
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
4 strTestingDirectory = file.path(c_strDir,c_strTestingDirectory)
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
5 sScriptMaaslin = file.path( c_strDir, "Maaslin.R" )
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
6
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
7 context("Test Run From Commandline")
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
8
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
9 #Input Files
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
10 sTestReadConfig = file.path(strTestingDirectory, c_strTestingInput, "TestMaaslin.read.config")
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
11 sTestCustomR = file.path(strTestingDirectory, c_strTestingInput, "TestMaaslin.R")
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
12 sTestMaaslinDirectory = file.path(strTestingDirectory, c_strTemporaryFiles, "testMaaslin")
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
13 sTestOutput = file.path(sTestMaaslinDirectory,"TestMaaslin_Summary.txt")
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
14 sTestTSV = file.path(strTestingDirectory, c_strTestingInput, "TestMaaslin.tsv")
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
15 #Test file answers
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
16 sTestOutputAnswer = file.path(strTestingDirectory, c_strCorrectAnswers, "TestMaaslin.tsv")
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
17
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
18 #Delete Test MaAsLin output
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
19 unlink(sTestMaaslinDirectory, recursive=TRUE)
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
20 #Make neccessary directories
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
21 dir.create(sTestMaaslinDirectory)
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
22 dir.create(file.path(sTestMaaslinDirectory,"QC"))
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
23
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
24 sCommand = paste(sScriptMaaslin, "-v", "ERROR", "-d", "0.25", "-r", "0.0001", "-p", "0.1", sTestOutput, sTestTSV, sTestReadConfig, sTestCustomR, sep=" ")
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
25 print(sCommand)
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
26 system(sCommand)
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
27
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
28 sExpectedTitle = "\tVariable\tFeature\tValue\tCoefficient\tN\tN.not.0\tP.value\tQ.value"
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
29 iExpectedNumberOfLines = 3
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
30 lsOutputSummaryFile = readLines(sTestOutput)
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
31
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
32 test_that("Make sure that the summary output file is what is expected (generally).",{
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
33 expect_equal(lsOutputSummaryFile[1], sExpectedTitle)
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
34 expect_equal(length(lsOutputSummaryFile),iExpectedNumberOfLines)
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
35 })
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
36
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
37 lsDirectoryStructure = list.files(sTestMaaslinDirectory)
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
38 lsDirectoryStructureAnswer = c(basename(sTestOutput),"QC","TestMaaslin-age.pdf","TestMaaslin-age.txt","TestMaaslin-dx.txt","TestMaaslin.pdf","TestMaaslin.txt")
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
39 test_that("Make sure the expected directory structure is created.",{
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
40 expect_equal(sort(lsDirectoryStructure), sort(lsDirectoryStructureAnswer))
ca61989bc3b4 Uploaded
george-weingart
parents:
diff changeset
41 })