|酒精| 0.082 || D |酮| 0.019 |+ ———– + ———- + ———— +
我已经在线性上运行了anova 模型 </跨度>
lm(比例〜治疗)
为每个人 复合 </跨度> 使用data.table 成功:
mytest&lt; - function(x){ 模型 </跨度> &lt; -lm(比例〜治疗,数据= t.df) aovmodel&LT; -aov( 模型 </跨度> ) tuks&LT; -TukeyHSD(aovmodel) }tapply((t.df [独特(t.df $ 复合 </跨度> )]),mytest的)
这将返回
要按照您的指定获得每种化合物的Tukey HSD,请尝试以下方法:
lapply(unique(t.df$Compound), function(x, df) TukeyHSD(aov(glm(Proportion ~ Treatment, data = df, subset = Compound == x)))[[1]], df = t.df)
对于每个独特的化合物,这个叫 TukeyHSD() 在ANOVA上,一般线性模型拟合对应于化合物的数据子集。它返回一个列表,其中每个元素对应一个化合物。
TukeyHSD()