2

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

我有以下代码：

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.

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

4

2 回答 2

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))
于 2013-02-28T05:19:19.277 回答
0

我认为这应该有效：

      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
于 2013-02-28T05:25:14.163 回答