Dear Dr. González,

I am interested in testing tweeDEseq on my own RNA-seq data. I have 
noticed that tweeDEseq supports glm fitting for more complex designs 
than simple pairwise models. However, it seems that the tweeDEglm 
function only supports a null model of "count ~ 1", which makes it 
impossible to do do things like test for an interaction term by 
comparing "count ~ var1 * var2" against "count ~ var1 + var2". Looking 
at the code for tweeDEseq:::anova.glmPT, it looks like the underlying 
mechanisms have support fur such comparisons. Are there any future plans 
to support null hypotheses other than "count ~ 1" in the tweeDEglm function?

Thank you,

-Ryan Thompson

	[[alternative HTML version deleted]]

