Question

下面的代码和输出给出了我想要的解释 - 即底部的两个对比没有显着不同，但第三个与另外两个不同。如何找到用于此比较的比较和p值？

即。我的目标是测试对比度（1-2）@目标3明显小于其他两个对比。

mod<-lmer(Mind_avg ~ MindType*Target + 1 +  (1|Subject) + (1|Category),data=dat)
m3<-lsmeans(mod,~MindType|Target)
m3u<-update(pairs(m3,by="target"),by=NULL)
cld(m3u)

contrast Target  estimate       SE      df t.ratio p.value .group  
1 - 2     3       3.502885 1.020477 6521.04   3.433  0.0006  1      
1 - 2     2      11.559040 1.022071 6521.00  11.309  <.0001   2     
1 - 2     1      11.657720 1.016615 6521.00  11.467  <.0001   2

Answer 1

知道了。

   > pairs(m3u)
     contrast            estimate       SE      df t.ratio p.value
     1 - 2,1 - 1 - 2,2 0.09867925 1.441574 6521.00   0.068  0.9974
     1 - 2,1 - 1 - 2,3 8.15483466 1.440444 6521.02   5.661  <.0001
     1 - 2,2 - 1 - 2,3 8.05615541 1.444300 6521.02   5.578  <.0001

通过R中的cld获得比较的p值

1 个答案: