Hi,
I am hoping that someone might be able to give me some advice on adjusting or controlling for batch effects in RNA-Seq data using Ballgown?
I have got RNA-seq data for different batches of sequencing for the replicates of the sample. I followed new tuxedo pipeline for my preliminary analysis of the data and found that there were difference in total number of PE reads between the two batches of the sequencing (batch 1: 27 million reads; batch 2: 20 million reads)- therefore I think this difference could be due to batch effects within my dataset.
Is it possible to control for batch effects within Ballgown pipeline using the stattest function?
Any advice on this would be gratefully appreciated!
Many thanks,
Amy
Hi everyone,
Also, to add to my question above, will this difference in number of reads cause a batch effect while estimating the transcripts for each samples. The reason, I am asking this question is that because using stringtie merge function, I merged the transcripts from all samples and then used it that merge file to estimate the transcript for each sample.
-Amy