如何在 rjags 中为非负因子加载指定尖峰和平板先验

编程问答 2022-05-18

问题描述

我正在尝试使用适用于非负因子载荷 (https://jmlr.org/papers/volume15/lan14a/lan14a.pdf,eq 10) 的尖峰板模型，这是我的代码：

for (k in 1:K){
  for (m in 1:M){
  psi_raw[m,k] ~ dsum(a[k]*dexp(lambda[k]) +(1-a[k])*dnorm(0,10000))
}
}

但是，它给了我错误消息，如果我指定它会起作用：

psi_raw[m,k] = dsum(a[k]*dexp(lambda[k]) +(1-a[k])*dnorm(0,10000))

经过一番搜索，似乎错误是我无法为 psi_raw[m,k] 指定两个随机部分的总和，但是如果我指定 = 而不是 ~，它将是确定性的，并且这不是我想要的。我真的很感激任何想法或想法！提前致谢！

（我也在这里问了同样的问题：https://sourceforge.net/p/mcmc-jags/discussion/610037/thread/2494a8f6fc/）

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）

rjags