Legofit
infers population history from nucleotide site patterns.
|
Calculate confidence interval from a flat file.
usage: bootci.py [options] <foo.flat> where <foo.flat> is a file with parameters in columns and data sets in rows. Such files are produced by flatfile.py and booma. The first row should contain a column header listing parameter names. The second should contain parameter estimates for the real data. Each succeeding row should contain a parameter estimate for a bootstrap replicate. Options may include: -c <x> or --conf <x> Set confidence level. -h or --help Print this message. -l <x> or --label <x> Set label. The program writes to standard output. # bootci.py run at: 2017-01-06 22:32:58.916983 # real data: lancre.legofit # bootstrap replicates: boot0.legofit boot1.legofit boot10.legofit ... # confidence: 0.950 par est low high Tnd 25882.18097500 25783.98581875 25918.98318675 Txy 1910.32458600 1898.67020375 2020.23384350 mD 0.00030900 0.00001650 0.00133550 mN 0.01769600 0.01396675 0.01715075 twoNnd 1204.62497100 1168.42503700 1327.14358450 twoNxy 59669.81112500 59120.53155525 61060.16095975 twoNxynd 34355.15657100 34093.55639600 34552.91827875
Here, "par" lists the parameter names, "est" gives the estimates, and "low" and "high" give the lower and upper confidence bounds.
By default, the confidence level is 95%, but this can be changed with the –conf command-line option.
The –label <x> option will add a column in which every row contains the character string <x>. This is useful in data analysis, if one concatenates the output from several runs of bootci.py. The column of labels makes it possible to distinguish the estimates from different runs of bootci.py.
For example,
bootci.py --label lancre lancre.legofit boot*.legofit
would generate
# bootci.py run at: 2017-01-06 22:32:58.916983 # real data: lancre.legofit # bootstrap replicates: boot0.legofit boot1.legofit boot10.legofit ... # confidence: 0.950 par est low high lbl Tnd 25882.18097500 25783.98581875 25918.98318675 lancre Txy 1910.32458600 1898.67020375 2020.23384350 lancre mD 0.00030900 0.00001650 0.00133550 lancre mN 0.01769600 0.01396675 0.01715075 lancre twoNnd 1204.62497100 1168.42503700 1327.14358450 lancre twoNxy 59669.81112500 59120.53155525 61060.16095975 lancre twoNxynd 34355.15657100 34093.55639600 34552.91827875 lancre