0

我的数据 (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函数吗?

4

0 回答 0