winbugs - 用 BUGS/JAGS/STAN 语言定义新函数和发行版

Question

我是统计分析领域的新手，最近对 BUGS/JAGS/STAN 建模语言产生了兴趣。真正令我惊讶的是，我没有看到任何定义新函数或发行版以避免代码重复的示例。例如，假设我经常使用泊松分布的平方，是否可以执行以下操作？

dsqpo <- function(lambda) {
  tmp ~ dpois(lambda)
  tmp2 <- tmp * tmp
  return(tmp2)
}

然后稍后

model{
  ...
  x ~ dsqpo(alpha)
  y ~ dsqpo(beta)
}

无需每次都定义新的临时变量。

Answer 1

对于 Stan，功能将在下一版本中提供。当前版本 v2.2.0 没有将用户定义函数作为语言的一部分。

有关建议的语法，请参阅:https://github.com/stan-dev/stan/wiki/Function-Syntax-and-Semantics-Design

有关 Stan 的其他帮助，请查看 stan-users google 群组:https://groups.google.com/forum/#!forum/stan-users