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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKYG


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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2559.9240	5151.8479	0.4413	0.4413
HKYI	16	-2560.6739	5153.3479	0.2085	0.6498
HKYIG	17	-2559.8089	5153.6178	0.1822	0.8320
GTRG	20	-2557.5118	5155.0237	0.0902	0.9222
GTRI	20	-2558.3225	5156.6449	0.0401	0.9623
GTRIG	21	-2557.3831	5156.7661	0.0377	1.0000
HKY	15	-2572.5568	5175.1136	0.0000	1.0000
GTR	19	-2570.2923	5178.5845	0.0000	1.0000
K2PG	13	-2594.6534	5215.3068	0.0000	1.0000
K2PI	13	-2595.3205	5216.6410	0.0000	1.0000
K2PIG	14	-2594.5486	5217.0972	0.0000	1.0000
SYMG	17	-2592.8977	5219.7954	0.0000	1.0000
SYMI	17	-2593.6221	5221.2441	0.0000	1.0000
SYMIG	18	-2592.7885	5221.5769	0.0000	1.0000
F81G	15	-2598.8333	5227.6666	0.0000	1.0000
F81I	15	-2599.4967	5228.9934	0.0000	1.0000
F81IG	16	-2598.7395	5229.4791	0.0000	1.0000
K2P	12	-2607.3256	5238.6512	0.0000	1.0000
SYM	16	-2605.4768	5242.9537	0.0000	1.0000
F81	14	-2610.9106	5249.8212	0.0000	1.0000
JC69G	12	-2629.7748	5283.5496	0.0000	1.0000
JC69I	12	-2630.4262	5284.8524	0.0000	1.0000
JC69IG	13	-2629.6847	5285.3694	0.0000	1.0000
JC69	11	-2641.9870	5305.9739	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2559.9240	5231.3495	0.6636	0.6636
HKYI	16	-2560.6739	5232.8495	0.3135	0.9771
HKYIG	17	-2559.8089	5238.0883	0.0228	0.9999
HKY	15	-2572.5568	5249.6463	0.0001	1.0000
GTRG	20	-2557.5118	5254.4007	0.0000	1.0000
GTRI	20	-2558.3225	5256.0219	0.0000	1.0000
GTRIG	21	-2557.3831	5261.1120	0.0000	1.0000
GTR	19	-2570.2923	5272.9927	0.0000	1.0000
K2PG	13	-2594.6534	5279.9019	0.0000	1.0000
K2PI	13	-2595.3205	5281.2361	0.0000	1.0000
K2PIG	14	-2594.5486	5286.6611	0.0000	1.0000
K2P	12	-2607.3256	5298.2774	0.0000	1.0000
F81G	15	-2598.8333	5302.1993	0.0000	1.0000
F81I	15	-2599.4967	5303.5262	0.0000	1.0000
SYMG	17	-2592.8977	5304.2659	0.0000	1.0000
SYMI	17	-2593.6221	5305.7146	0.0000	1.0000
F81IG	16	-2598.7395	5308.9807	0.0000	1.0000
SYMIG	18	-2592.7885	5311.0162	0.0000	1.0000
F81	14	-2610.9106	5319.3851	0.0000	1.0000
SYM	16	-2605.4768	5322.4553	0.0000	1.0000
JC69G	12	-2629.7748	5343.1758	0.0000	1.0000
JC69I	12	-2630.4262	5344.4786	0.0000	1.0000
JC69IG	13	-2629.6847	5349.9644	0.0000	1.0000
JC69	11	-2641.9870	5360.6313	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr9_5250.phylip = 1 - 1063;
 Partition Dummy = 1:591-tomodel-phylip/chr9_5250.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr9_5250.phylip = 1 - 1063;
 Partition Dummy = 1:591-tomodel-phylip/chr9_5250.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr9_5250.phylip = 1 - 1063;
 Partition Dummy = 1:591-tomodel-phylip/chr9_5250.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

