0

我在 R 中运行一个混合效果模型，我想知道是否有代码我可以运行来获得对比对比，以进行事后测试的交互。

答：即。我目前正在运行的 tukey 测试是这样的：

有 2 个区域（ML 和 HG）和 3 个处理（Low_control、Frog 和 High_control）。

L.S <- pairs(lsmeans(juvweightmodelfeb6th, ~ Region*Treatment_three))
test(L.S, adjust = "tukey")

B. 这让我在输出中得到了这些比较：

HG,Low_control - ML,Low_control   
 HG,Low_control - HG,Frog           
 HG,Low_control - ML,Frog          
 HG,Low_control - HG,High_control   
 HG,Low_control - ML,High_control   
 ML,Low_control - HG,Frog           
 ML,Low_control - ML,Frog           
 ML,Low_control - HG,High_control   
 ML,Low_control - ML,High_control   
 HG,Frog - ML,Frog                 
 HG,Frog - HG,High_control          
 HG,Frog - ML,High_control          
 ML,Frog - HG,High_control          
 ML,Frog - ML,High_control          
 HG,High_control - ML,High_control

C. 但我真正想要的是这三个对比，本质上是对比的对比：

1. HG, Low_control - HG, Frog    -                 ML, Low_control - ML, Frog

2. HG, High_control - HG, Frog     -               ML, High_control - ML, Frog

3. HG, High_control - HG, Low_control     -   ML, High_control - ML, Low_control

D. 如何在 R 中的混合效应模型的事后测试中为上述三个对比编码？

谢谢！！

4

0 回答 0