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

Input data from file "261-aggressive-phylip/3846_hg18.phylip" (29 taxa, 480 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	-1378.3965	2894.2631	0.4256	0.4256
HKYI	60	-1378.6365	2894.7432	0.3348	0.7604
HKYIG	61	-1378.3735	2896.8428	0.1172	0.8775
GTRG	64	-1375.0414	2898.1309	0.0615	0.9391
GTRI	64	-1375.3923	2898.8327	0.0433	0.9824
GTRIG	65	-1374.9540	2900.6327	0.0176	1.0000
K2PG	57	-1394.9903	2919.6488	0.0000	1.0000
K2PI	57	-1395.3232	2920.3146	0.0000	1.0000
HKY	59	-1393.0746	2921.0064	0.0000	1.0000
K2PIG	58	-1394.9384	2922.1334	0.0000	1.0000
GTR	63	-1390.1232	2925.6310	0.0000	1.0000
SYMG	61	-1393.5637	2927.2231	0.0000	1.0000
SYMI	61	-1394.0709	2928.2376	0.0000	1.0000
SYMIG	62	-1393.4134	2929.5607	0.0000	1.0000
F81G	59	-1400.0706	2934.9984	0.0000	1.0000
F81I	59	-1400.4368	2935.7306	0.0000	1.0000
F81IG	60	-1399.9815	2937.4332	0.0000	1.0000
K2P	56	-1409.4678	2946.0279	0.0000	1.0000
SYM	60	-1407.8883	2953.2467	0.0000	1.0000
JC69G	56	-1414.9031	2956.8983	0.0000	1.0000
JC69I	56	-1415.3292	2957.7507	0.0000	1.0000
JC69IG	57	-1414.7875	2959.2431	0.0000	1.0000
F81	58	-1414.0122	2960.2809	0.0000	1.0000
JC69	55	-1428.7623	2982.0529	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1378.3965	2876.7929	0.3060	0.3060
HKYI	60	-1378.6365	2877.2731	0.2407	0.5467
GTRG	64	-1375.0414	2878.0827	0.1606	0.7073
HKYIG	61	-1378.3735	2878.7471	0.1152	0.8225
GTRI	64	-1375.3923	2878.7845	0.1131	0.9355
GTRIG	65	-1374.9540	2879.9081	0.0645	1.0000
K2PG	57	-1394.9903	2903.9806	0.0000	1.0000
HKY	59	-1393.0746	2904.1492	0.0000	1.0000
K2PI	57	-1395.3232	2904.6463	0.0000	1.0000
K2PIG	58	-1394.9384	2905.8768	0.0000	1.0000
GTR	63	-1390.1232	2906.2464	0.0000	1.0000
SYMG	61	-1393.5637	2909.1274	0.0000	1.0000
SYMI	61	-1394.0709	2910.1419	0.0000	1.0000
SYMIG	62	-1393.4134	2910.8268	0.0000	1.0000
F81G	59	-1400.0706	2918.1413	0.0000	1.0000
F81I	59	-1400.4368	2918.8735	0.0000	1.0000
F81IG	60	-1399.9815	2919.9630	0.0000	1.0000
K2P	56	-1409.4678	2930.9357	0.0000	1.0000
SYM	60	-1407.8883	2935.7766	0.0000	1.0000
JC69G	56	-1414.9031	2941.8061	0.0000	1.0000
JC69I	56	-1415.3292	2942.6585	0.0000	1.0000
JC69IG	57	-1414.7875	2943.5749	0.0000	1.0000
F81	58	-1414.0122	2944.0244	0.0000	1.0000
JC69	55	-1428.7623	2967.5246	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1378.3965	3127.2201	0.5451	0.5451
HKYI	60	-1378.6365	3127.7002	0.4287	0.9738
HKYIG	61	-1378.3735	3133.3480	0.0255	0.9992
K2PG	57	-1394.9903	3141.8864	0.0004	0.9996
K2PI	57	-1395.3232	3142.5521	0.0003	0.9999
GTRG	64	-1375.0414	3145.2051	0.0001	0.9999
GTRI	64	-1375.3923	3145.9068	0.0000	1.0000
K2PIG	58	-1394.9384	3147.9564	0.0000	1.0000
HKY	59	-1393.0746	3150.4026	0.0000	1.0000
GTRIG	65	-1374.9540	3151.2042	0.0000	1.0000
SYMG	61	-1393.5637	3163.7284	0.0000	1.0000
F81G	59	-1400.0706	3164.3947	0.0000	1.0000
K2P	56	-1409.4678	3164.6677	0.0000	1.0000
SYMI	61	-1394.0709	3164.7428	0.0000	1.0000
F81I	59	-1400.4368	3165.1269	0.0000	1.0000
GTR	63	-1390.1232	3169.1949	0.0000	1.0000
SYMIG	62	-1393.4134	3169.6016	0.0000	1.0000
F81IG	60	-1399.9815	3170.3902	0.0000	1.0000
JC69G	56	-1414.9031	3175.5382	0.0000	1.0000
JC69I	56	-1415.3292	3176.3905	0.0000	1.0000
JC69IG	57	-1414.7875	3181.4807	0.0000	1.0000
F81	58	-1414.0122	3186.1040	0.0000	1.0000
SYM	60	-1407.8883	3186.2037	0.0000	1.0000
JC69	55	-1428.7623	3197.0828	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/3846_hg18.phylip = 1 - 480;
 Partition Dummy = 1:261-aggressive-phylip/3846_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/3846_hg18.phylip = 1 - 480;
 Partition Dummy = 1:261-aggressive-phylip/3846_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/3846_hg18.phylip = 1 - 480;
 Partition Dummy = 1:261-aggressive-phylip/3846_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

