r - GMNL 估计量不会包含观察到的异质性

Question

实用程序部分中的“观察到的异质性”似乎在 gmnl 估计器中不起作用。我已经用 model = "mixl" 测试了数据，它工作得完美无缺。

GMNL 估计器能否不包含观察到的异质性？这是一个错误？有解决方法吗？

我特别指出公式“| rp-1”的第四部分给我带来了麻烦。如前所述，如果我选择“mixl”，那么包含 rp-1 就可以了。

它返回此错误消息：

Error in dUds[[j]][theRows, ] <- tcrossprod(eta[[j]][theRows, , drop = FALSE],  : 
  number of items to replace is not a multiple of replacement length

下面是代码：

rpgmnl <- gmnl(choice ~ price + l| 0 | 0 | rp-1 | 0, 
       data = dat, 
       model = "gmnl", 
       R = 30, 
       ranp = c(price = "n", l = "n"), 
       mvar= list(price= c("rp"), l = c("rp")), 
       panel= TRUE,
       correlation = TRUE, 
       method="bfgs", 
       haltons=NA)
summary(rpgmnl)

head(dat)
id j t cs choice price l l160 l320 prov can organic grass btest bfree cn rp idx
1  1 1 1  1  FALSE    17 0    0    0    0   0       1     1     1     0  0  1 1:1
2  1 1 2  1  FALSE    21 1    0    0    0   0       1     1     0     1  0  1 1:2
3  1 1 3  1  FALSE     9 0    1    0    0   0       0     0     0     0  0  1 1:3

~~~ indexes ~~~~

  chid id alt
1    1  1   1
2    1  1   2
3    1  1   3
indexes:  1, 1, 2 

谢谢！