Entering edit mode
Jeremy Leipzig
▴
70
@jeremy-leipzig-4924
Last seen 7.0 years ago
BiocInstaller::biocLite(("RNAseqData.HNRNPC.bam.chr14"))
BiocInstaller::biocLite("Rsamtools")
library(Rsamtools)
library("RNAseqData.HNRNPC.bam.chr14")
fl <- RNAseqData.HNRNPC.bam.chr14_BAMFILES[1]
res <- pileup(fl)
head(res)
> res <- pileup(fl)
>
> head(res)
[1] seqnames pos strand nucleotide count
<0 rows> (or 0-length row.names)
> table(res$strand, res$nucleotide)
A C G T N = - +
+ 0 0 0 0 0 0 0 0
- 0 0 0 0 0 0 0 0
* 0 0 0 0 0 0 0 0
> sessionInfo()
R version 3.4.1 (2017-06-30)
Platform: x86_64-redhat-linux-gnu (64-bit)
Running under: Amazon Linux AMI 2017.09
Matrix products: default
BLAS/LAPACK: /efs/R/lib/libRblas.so
locale:
[1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C LC_TIME=en_US.UTF-8 LC_COLLATE=en_US.UTF-8
[5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8 LC_PAPER=en_US.UTF-8 LC_NAME=C
[9] LC_ADDRESS=C LC_TELEPHONE=C LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
attached base packages:
[1] stats4 parallel stats graphics grDevices utils datasets methods base
other attached packages:
[1] RNAseqData.HNRNPC.bam.chr14_0.12.0 Rsamtools_1.26.2 Biostrings_2.42.1
[4] XVector_0.14.1 GenomicRanges_1.26.4 GenomeInfoDb_1.10.3
[7] IRanges_2.8.2 S4Vectors_0.12.2 BiocGenerics_0.20.0
loaded via a namespace (and not attached):
[1] pcaPP_1.9-72 Rcpp_0.12.11 DEoptimR_1.0-8 compiler_3.4.1 BiocInstaller_1.24.0 RColorBrewer_1.1-2
[7] zlibbioc_1.20.0 bitops_1.0-6 tools_3.4.1 dotCall64_0.9-5 lattice_0.20-35 graph_1.52.0
[13] mvtnorm_1.0-6 spam_2.1-2 hexbin_1.27.1 cluster_2.0.6 IDPmisc_1.1.17 fields_9.0
[19] maps_3.2.0 grid_3.4.1 robustbase_0.92-8 Biobase_2.34.0 rrcov_1.4-3 R6_2.2.2
[25] BiocParallel_1.8.2 latticeExtra_0.6-28 corpcor_1.6.9 matrixStats_0.52.2 flowFramePlus_0.1.0 MASS_7.3-48
[31] assertthat_0.2.0 flowCore_1.40.6 KernSmooth_2.23-15 flowViz_1.38.0 RCurl_1.95-4.9

same after update to Rsamtools_1.30.0 and RNAseqData.HNRNPC.bam.chr14_0.16.0