Mlogit 错误：solve.default(H, g[!fixed]) 中的错误：Lapack 例程 dgesv：系统完全是奇异的：U[6,6] = 0

Question

我正在尝试对我的数据集运行多项回归，以查看渠道和接触点对选择价格和设备作为控件的影响，但不幸的是收到错误消息。

在我的 data.frame 上运行 mlogit.data 后的前几行如下所示：

ORDER_ID    PRODUCT_ID    DEVICE    PRICE    TOUCHPOINT    CHANNEL    1_or_2    CHOICE    chid    alt
123         566           laptop    99       paid          offline    1         TRUE      33      1
123         566           laptop    99       paid          offline    1         FALSE     33      2
123         534           phone     56       paid          offline    2         FALSE     45      1
123         534           phone     56       paid          offline    2         TRUE      45      2
124         876           laptop    85       unpaid        online     1         TRUE      111     1
124         876           laptop    85       unpaid        online     1         FALSE     111     2

我要运行的代码是：

Choice_mlg <- mlogit(Choice_A_or_B ~ 1 | Channel + Touchpoint + Price + Device, 
                                    data = ml_choice_1, reflevel = 1, na.action = na.exclude)

然后我收到以下错误消息：

Error in solve.default(H, g[!fixed]) : Lapack routine dgesv: system is exactly singular: U[6,6] = 0

谁能帮助我在这里做错了什么？

谢谢你，最好的

Answer 1

错误意味着Hessian矩阵是奇异的，即行列式等于0。实际上，您无法获得方差矩阵，它等于 Hessian 的负逆矩阵。

查看您的模型和数据，可能有几件事导致了它。您没有提供 MWE，所以我只能取消您提供的信息。

1. 查看您的 alt 变量，它在 1 和 2 之间变化，但两种选择是相同的。这意味着解释选择的备选方案之间绝对没有变化。这可能会导致计算上的奇异 Hessian（另一个原因是备选方案之间的相关性非常强）。
2. 您的选择变量称为Choice_A_or_B。它不是数据的一部分，所以很难看到，但在长格式中它应该取值TRUE/FALSE。它应该是TRUE 用于选择的替代品和FALSE 用于每个选择场合中的所有非选择的替代品。查看您的 CHOICE 变量，它看起来是在这里使用的正确变量。
3. 看ORDER_ID变量，两个chid是一样的。这是否意味着同一客户购买了两件商品？
4. 您指定Channel、Touchpoint、Price 和Device 是替代特定的原因吗？这对于上面的示例无关紧要（参见第 1 点），但应在最终模型中仔细考虑。

一般来说，当您设置数据（长格式）时，您希望每个备选方案有一行，其中包含个人、选择场合和备选方案的索引。