我在R中编程使用分层bayes作为基于选择的联合任务,并想知道如何编码设计矩阵中的“无上述”选项?
我正在使用的R代码如下,其中包含来自bayesm包的rhierMnlRwMixture()函数。响应选项来自两个产品包加上没有上述选项(即总共3个响应选项)。
在设计矩阵的第三行中使用全部0是否正确,这与上述选项无关(答案1首选的例子)?
R代码:
mcmc=list(R=20000,keep=10)
library(bayesm)
out=rhierMnlRwMixture(Data=list(p=3,lgtdata=lgtdata),Prior=list(ncomp=1),Mcmc=mcmc)
如果我要查看受访者1的数据,他们会喜欢这个,回答选择3,x是设计矩阵。
lgtdata[[1]]
$y
[1] 3
$X
[,1] [,2] [,3] [,4] [,5]
[1,] 1 0 0 0 1
[2,] 0 0 1 0 0
[3,] 0 0 0 0 0