我发现用于多项logit模型的mlogit-
package用于寻找估计多项式混合logit模型.在阅读了优秀的插图后,我发现我无法在任何描述的示例中应用我的数据.
我现在写信希望对我的问题有所帮助,并创建了一个简单的例子来说明我的情况.
问题如下:
在某处有辅音’Q’的字样.现在,我们进行了一项实验,他们的任务是听取这些话,并说他们是否听过Q,U或其他辅音.这必须依赖于诸如音节位置或真实/非真实单词之类的一些因素来建模.
在最小的例子中,我用音节位置创建了4个人和他们的答案.
library(mlogit)
library(nnet)
set.seed(1234)
data
decision = as.factor(sample(c("Q","U", "other"), 40, replace=TRUE)),
syllable = as.factor(sample(1:4, 40, replace=TRUE)))
summary(data)
personID decision syllable
1:11 other:10 1:18
2:10 Q :18 2: 9
3:10 U :12 3: 5
4: 9 4: 8
据我所知,nnet的multinomfunction不包括混合模型.
modNnet1
首先,我使用mlogit.data-function来重塑文件.在与同事讨论后,我们得出结论,没有替代品.特定.变量.
dataMod
mod1
data = dataMod,
reflevel="Q", rpar=c(personID="n"), panel=TRUE)
Error in names(sup.coef)
'names' attribute [1] must be the same length as the vector [0]
mod2
data = dataMod,
reflevel="Q", rpar=c(personID="n"), panel=TRUE)
Error in solve.default(H, g[!fixed]) :
Lapack routine dgesv: system is exactly singular: U[3,3] = 0
不,我不知道该怎么做,所以我在这里寻求帮助.但我相信这种问题可以通过mlogit解决,我只是看不到它;)