多因素方差分析
1、准备数据
至少需要4列，第一列是蛋白或者基因的id，第二列是相应的表达量，第三列是分组信息，第四列是形状信息。三四两列就是要分析的因素，这里我用R包里面的数据来代替
attach(ToothGrowth)
table(supp, dose)
ToothGrowth的数据格式如下，如果配置自己的数据，则只需要配置成下面的形式
fit <- aov(len ~ supp*dose,data=ToothGrowth)
中间的 * 代表是考虑supp和dose之间的交互作用，如果将 * 替换成+则不考虑交互作用
summary(fit)
supp：dose为两者之间的交互作用，后面的*代表显著
df为自由度，sum Sq为总的离差平方和，Mean Sq为平均离差平方和，F value为F检验的F值，Pr即为p,如果p<0.05,则分组对其是有影响的
