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

Input data from file "261-aggressive-phylip/6804_hg18.phylip" (29 taxa, 403 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	-1300.8307	2754.2767	0.3123	0.3123
HKYG	60	-1306.8178	2755.0390	0.2133	0.5256
GTRIG	65	-1300.0076	2755.4752	0.1715	0.6971
HKYIG	61	-1305.9980	2756.1779	0.1207	0.8178
GTRI	64	-1301.9426	2756.5005	0.1027	0.9205
HKYI	60	-1307.8059	2757.0154	0.0794	0.9999
K2PG	57	-1319.5310	2772.2272	0.0000	1.0000
K2PIG	58	-1318.6612	2773.2178	0.0000	1.0000
K2PI	57	-1320.5134	2774.1921	0.0000	1.0000
SYMG	61	-1316.7917	2777.7653	0.0000	1.0000
SYMIG	62	-1315.9738	2778.9241	0.0000	1.0000
SYMI	61	-1317.8255	2779.8328	0.0000	1.0000
GTR	63	-1329.5800	2808.9475	0.0000	1.0000
HKY	59	-1337.7712	2814.1838	0.0000	1.0000
F81G	59	-1342.6538	2823.9490	0.0000	1.0000
F81IG	60	-1341.9198	2825.2431	0.0000	1.0000
F81I	59	-1343.7089	2826.0592	0.0000	1.0000
K2P	56	-1351.5166	2833.4840	0.0000	1.0000
SYM	60	-1346.8702	2835.1440	0.0000	1.0000
JC69G	56	-1354.1479	2838.7467	0.0000	1.0000
JC69IG	57	-1353.4289	2840.0229	0.0000	1.0000
JC69I	56	-1355.1468	2840.7444	0.0000	1.0000
F81	58	-1372.1036	2880.1024	0.0000	1.0000
JC69	55	-1384.7898	2897.3319	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1300.8307	2729.6613	0.4050	0.4050
GTRIG	65	-1300.0076	2730.0152	0.3393	0.7442
GTRI	64	-1301.9426	2731.8852	0.1332	0.8774
HKYG	60	-1306.8178	2733.6355	0.0555	0.9330
HKYIG	61	-1305.9980	2733.9960	0.0464	0.9793
HKYI	60	-1307.8059	2735.6119	0.0207	1.0000
K2PG	57	-1319.5310	2753.0619	0.0000	1.0000
K2PIG	58	-1318.6612	2753.3224	0.0000	1.0000
K2PI	57	-1320.5134	2755.0268	0.0000	1.0000
SYMG	61	-1316.7917	2755.5835	0.0000	1.0000
SYMIG	62	-1315.9738	2755.9476	0.0000	1.0000
SYMI	61	-1317.8255	2757.6510	0.0000	1.0000
GTR	63	-1329.5800	2785.1599	0.0000	1.0000
HKY	59	-1337.7712	2793.5424	0.0000	1.0000
F81G	59	-1342.6538	2803.3076	0.0000	1.0000
F81IG	60	-1341.9198	2803.8396	0.0000	1.0000
F81I	59	-1343.7089	2805.4178	0.0000	1.0000
SYM	60	-1346.8702	2813.7405	0.0000	1.0000
K2P	56	-1351.5166	2815.0331	0.0000	1.0000
JC69G	56	-1354.1479	2820.2959	0.0000	1.0000
JC69IG	57	-1353.4289	2820.8577	0.0000	1.0000
JC69I	56	-1355.1468	2822.2935	0.0000	1.0000
F81	58	-1372.1036	2860.2071	0.0000	1.0000
JC69	55	-1384.7898	2879.5797	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1306.8178	2973.5717	0.6556	0.6556
HKYI	60	-1307.8059	2975.5481	0.2441	0.8997
HKYIG	61	-1305.9980	2977.9312	0.0741	0.9738
K2PG	57	-1319.5310	2981.0013	0.0160	0.9898
K2PI	57	-1320.5134	2982.9662	0.0060	0.9958
K2PIG	58	-1318.6612	2985.2608	0.0019	0.9977
GTRG	64	-1300.8307	2985.5933	0.0016	0.9993
GTRI	64	-1301.9426	2987.8171	0.0005	0.9998
GTRIG	65	-1300.0076	2989.9461	0.0002	1.0000
SYMG	61	-1316.7917	2999.5186	0.0000	1.0000
SYMI	61	-1317.8255	3001.5861	0.0000	1.0000
SYMIG	62	-1315.9738	3003.8816	0.0000	1.0000
HKY	59	-1337.7712	3029.4796	0.0000	1.0000
GTR	63	-1329.5800	3037.0929	0.0000	1.0000
K2P	56	-1351.5166	3038.9736	0.0000	1.0000
F81G	59	-1342.6538	3039.2449	0.0000	1.0000
F81I	59	-1343.7089	3041.3551	0.0000	1.0000
F81IG	60	-1341.9198	3043.7758	0.0000	1.0000
JC69G	56	-1354.1479	3044.2363	0.0000	1.0000
JC69I	56	-1355.1468	3046.2340	0.0000	1.0000
JC69IG	57	-1353.4289	3048.7971	0.0000	1.0000
SYM	60	-1346.8702	3053.6767	0.0000	1.0000
F81	58	-1372.1036	3092.1454	0.0000	1.0000
JC69	55	-1384.7898	3099.5212	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/6804_hg18.phylip = 1 - 403;
 Partition Dummy = 1:261-aggressive-phylip/6804_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/6804_hg18.phylip = 1 - 403;
 Partition Dummy = 1:261-aggressive-phylip/6804_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/6804_hg18.phylip = 1 - 403;
 Partition Dummy = 1:261-aggressive-phylip/6804_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

