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

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

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


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	-1474.9969	2991.2148	0.3475	0.3475
GTRI	20	-1475.2389	2991.6988	0.2728	0.6204
GTR	19	-1476.4851	2992.0733	0.2262	0.8466
GTRIG	21	-1475.0094	2993.3638	0.1187	0.9653
HKYG	16	-1482.5822	2997.9506	0.0120	0.9773
HKYI	16	-1482.8138	2998.4137	0.0095	0.9868
HKY	15	-1483.9032	2998.4991	0.0091	0.9959
HKYIG	17	-1482.5999	3000.0854	0.0041	1.0000
SYMG	17	-1495.4753	3025.8363	0.0000	1.0000
SYMI	17	-1495.7193	3026.3242	0.0000	1.0000
SYM	16	-1496.8533	3026.4927	0.0000	1.0000
K2PG	13	-1500.3455	3027.2147	0.0000	1.0000
K2PI	13	-1500.5818	3027.6873	0.0000	1.0000
K2P	12	-1501.7567	3027.9616	0.0000	1.0000
SYMIG	18	-1495.4903	3027.9718	0.0000	1.0000
K2PIG	14	-1500.3586	3029.3224	0.0000	1.0000
F81	14	-1506.4459	3041.4969	0.0000	1.0000
F81G	15	-1505.4076	3041.5078	0.0000	1.0000
F81I	15	-1505.6011	3041.8949	0.0000	1.0000
F81IG	16	-1505.4076	3043.6013	0.0000	1.0000
JC69G	12	-1521.5997	3067.6477	0.0000	1.0000
JC69	11	-1522.6994	3067.7776	0.0000	1.0000
JC69I	12	-1521.7996	3068.0474	0.0000	1.0000
JC69IG	13	-1521.5998	3069.7233	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1474.9969	2989.9938	0.3519	0.3519
GTRI	20	-1475.2389	2990.4779	0.2762	0.6281
GTR	19	-1476.4851	2990.9702	0.2160	0.8441
GTRIG	21	-1475.0094	2992.0189	0.1278	0.9719
HKYG	16	-1482.5822	2997.1645	0.0098	0.9817
HKYI	16	-1482.8138	2997.6276	0.0077	0.9894
HKY	15	-1483.9032	2997.8064	0.0071	0.9965
HKYIG	17	-1482.5999	2999.1997	0.0035	1.0000
SYMG	17	-1495.4753	3024.9506	0.0000	1.0000
SYMI	17	-1495.7193	3025.4386	0.0000	1.0000
SYM	16	-1496.8533	3025.7066	0.0000	1.0000
K2PG	13	-1500.3455	3026.6910	0.0000	1.0000
SYMIG	18	-1495.4903	3026.9805	0.0000	1.0000
K2PI	13	-1500.5818	3027.1636	0.0000	1.0000
K2P	12	-1501.7567	3027.5134	0.0000	1.0000
K2PIG	14	-1500.3586	3028.7172	0.0000	1.0000
F81G	15	-1505.4076	3040.8151	0.0000	1.0000
F81	14	-1506.4459	3040.8917	0.0000	1.0000
F81I	15	-1505.6011	3041.2023	0.0000	1.0000
F81IG	16	-1505.4076	3042.8152	0.0000	1.0000
JC69G	12	-1521.5997	3067.1994	0.0000	1.0000
JC69	11	-1522.6994	3067.3988	0.0000	1.0000
JC69I	12	-1521.7996	3067.5992	0.0000	1.0000
JC69IG	13	-1521.5998	3069.1995	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKY	15	-1483.9032	3066.2643	0.7922	0.7922
HKYG	16	-1482.5822	3070.1862	0.1115	0.9037
HKYI	16	-1482.8138	3070.6492	0.0884	0.9921
HKYIG	17	-1482.5999	3076.7853	0.0041	0.9962
GTR	19	-1476.4851	3077.6835	0.0026	0.9989
GTRG	20	-1474.9969	3081.2709	0.0004	0.9993
GTRI	20	-1475.2389	3081.7550	0.0003	0.9996
K2P	12	-1501.7567	3082.2796	0.0003	0.9999
K2PG	13	-1500.3455	3086.0211	0.0000	1.0000
K2PI	13	-1500.5818	3086.4937	0.0000	1.0000
GTRIG	21	-1475.0094	3087.8598	0.0000	1.0000
K2PIG	14	-1500.3586	3092.6112	0.0000	1.0000
SYM	16	-1496.8533	3098.7282	0.0000	1.0000
SYMG	17	-1495.4753	3102.5362	0.0000	1.0000
SYMI	17	-1495.7193	3103.0241	0.0000	1.0000
F81	14	-1506.4459	3104.7857	0.0000	1.0000
SYMIG	18	-1495.4903	3109.1299	0.0000	1.0000
F81G	15	-1505.4076	3109.2730	0.0000	1.0000
F81I	15	-1505.6011	3109.6601	0.0000	1.0000
F81IG	16	-1505.4076	3115.8369	0.0000	1.0000
JC69	11	-1522.6994	3117.6013	0.0000	1.0000
JC69G	12	-1521.5997	3121.9657	0.0000	1.0000
JC69I	12	-1521.7996	3122.3654	0.0000	1.0000
JC69IG	13	-1521.5998	3128.5296	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_32429.phylip = 1 - 709;
 Partition Dummy = 1:591-tomodel-phylip/chr1_32429.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_32429.phylip = 1 - 709;
 Partition Dummy = 1:591-tomodel-phylip/chr1_32429.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 (HKY)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_32429.phylip = 1 - 709;
 Partition Dummy = 1:591-tomodel-phylip/chr1_32429.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

