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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: HKYI


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
GTRI	20	-1411.8524	2865.3079	0.4680	0.4680
GTRG	20	-1412.0268	2865.6567	0.3931	0.8610
GTRIG	21	-1412.0403	2867.8474	0.1315	0.9925
GTR	19	-1418.1244	2875.6963	0.0026	0.9951
HKYI	16	-1421.5074	2876.0451	0.0022	0.9973
HKYG	16	-1421.5734	2876.1772	0.0020	0.9993
HKYIG	17	-1421.5868	2878.3349	0.0007	1.0000
HKY	15	-1427.6790	2886.2655	0.0000	1.0000
SYMI	17	-1432.8440	2900.8494	0.0000	1.0000
SYMG	17	-1433.0539	2901.2690	0.0000	1.0000
SYMIG	18	-1433.0675	2903.4353	0.0000	1.0000
K2PI	13	-1441.3242	2909.3339	0.0000	1.0000
K2PG	13	-1441.3763	2909.4381	0.0000	1.0000
K2PIG	14	-1441.3866	2911.5656	0.0000	1.0000
SYM	16	-1439.4932	2912.0168	0.0000	1.0000
F81I	15	-1442.2943	2915.4961	0.0000	1.0000
F81G	15	-1442.4190	2915.7455	0.0000	1.0000
F81IG	16	-1442.4253	2917.8809	0.0000	1.0000
K2P	12	-1447.7552	2920.0969	0.0000	1.0000
F81	14	-1448.0447	2924.8818	0.0000	1.0000
JC69I	12	-1459.9877	2944.5618	0.0000	1.0000
JC69G	12	-1460.1118	2944.8101	0.0000	1.0000
JC69IG	13	-1460.1178	2946.9211	0.0000	1.0000
JC69	11	-1465.9126	2954.3205	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-1411.8524	2863.7048	0.4634	0.4634
GTRG	20	-1412.0268	2864.0537	0.3892	0.8526
GTRIG	21	-1412.0403	2866.0806	0.1413	0.9939
GTR	19	-1418.1244	2874.2487	0.0024	0.9963
HKYI	16	-1421.5074	2875.0148	0.0016	0.9979
HKYG	16	-1421.5734	2875.1469	0.0015	0.9994
HKYIG	17	-1421.5868	2877.1736	0.0006	1.0000
HKY	15	-1427.6790	2885.3581	0.0000	1.0000
SYMI	17	-1432.8440	2899.6881	0.0000	1.0000
SYMG	17	-1433.0539	2900.1077	0.0000	1.0000
SYMIG	18	-1433.0675	2902.1349	0.0000	1.0000
K2PI	13	-1441.3242	2908.6484	0.0000	1.0000
K2PG	13	-1441.3763	2908.7526	0.0000	1.0000
K2PIG	14	-1441.3866	2910.7732	0.0000	1.0000
SYM	16	-1439.4932	2910.9865	0.0000	1.0000
F81I	15	-1442.2943	2914.5887	0.0000	1.0000
F81G	15	-1442.4190	2914.8381	0.0000	1.0000
F81IG	16	-1442.4253	2916.8506	0.0000	1.0000
K2P	12	-1447.7552	2919.5105	0.0000	1.0000
F81	14	-1448.0447	2924.0893	0.0000	1.0000
JC69I	12	-1459.9877	2943.9754	0.0000	1.0000
JC69G	12	-1460.1118	2944.2236	0.0000	1.0000
JC69IG	13	-1460.1178	2946.2356	0.0000	1.0000
JC69	11	-1465.9126	2953.8252	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1421.5074	2943.8274	0.4705	0.4705
HKYG	16	-1421.5734	2943.9594	0.4405	0.9110
GTRI	20	-1411.8524	2949.7205	0.0247	0.9357
HKY	15	-1427.6790	2949.8699	0.0229	0.9586
GTRG	20	-1412.0268	2950.0694	0.0208	0.9794
HKYIG	17	-1421.5868	2950.2870	0.0186	0.9980
GTR	19	-1418.1244	2955.9637	0.0011	0.9991
GTRIG	21	-1412.0403	2956.3971	0.0009	1.0000
K2PI	13	-1441.3242	2964.5586	0.0000	1.0000
K2PG	13	-1441.3763	2964.6629	0.0000	1.0000
K2PIG	14	-1441.3866	2970.9842	0.0000	1.0000
K2P	12	-1447.7552	2971.1199	0.0000	1.0000
SYMI	17	-1432.8440	2972.8014	0.0000	1.0000
SYMG	17	-1433.0539	2973.2211	0.0000	1.0000
F81I	15	-1442.2943	2979.1005	0.0000	1.0000
F81G	15	-1442.4190	2979.3499	0.0000	1.0000
SYMIG	18	-1433.0675	2979.5490	0.0000	1.0000
SYM	16	-1439.4932	2979.7990	0.0000	1.0000
F81	14	-1448.0447	2984.3003	0.0000	1.0000
F81IG	16	-1442.4253	2985.6631	0.0000	1.0000
JC69I	12	-1459.9877	2995.5848	0.0000	1.0000
JC69G	12	-1460.1118	2995.8330	0.0000	1.0000
JC69	11	-1465.9126	3001.1339	0.0000	1.0000
JC69IG	13	-1460.1178	3002.1459	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_23548.phylip = 1 - 545;
 Partition Dummy = 1:591-tomodel-phylip/chr2_23548.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_23548.phylip = 1 - 545;
 Partition Dummy = 1:591-tomodel-phylip/chr2_23548.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_23548.phylip = 1 - 545;
 Partition Dummy = 1:591-tomodel-phylip/chr2_23548.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

