Thank you for your answer, it's been very useful. I was hoping to ask you for some additional help. I also have a RIP-seq experiment with three assays - IP, input and IgG non-specific antibody, and two conditions - A and B.
So far I have compared input-A vs IP-A (filtering q<0.05 and FC 1.2, let's call this real-IP-A) and input-A vs IgG-A (filtering q<0.05 and FC 1.2, let's call this real-IgG-A). Then I compare real-IP-A vs. real-IP-B (filtering q<0.05 and FC 1.5), and I compare real-IgG-A vs. real-IgG-B. Any genes in the IgG list I remove from the IP list (there's quite a few to remove).
Do you think this is the best way to handle the data, or should I be doing a differences-in-differences as you mentioned above? In my case, the non-specific binding is not the same in A and B. In B, I appear to pick up more non-specific RNA.
Also, do you think that a FC cut off of 1.2 is appropriate when comparing input vs. IP or would you go higher?
I really look forward to your reply, and appreciate your help in advance.
Best wishes,
Jo
Thank you so much for your answer! I apologize for such a late reply but the I was finally able to go back to the lab so the analysis part of the project got a bit delayed.
That's very helpful and confirms what I thought about possibly dropping the IgG when comparing my treated vs untreated samples! I assumed that any non-specific binding would be comparable across the samples so I thought it could be omitted from the analysis. However, I was told by a fellow DIY bioinformatician that I should still include it and that's where my confusion came from.
Once again, thank you for your input (pun intended).
Please use
ADD COMMENT/ADD REPLYwhen responding to existing posts to keep threads logically organized. This comment belongs under @Ian's answer.SUBMIT ANSWERis for new answers to original question.