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 值的方式。