模型框架检错中的R Hosmer-Lemeshow错误

Question

我有一个数据框架，其中包含来自患者的119个观察值和565个变量。对于回归线性模型，我感兴趣的变量中有一些NA。 我首先做了线性回归模型：

RL1 <- glm(x~Var1+
         Var2+
         Var3+
         Var4, Dataframe, family = "binomial")
summary(RL1)
enteexp(coefficients(RL1))

然后：

library(ResourceSelection)
hoslem.test(Dataframe$x, fitted(RL1), g=4)

但是我有这个答案：

Error in model.frame.default(formula = cbind(y0 = 1 - y, y1 = y) ~ cutyhat)

我已经在hoslem测试中尝试使用na.omit了，但是还是行不通...有人可以帮助我吗？