R中的ols命令（rms软件包）的linearHypothesis等效项

Question

我正在尝试使用“ car”包中的“ linearHypothesis”函数来测试用“ rms”包中的“ ols”估算的模型的系数。该函数适用于“ lrm”对象，但不适用于“ ols”对象。您有其他选择吗？我知道使用“ lm”可以解决问题，但是我想使用“ ols”，因为在这里更容易得到群集的标准错误。

Answer 1

您可以使用glht软件包中的multcomp。

library(rms)
library(multcomp)

d <- datadist(swiss); options(datadist="d")
fit <- ols(Fertility ~ ., data = swiss)
summary(fit)

test <- glht(fit, linfct = "Agriculture = 0")
summary(test)
# Fit: ols(formula = Fertility ~ ., data = swiss, x = TRUE)
# 
# Linear Hypotheses:
#                  Estimate Std. Error z value Pr(>|z|)  
# Agriculture == 0  -0.1721     0.0703  -2.448   0.0144 *
#   ---
#   Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1