Browse code

updated a few counts to accomodate jaspar2014

git-svn-id: https://hedgehog.fhcrc.org/bioconductor/trunk/madman/Rpacks/MotifDb@88739 bc3139a8-67e5-0310-9ffc-ced21a209358

p.shannon authored on 10/04/2014 17:21:34
Showing 1 changed files

... ...
@@ -156,21 +156,13 @@ test.geneIdsAndTypes = function ()
156 156
   geneIdTypes = tbl$geneIdType
157 157
   typeCounts = as.list (table (geneIdTypes))
158 158
 
159
-  #browser()
160
-  #fout <- file("/tmp/pshannon.txt", "w")
161
-  #write("ENTREZ", file=fout)
162
-  #write(typeCounts$ENTREZ, file=fout)
163 159
   
164
-  checkTrue(typeCounts$ENTREZ == 2180)
160
+  checkTrue(typeCounts$ENTREZ == 2347)
165 161
   checkTrue(typeCounts$FLYBASE >= 47)
166
-  checkTrue(typeCounts$SGD >= 453)
167
-  checkEquals(nrow(subset(tbl, is.na(geneIdType))), 932)
162
+  checkTrue(typeCounts$SGD >= 629)
163
+  checkEquals(nrow(subset(tbl, is.na(geneIdType))), 1181)
168 164
   
169 165
 
170
-  #checkEquals(typeCounts, list(ENTREZ=2197, FLYBASE=30, SGD=453, comment=683))
171
-
172
-  #na.count = length (which (is.na (geneIds)))
173
-  #checkEquals (na.count, 304) 
174 166
   empty.count = length (which (geneIds == ''))
175 167
   checkEquals (empty.count, 0)
176 168
 
... ...
@@ -278,7 +270,7 @@ test.flyBindingDomains = function ()
278 270
 
279 271
   checkEquals (tmp$Homeobox, 212)
280 272
   checkEquals (tmp[['zf-C2H2']], 160)
281
-  checkEquals (tmp$HLH, 109)
273
+  checkEquals (tmp[["Helix-Turn-Helix"]], 182)
282 274
   checkEquals (length (which (is.na (subset (x, organism=='Dmelanogaster')$bindingDomain))), 24)
283 275
 
284 276
 } # test.flyBindingDomains
... ...
@@ -514,7 +506,7 @@ test.export_memeFormatStdOut = function ()
514 506
   print ('--- test.export_memeFormatStdOut')
515 507
   mdb = MotifDb # ()
516 508
   mdb.chicken = subset (mdb, organism=='Gallus')
517
-  checkEquals (length (mdb.chicken), 2)
509
+  checkEquals (length (mdb.chicken), 3)
518 510
     # text is cat-ed to stdout, so not avaialable here to check.
519 511
     # but just like print, export also returns the text invisibly.
520 512
     # so that CAN be checked.
... ...
@@ -532,7 +524,7 @@ test.export_memeFormatToFile = function ()
532 524
   print ('--- test.export_memeFormatToFile')
533 525
   mdb = MotifDb # ()
534 526
   mdb.chicken = subset (mdb, organism=='Gallus')
535
-  checkEquals (length (mdb.chicken), 2)
527
+  checkEquals (length (mdb.chicken), 3)
536 528
   output.file = tempfile ()
537 529
   meme.text = export (mdb.chicken, output.file, 'meme')
538 530
   retrieved = scan (output.file, what=character (0), sep='\n', quiet=TRUE)
... ...
@@ -545,7 +537,7 @@ test.export_memeFormatToFileDuplication = function ()
545 537
   print ('--- test.export_memeFormatToFileDuplication')
546 538
   mdb = MotifDb # ()
547 539
   mdb.mouse = subset (mdb, organism=='Mmusculus')
548
-  checkEquals (length (mdb.mouse), 462)
540
+  checkEquals (length (mdb.mouse), 528)
549 541
   output.file = 'mouse.txt' # tempfile ()
550 542
   max = 3
551 543
   meme.text = export (mdb.mouse [1:max], output.file, 'meme')