1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
|
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;
|
