Ошибка в двумерной логистической модели в R

Question

У меня неожиданная ошибка в моем исследовании.Позвольте мне показать вам несколько фрагментов кода из моего исследования.Надеюсь, вы мне поможете.

У меня есть две двоичные переменные: 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

Iне удалось получить его, поскольку переменные в cbind являются биноминальными.

Answer 1

Было бы хорошо узнать, что вы пытаетесь подогнать и как выглядят ваши данные. Здесь есть несколько рекомендаций, чтобы задать хорошие вопросы.

Я предполагаю, что вы пытаетесь запустить (двоичный) логит.Если да, то glm() может оценить такую ​​модель

Например:

df = data.frame(x = factor(sample(0:1, 25, replace = TRUE)), 
                y = factor(sample(1:4, 25, replace = TRUE)), 
                z = sample(18:65, 25, replace = TRUE))


summary(glm(x ~ y + z, family = binomial(link = "logit"), 
            data = df))

Если у вас есть более двух категорий в вашей переменной результата, и они упорядочены.clm() из ordinal пакет может быть опцией:

 library(ordinal)

 summary(clm(y ~ x + z, data = df, link = "logit")

Надеюсь, это поможет