在用户定义的函数内使用STAN函数

时间:2016-11-23 22:04:55

标签: r bayesian stan rstan

我的数据(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;在用户定义的函数内部?

0 个答案:

没有答案