
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
01:00 02/17/2011

Input data from file "183-phylip/chr9_5190.phylip" (29 taxa, 719 characters)

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


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	64	-2556.7770	5254.2757	0.6198	0.6198
GTRIG	65	-2556.0607	5255.2607	0.3787	0.9985
GTRI	64	-2562.9495	5266.6207	0.0013	0.9998
HKYG	60	-2570.0395	5271.2037	0.0001	0.9999
HKYIG	61	-2569.4463	5272.4056	0.0001	1.0000
HKYI	60	-2575.5912	5282.3071	0.0000	1.0000
K2PG	57	-2601.8395	5327.6820	0.0000	1.0000
K2PIG	58	-2601.6211	5329.6119	0.0000	1.0000
SYMG	61	-2598.5671	5330.6472	0.0000	1.0000
K2PI	57	-2603.9932	5331.9894	0.0000	1.0000
SYMIG	62	-2598.2671	5332.4427	0.0000	1.0000
SYMI	61	-2603.1357	5339.7844	0.0000	1.0000
GTR	63	-2617.9868	5374.2850	0.0000	1.0000
HKY	59	-2629.3463	5387.4362	0.0000	1.0000
K2P	56	-2655.3954	5432.4344	0.0000	1.0000
SYM	60	-2651.0674	5433.2595	0.0000	1.0000
F81G	59	-2663.9085	5456.5605	0.0000	1.0000
F81IG	60	-2663.2901	5457.7049	0.0000	1.0000
F81I	59	-2668.6345	5466.0125	0.0000	1.0000
JC69G	56	-2687.5410	5496.7255	0.0000	1.0000
JC69IG	57	-2686.8179	5497.6389	0.0000	1.0000
JC69I	56	-2688.6916	5499.0267	0.0000	1.0000
F81	58	-2718.1425	5562.6547	0.0000	1.0000
JC69	55	-2738.1333	5595.5577	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2556.7770	5241.5540	0.5697	0.5697
GTRIG	65	-2556.0607	5242.1213	0.4290	0.9987
GTRI	64	-2562.9495	5253.8990	0.0012	0.9999
HKYG	60	-2570.0395	5260.0790	0.0001	1.0000
HKYIG	61	-2569.4463	5260.8926	0.0000	1.0000
HKYI	60	-2575.5912	5271.1825	0.0000	1.0000
K2PG	57	-2601.8395	5317.6790	0.0000	1.0000
SYMG	61	-2598.5671	5319.1342	0.0000	1.0000
K2PIG	58	-2601.6211	5319.2422	0.0000	1.0000
SYMIG	62	-2598.2671	5320.5342	0.0000	1.0000
K2PI	57	-2603.9932	5321.9864	0.0000	1.0000
SYMI	61	-2603.1357	5328.2714	0.0000	1.0000
GTR	63	-2617.9868	5361.9736	0.0000	1.0000
HKY	59	-2629.3463	5376.6927	0.0000	1.0000
SYM	60	-2651.0674	5422.1349	0.0000	1.0000
K2P	56	-2655.3954	5422.7909	0.0000	1.0000
F81G	59	-2663.9085	5445.8170	0.0000	1.0000
F81IG	60	-2663.2901	5446.5803	0.0000	1.0000
F81I	59	-2668.6345	5455.2689	0.0000	1.0000
JC69G	56	-2687.5410	5487.0820	0.0000	1.0000
JC69IG	57	-2686.8179	5487.6359	0.0000	1.0000
JC69I	56	-2688.6916	5489.3832	0.0000	1.0000
F81	58	-2718.1425	5552.2850	0.0000	1.0000
JC69	55	-2738.1333	5586.2666	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-2556.7770	5534.5371	0.4899	0.4899
HKYG	60	-2570.0395	5534.7507	0.4403	0.9302
GTRIG	65	-2556.0607	5539.6823	0.0374	0.9676
HKYIG	61	-2569.4463	5540.1422	0.0297	0.9973
HKYI	60	-2575.5912	5545.8542	0.0017	0.9990
GTRI	64	-2562.9495	5546.8821	0.0010	1.0000
K2PG	57	-2601.8395	5578.6171	0.0000	1.0000
K2PI	57	-2603.9932	5582.9245	0.0000	1.0000
K2PIG	58	-2601.6211	5584.7582	0.0000	1.0000
SYMG	61	-2598.5671	5598.3838	0.0000	1.0000
SYMIG	62	-2598.2671	5604.3616	0.0000	1.0000
SYMI	61	-2603.1357	5607.5210	0.0000	1.0000
HKY	59	-2629.3463	5646.7865	0.0000	1.0000
GTR	63	-2617.9868	5650.3788	0.0000	1.0000
K2P	56	-2655.3954	5679.1511	0.0000	1.0000
SYM	60	-2651.0674	5696.8066	0.0000	1.0000
F81G	59	-2663.9085	5715.9108	0.0000	1.0000
F81IG	60	-2663.2901	5721.2520	0.0000	1.0000
F81I	59	-2668.6345	5725.3627	0.0000	1.0000
JC69G	56	-2687.5410	5743.4422	0.0000	1.0000
JC69I	56	-2688.6916	5745.7435	0.0000	1.0000
JC69IG	57	-2686.8179	5748.5740	0.0000	1.0000
F81	58	-2718.1425	5817.8009	0.0000	1.0000
JC69	55	-2738.1333	5838.0490	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 183-phylip/chr9_5190.phylip = 1 - 719;
 Partition Dummy = 1:183-phylip/chr9_5190.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 183-phylip/chr9_5190.phylip = 1 - 719;
 Partition Dummy = 1:183-phylip/chr9_5190.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 (GTRG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr9_5190.phylip = 1 - 719;
 Partition Dummy = 1:183-phylip/chr9_5190.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;


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

