我的数据(m)的似然函数如下所示:
"dexp(m,beta) * pnorm(m,mu,sigma)".
我在STAN代码中写的用户定义的loglikelihood函数就像:
functions{
real loglikelihood(int N,
real mu,
real sigma,
real beta,
real[] m
){
real a[N];
real b[N];
real c[N];
real final;
for(i in 1:N){
a[i]<- exponential_log(m[i],beta);
b[i]<- normal_cdf_log( m[i], mu, sigma);
c[i]<- a[i]+b[i];
}
final<- sum(c);
return(final);
}
}
。我想知道我做得对吗?我可以使用STAN函数,例如&#34; exponential_log&#34;和&#34; normal_cdf_log&#34;在用户定义的函数内部?