2016-10-28 6 views
0

以下のコードと出力は、私が望む解釈を与えます。つまり、下の2つのコントラストは大きく異なるわけではありませんが、3番目のコントラストは他の2つと異なります。この比較に使用される比較値とp値は、どのようにして見つけることができますか?Rでのcldによる比較のためのp値の取得

私の目標は、コントラスト(1-2)@ターゲット3が他の2つのコントラストよりも大幅に小さいことをテストすることです。

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 

答えて

0

これが判明しました。

> 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 
関連する問題