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

Input data from file "591-tomodel-phylip/chr3_5552.phylip" (7 taxa, 496 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	-1369.7034	2781.1752	0.4475	0.4475
GTRIG	21	-1369.5380	2783.0254	0.1774	0.6249
GTRI	20	-1370.6409	2783.0501	0.1752	0.8001
HKYG	16	-1375.4242	2783.9842	0.1099	0.9100
HKYIG	17	-1375.1768	2785.6340	0.0481	0.9581
HKYI	16	-1376.3885	2785.9126	0.0419	1.0000
GTR	19	-1387.8363	2815.2692	0.0000	1.0000
K2PG	13	-1394.7083	2816.1718	0.0000	1.0000
HKY	15	-1392.8410	2816.6819	0.0000	1.0000
K2PIG	14	-1394.4693	2817.8118	0.0000	1.0000
K2PI	13	-1395.6483	2818.0518	0.0000	1.0000
SYMG	17	-1392.1967	2819.6737	0.0000	1.0000
SYMIG	18	-1391.9872	2821.4084	0.0000	1.0000
SYMI	17	-1393.0952	2821.4707	0.0000	1.0000
K2P	12	-1410.9158	2846.4775	0.0000	1.0000
SYM	16	-1408.3452	2849.8262	0.0000	1.0000
F81G	15	-1409.7703	2850.5407	0.0000	1.0000
F81I	15	-1410.3799	2851.7599	0.0000	1.0000
F81IG	16	-1409.6353	2852.4063	0.0000	1.0000
JC69G	12	-1424.9500	2874.5460	0.0000	1.0000
JC69I	12	-1425.4866	2875.6191	0.0000	1.0000
JC69IG	13	-1424.8044	2876.3641	0.0000	1.0000
F81	14	-1424.6513	2878.1758	0.0000	1.0000
JC69	11	-1439.3483	2901.2421	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1369.7034	2779.4068	0.4636	0.4636
GTRIG	21	-1369.5380	2781.0760	0.2012	0.6648
GTRI	20	-1370.6409	2781.2817	0.1815	0.8464
HKYG	16	-1375.4242	2782.8485	0.0829	0.9293
HKYIG	17	-1375.1768	2784.3537	0.0391	0.9684
HKYI	16	-1376.3885	2784.7769	0.0316	1.0000
GTR	19	-1387.8363	2813.6726	0.0000	1.0000
K2PG	13	-1394.7083	2815.4166	0.0000	1.0000
HKY	15	-1392.8410	2815.6819	0.0000	1.0000
K2PIG	14	-1394.4693	2816.9387	0.0000	1.0000
K2PI	13	-1395.6483	2817.2966	0.0000	1.0000
SYMG	17	-1392.1967	2818.3934	0.0000	1.0000
SYMIG	18	-1391.9872	2819.9744	0.0000	1.0000
SYMI	17	-1393.0952	2820.1904	0.0000	1.0000
K2P	12	-1410.9158	2845.8316	0.0000	1.0000
SYM	16	-1408.3452	2848.6905	0.0000	1.0000
F81G	15	-1409.7703	2849.5407	0.0000	1.0000
F81I	15	-1410.3799	2850.7599	0.0000	1.0000
F81IG	16	-1409.6353	2851.2706	0.0000	1.0000
JC69G	12	-1424.9500	2873.9001	0.0000	1.0000
JC69I	12	-1425.4866	2874.9731	0.0000	1.0000
JC69IG	13	-1424.8044	2875.6089	0.0000	1.0000
F81	14	-1424.6513	2877.3027	0.0000	1.0000
JC69	11	-1439.3483	2900.6966	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1375.4242	2850.1537	0.6941	0.6941
HKYI	16	-1376.3885	2852.0822	0.2647	0.9588
HKYIG	17	-1375.1768	2855.8655	0.0399	0.9987
GTRG	20	-1369.7034	2863.5383	0.0009	0.9996
GTRI	20	-1370.6409	2865.4132	0.0003	0.9999
GTRIG	21	-1369.5380	2869.4141	0.0000	1.0000
K2PG	13	-1394.7083	2870.1021	0.0000	1.0000
K2PI	13	-1395.6483	2871.9821	0.0000	1.0000
K2PIG	14	-1394.4693	2875.8307	0.0000	1.0000
HKY	15	-1392.8410	2878.7806	0.0000	1.0000
SYMG	17	-1392.1967	2889.9052	0.0000	1.0000
SYMI	17	-1393.0952	2891.7022	0.0000	1.0000
GTR	19	-1387.8363	2893.5975	0.0000	1.0000
SYMIG	18	-1391.9872	2895.6928	0.0000	1.0000
K2P	12	-1410.9158	2896.3105	0.0000	1.0000
F81G	15	-1409.7703	2912.6393	0.0000	1.0000
F81I	15	-1410.3799	2913.8585	0.0000	1.0000
SYM	16	-1408.3452	2915.9957	0.0000	1.0000
F81IG	16	-1409.6353	2918.5758	0.0000	1.0000
JC69G	12	-1424.9500	2924.3790	0.0000	1.0000
JC69I	12	-1425.4866	2925.4520	0.0000	1.0000
JC69IG	13	-1424.8044	2930.2944	0.0000	1.0000
F81	14	-1424.6513	2936.1947	0.0000	1.0000
JC69	11	-1439.3483	2946.9690	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_5552.phylip = 1 - 496;
 Partition Dummy = 1:591-tomodel-phylip/chr3_5552.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_5552.phylip = 1 - 496;
 Partition Dummy = 1:591-tomodel-phylip/chr3_5552.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_5552.phylip = 1 - 496;
 Partition Dummy = 1:591-tomodel-phylip/chr3_5552.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

