分层贝叶斯 - 如何编码“以上都不是”选项

时间:2017-01-11 17:46:06

标签: r bayesian

我在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

0 个答案:

没有答案