У меня неожиданная ошибка в моем исследовании. Позвольте мне показать вам несколько фрагментов кода из моего исследования. Надеюсь, вы мне поможете.
У меня есть две бинарные переменные: alco и smoke, которые были сгенерированы так:
smoke<- factor(with(df, ifelse((q34<2),1,0)))
alco<-factor(with(df, ifelse((q47==1), 1,0)))
df<- cbind(df, smoke, alco, educ_3, smoke_14)
Я попытался проанализировать модель с помощью пакета zeligverse.
m3<-zelig(cbind(smoke,alco) ~ fem+age+age2+smoke_14+ninc, model = "blogit", data = df)
которые приводят к ошибке
Ошибка в eval(process.binomial2.data.VGAM): ответ должен содержать только 0 и 1
Я не мог получить это, поскольку переменные в cbind являются биноминальными.
cbind(smoke,alco)
допустимо в качестве переменной ответа в соответствии с вашей моделью. тогда я думаю, чтоsmoke
илиalco
могут иметь NA. - person A. Suliman   schedule 10.06.2018