多项/条件 Logit 回归,为什么 StatsModel 在 mlogit 包示例上失败?

2023-12-29

我正在尝试重现 R 中 mlogit 包的多项 logit 回归的示例。

data("Fishing", package = "mlogit")
Fish <- mlogit.data(Fishing, varying = c(2:9), shape = "wide", choice = "mode")
#a pure "conditional" model
summary(mlogit(mode ~ price + catch, data = Fish))

为了使用 statsmodel 函数 MNLogit 重现此示例,我将 Fishing 数据集导出为 csv 文件并执行以下操作

import pandas
import statsmodels.api as st
#load data
df = pandas.read_csv("Fishing.csv")
x = df.drop('mode', axis = 1)
y = df['mode']
mdl = st.MNLogit(y, x)
mdl_fit = mdl.fit()  

我收到以下错误

LinAlgError: Singular matrix

我试图弄清楚如何重新组织原始数据集 Fishing,因为我知道 mlogit 包在拟合之前重新组织数据,但无法弄清楚如何在 statsmodel 中更改它。任何帮助将非常感激。


statsmodels 中的 MNLogit 实现了多项 Logit 的不同版本。 AFAICS,对应R中的nnet multinomhttps://stats.stackexchange.com/questions/186344/r-interpreting-the-multinom-output-using-the-iris-dataset/188426 https://stats.stackexchange.com/questions/186344/r-interpreting-the-multinom-output-using-the-iris-dataset/188426

在这种情况下,参数因选择而异,但解释变量不变。在多选 CLogit 版本或 R 中的 mlogit 版本中,解释变量因选择而异,但参数与选择无关。

CLogit 和其他多项式 logit 版本正在等待 statsmodels 的拉取请求,目前在主分支中不可用。

本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)

多项/条件 Logit 回归,为什么 StatsModel 在 mlogit 包示例上失败? 的相关文章

随机推荐