double thp1 = 1.0, thp2 = 5.0; /* params for theta */
double thxp1 = 1.0, thxp2 = 10.0; /* params for theta X  */
double thxp0 = 40; /* cross term */
/* -1 is default value (=0 logically)  */

double lp1 = 10.2, lp2 = 2;
double lxp1 = 10.2, lxp2 = 2;

double priorlmean = 6;
double priorlmsig = 5;
/* hyperprior on mean, s.dev for gamma prior on lambda */

double qtrbase = 0.0;

double loclip = -20.0;
double hiclip = 15.0;

double a1 = 2, b1 = 8;
double aa2 = 2, bb2 = 14, cc2 = 85;
double p1 = 18, psi1 = 3;
/* for toys */

double muval = 0.0, tmumean = 0.2, muval1;
/* for bridge sampler */

int pubxindiv = -1;
int alkesmode = NO;
int malexhet = NO;
int familynames = YES;

int decim = 0, dmindis = 200000, dmaxdis = 500000; // decimation parameters
int hashcheck = YES;
int outputall = NO;
int sevencolumnped = NO;

FILE *fstdetails = NULL;

