Получение p-значений для сравнения по cld в R
Приведенный ниже код и выходные данные дают мне требуемую интерпретацию - то есть два нижних контраста существенно не отличаются, но третий отличается от двух других. Как я могу найти сравнения и 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
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