Browse code

latest available flybase identifers no longer map FBgn029750 to ab

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

p.shannon authored on 04/02/2014 19:29:41
Showing 1 changed files

... ...
@@ -151,12 +151,21 @@ test.geneIdsAndTypes = function ()
151 151
 {
152 152
   print ('--- test.getGeneIdsAndTypes')
153 153
   mdb = MotifDb
154
-  geneIds = mcols(mdb)$geneId
155
-  geneIdTypes = mcols(mdb)$geneIdType
154
+  tbl <- mcols(mdb)
155
+  geneIds = tbl$geneId
156
+  geneIdTypes = tbl$geneIdType
156 157
   typeCounts = as.list (table (geneIdTypes))
157
-  checkEquals(typeCounts$ENTREZ, 2197)
158
-  checkEquals(typCounts$FLYBASE, 30)
159
-  checkEquals(typeCounts$SGD, 453)
158
+
159
+  #browser()
160
+  #fout <- file("/tmp/pshannon.txt", "w")
161
+  #write("ENTREZ", file=fout)
162
+  #write(typeCounts$ENTREZ, file=fout)
163
+  
164
+  checkTrue(typeCounts$ENTREZ == 2180)
165
+  checkTrue(typeCounts$FLYBASE >= 47)
166
+  checkTrue(typeCounts$SGD >= 453)
167
+  checkEquals(nrow(subset(tbl, is.na(geneIdType))), 932)
168
+  
160 169
 
161 170
   #checkEquals(typeCounts, list(ENTREZ=2197, FLYBASE=30, SGD=453, comment=683))
162 171
 
... ...
@@ -685,7 +694,8 @@ test.flyFactorGeneSymbols <- function()
685 694
 {
686 695
     print ("--- test.flyFactorGeneSymbols")
687 696
     mdb = MotifDb
688
-    checkEquals(mcols(query(mdb, "FBgn0259750"))$geneSymbol, c("ab", "ab"))
697
+    checkEquals(mcols(query(mdb, "FBgn0259750"))$geneSymbol,
698
+                c("FBgn0259750", "FBgn0259750"))
689 699
     checkEquals(mcols(query(mdb, "FBgn0000014"))$geneSymbol, rep("abd-A", 3))
690 700
     checkEquals(mcols(query(mdb, "FBgn0000015"))$geneSymbol, rep("Abd-B", 3))
691 701