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

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

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


r regression mlogit
person Galina Pazhitnova    schedule 09.06.2018    source источник
comment
Ошибка говорит вам, что две переменные ответа должны быть равны нулю и единицам, поэтому попробуйте не преобразовывать их в множители, а оставить их как 0/1.   -  person user20650    schedule 09.06.2018
comment
если выполнение cbind(smoke,alco) допустимо в качестве переменной ответа в соответствии с вашей моделью. тогда я думаю, что smoke или alco могут иметь NA.   -  person A. Suliman    schedule 10.06.2018

Ответы (1)


arrow_upward
0
arrow_downward

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

Я предполагаю, что вы пытаетесь запустить (двоичный) логит. Если да, то 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")

Я надеюсь, что это помогает

person Edu    schedule 09.06.2018
comment
Эду; похоже, это модель, которой соответствует оператор docs.zeligproject.org/articles/zeligchoice_blogit.html< /а> - person user20650; 09.06.2018