Линейный эквивалент гипотезы для команды ols (rms package) в R

Question

Я пытаюсь использовать функцию «linearHypothesis» из пакета «car» для проверки коэффициентов модели, оцененной с помощью «ols» из пакета «rms».Функция работает с объектами "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