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)

Я пробовал использовать na.omit на тесте Hoslem, но он тоже не работает ... Кто-то может мне помочь?