Question

是否有一种简单的方法可以对摘要的元素进行后续数学计算？我有通过anova分析运行的日志转换数据。我想计算摘要输出的反对数。

我有以下代码：

require(multcomp)
inc <- log(Inc)
myanova <- aov(inc ~ educ)    
tukey <- glht(myanova, linfct = mcp(educ = "Tukey"))
summary(tukey)

产生如下输出：

                      Estimate Std. Error t value Pr(>|t|)    
12 - under12 == 0      0.32787    0.08493   3.861  0.00104 ** 
13to15 - under12 == 0  0.49187    0.08775   5.606  < 0.001 ***
16 - under12 == 0      0.89775    0.09217   9.740  < 0.001 ***
over16 - under12 == 0  0.99856    0.09316  10.719  < 0.001 ***
13to15 - 12 == 0       0.16400    0.04674   3.509  0.00394 ** 
etc.

如何在估算值上轻松执行反对数计算？

Answer 1

这有点像黑客，所以我建议进一步检查，但如果您只想查看指数估计和标准错误，我认为类似于以下内容将起作用（我使用了不同的数据）。

> amod <- aov(breaks ~ tension, data = warpbreaks)
> tukey = glht(amod, linfct = mcp(tension = "Tukey"))

> tsum = summary(tukey)
> tsum[[10]]$coefficients = exp(tsum[[10]]$coefficients)
> tsum[[10]]$sigma = exp(tsum[[10]]$sigma)
> tsum

如果您想使用coef（tukey）为您提供估算值，那么您将使用以下方法进行逆转换：

exp(coef(tukey))

Answer 2

我认为这应该有效：

      coef(tukey)

获取估算值。这里有一个例子：

  amod <- aov(breaks ~ tension, data = warpbreaks)
  tukey <- glht(amod, linfct = mcp(tension = "Tukey"))

现在，如果想要获取您键入的所有tukey摘要元素，请应用head或tail以获取包含摘要元素的命名列表。

head(summary(tukey))
$model
Call:
   aov(formula = breaks ~ tension, data = warpbreaks)

Terms:
                 tension Residuals
Sum of Squares  2034.259  7198.556
Deg. of Freedom        2        51

Residual standard error: 11.88058 
Estimated effects may be unbalanced

$linfct
      (Intercept) tensionM tensionH
M - L           0        1        0
H - L           0        0        1
H - M           0       -1        1
attr(,"type")
[1] "Tukey"

$rhs
[1] 0 0 0

$coef
(Intercept)    tensionM    tensionH 
   36.38889   -10.00000   -14.72222 

$vcov
            (Intercept)  tensionM  tensionH
(Intercept)    7.841564 -7.841564 -7.841564
tensionM      -7.841564 15.683128  7.841564
tensionH      -7.841564  7.841564 15.683128

$df
[1] 51

对摘要元素进行计算

2 个答案: