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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-3547.8977	7135.7954	0.5696	0.5696
GTRIG	21	-3547.9374	7137.8747	0.2014	0.7710
HKYG	16	-3553.5337	7139.0674	0.1109	0.8819
GTRI	20	-3550.1693	7140.3386	0.0588	0.9407
HKYIG	17	-3553.5763	7141.1527	0.0391	0.9798
HKYI	16	-3555.2363	7142.4726	0.0202	1.0000
HKY	15	-3588.0580	7206.1159	0.0000	1.0000
GTR	19	-3584.1408	7206.2815	0.0000	1.0000
F81G	15	-3659.8195	7349.6391	0.0000	1.0000
F81I	15	-3660.7994	7351.5987	0.0000	1.0000
F81IG	16	-3659.8437	7351.6874	0.0000	1.0000
F81	14	-3690.0078	7408.0156	0.0000	1.0000
K2PG	13	-3704.6829	7435.3659	0.0000	1.0000
SYMG	17	-3701.1842	7436.3685	0.0000	1.0000
K2PIG	14	-3704.7096	7437.4193	0.0000	1.0000
SYMIG	18	-3701.1983	7438.3966	0.0000	1.0000
K2PI	13	-3706.3649	7438.7298	0.0000	1.0000
SYMI	17	-3703.0440	7440.0880	0.0000	1.0000
K2P	12	-3740.0226	7504.0453	0.0000	1.0000
SYM	16	-3737.4293	7506.8585	0.0000	1.0000
JC69G	12	-3779.9118	7583.8235	0.0000	1.0000
JC69IG	13	-3779.9272	7585.8543	0.0000	1.0000
JC69I	12	-3781.1551	7586.3101	0.0000	1.0000
JC69	11	-3811.7159	7645.4318	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3553.5337	7219.2468	0.8256	0.8256
HKYI	16	-3555.2363	7222.6520	0.1504	0.9760
HKYIG	17	-3553.5763	7226.3433	0.0238	0.9998
GTRG	20	-3547.8977	7236.0197	0.0002	1.0000
GTRI	20	-3550.1693	7240.5629	0.0000	1.0000
GTRIG	21	-3547.9374	7243.1102	0.0000	1.0000
HKY	15	-3588.0580	7281.2841	0.0000	1.0000
GTR	19	-3584.1408	7301.4946	0.0000	1.0000
F81G	15	-3659.8195	7424.8073	0.0000	1.0000
F81I	15	-3660.7994	7426.7669	0.0000	1.0000
F81IG	16	-3659.8437	7431.8668	0.0000	1.0000
F81	14	-3690.0078	7478.1726	0.0000	1.0000
K2PG	13	-3704.6829	7500.5117	0.0000	1.0000
K2PI	13	-3706.3649	7503.8756	0.0000	1.0000
K2PIG	14	-3704.7096	7507.5763	0.0000	1.0000
SYMG	17	-3701.1842	7521.5591	0.0000	1.0000
SYMI	17	-3703.0440	7525.2787	0.0000	1.0000
SYMIG	18	-3701.1983	7528.5985	0.0000	1.0000
K2P	12	-3740.0226	7564.1798	0.0000	1.0000
SYM	16	-3737.4293	7587.0379	0.0000	1.0000
JC69G	12	-3779.9118	7643.9581	0.0000	1.0000
JC69I	12	-3781.1551	7646.4447	0.0000	1.0000
JC69IG	13	-3779.9272	7651.0001	0.0000	1.0000
JC69	11	-3811.7159	7700.5552	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_5492.phylip = 1 - 1109;
 Partition Dummy = 1:591-tomodel-phylip/chr1_5492.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_5492.phylip = 1 - 1109;
 Partition Dummy = 1:591-tomodel-phylip/chr1_5492.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_5492.phylip = 1 - 1109;
 Partition Dummy = 1:591-tomodel-phylip/chr1_5492.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

