int main_IHM(int *arga,int *argb,int *argc,int *argd,double *OUT,
char **HEADER,int *QFAIA,double *QFADyA,int *QFADNoORFA,int *QFAIB,
double *QFADyB,int *QFADNoORFB,double *PRIORS,double *TUNING)
{
GetRNGstate();
struct_data_IHM *data= malloc(sizeof(struct_data_IHM));
struct_para_IHM *para= malloc(sizeof(struct_para_IHM));
struct_priors_IHM *priors= malloc(sizeof(struct_priors_IHM));
struct_tuning_IHM *tuning = malloc(sizeof(struct_tuning_IHM));
struct_adaptive_IHM *adaptive = malloc(sizeof(struct_tuning_IHM));
int burn,iters,thin,adaptive_phase;
burn=*arga; /*burn in*/
iters=*argb; /*iterations*/
thin=*argc; /*thinning*/
adaptive_phase=*argd; /*adaptive phase length within burn in*/
inzstruct_data_IHM(data,QFAIA,QFADyA,QFADNoORFA,QFAIB,QFADyB,QFADNoORFB);
inzstruct_tuning_IHM(tuning,TUNING);
inzstruct_adaptive_IHM(adaptive);
inzstruct_priors_IHM(priors,PRIORS);
inzstruct_para_IHM(para,data,priors);
gibbsandMHloop_IHM(burn,1,data,para,priors,tuning,adaptive,0,adaptive_phase,OUT,HEADER);
gibbsandMHloop_IHM(iters,thin,data,para,priors,tuning,adaptive,1,0,OUT,HEADER);
PutRNGstate();
return 0;
}
int
main_IHM(int *arga,int *argb,int *argc,double *OUT,char **HEADER,int *QFAIA,double *QFADyA,int *QFADNoORFA,int *QFAIB,double *QFADyB,int *QFADNoORFB,double *PRIORS)/****/
{
GetRNGstate();
struct_data_IHM *data= malloc(sizeof(struct_data_IHM));
struct_para_IHM *para= malloc(sizeof(struct_para_IHM));
struct_priors_IHM *priors= malloc(sizeof(struct_priors_IHM));
struct_MH_IHM *MH = malloc(sizeof(struct_MH_IHM));
int burn,iters,thin;
burn=*arga; /*Burn in*/
iters=*argb; /*iterations*/
thin=*argc; /*thining*/
inzstruct_data_IHM(data,QFAIA,QFADyA,QFADNoORFA,QFAIB,QFADyB,QFADNoORFB);
inzstruct_priors_IHM(priors,PRIORS);
inzstruct_para_IHM(para,data,priors);
inzstruct_MH_IHM(MH);
gibbsandMHloop_IHM(burn,1,data,para,priors,MH,0,OUT,HEADER);
gibbsandMHloop_IHM(iters,thin,data,para,priors,MH,1,OUT,HEADER);
PutRNGstate();
return 0;
}