Displaying 1 result from an estimated 1 matches for "app1_3".
Did you mean:
app1_1
2013 Aug 27
1
Error in simulation. NAN
...with a small
script and then i post the original one.
The problem is after the line
for(MCMC_iter2=0;MCMC_iter2<thin;MCMC_iter2++)
For the first 34 iterations all work fine, after, all the simulations
of mu_acc_P return an "nan". If i delete the line
alpha_acc_P[0] = rnorm(app1_3,sqrt(app1_4) );
the simulations of mu_acc_P work fine. For me it is really strange,
Thanks!.
#include <iostream>
#include <string>
using namespace std;
#include <R.h>
#include <Rinternals.h>
#include <R_ext/Linpack.h>
#include <R_ext/Lapack.h>
#include <...