[BioC] GenomicRanges:::.filterMatchMatrix blows up nearest()
Harris A. Jaffee
hj at jhu.edu
Thu Oct 25 21:09:26 CEST 2012
The statement
m[, 1L] <- map[m[, 1L]]
seems pretty likely to cause
Error in map[m[, 1L]] : object of type 'closure' is not subsettable
This occurs with nearest(select="all") when there is no overlap.
Slightly sanitized real-life example:
> library(TxDb.Hsapiens.UCSC.hg19.knownGene)
Loading required package: GenomicFeatures
Loading required package: BiocGenerics
...
> TT <- transcriptsBy(TxDb.Hsapiens.UCSC.hg19.knownGene, by="gene")
> tt <- TT[[1]]
> tt
GRanges with 2 ranges and 2 metadata columns:
seqnames ranges strand | tx_id tx_name
<Rle> <IRanges> <Rle> | <integer> <character>
[1] chr19 [58858172, 58864865] - | 68796 uc002qsd.4
[2] chr19 [58859832, 58874214] - | 68797 uc002qsf.2
---
seqlengths:
chr1 chr2 ... chrUn_gl000249
249250621 243199373 ... 38502
> q <- GRanges(IRanges(1,1), seqnames="chr19", strand="*")
> nearest(q, tt)
[1] 1
> nearest(q, tt, select="all")
Error in map[m[, 1L]] : object of type 'closure' is not subsettable
> sessionInfo()
R version 2.15.1 Patched (2012-07-01 r59713)
Platform: x86_64-unknown-linux-gnu (64-bit)
locale:
[1] LC_CTYPE=en_US.iso885915 LC_NUMERIC=C
[3] LC_TIME=en_US.iso885915 LC_COLLATE=en_US.iso885915
[5] LC_MONETARY=en_US.iso885915 LC_MESSAGES=en_US.iso885915
[7] LC_PAPER=C LC_NAME=C
[9] LC_ADDRESS=C LC_TELEPHONE=C
[11] LC_MEASUREMENT=en_US.iso885915 LC_IDENTIFICATION=C
attached base packages:
[1] stats graphics grDevices datasets utils methods base
other attached packages:
[1] TxDb.Hsapiens.UCSC.hg19.knownGene_2.8.0
[2] GenomicFeatures_1.10.0
[3] AnnotationDbi_1.20.2
[4] Biobase_2.18.0
[5] GenomicRanges_1.10.2
[6] IRanges_1.16.3
[7] BiocGenerics_0.4.0
loaded via a namespace (and not attached):
[1] biomaRt_2.14.0 Biostrings_2.26.2 bitops_1.0-4.1 BSgenome_1.26.1
[5] DBI_0.2-5 parallel_2.15.1 RCurl_1.95-1.1 Rsamtools_1.10.1
[9] RSQLite_0.11.2 rtracklayer_1.18.0 stats4_2.15.1 tools_2.15.1
[13] XML_3.95-0.1 zlibbioc_1.4.0
More information about the Bioconductor
mailing list