r - 如何在 R 中的 TMB 中包含不同的概率分布？

Question

我在 R 中使用 TMB 包，我想使用R 风格分布中不可用的密度分布。

我应该如何包含它？我在文档中没有找到任何东西。

例如，我想使用我已经用 R 编写的 Simplex 分布：

dsimplex <- function(y, mu, sigma, log=TRUE){
  dis <- (y-mu)^2 / (y*(1-y)*mu^2 *(1-mu)^2)
  dsim <- -0.5*log(2*pi) -0.5*log(sigma) -
    (3/2)*log(y*(1-y)) - (1/(2*sigma))*dis
  if(log == FALSE){dsim <- exp(dsim)}
  return(dsim)
}

它在simplexreg和Stan包中可用。

• 处理它的最佳方法是什么？（我想要一个包含此密度的外部文件，并从我正在编写日志可能性的文件中调用它）