我的数据 (m) 的似然函数如下所示:
"dexp(m,beta) * pnorm(m,mu,sigma)".
我在我的 STAN 代码中编写的用户定义的对数似然函数如下:
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);
}
}
. 我想知道我做得对吗?我可以在用户定义的函数中使用“exponential_log”和“normal_cdf_log”等STAN函数吗?