预测拟合gmnl回归的概率

Question

我正在使用gmnl函数来拟合混合多项logit模型。 由于我对该模型的预测生育能力更感兴趣，我想通过应用predic函数之类的东西来获得它们。

 m4=gmnl(int_choice ~ 1+fico+annual_inc+int_emp_length+| time +grade+ last_fico |0, data = mldata, model="mixl",R=50,panel=TRUE,correlation = TRUE,ranp=c(annual_inc="n",int_emp_length="n"))

## how to mimic predict??
p_hat=predict(m4,type="probs")

有什么建议吗？

Answer 1

您正在寻找的是一个简单的转换规则：

• 通过计算系数的exp（）将logit（glm输出）转换为概率，这将为您提供赔率。
• 使用以下公式将赔率转换为概率：prob = odds / (1 + odds)。

可以在此处找到有关示例的非常好的解释：https://sebastiansauer.github.io/convert_logit2prob/

Answer 2

无法在插图中找到解决方案，但会在帮助文件中记录。键入help("fitted.gmnl")会产生以下结果：

fitted(object, outcome = TRUE, ...)

如果为TRUE，则拟合和残差方法返回对应于所选替代的向量，否则返回矩阵，其中每列对应于每个替代。

Answer 3

我在gmnl模型上键入str（） 我发现一个名为prob.alt的内部属性可以提供选择-残差 因此，在您的情况下，m4 $ prob.alt给出了一些有用的值（逐行查找最大值给出了预测的选择）

（在我的情况下（潜在类别mnl），这无济于事，因为它具有预测的潜在类别概率）