我发现用于多项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解决,我只是看不到它;)