View Single Post
Old 04-27-2016, 11:03 PM   #6
Devon Ryan
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480

By intersecting lists you're not measuring anything, just asking, "if I take two sets of results and intersect them, what do I get?" That's not answering the question, "what's the additional effect of condition on genotype B?".

Regarding the code, you can skip the LRT stuff and just use "results(ddsM)" and specify the interaction coefficient. For the actual LRT test, you're using the same formula for the full and reduced design. You want just "reduced = ~condition+genotype". "resMTolFact" is never defined. You can skip all of the "ddsMN" definition stuff.
dpryan is offline   Reply With Quote