优势比与三级因变量

Question

我试图计算3级因变量的优势比。

通过这样做，我已经找到了如何用二元抑郁评分（测试分数的最低25％= 1，所有其他= 0）进行比值比：

ODDSR <-glm(DEPRESSION~AGE+MALE+LITERACY+COLLEGE+MARRIED, data = df, family = "binomial")

summary(ODDSR)

exp(cbind(OR = coef(ODDSR), confint(ODDSR)))

但是当我尝试将因变量变为3个级别时： - 最低25％的抑郁评分= 0 - 中间50％抑郁评分= 1 - 前25％抑郁评分= 2

glm报告错误并且不会运行。

  

eval（expr，envir，enclos）中的错误：y值必须为0＆lt; = y＆lt; = 1

甚至可能值得查看4个级别的因变量（每个四分位数编号为0,1,2,3），但我不知道是否可能。

Answer 1

安装包nnet，

library(nnet)
ODDSR <-multinom(DEPRESSION~AGE+MALE+LITERACY+COLLEGE+MARRIED, data = df)
exp(confint(ODDSR)) 
exp(coef(ODDSR))

让它发挥作用