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

如何解决GMNL 估计量不会包含观察到的异质性

效用部分中的“观察到的异质性”似乎在 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 

谢谢！