r - 「summary.glht」からp値を抽出する方法

Question

Pr(>|t|) 列から値を抽出したい

library(lsmeans)    
warp.lm = lm(breaks ~ wool * tension, data = warpbreaks)
toP<-lsmeans(warp.lm, pairwise ~ wool | tension, glhargs=list())
toP[[2]]

         Simultaneous Tests for General Linear Hypotheses

Fit: lm(formula = breaks ~ wool * tension, data = warpbreaks)

Linear Hypotheses:
               Estimate Std. Error t value Pr(>|t|)   
A - B | L == 0   16.333      5.157   3.167  0.00797 **
A - B | M == 0   -4.778      5.157  -0.926  0.73187   
A - B | H == 0    5.778      5.157   1.120  0.60282   
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 
(Adjusted p values reported -- single-step method)

これを行う方法？class(toP[[2]])それを言う"summary.glht" "glht"。toP[[2]][9] には $pvalues がありますが、toP[[2]][9]$pvalues は NULL を返します

Answer 1

特定のオブジェクトにアクセスできる「要素」を知りたい場合はclass、使用せず、次を使用してnamesください。

R> names(toP[[2]])
[1] "model"       "linfct"      "rhs"         "coef"        "vcov"       
[6] "df"          "alternative" "type"        "test"       

ここでは、という要素があることがわかりますtest。それを見てみましょう：

R> names(toP[[2]]$test)
[1] "pfunction"    "qfunction"    "coefficients" "sigma"       
[5] "tstat"        "pvalues"      "type"        

うーん、という要素がありpvaluesます。いいですね。あなたはそれにアクセスすることができます：

R> toP[[2]]$test$pvalues
[1] 0.007954354 0.731843623 0.602840958
attr(,"error")
[1] 7.579683e-05

そして、ここであなたはあなたのp値を取得します...

オブジェクトの構造を取得する別の方法は、を使用することstr()です。ケース（str(toP[[2]])）に適用すると、出力が少し長くなりますが、p値にアクセスする方法を直接決定できる可能性があります。

Answer 2

試しましたか：

 toP[[2]]$test$pvalues