
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:18 04/04/2011

Input data from file "683-tomodel-phylip/chr7_4715.phylip" (7 taxa, 440 characters)

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


Note: number of parameters compared to sample size
is low (Nchar/Nparams < 40). Consider the use of AICc.


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

Model	df	lnL		AICc		wAICc	cwAICc
GTRG	20	-1596.9022	3235.8093	0.6586	0.6586
GTRIG	21	-1596.9024	3238.0153	0.2186	0.8772
GTRI	20	-1598.6299	3239.2646	0.1170	0.9942
GTR	19	-1603.0922	3245.9939	0.0040	0.9983
SYMG	17	-1606.5935	3248.6372	0.0011	0.9994
SYMIG	18	-1606.5936	3250.8119	0.0004	0.9997
SYMI	17	-1608.1052	3251.6605	0.0002	1.0000
HKYG	16	-1611.9964	3257.2789	0.0000	1.0000
SYM	16	-1612.6664	3258.6189	0.0000	1.0000
HKYIG	17	-1611.9966	3259.4434	0.0000	1.0000
HKYI	16	-1613.6007	3260.4874	0.0000	1.0000
HKY	15	-1618.0568	3267.2457	0.0000	1.0000
K2PG	13	-1622.4753	3271.8050	0.0000	1.0000
K2PIG	14	-1622.4754	3273.9390	0.0000	1.0000
K2PI	13	-1623.8491	3274.5526	0.0000	1.0000
K2P	12	-1628.7322	3282.1950	0.0000	1.0000
F81G	15	-1662.6913	3356.5146	0.0000	1.0000
F81IG	16	-1662.6914	3358.6688	0.0000	1.0000
F81I	15	-1663.8710	3358.8740	0.0000	1.0000
F81	14	-1668.5406	3366.0695	0.0000	1.0000
JC69G	12	-1676.4486	3377.6278	0.0000	1.0000
JC69I	12	-1677.5002	3379.7310	0.0000	1.0000
JC69IG	13	-1676.4486	3379.7517	0.0000	1.0000
JC69	11	-1681.7021	3386.0210	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1596.9022	3233.8045	0.6439	0.6439
GTRIG	21	-1596.9024	3235.8047	0.2368	0.8807
GTRI	20	-1598.6299	3237.2598	0.1144	0.9951
GTR	19	-1603.0922	3244.1844	0.0036	0.9987
SYMG	17	-1606.5935	3247.1869	0.0008	0.9995
SYMIG	18	-1606.5936	3249.1872	0.0003	0.9998
SYMI	17	-1608.1052	3250.2103	0.0002	1.0000
HKYG	16	-1611.9964	3255.9928	0.0000	1.0000
SYM	16	-1612.6664	3257.3328	0.0000	1.0000
HKYIG	17	-1611.9966	3257.9932	0.0000	1.0000
HKYI	16	-1613.6007	3259.2014	0.0000	1.0000
HKY	15	-1618.0568	3266.1136	0.0000	1.0000
K2PG	13	-1622.4753	3270.9505	0.0000	1.0000
K2PIG	14	-1622.4754	3272.9508	0.0000	1.0000
K2PI	13	-1623.8491	3273.6982	0.0000	1.0000
K2P	12	-1628.7322	3281.4643	0.0000	1.0000
F81G	15	-1662.6913	3355.3826	0.0000	1.0000
F81IG	16	-1662.6914	3357.3828	0.0000	1.0000
F81I	15	-1663.8710	3357.7420	0.0000	1.0000
F81	14	-1668.5406	3365.0812	0.0000	1.0000
JC69G	12	-1676.4486	3376.8971	0.0000	1.0000
JC69IG	13	-1676.4486	3378.8973	0.0000	1.0000
JC69I	12	-1677.5002	3379.0003	0.0000	1.0000
JC69	11	-1681.7021	3385.4041	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	20	-1596.9022	3315.5400	0.4743	0.4743
SYMG	17	-1606.5935	3316.6621	0.2706	0.7450
GTRI	20	-1598.6299	3318.9953	0.0843	0.8292
SYMI	17	-1608.1052	3319.6855	0.0597	0.8889
HKYG	16	-1611.9964	3321.3812	0.0256	0.9145
GTRIG	21	-1596.9024	3321.6270	0.0226	0.9371
GTR	19	-1603.0922	3321.8331	0.0204	0.9575
SYM	16	-1612.6664	3322.7212	0.0131	0.9706
SYMIG	18	-1606.5936	3322.7492	0.0129	0.9835
K2PG	13	-1622.4753	3324.0786	0.0066	0.9901
HKYI	16	-1613.6007	3324.5898	0.0051	0.9953
K2PI	13	-1623.8491	3326.8263	0.0017	0.9969
HKY	15	-1618.0568	3327.4152	0.0013	0.9982
HKYIG	17	-1611.9966	3327.4683	0.0012	0.9994
K2PIG	14	-1622.4754	3330.1656	0.0003	0.9997
K2P	12	-1628.7322	3330.5056	0.0003	1.0000
F81G	15	-1662.6913	3416.6842	0.0000	1.0000
F81I	15	-1663.8710	3419.0436	0.0000	1.0000
F81	14	-1668.5406	3422.2961	0.0000	1.0000
F81IG	16	-1662.6914	3422.7712	0.0000	1.0000
JC69G	12	-1676.4486	3425.9384	0.0000	1.0000
JC69I	12	-1677.5002	3428.0416	0.0000	1.0000
JC69	11	-1681.7021	3430.3587	0.0000	1.0000
JC69IG	13	-1676.4486	3432.0254	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 683-tomodel-phylip/chr7_4715.phylip = 1 - 440;
 Partition Dummy = 1:683-tomodel-phylip/chr7_4715.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 683-tomodel-phylip/chr7_4715.phylip = 1 - 440;
 Partition Dummy = 1:683-tomodel-phylip/chr7_4715.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 (GTRG)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr7_4715.phylip = 1 - 440;
 Partition Dummy = 1:683-tomodel-phylip/chr7_4715.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;


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

