
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
21:50 04/04/2011

Input data from file "591-tomodel-phylip/chr9_6300.phylip" (7 taxa, 961 characters)

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: HKYI


(Output sorted by AIC. w: Akaike weight, cw: cumulative w.)

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-2496.5480	5033.0960	0.6365	0.6365
GTRG	20	-2497.4220	5034.8441	0.2656	0.9020
GTRIG	21	-2497.4220	5036.8441	0.0977	0.9997
HKYI	16	-2508.9400	5049.8801	0.0001	0.9999
HKYG	16	-2509.5500	5051.0999	0.0001	1.0000
HKYIG	17	-2509.5500	5053.1001	0.0000	1.0000
GTR	19	-2521.8986	5081.7973	0.0000	1.0000
HKY	15	-2532.1832	5094.3664	0.0000	1.0000
SYMI	17	-2532.2934	5098.5868	0.0000	1.0000
SYMG	17	-2532.7825	5099.5649	0.0000	1.0000
SYMIG	18	-2532.7824	5101.5649	0.0000	1.0000
K2PI	13	-2541.4294	5108.8587	0.0000	1.0000
K2PG	13	-2541.8113	5109.6226	0.0000	1.0000
K2PIG	14	-2541.8114	5111.6227	0.0000	1.0000
SYM	16	-2553.5724	5139.1448	0.0000	1.0000
F81I	15	-2556.3047	5142.6093	0.0000	1.0000
F81G	15	-2556.5281	5143.0562	0.0000	1.0000
F81IG	16	-2556.6190	5145.2380	0.0000	1.0000
K2P	12	-2562.4619	5148.9237	0.0000	1.0000
F81	14	-2574.9330	5177.8660	0.0000	1.0000
JC69I	12	-2585.0865	5194.1731	0.0000	1.0000
JC69G	12	-2585.3125	5194.6250	0.0000	1.0000
JC69IG	13	-2585.3853	5196.7706	0.0000	1.0000
JC69	11	-2603.2455	5228.4911	0.0000	1.0000



(Output sorted by BIC. w: Akaike weight, cw: cumulative w.)

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2508.9400	5127.7677	0.5170	0.5170
HKYG	16	-2509.5500	5128.9875	0.2810	0.7980
GTRI	20	-2496.5480	5130.4555	0.1349	0.9329
GTRG	20	-2497.4220	5132.2035	0.0563	0.9891
HKYIG	17	-2509.5500	5135.8557	0.0091	0.9982
GTRIG	21	-2497.4220	5139.0715	0.0018	1.0000
HKY	15	-2532.1832	5167.3860	0.0000	1.0000
K2PI	13	-2541.4294	5172.1424	0.0000	1.0000
K2PG	13	-2541.8113	5172.9063	0.0000	1.0000
GTR	19	-2521.8986	5174.2888	0.0000	1.0000
K2PIG	14	-2541.8114	5179.7744	0.0000	1.0000
SYMI	17	-2532.2934	5181.3423	0.0000	1.0000
SYMG	17	-2532.7825	5182.3205	0.0000	1.0000
SYMIG	18	-2532.7824	5189.1884	0.0000	1.0000
K2P	12	-2562.4619	5207.3394	0.0000	1.0000
F81I	15	-2556.3047	5215.6289	0.0000	1.0000
F81G	15	-2556.5281	5216.0759	0.0000	1.0000
SYM	16	-2553.5724	5217.0324	0.0000	1.0000
F81IG	16	-2556.6190	5223.1256	0.0000	1.0000
F81	14	-2574.9330	5246.0176	0.0000	1.0000
JC69I	12	-2585.0865	5252.5888	0.0000	1.0000
JC69G	12	-2585.3125	5253.0407	0.0000	1.0000
JC69IG	13	-2585.3853	5260.0543	0.0000	1.0000
JC69	11	-2603.2455	5282.0388	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (GTRI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr9_6300.phylip = 1 - 961;
 Partition Dummy = 1:591-tomodel-phylip/chr9_6300.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (GTRI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr9_6300.phylip = 1 - 961;
 Partition Dummy = 1:591-tomodel-phylip/chr9_6300.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best BIC model (HKYI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr9_6300.phylip = 1 - 961;
 Partition Dummy = 1:591-tomodel-phylip/chr9_6300.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


-------------------------------------------------------------
End of Output

