[BioC] problem with TXNAME -> SYMBOL mapping in Homo.sapiens library

Aleksandra Pfeifer [guest] guest at bioconductor.org
Fri Sep 27 11:28:44 CEST 2013


Hello,
  I have a problem with the maping from txname to symbol of the gene. For most transcripts it works ok, but for some it doesn't:

> library(Homo.sapiens)
> select(Homo.sapiens, cols="SYMBOL", keys= "uc021wml.1", keytype="TXNAME")
Error in .testIfKeysAreOfProposedKeytype(x, keys, keytype) :
  None of the keys entered are valid keys for the keytype specified.

The traceback is as follows:
> traceback()
10: stop("None of the keys entered are valid keys for the keytype specified.")
9: .testIfKeysAreOfProposedKeytype(x, keys, keytype)
8: .select(x, keys, cols, keytype, jointype = jointype)
7: .local(x, keys, cols, keytype, ...)
6: select(.makeReal(nodeName), keys = fromKeys, cols = needCols[[nodeName]],
       keytype = toKey)
5: select(.makeReal(nodeName), keys = fromKeys, cols = needCols[[nodeName]],
       keytype = toKey)
4: .getSelects(x, keytype, keys, needCols, visitNodes)
3: .select(x, keys, cols, keytype, ...)
2: select(Homo.sapiens, cols = "SYMBOL", keys = "uc021wml.1", keytype = "TXNAME")
1: select(Homo.sapiens, cols = "SYMBOL", keys = "uc021wml.1", keytype = "TXNAME")


However, When I try to check whether the problematic txname is present in Homo.sapiens database, it occurs that it is there. I can also find some other information about this transcript:
> "uc021wml.1" %in% keys(Homo.sapiens, keytype="TXNAME")
[1] TRUE
> select(Homo.sapiens, cols="TXSTART", keys= "uc021wml.1", keytype="TXNAME")
      TXNAME  TXSTART
1 uc021wml.1 22385572

Is there a way to solve that problem? I would be appreciated for your help.

Best regards,
Aleksandra Pfeifer



 -- output of sessionInfo(): 

> sessionInfo()
R version 3.0.1 (2013-05-16)
Platform: x86_64-unknown-linux-gnu (64-bit)

locale:
[1] C

attached base packages:
[1] parallel  stats     graphics  grDevices utils     datasets  methods
[8] base

other attached packages:
 [1] Homo.sapiens_1.1.1
 [2] TxDb.Hsapiens.UCSC.hg19.knownGene_2.9.2
 [3] org.Hs.eg.db_2.9.0
 [4] GO.db_2.9.0
 [5] RSQLite_0.11.4
 [6] DBI_0.2-7
 [7] OrganismDbi_1.2.0
 [8] GenomicFeatures_1.12.4
 [9] GenomicRanges_1.12.5
[10] IRanges_1.18.4
[11] AnnotationDbi_1.22.6
[12] Biobase_2.20.1
[13] BiocGenerics_0.6.0

loaded via a namespace (and not attached):
 [1] BSgenome_1.28.0    Biostrings_2.28.0  RBGL_1.36.2        RCurl_1.95-4.1
 [5] Rsamtools_1.12.4   XML_3.98-1.1       biomaRt_2.16.0     bitops_1.0-6
 [9] graph_1.38.3       rtracklayer_1.20.4 stats4_3.0.1       tools_3.0.1
[13] zlibbioc_1.6.0
>


--
Sent via the guest posting facility at bioconductor.org.



More information about the Bioconductor mailing list