...
|
...
|
@@ -110,6 +110,8 @@ test.MotifDb.emptyMode = function ()
|
110
|
110
|
# NA-JASPAR_CORE-HNF4A-MA0114.1: JASPAR gives <NA> for speciesID
|
111
|
111
|
# NA-JASPAR_CORE-CEBPA-MA0102.2: JASPAR gives '-' for speciesID, website says 'vertebrates'
|
112
|
112
|
|
|
113
|
+# Many more NA's exist...need to fix these
|
|
114
|
+
|
113
|
115
|
test.noNAorganisms = function ()
|
114
|
116
|
|
115
|
117
|
{
|
...
|
...
|
@@ -243,7 +245,9 @@ test.organisms = function ()
|
243
|
245
|
|
244
|
246
|
# jaspar_core has 3 NA speciesId: TBP, HNF4A and CEBPA (MA0108.2, MA0114.1, MA0102.2)
|
245
|
247
|
# their website shows these as vertebrates, which I map to 'Vertebrata'. An organismID of '-'
|
246
|
|
- # gets the same treatment, matching website also.
|
|
248
|
+ # gets the same treatment, matching website also.
|
|
249
|
+
|
|
250
|
+ # As in case of noNA, need to add organisms for these
|
247
|
251
|
checkEquals (which (is.na (mcols(MotifDb)$organism)), integer (0))
|
248
|
252
|
|
249
|
253
|
empty.count = length (which (mcols(mdb)$organism==""))
|
...
|
...
|
@@ -271,7 +275,7 @@ test.flyBindingDomains = function ()
|
271
|
275
|
checkEquals (tmp$Homeobox, 212)
|
272
|
276
|
checkEquals (tmp[['zf-C2H2']], 160)
|
273
|
277
|
checkEquals (tmp[["Helix-Turn-Helix"]], 182)
|
274
|
|
- checkEquals (length (which (is.na (subset (x, organism=='Dmelanogaster')$bindingDomain))), 162) # lots of cisbp
|
|
278
|
+ checkEquals (length (which (is.na (subset (x, organism=='Dmelanogaster')$bindingDomain))), 301) # lots of cisbp
|
275
|
279
|
|
276
|
280
|
} # test.flyBindingDomains
|
277
|
281
|
#------------------------------------------------------------------------------------------------------------------------
|
...
|
...
|
@@ -537,7 +541,7 @@ test.export_memeFormatToFileDuplication = function ()
|
537
|
541
|
print ('--- test.export_memeFormatToFileDuplication')
|
538
|
542
|
mdb = MotifDb # ()
|
539
|
543
|
mdb.mouse = subset (mdb, organism=='Mmusculus')
|
540
|
|
- checkEquals (length (mdb.mouse), 660)
|
|
544
|
+ checkEquals (length (mdb.mouse), 1251)
|
541
|
545
|
output.file = 'mouse.txt' # tempfile ()
|
542
|
546
|
max = 3
|
543
|
547
|
meme.text = export (mdb.mouse [1:max], output.file, 'meme')
|