/* the following lines are all needed in divcntrl.dat ...
do not add or delete lines but change entry on left of each
line as you see fit ... */
Coe29.tree
12 ... number of genes ... FOLLOWING LINES CONTAIN ONLY NAMES OF DATA FILES
oest.Coe29ATP6
oest.Coe29ATP8
oest.Coe29CO1
oest.Coe29CO2
oest.Coe29CO3
oest.Coe29CYB
oest.Coe29ND1
oest.Coe29ND2
oest.Coe29ND3
oest.Coe29ND4
oest.Coe29ND4L
oest.Coe29ND5
10000 ... numsamps: How many times should the Markov chain be sampled?
100 ... sampfreq: How many cycles between samples of the Markov chain?
100000 ... burnin: How many cycles before the first sample of Markov chain?
0 ... starttype: 1 means start at not-so-great state,
0 means start at hopefully better state.
4.5 ... rttm: a priori expected number of time units between tip and root
2.25 ... rttmsd: standard deviation of prior for time between tip and root
6.0 ... rtrate: mean of prior distribution for rate at root node
3.0 ... rtratesd: standard deviation of prior for rate at root node
0.22 ... brownmean: mean of prior for brownian motion constant "nu"
0.22 ... brownsd: std. deviation of prior for brownian motion constant "nu"
/* the following lines are all needed (i.e., do not delete them) but you may
not want to alter entries unless you are familiar with the computer code */
1.0 ... minab: parameter for beta prior on proportional node depth
0.1 ... newk: parameter in Markov chain proposal step
0.5 ... othk: parameter in Markov chain proposal step
0.5 ... thek: parameter in Markov chain proposal step
5.0 ... bigtime: number higher than time units between tip and root could
be in your wildest imagination
/* the program will expect the entry below to be the number of constraints
and then the specified number of constraints should follow on
subsequent lines */
16 ... number of constraints on node times
L 28 0.9
L 32 1.46
L 33 0.5
L 34 0.57
L 35 1.46
L 37 0.5
L 39 0.9
L 40 1.52
L 41 1.12
L 46 0.89
L 47 1.41
L 48 1.55
L 49 2.23
U 52 4.5
L 54 4.11
U 54 5.28
0 ... number of tips which are not collected at time 0
0 ... nodata: 1 means approximate prior, 0 means approximate posterior
0 ...commonbrown: 1 if all genes have same tendency to change rate, 0 otherwise

//end of file

トップページ