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

Input data from file "591-tomodel-phylip/chr6_3033.phylip" (7 taxa, 744 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	-2594.1324	5229.4266	0.5321	0.5321
GTRI	20	-2594.8122	5230.7862	0.2696	0.8017
GTRIG	21	-2594.1409	5231.5617	0.1830	0.9847
HKYG	16	-2602.6038	5237.9560	0.0075	0.9922
HKYI	16	-2602.9782	5238.7046	0.0051	0.9973
HKYIG	17	-2602.5800	5240.0030	0.0027	1.0000
GTR	19	-2605.8392	5250.7281	0.0000	1.0000
HKY	15	-2613.7246	5258.1086	0.0000	1.0000
SYMI	17	-2617.9769	5270.7969	0.0000	1.0000
SYMG	17	-2618.1916	5271.2261	0.0000	1.0000
K2PI	13	-2622.7016	5271.9019	0.0000	1.0000
K2PG	13	-2622.7582	5272.0150	0.0000	1.0000
SYMIG	18	-2617.9417	5272.8269	0.0000	1.0000
K2PIG	14	-2622.5823	5273.7407	0.0000	1.0000
SYM	16	-2628.7266	5290.2015	0.0000	1.0000
K2P	12	-2633.3073	5291.0415	0.0000	1.0000
F81I	15	-2668.5044	5367.6681	0.0000	1.0000
F81G	15	-2668.7601	5368.1795	0.0000	1.0000
F81IG	16	-2668.4798	5369.7079	0.0000	1.0000
F81	14	-2678.1890	5384.9542	0.0000	1.0000
JC69I	12	-2683.5366	5391.5000	0.0000	1.0000
JC69G	12	-2684.0195	5392.4658	0.0000	1.0000
JC69IG	13	-2683.5368	5393.5721	0.0000	1.0000
JC69	11	-2693.3069	5408.9744	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2594.1324	5228.2648	0.5277	0.5277
GTRI	20	-2594.8122	5229.6244	0.2674	0.7951
GTRIG	21	-2594.1409	5230.2819	0.1925	0.9875
HKYG	16	-2602.6038	5237.2077	0.0060	0.9936
HKYI	16	-2602.9782	5237.9564	0.0041	0.9977
HKYIG	17	-2602.5800	5239.1600	0.0023	1.0000
GTR	19	-2605.8392	5249.6784	0.0000	1.0000
HKY	15	-2613.7246	5257.4492	0.0000	1.0000
SYMI	17	-2617.9769	5269.9539	0.0000	1.0000
SYMG	17	-2618.1916	5270.3831	0.0000	1.0000
K2PI	13	-2622.7016	5271.4033	0.0000	1.0000
K2PG	13	-2622.7582	5271.5163	0.0000	1.0000
SYMIG	18	-2617.9417	5271.8835	0.0000	1.0000
K2PIG	14	-2622.5823	5273.1645	0.0000	1.0000
SYM	16	-2628.7266	5289.4532	0.0000	1.0000
K2P	12	-2633.3073	5290.6147	0.0000	1.0000
F81I	15	-2668.5044	5367.0088	0.0000	1.0000
F81G	15	-2668.7601	5367.5202	0.0000	1.0000
F81IG	16	-2668.4798	5368.9596	0.0000	1.0000
F81	14	-2678.1890	5384.3781	0.0000	1.0000
JC69I	12	-2683.5366	5391.0731	0.0000	1.0000
JC69G	12	-2684.0195	5392.0389	0.0000	1.0000
JC69IG	13	-2683.5368	5393.0735	0.0000	1.0000
JC69	11	-2693.3069	5408.6137	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2602.6038	5311.0003	0.5750	0.5750
HKYI	16	-2602.9782	5311.7490	0.3955	0.9705
HKYIG	17	-2602.5800	5317.5647	0.0216	0.9921
GTRG	20	-2594.1324	5320.5056	0.0050	0.9970
GTRI	20	-2594.8122	5321.8652	0.0025	0.9995
HKY	15	-2613.7246	5326.6298	0.0002	0.9998
GTRIG	21	-2594.1409	5327.1348	0.0002	1.0000
K2PI	13	-2622.7016	5331.3598	0.0000	1.0000
K2PG	13	-2622.7582	5331.4729	0.0000	1.0000
GTR	19	-2605.8392	5337.3072	0.0000	1.0000
K2PIG	14	-2622.5823	5337.7331	0.0000	1.0000
K2P	12	-2633.3073	5345.9592	0.0000	1.0000
SYMI	17	-2617.9769	5348.3586	0.0000	1.0000
SYMG	17	-2618.1916	5348.7878	0.0000	1.0000
SYMIG	18	-2617.9417	5354.9002	0.0000	1.0000
SYM	16	-2628.7266	5363.2459	0.0000	1.0000
F81I	15	-2668.5044	5436.1894	0.0000	1.0000
F81G	15	-2668.7601	5436.7008	0.0000	1.0000
F81IG	16	-2668.4798	5442.7523	0.0000	1.0000
JC69I	12	-2683.5366	5446.4176	0.0000	1.0000
JC69G	12	-2684.0195	5447.3834	0.0000	1.0000
F81	14	-2678.1890	5448.9466	0.0000	1.0000
JC69IG	13	-2683.5368	5453.0300	0.0000	1.0000
JC69	11	-2693.3069	5459.3462	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/chr6_3033.phylip = 1 - 744;
 Partition Dummy = 1:591-tomodel-phylip/chr6_3033.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/chr6_3033.phylip = 1 - 744;
 Partition Dummy = 1:591-tomodel-phylip/chr6_3033.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/chr6_3033.phylip = 1 - 744;
 Partition Dummy = 1:591-tomodel-phylip/chr6_3033.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

