使用rstanarm的多项logit

Question

我知道您可以使用rstanarm并使用stan_glm来设置family = binomial(link="logit")包的二项式logit模型，但是rstanarm可以让您选择适合多项logit模型？

我发现了https://stats.stackexchange.com/questions/24705/can-i-use-glm-algorithms-to-do-a-multinomial-logistic-regression

但我对如何将模型Y=A+BX翻译成那里提到的表格感到困惑。

Answer 1

目前无法使用 rstanarm R软件包估算多项logit模型。有一个长期的issue来实现它，这不会太困难，但我们更关注实现多项式概率模型的更难的问题。我相信你可以在 brms R包中使用brm函数做一个多项logit模型，它也使用Stan从后验分布中绘制。

Answer 2

brms 支持多项模型。这里有两种方法：

library(foreign) 
library(brms)

#data source: https://stats.idre.ucla.edu/r/dae/multinomial-logistic-regression/
ml <- read.dta("https://stats.idre.ucla.edu/stat/data/hsbdemo.dta")
ml$prog2 <- relevel(ml$prog, ref = "academic")

#Method 1:
m1 <- brm(prog2 ~ ses + write, data = ml, family = multinomial(link = "logit"))
summary(m1)

#Method 2:
m2 <- brm(prog ~ ses + write, data = ml, family = categorical(link = "logit"))
summary(m2)