[BioC] zero-row result breaks select() on PolyPhen.Hsapiens.* and SIFT.Hsapiens.*

Robert Castelo robert.castelo at upf.edu
Sat Sep 21 02:25:21 CEST 2013


Dear list,

interrogating the TxDb.Hsapiens.UCSC.hg19.knownGene package with no 
result gives the following expected result:

library(TxDb.Hsapiens.UCSC.hg19.knownGene)
select(TxDb.Hsapiens.UCSC.hg19.knownGene, keys="dummy", 
keytype="GENEID", cols="SYMBOL")
[1] GENEID
<0 rows> (or 0-length row.names)

however, when i try the same with the annotation packages 
PolyPhen.Hsapiens.dbSNP131 and SIFT.Hsapiens.dbSNP132, the select 
instruction breaks into an error:

library(SIFT.Hsapiens.dbSNP132)
library(PolyPhen.Hsapiens.dbSNP131)

select(SIFT.Hsapiens.dbSNP132, keys=c("dummy"))
Error in data.frame(RSID = unlist(rsid), PROTEINID = unlist(protein_id),  :
   arguments imply differing number of rows: 1, 0

select(PolyPhen.Hsapiens.dbSNP131, keys="dummy")
Error in `*tmp*`$RSID : $ operator is invalid for atomic vectors

i guess these two annotation packages should work analogously to 
TxDb.Hsapiens.UCSC.hg19.knownGene, and give just a 0-row data.frame 
object, right?

these errors reproduce also with the current devel version of BioC, 
please find below both sessionInfo() outputs.

cheers,
robert.

=====RELEASE====
R version 3.0.1 (2013-05-16)
Platform: x86_64-apple-darwin10.8.0 (64-bit)

locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8

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

other attached packages:
  [1] TxDb.Hsapiens.UCSC.hg19.knownGene_2.9.2 GenomicFeatures_1.12.3
  [3] AnnotationDbi_1.22.6 Biobase_2.20.1
  [5] PolyPhen.Hsapiens.dbSNP131_1.0.2 SIFT.Hsapiens.dbSNP132_1.0.2
  [7] RSQLite_0.11.4 DBI_0.2-7
  [9] VariantAnnotation_1.6.7 Rsamtools_1.12.4
[11] Biostrings_2.28.0 GenomicRanges_1.12.5
[13] IRanges_1.18.3 BiocGenerics_0.6.0
[15] vimcom_0.9-8 setwidth_1.0-3
[17] colorout_1.0-0

loaded via a namespace (and not attached):
[1] biomaRt_2.16.0     bitops_1.0-6       BSgenome_1.28.0 
RCurl_1.95-4.1     rtracklayer_1.20.4
[6] stats4_3.0.1       tools_3.0.1        XML_3.95-0.2 zlibbioc_1.6.0



=====DEVEL=====
R version 3.0.1 (2013-05-16)
Platform: x86_64-apple-darwin10.8.0 (64-bit)

locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8

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

other attached packages:
  [1] TxDb.Hsapiens.UCSC.hg19.knownGene_2.9.2 GenomicFeatures_1.13.40
  [3] AnnotationDbi_1.23.23 Biobase_2.21.7
  [5] PolyPhen.Hsapiens.dbSNP131_1.0.2 SIFT.Hsapiens.dbSNP132_1.0.2
  [7] RSQLite_0.11.4 DBI_0.2-7
  [9] VariantAnnotation_1.7.46 Rsamtools_1.13.41
[11] Biostrings_2.29.19 GenomicRanges_1.13.44
[13] XVector_0.1.4 IRanges_1.19.37
[15] BiocGenerics_0.7.5 vimcom_0.9-8
[17] setwidth_1.0-3 colorout_1.0-0

loaded via a namespace (and not attached):
[1] biomaRt_2.17.3      bitops_1.0-6        BSgenome_1.29.1 
RCurl_1.95-4.1      rtracklayer_1.21.12
[6] stats4_3.0.1        tools_3.0.1         XML_3.95-0.2 zlibbioc_1.7.0



More information about the Bioconductor mailing list