Провести проверку линейной гипотезы о предполагаемых коэффициентах логистической регрессии в R

Question

У меня логистическая регрессия формы

, где

Соответствующий код R следующий.

Regression_hat<-glm(Y~X,family=binomial(link='logit'))

Я знаю, что для проверки совместной гипотезы

тест linearHypothesis может использоваться со следующей формой:

linearHypothesis(Regression_hat,c("(Intercept)=0","X=0"),test=c("F"))

Я хотел бы проверить (совместно) гипотезу о том, что оба оценочных параметра равны двум произвольным значениям соответственно. Эти два значения хранятся под разными именами переменных (например, V1 и V2) в R, хотя использование имен переменных в приведенном выше коде (linearHypothesis(Regression_hat,c("(Intercept)=V1","X=V2"),test=c("F"))) не работает.

Answer 1

Попробуйте

linearHypothesis(Regression_hat, paste(c("(Intercept)", "X"), "=", c(V1, V2)), test = "F")

Проблема в том, что запись, например, "X=V2" не делает X2 похожим на переменную; это только часть этого персонажа. Использование paste помогает вам построить, например, "X=3", когда V1 принимает значение 3.