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

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

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


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
HKYG	16	-1607.5762	3248.0698	0.3484	0.3484
HKYI	16	-1607.8049	3248.5272	0.2772	0.6255
HKYIG	17	-1607.5603	3250.1544	0.1229	0.7484
GTRG	20	-1604.5002	3250.4265	0.1072	0.8556
GTRI	20	-1604.6229	3250.6719	0.0948	0.9504
GTRIG	21	-1604.4777	3252.5269	0.0375	0.9879
HKY	15	-1612.2076	3255.2232	0.0097	0.9977
GTR	19	-1609.4054	3258.0989	0.0023	1.0000
K2PG	13	-1625.1376	3276.8859	0.0000	1.0000
K2PI	13	-1625.4634	3277.5376	0.0000	1.0000
K2PIG	14	-1625.1315	3278.9689	0.0000	1.0000
SYMG	17	-1623.2953	3281.6243	0.0000	1.0000
SYMI	17	-1623.5572	3282.1482	0.0000	1.0000
K2P	12	-1629.4679	3283.4584	0.0000	1.0000
SYMIG	18	-1623.2843	3283.7260	0.0000	1.0000
SYM	16	-1627.5988	3288.1150	0.0000	1.0000
F81G	15	-1634.7449	3300.2978	0.0000	1.0000
F81I	15	-1635.0281	3300.8643	0.0000	1.0000
F81IG	16	-1634.7418	3302.4010	0.0000	1.0000
F81	14	-1638.8621	3306.4301	0.0000	1.0000
JC69G	12	-1650.0339	3324.5904	0.0000	1.0000
JC69I	12	-1650.3494	3325.2213	0.0000	1.0000
JC69IG	13	-1650.0344	3326.6795	0.0000	1.0000
JC69	11	-1654.0443	3330.5301	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1607.5762	3247.1524	0.3225	0.3225
HKYI	16	-1607.8049	3247.6098	0.2565	0.5790
GTRG	20	-1604.5002	3249.0003	0.1280	0.7070
HKYIG	17	-1607.5603	3249.1206	0.1205	0.8275
GTRI	20	-1604.6229	3249.2457	0.1132	0.9407
GTRIG	21	-1604.4777	3250.9554	0.0482	0.9889
HKY	15	-1612.2076	3254.4151	0.0085	0.9974
GTR	19	-1609.4054	3256.8108	0.0026	1.0000
K2PG	13	-1625.1376	3276.2752	0.0000	1.0000
K2PI	13	-1625.4634	3276.9268	0.0000	1.0000
K2PIG	14	-1625.1315	3278.2631	0.0000	1.0000
SYMG	17	-1623.2953	3280.5905	0.0000	1.0000
SYMI	17	-1623.5572	3281.1145	0.0000	1.0000
SYMIG	18	-1623.2843	3282.5687	0.0000	1.0000
K2P	12	-1629.4679	3282.9357	0.0000	1.0000
SYM	16	-1627.5988	3287.1976	0.0000	1.0000
F81G	15	-1634.7449	3299.4897	0.0000	1.0000
F81I	15	-1635.0281	3300.0562	0.0000	1.0000
F81IG	16	-1634.7418	3301.4836	0.0000	1.0000
F81	14	-1638.8621	3305.7242	0.0000	1.0000
JC69G	12	-1650.0339	3324.0678	0.0000	1.0000
JC69I	12	-1650.3494	3324.6987	0.0000	1.0000
JC69IG	13	-1650.0344	3326.0688	0.0000	1.0000
JC69	11	-1654.0443	3330.0886	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1607.5762	3317.7678	0.4812	0.4812
HKYI	16	-1607.8049	3318.2251	0.3828	0.8640
HKY	15	-1612.2076	3320.6170	0.1158	0.9798
HKYIG	17	-1607.5603	3324.1494	0.0198	0.9996
K2PG	13	-1625.1376	3333.6501	0.0002	0.9998
K2PI	13	-1625.4634	3334.3018	0.0001	0.9999
K2P	12	-1629.4679	3335.8972	0.0001	0.9999
GTRG	20	-1604.5002	3337.2695	0.0000	1.0000
GTRI	20	-1604.6229	3337.5149	0.0000	1.0000
K2PIG	14	-1625.1315	3340.0515	0.0000	1.0000
GTR	19	-1609.4054	3340.6665	0.0000	1.0000
GTRIG	21	-1604.4777	3343.6381	0.0000	1.0000
SYMG	17	-1623.2953	3355.6193	0.0000	1.0000
SYMI	17	-1623.5572	3356.1433	0.0000	1.0000
SYM	16	-1627.5988	3357.8129	0.0000	1.0000
SYMIG	18	-1623.2843	3362.0109	0.0000	1.0000
F81G	15	-1634.7449	3365.6916	0.0000	1.0000
F81I	15	-1635.0281	3366.2581	0.0000	1.0000
F81	14	-1638.8621	3367.5127	0.0000	1.0000
F81IG	16	-1634.7418	3372.0990	0.0000	1.0000
JC69G	12	-1650.0339	3377.0293	0.0000	1.0000
JC69I	12	-1650.3494	3377.6602	0.0000	1.0000
JC69	11	-1654.0443	3378.6367	0.0000	1.0000
JC69IG	13	-1650.0344	3383.4437	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/chr1_32310.phylip = 1 - 610;
 Partition Dummy = 1:591-tomodel-phylip/chr1_32310.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/chr1_32310.phylip = 1 - 610;
 Partition Dummy = 1:591-tomodel-phylip/chr1_32310.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/chr1_32310.phylip = 1 - 610;
 Partition Dummy = 1:591-tomodel-phylip/chr1_32310.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

