1

R で glht を使用して仮説検定を実行しましたが、検定の t-stat を抽出したいと考えています。受け取ったクラス「glht」に要素「test」があるはずだと読みましたが、glht を実行して取得した要素を調べると表示されません。

コードは非常に単純で、以下に示します。

reg = lm(dep ~ indep)
htest = glht(reg,linfct = c("indep  = 0.5"))
names(htest)

最後の行を実行すると、次のようになります: [1] "model" "linfct" "rhs" "coef" "vcov" "df" "alternative" "type"

誰もそれに対する答えを持っていますか？

ありがとうございました。

4

1 に答える 1

0

使用summary:

indep <- 1:10
set.seed(42)
dep <- indep/2+5+rnorm(10)

reg = lm(dep ~ indep)

library(multcomp)
htest = glht(reg,linfct = c("indep  = 0.5"))
summary(htest)
#Simultaneous Tests for General Linear Hypotheses
#
#Fit: lm(formula = dep ~ indep)
#
#Linear Hypotheses:
#             Estimate Std. Error t value Pr(>|t|)
#indep == 0.5  0.53040    0.09697   0.313    0.762
#(Adjusted p values reported -- single-step method)

次のように値を抽出できます。

res <- summary(htest)
res$test[-(1:2)]
# $coefficients
# indep 
# 0.5303967 
# 
# $sigma
# indep 
# 0.09696568 
# 
# $tstat
# indep 
# 0.3134785 
# 
# $pvalues
# [1] 0.7619346
# attr(,"error")
# [1] 0
# 
# $type
# [1] "single-step"
于 2013-10-29T16:35:55.027 に答える