Browse code

debug for uniTest

git-svn-id: file:///home/git/hedgehog.fhcrc.org/bioconductor/trunk/madman/Rpacks/SICtools@115482 bc3139a8-67e5-0310-9ffc-ced21a209358

Xiaobin Xing authored on 31/03/2016 07:31:26
Showing 2 changed files

... ...
@@ -7,19 +7,19 @@ test_indelDiff <- function(){
7 7
 	refFsa <- system.file(package='SICtools','extdata','example.ref.fasta')
8 8
 	
9 9
 	## test deletion
10
-	deletionTest <- indelDiff(bam1,bam2,refFsa,'chr04',962501,1026983,pValueCutOff=1,gtDistCutOff=0,verbose=FALSE)
10
+	deletionTest <- indelDiff(bam1,bam2,refFsa,'chr04',962501,1026983,pValueCutOff=1,gtDistCutOff=0,verbose=FALSE,nCores=2)
11 11
 	checkEquals(deletionTest$pos, 1026860)
12 12
 	checkEquals(deletionTest$ref, 'TCG')
13 13
 	checkEquals(deletionTest$altGt1, 'T')
14 14
 
15 15
 	## test insertion
16
-	insertionTest <- indelDiff(bam1,bam2,refFsa,'chr07',828514,828914,pValueCutOff=1,gtDistCutOff=0) 
17
-	checkEquals(insertionTest$pos, 828714)
18
-	checkEquals(insertionTest$ref, 'CAAAAAAA')
19
-	checkEquals(insertionTest$altGt1, 'CAAAAAAAA')
16
+	insertionTest <- indelDiff(bam1,bam2,refFsa,'chr07',828514,828914,pValueCutOff=1,gtDistCutOff=0,nCores=2) 
17
+	checkEquals(insertionTest$pos, c(828636,828714))
18
+	checkEquals(insertionTest$ref, c('GCCA','CAAAAAAA'))
19
+	checkEquals(insertionTest$altGt1, c('GCCACCA','CAAAAAAAA'))
20 20
 	
21 21
 	## test NULL
22
-	checkTrue(is.null(indelDiff(bam1,bam2,refFsa,'chr04',962501,1026983,pValueCutOff=1,gtDistCutOff=2,verbose=FALSE)))
22
+	checkTrue(is.null(indelDiff(bam1,bam2,refFsa,'chr04',962501,1026983,pValueCutOff=1,gtDistCutOff=2,verbose=FALSE,nCores=2)))
23 23
 	
24 24
 }
25 25
 
... ...
@@ -8,12 +8,12 @@ test_snpDiff <- function(){
8 8
 	refFsa <- system.file(package='SICtools','extdata','example.ref.fasta')
9 9
 	
10 10
 	## check
11
-	snpTest <- snpDiff(bam1,bam2,refFsa,'chr04',962501,1026983,pValueCutOff=1,baseDistCutOff=0)
11
+	snpTest <- snpDiff(bam1,bam2,refFsa,'chr04',962501,1026983,pValueCutOff=1,baseDistCutOff=0,nCores=2)
12 12
 	checkEquals(nrow(snpTest),16)
13 13
 	checkEquals(subset(snpTest,p.value < 1e-10)$pos,c(962801,1026683))
14 14
 	
15 15
 	## check NULL
16
-	checkTrue(is.null(snpDiff(bam1,bam2,refFsa,'chr04',962501,1026983,pValueCutOff=0,baseDistCutOff=0)))
16
+	checkTrue(is.null(snpDiff(bam1,bam2,refFsa,'chr04',962501,1026983,pValueCutOff=0,baseDistCutOff=0,nCores=2)))
17 17
 	
18 18
 }
19 19