Using external C++ function in the model block

[bgoodri]

Well, you just call it in the model block like any other function. There was an example of that at

Problem implementing Zipf (power law) distribution -- external C++

OK. I have the problem surrounded. This works, i.e. just declare real zeta(real alpha); but don’t define a power_law function that calls it: functions{ real zeta(real s); /*TODO: implement rng for prior and posterior predictive checks*/ } data{ int<lower=0> K; // number of unique values int values[K]; int<lower=0> frequencies[K]; } parameters{ real <lower=1> alpha; } model{ real constant = log(zeta(alpha)); for (k in 1:K) { target += frequencies[k] * (-alpha * log(values…