我正在尝试使用 r 函数 mlogit 预测选择实验的选项。我有替代的特定变量和单独的特定变量。现在，我想介绍一个替代特定变量和一个单独的特定变量之间的交互作用项。因此，我的问题是：这可能吗？如果是，必须如何指定公式？如果没有，是否有另一个软件包支持此功能？
在介绍mlogit(http://www2.uaem.mx/r-mirror/web/packages/mlogit/vignettes/mlogit.pdf(的论文中，替代变量和单个特定变量的使用指定如下(x1是替代特定变量，x2是单个特定变量(：
f1 <- mFormula(y ~ x1 | x2)
所以我尝试了几种变体，但没有一个产生预期的结果：
f1 <- mFormula(y ~ x1 | x1*x2 | x2)
f1 <- mFormula(y ~ x1*x2 | x2)
f1 <- mFormula(y ~ x1 | x1*x2)
或由于奇点，我总是收到错误(倒数条件数不同，但错误保持不变(：
Error in solve.default(H, g[!fixed]) :
system is computationally singular: reciprocal condition number = 3.99539e-21
用替代和单个交互变量的相互作用来估计系统在数学上是不可能的，因为系统变得奇异。直观地说，通过引入这种交互，你想说某些人更有可能购买某些替代品 - 但这种情况已经被个人特定的变量"覆盖"。
见论文第7页：
必须考虑三种变量：
- 具有通用系数 beta x_ij的替代特定变量，
- 单个特定变量z_i具有替代特定系数gamma_j，
具有替代特定- 系数delta_j w_ij的替代特定变量
我想情况 2 是你想要的。