r - 使用 rstanarm 的多项式 logit

Question

我知道您可以rstanarm通过使用stan_glm和设置将二项式 logit 模型与包拟合family = binomial(link="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

多项式 logit 模型目前无法使用rstanarm R 包进行估计。实现它是一个长期存在的问题，这不会太困难，但我们更关注实现多项概率模型的更困难的问题。我相信你可以用brmsbrm R 包中的函数做一个多项式 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)