问题描述
我尝试从gmnl参考手册示例中估算潜在类多项式logit模型，它导致错误"Error in t.default(x) : argument is not a matrix"
其他模型（例如“ mixl”）有效
R版本4.0.2（2020-06-22）-“再次起飞”
平台：x86_64-w64-mingw32 / x64（64位）
data("Electricity",package = "mlogit")
Electr <- mlogit.data(Electricity,id.var = "id",choice = "choice",varying = 3:26,shape = "wide",sep = "")
Elec.lc <- gmnl(choice ~ pf + cl + loc + wk + tod + seas| 0 | 0 | 0 | 1,data = Electr,subset = 1:3000,model = 'lc',panel = TRUE,Q = 2)
#Estimating LC model
#Error in t.default(x) : argument is not a matrix
解决方法
我收到了相同的错误消息。
问题似乎与mlogit-package有关。
回到以前的mlogit版本（1.0-2）之后，我得以使代码运行。您可以使用以下代码安装该版本：
install.packages("https://cran.r-project.org/src/contrib/Archive/mlogit/mlogit_1.0-2.tar.gz",repos=NULL,type="source")
，然后尝试重新运行潜在类模型：
data("Electricity",package = "mlogit")
Electr <- mlogit.data(Electricity,id.var = "id",choice = "choice",varying = 3:26,shape = "wide",sep = "")
Elec.lc <- gmnl(choice ~ pf + cl + loc + wk + tod + seas | 0 | 0 | 0 | 1,data = Electr,subset = 1:3000,model = "lc",panel = TRUE,Q = 2)
summary(Elec.lc)
但是，这只是一个临时解决方法。