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

Input data from file "261-aggressive-phylip/4111_hg18.phylip" (29 taxa, 311 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	64	-1166.9848	2495.7908	0.3253	0.3253
HKYG	60	-1173.2846	2495.8492	0.3160	0.6413
HKYIG	61	-1172.0934	2496.5642	0.2210	0.8623
GTRIG	65	-1166.4331	2497.8865	0.1141	0.9763
HKYI	60	-1176.2289	2501.7379	0.0166	0.9930
GTRI	64	-1170.8190	2503.4591	0.0070	1.0000
SYMG	61	-1187.5024	2527.3823	0.0000	1.0000
SYMIG	62	-1186.8918	2529.2835	0.0000	1.0000
K2PG	57	-1196.3379	2532.8101	0.0000	1.0000
SYMI	61	-1190.6145	2533.6066	0.0000	1.0000
K2PIG	58	-1195.2518	2533.6624	0.0000	1.0000
F81G	59	-1193.9312	2534.0695	0.0000	1.0000
F81IG	60	-1193.0420	2535.3640	0.0000	1.0000
K2PI	57	-1199.1785	2538.4914	0.0000	1.0000
F81I	59	-1196.4278	2539.0627	0.0000	1.0000
JC69G	56	-1212.4158	2561.9655	0.0000	1.0000
JC69IG	57	-1211.4103	2562.9550	0.0000	1.0000
JC69I	56	-1215.0383	2567.2104	0.0000	1.0000
GTR	63	-1207.4105	2573.4687	0.0000	1.0000
HKY	59	-1215.3941	2576.9954	0.0000	1.0000
SYM	60	-1225.2865	2599.8531	0.0000	1.0000
K2P	56	-1235.7129	2608.5597	0.0000	1.0000
F81	58	-1233.3955	2609.9498	0.0000	1.0000
JC69	55	-1251.1802	2636.5173	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1166.9848	2461.9696	0.5298	0.5298
GTRIG	65	-1166.4331	2462.8661	0.3384	0.8683
HKYIG	61	-1172.0934	2466.1867	0.0643	0.9326
HKYG	60	-1173.2846	2466.5692	0.0531	0.9857
GTRI	64	-1170.8190	2469.6380	0.0115	0.9972
HKYI	60	-1176.2289	2472.4579	0.0028	1.0000
SYMG	61	-1187.5024	2497.0047	0.0000	1.0000
SYMIG	62	-1186.8918	2497.7835	0.0000	1.0000
SYMI	61	-1190.6145	2503.2291	0.0000	1.0000
F81G	59	-1193.9312	2505.8623	0.0000	1.0000
F81IG	60	-1193.0420	2506.0840	0.0000	1.0000
K2PIG	58	-1195.2518	2506.5037	0.0000	1.0000
K2PG	57	-1196.3379	2506.6757	0.0000	1.0000
F81I	59	-1196.4278	2510.8556	0.0000	1.0000
K2PI	57	-1199.1785	2512.3571	0.0000	1.0000
JC69IG	57	-1211.4103	2536.8206	0.0000	1.0000
JC69G	56	-1212.4158	2536.8316	0.0000	1.0000
GTR	63	-1207.4105	2540.8209	0.0000	1.0000
JC69I	56	-1215.0383	2542.0765	0.0000	1.0000
HKY	59	-1215.3941	2548.7882	0.0000	1.0000
SYM	60	-1225.2865	2570.5731	0.0000	1.0000
F81	58	-1233.3955	2582.7911	0.0000	1.0000
K2P	56	-1235.7129	2583.4258	0.0000	1.0000
JC69	55	-1251.1802	2612.3605	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1173.2846	2690.9568	0.8028	0.8028
HKYIG	61	-1172.0934	2694.3141	0.1498	0.9527
HKYI	60	-1176.2289	2696.8455	0.0423	0.9949
GTRG	64	-1166.9848	2701.3164	0.0045	0.9995
GTRIG	65	-1166.4331	2705.9527	0.0004	0.9999
GTRI	64	-1170.8190	2708.9847	0.0001	1.0000
K2PG	57	-1196.3379	2719.8439	0.0000	1.0000
K2PIG	58	-1195.2518	2723.4117	0.0000	1.0000
SYMG	61	-1187.5024	2725.1321	0.0000	1.0000
K2PI	57	-1199.1785	2725.5253	0.0000	1.0000
F81G	59	-1193.9312	2726.5101	0.0000	1.0000
SYMIG	62	-1186.8918	2729.6507	0.0000	1.0000
F81IG	60	-1193.0420	2730.4716	0.0000	1.0000
SYMI	61	-1190.6145	2731.3564	0.0000	1.0000
F81I	59	-1196.4278	2731.5033	0.0000	1.0000
JC69G	56	-1212.4158	2746.2600	0.0000	1.0000
JC69IG	57	-1211.4103	2749.9888	0.0000	1.0000
JC69I	56	-1215.0383	2751.5049	0.0000	1.0000
HKY	59	-1215.3941	2769.4360	0.0000	1.0000
GTR	63	-1207.4105	2776.4279	0.0000	1.0000
K2P	56	-1235.7129	2792.8542	0.0000	1.0000
SYM	60	-1225.2865	2794.9606	0.0000	1.0000
F81	58	-1233.3955	2799.6991	0.0000	1.0000
JC69	55	-1251.1802	2818.0491	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 261-aggressive-phylip/4111_hg18.phylip = 1 - 311;
 Partition Dummy = 1:261-aggressive-phylip/4111_hg18.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 261-aggressive-phylip/4111_hg18.phylip = 1 - 311;
 Partition Dummy = 1:261-aggressive-phylip/4111_hg18.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 261-aggressive-phylip/4111_hg18.phylip = 1 - 311;
 Partition Dummy = 1:261-aggressive-phylip/4111_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

