Question: error in "count" for SummarizedExperiment in easyRNASeq
gravatar for Lescai, Francesco
6.2 years ago by
Lescai, Francesco380 wrote:
Hi guys, I'm trying to use the function "count" in easyRNASeq to get a SummarizedExp to be used with DESeq2. I'm aware it's in active development, so just flagging this up in case I'm just doing something wrong Here's what I've done bamfiles<-dir("RNAseq/BAM files/", pattern="bam$") conditions <- c(rep("treated",8),rep("normal",8)) names(conditions) <- bamfiles rawCounts <- count(filesDirectory="RNAseq/BAM files/", pattern="bam$", organism="Mmusculus", annotationMethod="gtf", annotationFile="RNAseq/BAM files/Mus_musculus.GRCm38.72.gtf", count="genes", summarization="geneModels", conditions=conditions ) Checking arguments... Fetching annotations... Read 1142508 records Computing gene models... Summarizing counts... Processing Tophat2R1BAM.bam Updating the read length information. The reads are of 49 bp. [ ] Processing Tophat2W8BAM.bam Updating the read length information. The reads are of 49 bp. Preparing output Error in validObject(.Object) : invalid class “SummarizedExperiment” object: 'rowData' length differs from 'assays' nrow Any clues about this error? My session info below. thanks very much, Francesco > sessionInfo() R version 3.0.1 (2013-05-16) Platform: x86_64-apple-darwin10.8.0 (64-bit) locale: [1] en_GB.UTF-8/en_GB.UTF-8/en_GB.UTF-8/C/en_GB.UTF-8/en_GB.UTF-8 attached base packages: [1] parallel stats graphics grDevices utils datasets methods base other attached packages: [1] DESeq2_1.0.19 RcppArmadillo_0.3.910.0 Rcpp_0.10.4 easyRNASeq_1.6.0 ShortRead_1.18.0 latticeExtra_0.6-26 RColorBrewer_1.0-5 [8] Rsamtools_1.12.4 DESeq_1.12.1 lattice_0.20-23 locfit_1.5-9.1 BSgenome_1.28.0 GenomicRanges_1.12.5 Biostrings_2.28.0 [15] IRanges_1.18.3 edgeR_3.2.4 limma_3.16.7 biomaRt_2.16.0 Biobase_2.20.1 genomeIntervals_1.16.0 BiocGenerics_0.6.0 [22] intervals_0.14.0 BiocInstaller_1.10.3 loaded via a namespace (and not attached): [1] annotate_1.38.0 AnnotationDbi_1.22.6 bitops_1.0-6 DBI_0.2-7 genefilter_1.42.0 geneplotter_1.38.0 grid_3.0.1 hwriter_1.3 [9] RCurl_1.95-4.1 RSQLite_0.11.4 splines_3.0.1 stats4_3.0.1 survival_2.37-4 tools_3.0.1 XML_3.95-0.2 xtable_1.7-1 [17] zlibbioc_1.6.0 ________________________________ [[alternative HTML version deleted]]
easyrnaseq deseq2 • 708 views
ADD COMMENTlink written 6.2 years ago by Lescai, Francesco380
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 178 users visited in the last hour