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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	20	-1324.2172	2690.6163	0.4684	0.4684
GTRI	20	-1325.2721	2692.7260	0.1631	0.6315
GTRIG	21	-1324.2153	2692.8368	0.1543	0.7859
HKYG	16	-1329.9807	2693.3598	0.1188	0.9047
HKYI	16	-1330.7602	2694.9189	0.0545	0.9592
HKYIG	17	-1329.9685	2695.5142	0.0405	0.9996
GTR	19	-1332.7319	2705.4327	0.0003	0.9999
HKY	15	-1338.4262	2708.0832	0.0001	1.0000
K2PG	13	-1355.1597	2737.2480	0.0000	1.0000
K2PI	13	-1355.5917	2738.1121	0.0000	1.0000
K2PIG	14	-1355.1145	2739.3031	0.0000	1.0000
SYMG	17	-1352.3945	2740.3664	0.0000	1.0000
SYMI	17	-1352.8637	2741.3047	0.0000	1.0000
SYMIG	18	-1352.3609	2742.4892	0.0000	1.0000
K2P	12	-1364.1097	2753.0133	0.0000	1.0000
SYM	16	-1361.4579	2756.3143	0.0000	1.0000
F81G	15	-1368.3792	2767.9892	0.0000	1.0000
F81I	15	-1368.7292	2768.6891	0.0000	1.0000
F81IG	16	-1368.3579	2770.1143	0.0000	1.0000
F81	14	-1375.1931	2779.4604	0.0000	1.0000
JC69G	12	-1389.7843	2804.3625	0.0000	1.0000
JC69I	12	-1390.0110	2804.8160	0.0000	1.0000
JC69IG	13	-1389.7349	2806.3984	0.0000	1.0000
JC69	11	-1396.9542	2816.5785	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1324.2172	2688.4344	0.4923	0.4923
GTRIG	21	-1324.2153	2690.4305	0.1814	0.6737
GTRI	20	-1325.2721	2690.5442	0.1714	0.8451
HKYG	16	-1329.9807	2691.9614	0.0844	0.9295
HKYI	16	-1330.7602	2693.5204	0.0387	0.9683
HKYIG	17	-1329.9685	2693.9369	0.0314	0.9997
GTR	19	-1332.7319	2703.4638	0.0003	1.0000
HKY	15	-1338.4262	2706.8524	0.0000	1.0000
K2PG	13	-1355.1597	2736.3195	0.0000	1.0000
K2PI	13	-1355.5917	2737.1835	0.0000	1.0000
K2PIG	14	-1355.1145	2738.2290	0.0000	1.0000
SYMG	17	-1352.3945	2738.7891	0.0000	1.0000
SYMI	17	-1352.8637	2739.7273	0.0000	1.0000
SYMIG	18	-1352.3609	2740.7218	0.0000	1.0000
K2P	12	-1364.1097	2752.2194	0.0000	1.0000
SYM	16	-1361.4579	2754.9158	0.0000	1.0000
F81G	15	-1368.3792	2766.7584	0.0000	1.0000
F81I	15	-1368.7292	2767.4583	0.0000	1.0000
F81IG	16	-1368.3579	2768.7159	0.0000	1.0000
F81	14	-1375.1931	2778.3863	0.0000	1.0000
JC69G	12	-1389.7843	2803.5686	0.0000	1.0000
JC69I	12	-1390.0110	2804.0221	0.0000	1.0000
JC69IG	13	-1389.7349	2805.4699	0.0000	1.0000
JC69	11	-1396.9542	2815.9084	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1329.9807	2756.0630	0.6597	0.6597
HKYI	16	-1330.7602	2757.6221	0.3025	0.9622
HKYIG	17	-1329.9685	2762.0449	0.0331	0.9954
HKY	15	-1338.4262	2766.9477	0.0029	0.9982
GTRG	20	-1324.2172	2768.5615	0.0013	0.9995
GTRI	20	-1325.2721	2770.6713	0.0004	0.9999
GTRIG	21	-1324.2153	2774.5640	0.0001	1.0000
GTR	19	-1332.7319	2779.5845	0.0000	1.0000
K2PG	13	-1355.1597	2788.4021	0.0000	1.0000
K2PI	13	-1355.5917	2789.2661	0.0000	1.0000
K2PIG	14	-1355.1145	2794.3179	0.0000	1.0000
K2P	12	-1364.1097	2800.2957	0.0000	1.0000
SYMG	17	-1352.3945	2806.8971	0.0000	1.0000
SYMI	17	-1352.8637	2807.8353	0.0000	1.0000
SYMIG	18	-1352.3609	2812.8361	0.0000	1.0000
SYM	16	-1361.4579	2819.0175	0.0000	1.0000
F81G	15	-1368.3792	2826.8537	0.0000	1.0000
F81I	15	-1368.7292	2827.5536	0.0000	1.0000
F81IG	16	-1368.3579	2832.8175	0.0000	1.0000
F81	14	-1375.1931	2834.4752	0.0000	1.0000
JC69G	12	-1389.7843	2851.6448	0.0000	1.0000
JC69I	12	-1390.0110	2852.0983	0.0000	1.0000
JC69IG	13	-1389.7349	2857.5525	0.0000	1.0000
JC69	11	-1396.9542	2859.9783	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/chr3_3107.phylip = 1 - 406;
 Partition Dummy = 1:591-tomodel-phylip/chr3_3107.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/chr3_3107.phylip = 1 - 406;
 Partition Dummy = 1:591-tomodel-phylip/chr3_3107.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/chr3_3107.phylip = 1 - 406;
 Partition Dummy = 1:591-tomodel-phylip/chr3_3107.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

