git-svn-id: https://hedgehog.fhcrc.org/bioconductor/trunk/madman/Rpacks/MotifDb@86014 bc3139a8-67e5-0310-9ffc-ced21a209358
... | ... |
@@ -195,7 +195,7 @@ createMetadataTable = function (matrices, novels) |
195 | 195 |
novels.ordered = novels [tbl.md$providerName] # make sure we follow the order in the tbl |
196 | 196 |
novelPFM [which (novels.ordered)] = 'novelMotif' |
197 | 197 |
tbl.md$geneId = novelPFM |
198 |
- tbl.md$geneIdType = rep ('comment', nrow (tbl.md)) |
|
198 |
+ tbl.md$geneIdType = rep (NA_character_, nrow (tbl.md)) |
|
199 | 199 |
|
200 | 200 |
invisible (tbl.md) |
201 | 201 |
|
... | ... |
@@ -93,13 +93,13 @@ test.createAnnotationTable = function () |
93 | 93 |
#------------------------------------------------------------------------------------------------------------------------ |
94 | 94 |
test.createMetadataTable = function (x.matrices, x.novels) |
95 | 95 |
{ |
96 |
- print ('--- test.createMetadataTable') |
|
97 |
- # try it first with just two matrices |
|
98 |
- tbl.md = createMetadataTable (x.matrices [1:12], x.novels [1:12]) |
|
99 |
- checkEquals (dim (tbl.md), c (12, 15)) |
|
100 |
- checkEquals (colnames (tbl.md), c ("providerName", "providerId", "dataSource", "geneSymbol", "geneId", "geneIdType", |
|
101 |
- "proteinId", "proteinIdType", "organism", "sequenceCount", "bindingSequence", |
|
102 |
- "bindingDomain", "tfFamily", "experimentType", "pubmedID")) |
|
96 |
+ print ('--- test.createMetadataTable') |
|
97 |
+ # try it first with just two matrices |
|
98 |
+ tbl.md = createMetadataTable (x.matrices [1:12], x.novels [1:12]) |
|
99 |
+ checkEquals (dim (tbl.md), c (12, 15)) |
|
100 |
+ checkEquals (colnames (tbl.md), c ("providerName", "providerId", "dataSource", "geneSymbol", "geneId", "geneIdType", |
|
101 |
+ "proteinId", "proteinIdType", "organism", "sequenceCount", "bindingSequence", |
|
102 |
+ "bindingDomain", "tfFamily", "experimentType", "pubmedID")) |
|
103 | 103 |
checkEquals (tbl.md$providerName [1:2], c ('UW.Motif.0001', 'UW.Motif.0002')) |
104 | 104 |
checkEquals (tbl.md$providerId [1:2], c ('UW.Motif.0001', 'UW.Motif.0002')) |
105 | 105 |
checkEquals (tbl.md$pubmedID [1:2], c ('22959076', '22959076')) |
... | ... |
@@ -107,8 +107,9 @@ test.createMetadataTable = function (x.matrices, x.novels) |
107 | 107 |
checkEquals (tbl.md$organism [1:2], c ('Hsapiens', 'Hsapiens')) |
108 | 108 |
checkEquals (tbl.md$experimentType [1:2], c ('digital genomic footprinting', 'digital genomic footprinting')) |
109 | 109 |
checkEquals (tbl.md$geneId, c (rep ('knownMotif', 11), 'novelMotif')) |
110 |
+ checkTrue(all(is.na(tbl.md$geneIdType))) |
|
110 | 111 |
|
111 |
- invisible (tbl.md) |
|
112 |
+ invisible (tbl.md) |
|
112 | 113 |
|
113 | 114 |
} # test.createMetadataTable |
114 | 115 |
#------------------------------------------------------------------------------------------------------------------------ |
... | ... |
@@ -154,10 +154,14 @@ test.geneIdsAndTypes = function () |
154 | 154 |
geneIds = mcols(mdb)$geneId |
155 | 155 |
geneIdTypes = mcols(mdb)$geneIdType |
156 | 156 |
typeCounts = as.list (table (geneIdTypes)) |
157 |
- checkEquals(typeCounts, list(ENTREZ=2197, FLYBASE=30, SGD=453, comment=683)) |
|
157 |
+ checkEquals(typeCounts$ENTREZ, 2197) |
|
158 |
+ checkEquals(typCounts$FLYBASE, 30) |
|
159 |
+ checkEquals(typeCounts$SGD, 453) |
|
158 | 160 |
|
159 |
- na.count = length (which (is.na (geneIds))) |
|
160 |
- checkEquals (na.count, 304) |
|
161 |
+ #checkEquals(typeCounts, list(ENTREZ=2197, FLYBASE=30, SGD=453, comment=683)) |
|
162 |
+ |
|
163 |
+ #na.count = length (which (is.na (geneIds))) |
|
164 |
+ #checkEquals (na.count, 304) |
|
161 | 165 |
empty.count = length (which (geneIds == '')) |
162 | 166 |
checkEquals (empty.count, 0) |
163 | 167 |
|