Question: read count normalization for spike
gravatar for Sara
3.9 years ago by
Sara150 wrote:

I have 5 spikes (short sequence) and 4 samples (RNA-SEQ), one of them is control and 3 different treatment. I used different concentration for each spike but concentration of each spike in each sample is equal. now I counted the spikes in each sample. in two files the counts are almost the same but the other two samples have different spike counts. meaning in one of them they are 9 times less than control (comparing all 5 spikes between two samples) and the other one is 3 times less than control. now I am going to normalize the read counts using spikes? do you know how I can do that?

rna-seq • 1.4k views
ADD COMMENTlink modified 8 months ago by Biostar ♦♦ 20 • written 3.9 years ago by Sara150

Trying to normalize with only 5 spike-ins sounds like a recipe for disaster. It there a reason that you actually need to normalize with the spike-ins? If so, which bioconductor package are you using for the analysis?

ADD REPLYlink written 3.9 years ago by Devon Ryan98k

do you think one is enough? I am planning to use edgeR. what do you think? is that fine?

ADD REPLYlink written 3.9 years ago by ashkan110

Without a good reason otherwise, I would aim for more like 100... Then again, one shouldn't use spike-ins without a good reason.

ADD REPLYlink written 3.9 years ago by Devon Ryan98k
Please log in to add an answer.


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