在R中运行mlogit函数后生成预测概率的置信区间

Question

我一直在努力解决以下问题，并且非常感谢任何帮助。我使用mlogit函数在R中运行logit模型，并且能够生成为预测变量的给定值选择每个备选项的预测概率，如下所示：

library(mlogit)
data("Fishing", package = "mlogit")
Fish <- mlogit.data(Fishing, varying = c(2:9), shape = "wide", choice = "mode")

Fish_fit<-Fish[-(1:4),]
Fish_test<-Fish[1:4,]
m <- mlogit(mode ~price+ catch | income, data = Fish_fit)
predict(m,newdata=Fish_test,)

但是，我不能弄清楚如何在预测概率估计中增加置信区间。我已经尝试为predict函数添加参数，但似乎都没有生成它们。关于如何实现它的任何想法将不胜感激。

Answer 1

这里的一种方法是蒙特卡罗模拟。您可以模拟多变量正态采样分布的重复绘制，其参数由模型结果给出。

对于每个模拟，估计您的预测概率，并使用他们的经验分布而不是模拟来获得您的置信区间。

"