從glht中提取t-stat在R

Question

我在R中運行了glht的假設測試，我想提取測試的t-stat。我讀過那裏應該有我接受的「glht」類中的一個元素「測試」，但是當我通過運行glht得到的元素不會出現。

的代碼非常簡單，下面給出：

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

運行的最後一行給我： [1] 「模範」， 「linfct」 「RHS」 「COEF」 「vcov」， 「DF」 「替代」「類型」

任何人都有一個答案呢？

謝謝。

Answer 1

使用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"