用rstanarm的多项式logit

Question

我知道，您可以通过使用rstanarm包和设置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

多项logit模型目前无法用rstanarm R包估计。有一个长期的问题来实现它，这将不会太困难，但我们更关注的是更困难的问题，得到一个多项概率模型的实现。我相信你可以用brm函数在brms 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)

Answer 3

您可以使用stan_polr，参见"回归和其他故事“作为示例(第15.5节)。例如。

library(rstanarm)
fit <- stan_polr(factor(vote) ~ value, data = data_401, prior = R2(0.3, "mean))

取自https://avehtari.github.io/ROS-Examples/Storable/storable.html