Как включить другое распределение вероятностей в TMB в R?

Question

Я использую пакет TMB в R и хочу использовать распределение плотности, которое недоступно в стиле R .

Как мне его включить? Я ничего не нашел в документации.

Например, я хочу использовать Simplex-дистрибутив, который я уже написал в R:

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 пакетов.

• Как лучше всего с этим справиться? (Я хотел бы иметь внешний файл, который содержит эту плотность, и вызывать его из файла, в который я записываю правдоподобие моего журнала)