是否有可能在mlogit模型中引入替代特定变量和单个特定变量之间的交互项

Question

我正在尝试使用r函数的口语来预测选择实验的选项。我有替代的特定变量和各个特定变量。现在，我想介绍一个替代特定变量之一与单个特定变量之一之间的交互项。因此，我的问题是：这可能吗？如果是，必须如何指定公式？如果没有，是否还有另一个支持此功能的软件包？

在介绍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

Answer 1

在数学上不可能用替代变量和各个相互作用变量的相互作用来估计系统，因为系统变得单一。直观地讲，通过引入这种交互作用，您想说某些人更有可能购买某些替代品，但是这种情况已经被特定于个人的变量“覆盖”了。

请参阅论文的第7页：

  

必须考虑三种变量：

     

  
• 具有一般系数beta的替代特定变量x_ij
  
• 具有其他特定系数gamma_j的各个特定变量z_i，
  
• 具有其他特定系数delta_j的替代特定变量w_ij
  

我想情况2是您想要的。