Hi, I have a quick question regarding the merits of modeling data using a multivariate linear model versus a univariate linear model with an interaction term. Here's the data set. Say I have gene expression from two tissues liver and kidney and I'm interested in identifying differentially expressed genes in liver only, in kidney only or in both tissues. I could fit the following model:
expression ~ condition*tissue but I could also theoretically treat the expression levels as a two dimensional response and fit it as: Expression ~ condition where Expression is now a matrix although in the latter case I'm not sure how to test for differential expression in only one tissue, any ideas? Thanks, ~Jimmie [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.