[BioC] easyRNASeq read counts varying when recounting same sample!

[Nicolas Delhomme]

Hi René,

This is indeed fairly strange and nothing I have seen before. Can you please post your full script as well as your sessionInfo()? The console log (warnings, etc) would be helpful as well. I might have to have a go at your data; do you think that this could be arranged (off-list obviously)?

Cheers,

Nico

---------------------------------------------------------------
Nicolas Delhomme

Genome Biology Computational Support

European Molecular Biology Laboratory

Tel: +49 6221 387 8310
Email: nicolas.delhomme at embl.de
Meyerhofstrasse 1 - Postfach 10.2209
69102 Heidelberg, Germany
---------------------------------------------------------------





On 14 Feb 2013, at 16:26, René Böttcher wrote:

> Dear Nicolas,
> 
> I noticed that my read counts are varying when I count the exact same sample twice. Specifically, genes that have counted reads in the first run are not represented at all (0 counts) in the second run, while others have exactly the same counts!
> 
> I attached a shortened list of what my count tables look like in a tsv file, please note that I recounted the samples with the same options and annotation was provided via biomaRt.
> To explain the attached .tsv-file, I recounted the samples in order to make a comparison of Group1 vs Group1 treated and Group1 vs Group2 using separate objects.
> 
> You can easily see the differences in reads counts in the first sample of both comparisons (Group1_1412.bam).
> 
> Here's how I used the easyRNASeq function:
> 
> DGElist = easyRNASeq(organism="Hsapiens",
> 		                annotationMethod="biomaRt",
> 		                gapped=TRUE,
> 				count="genes",
> 		                summarization="geneModels",
> 		                filesDirectory=getwd(),
> 				filenames=bamfiles,
> 				recursive=T,
> 		                normalize=F,
> 				outputFormat="edgeR",
> 				conditions=conditions)
> 
> I hope this problem can be solved easily, because I already handed the list of differentially expressed genes over to a colleague of mine who wants to validate these genes now.
> 
> Best regards,
> René
> <readCount1.tsv>