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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
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
GTRI	64	-613.5617	1394.3687	0.2623	0.2623
GTRG	64	-613.5918	1394.4289	0.2546	0.5169
HKYG	60	-620.5270	1394.9429	0.1969	0.7138
HKYI	60	-620.5900	1395.0690	0.1848	0.8986
GTRIG	65	-613.5545	1397.7724	0.0478	0.9464
HKYIG	61	-620.3319	1397.8452	0.0461	0.9926
HKY	59	-626.3436	1403.3139	0.0030	0.9956
GTR	63	-620.8779	1405.6149	0.0009	0.9965
F81G	59	-627.7229	1406.0725	0.0008	0.9973
K2PG	57	-630.9776	1406.1470	0.0007	0.9980
F81I	59	-627.8067	1406.2401	0.0007	0.9987
K2PI	57	-631.0707	1406.3331	0.0007	0.9993
F81IG	60	-627.5569	1409.0027	0.0002	0.9995
K2PIG	58	-630.8209	1409.0363	0.0002	0.9997
SYMG	61	-626.2877	1409.7568	0.0001	0.9998
SYMI	61	-626.3370	1409.8554	0.0001	0.9999
SYMIG	62	-626.2073	1412.9192	0.0000	0.9999
K2P	56	-636.2598	1413.5379	0.0000	1.0000
F81	58	-633.1647	1413.7239	0.0000	1.0000
JC69G	56	-636.9831	1414.9843	0.0000	1.0000
JC69I	56	-637.0783	1415.1747	0.0000	1.0000
SYM	60	-631.8129	1417.5146	0.0000	1.0000
JC69IG	57	-636.8330	1417.8578	0.0000	1.0000
JC69	55	-642.2130	1422.2993	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	64	-613.5617	1355.1235	0.4054	0.4054
GTRG	64	-613.5918	1355.1836	0.3934	0.7988
GTRIG	65	-613.5545	1357.1089	0.1502	0.9490
HKYG	60	-620.5270	1361.0540	0.0209	0.9699
HKYI	60	-620.5900	1361.1801	0.0196	0.9896
HKYIG	61	-620.3319	1362.6638	0.0093	0.9989
GTR	63	-620.8779	1367.7557	0.0007	0.9996
HKY	59	-626.3436	1370.6872	0.0002	0.9998
F81G	59	-627.7229	1373.4458	0.0000	0.9999
F81I	59	-627.8067	1373.6133	0.0000	0.9999
SYMG	61	-626.2877	1374.5754	0.0000	0.9999
SYMI	61	-626.3370	1374.6740	0.0000	0.9999
F81IG	60	-627.5569	1375.1138	0.0000	1.0000
K2PG	57	-630.9776	1375.9552	0.0000	1.0000
K2PI	57	-631.0707	1376.1413	0.0000	1.0000
SYMIG	62	-626.2073	1376.4145	0.0000	1.0000
K2PIG	58	-630.8209	1377.6418	0.0000	1.0000
F81	58	-633.1647	1382.3294	0.0000	1.0000
SYM	60	-631.8129	1383.6257	0.0000	1.0000
K2P	56	-636.2598	1384.5197	0.0000	1.0000
JC69G	56	-636.9831	1385.9661	0.0000	1.0000
JC69I	56	-637.0783	1386.1565	0.0000	1.0000
JC69IG	57	-636.8330	1387.6661	0.0000	1.0000
JC69	55	-642.2130	1394.4260	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-620.5270	1578.4951	0.4162	0.4162
HKYI	60	-620.5900	1578.6211	0.3908	0.8070
K2PG	57	-630.9776	1582.5242	0.0555	0.8625
K2PI	57	-631.0707	1582.7103	0.0506	0.9131
HKYIG	61	-620.3319	1583.7289	0.0304	0.9435
HKY	59	-626.3436	1584.5042	0.0206	0.9641
GTRI	64	-613.5617	1587.0606	0.0057	0.9698
GTRG	64	-613.5918	1587.1207	0.0056	0.9754
F81G	59	-627.7229	1587.2628	0.0052	0.9806
F81I	59	-627.8067	1587.4304	0.0048	0.9854
K2P	56	-636.2598	1587.4647	0.0047	0.9901
K2PIG	58	-630.8209	1587.8348	0.0039	0.9940
JC69G	56	-636.9831	1588.9111	0.0023	0.9962
JC69I	56	-637.0783	1589.1015	0.0021	0.9983
F81	58	-633.1647	1592.5225	0.0004	0.9987
F81IG	60	-627.5569	1592.5549	0.0004	0.9991
GTRIG	65	-613.5545	1592.6701	0.0003	0.9994
JC69	55	-642.2130	1593.7469	0.0002	0.9996
JC69IG	57	-636.8330	1594.2351	0.0002	0.9998
SYMG	61	-626.2877	1595.6405	0.0001	0.9999
SYMI	61	-626.3370	1595.7391	0.0001	0.9999
GTR	63	-620.8779	1596.0688	0.0001	1.0000
SYM	60	-631.8129	1601.0668	0.0000	1.0000
SYMIG	62	-626.2073	1601.1036	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/11500_hg18.phylip = 1 - 277;
 Partition Dummy = 1:261-aggressive-phylip/11500_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/11500_hg18.phylip = 1 - 277;
 Partition Dummy = 1:261-aggressive-phylip/11500_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

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

