
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:52 08/01/2011

Input data from file "261-aggressive-phylip/7723_hg18.phylip" (29 taxa, 507 characters)

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	60	-1226.3691	2589.1508	0.4717	0.4717
HKYI	60	-1226.9118	2590.2361	0.2741	0.7458
HKYIG	61	-1225.8083	2590.6144	0.2269	0.9727
GTRG	64	-1224.7063	2596.2362	0.0136	0.9863
GTRI	64	-1225.2807	2597.3850	0.0077	0.9940
GTRIG	65	-1224.2192	2597.8942	0.0060	1.0000
HKY	59	-1237.5374	2608.9137	0.0000	1.0000
K2PG	57	-1241.9126	2612.5512	0.0000	1.0000
K2PI	57	-1242.2845	2613.2950	0.0000	1.0000
K2PIG	58	-1241.4483	2614.1733	0.0000	1.0000
GTR	63	-1235.5413	2615.2859	0.0000	1.0000
SYMG	61	-1239.6960	2618.3898	0.0000	1.0000
SYMI	61	-1240.2245	2619.4468	0.0000	1.0000
SYMIG	62	-1239.1674	2619.9295	0.0000	1.0000
F81G	59	-1243.6553	2621.1494	0.0000	1.0000
F81I	59	-1244.1988	2622.2365	0.0000	1.0000
F81IG	60	-1243.1512	2622.7149	0.0000	1.0000
K2P	56	-1253.0781	2632.3428	0.0000	1.0000
SYM	60	-1250.7726	2637.9578	0.0000	1.0000
JC69G	56	-1256.5458	2639.2783	0.0000	1.0000
F81	58	-1254.4865	2640.2498	0.0000	1.0000
JC69I	56	-1257.0821	2640.3509	0.0000	1.0000
JC69IG	57	-1256.0127	2640.7514	0.0000	1.0000
JC69	55	-1267.5541	2658.7668	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1226.3691	2572.7382	0.4107	0.4107
HKYIG	61	-1225.8083	2573.6166	0.2647	0.6755
HKYI	60	-1226.9118	2573.8235	0.2387	0.9142
GTRG	64	-1224.7063	2577.4127	0.0397	0.9539
GTRIG	65	-1224.2192	2578.4384	0.0238	0.9776
GTRI	64	-1225.2807	2578.5615	0.0223	1.0000
HKY	59	-1237.5374	2593.0748	0.0000	1.0000
GTR	63	-1235.5413	2597.0827	0.0000	1.0000
K2PG	57	-1241.9126	2597.8252	0.0000	1.0000
K2PI	57	-1242.2845	2598.5689	0.0000	1.0000
K2PIG	58	-1241.4483	2598.8965	0.0000	1.0000
SYMG	61	-1239.6960	2601.3920	0.0000	1.0000
SYMIG	62	-1239.1674	2602.3349	0.0000	1.0000
SYMI	61	-1240.2245	2602.4490	0.0000	1.0000
F81G	59	-1243.6553	2605.3105	0.0000	1.0000
F81IG	60	-1243.1512	2606.3024	0.0000	1.0000
F81I	59	-1244.1988	2606.3976	0.0000	1.0000
K2P	56	-1253.0781	2618.1561	0.0000	1.0000
SYM	60	-1250.7726	2621.5453	0.0000	1.0000
F81	58	-1254.4865	2624.9730	0.0000	1.0000
JC69G	56	-1256.5458	2625.0917	0.0000	1.0000
JC69IG	57	-1256.0127	2626.0253	0.0000	1.0000
JC69I	56	-1257.0821	2626.1642	0.0000	1.0000
JC69	55	-1267.5541	2645.1083	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1226.3691	2826.4489	0.6013	0.6013
HKYI	60	-1226.9118	2827.5342	0.3495	0.9508
HKYIG	61	-1225.8083	2831.5558	0.0468	0.9976
K2PG	57	-1241.9126	2838.8503	0.0012	0.9989
K2PI	57	-1242.2845	2839.5940	0.0008	0.9997
HKY	59	-1237.5374	2842.5569	0.0002	0.9999
K2PIG	58	-1241.4483	2844.1502	0.0001	1.0000
GTRG	64	-1224.7063	2848.0374	0.0000	1.0000
GTRI	64	-1225.2807	2849.1862	0.0000	1.0000
GTRIG	65	-1224.2192	2853.2916	0.0000	1.0000
F81G	59	-1243.6553	2854.7926	0.0000	1.0000
K2P	56	-1253.0781	2854.9527	0.0000	1.0000
F81I	59	-1244.1988	2855.8797	0.0000	1.0000
SYMG	61	-1239.6960	2859.3312	0.0000	1.0000
F81IG	60	-1243.1512	2860.0130	0.0000	1.0000
SYMI	61	-1240.2245	2860.3882	0.0000	1.0000
JC69G	56	-1256.5458	2861.8883	0.0000	1.0000
JC69I	56	-1257.0821	2862.9608	0.0000	1.0000
GTR	63	-1235.5413	2863.4789	0.0000	1.0000
SYMIG	62	-1239.1674	2864.5026	0.0000	1.0000
JC69IG	57	-1256.0127	2867.0505	0.0000	1.0000
F81	58	-1254.4865	2870.2267	0.0000	1.0000
SYM	60	-1250.7726	2875.2559	0.0000	1.0000
JC69	55	-1267.5541	2877.6764	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/7723_hg18.phylip = 1 - 507;
 Partition Dummy = 1:261-aggressive-phylip/7723_hg18.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;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/7723_hg18.phylip = 1 - 507;
 Partition Dummy = 1:261-aggressive-phylip/7723_hg18.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;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/7723_hg18.phylip = 1 - 507;
 Partition Dummy = 1:261-aggressive-phylip/7723_hg18.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

