Пример #1
0
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;
}
Пример #2
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;
}