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

Input data from file "261-aggressive-phylip/7309_hg18.phylip" (29 taxa, 511 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	-2965.2466	6077.1479	0.6982	0.6982
GTRIG	65	-2965.1377	6079.5563	0.2094	0.9076
HKYG	60	-2972.7985	6081.8637	0.0661	0.9737
HKYIG	61	-2972.4287	6083.7038	0.0263	1.0000
HKYI	60	-2987.1647	6110.5960	0.0000	1.0000
GTRI	64	-2982.6756	6112.0060	0.0000	1.0000
K2PG	57	-3013.7871	6156.1703	0.0000	1.0000
K2PIG	58	-3014.8329	6160.8074	0.0000	1.0000
SYMG	61	-3013.1476	6165.1415	0.0000	1.0000
SYMIG	62	-3012.9148	6167.2672	0.0000	1.0000
K2PI	57	-3028.1488	6184.8936	0.0000	1.0000
SYMI	61	-3025.8274	6190.5010	0.0000	1.0000
HKY	59	-3042.3650	6218.4284	0.0000	1.0000
GTR	63	-3039.8059	6223.6520	0.0000	1.0000
F81G	59	-3069.0153	6271.7291	0.0000	1.0000
F81IG	60	-3068.5642	6273.3950	0.0000	1.0000
F81I	59	-3081.1761	6296.0505	0.0000	1.0000
K2P	56	-3089.9682	6305.9980	0.0000	1.0000
SYM	60	-3089.4370	6315.1407	0.0000	1.0000
JC69G	56	-3097.1462	6320.3540	0.0000	1.0000
JC69IG	57	-3096.6398	6321.8756	0.0000	1.0000
JC69I	56	-3109.3417	6344.7450	0.0000	1.0000
F81	58	-3143.5402	6418.2220	0.0000	1.0000
JC69	55	-3174.8851	6473.3086	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2965.2466	6058.4932	0.6877	0.6877
GTRIG	65	-2965.1377	6060.2754	0.2821	0.9698
HKYG	60	-2972.7985	6065.5970	0.0197	0.9895
HKYIG	61	-2972.4287	6066.8575	0.0105	1.0000
GTRI	64	-2982.6756	6093.3513	0.0000	1.0000
HKYI	60	-2987.1647	6094.3294	0.0000	1.0000
K2PG	57	-3013.7871	6141.5743	0.0000	1.0000
K2PIG	58	-3014.8329	6145.6658	0.0000	1.0000
SYMG	61	-3013.1476	6148.2952	0.0000	1.0000
SYMIG	62	-3012.9148	6149.8297	0.0000	1.0000
K2PI	57	-3028.1488	6170.2976	0.0000	1.0000
SYMI	61	-3025.8274	6173.6547	0.0000	1.0000
HKY	59	-3042.3650	6202.7299	0.0000	1.0000
GTR	63	-3039.8059	6205.6117	0.0000	1.0000
F81G	59	-3069.0153	6256.0306	0.0000	1.0000
F81IG	60	-3068.5642	6257.1284	0.0000	1.0000
F81I	59	-3081.1761	6280.3521	0.0000	1.0000
K2P	56	-3089.9682	6291.9364	0.0000	1.0000
SYM	60	-3089.4370	6298.8740	0.0000	1.0000
JC69G	56	-3097.1462	6306.2923	0.0000	1.0000
JC69IG	57	-3096.6398	6307.2796	0.0000	1.0000
JC69I	56	-3109.3417	6330.6833	0.0000	1.0000
F81	58	-3143.5402	6403.0804	0.0000	1.0000
JC69	55	-3174.8851	6459.7702	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2972.7985	6319.7792	0.9331	0.9331
HKYIG	61	-2972.4287	6325.2760	0.0597	0.9929
GTRG	64	-2965.2466	6329.6208	0.0068	0.9997
GTRIG	65	-2965.1377	6335.6395	0.0003	1.0000
HKYI	60	-2987.1647	6348.5115	0.0000	1.0000
GTRI	64	-2982.6756	6364.4789	0.0000	1.0000
K2PG	57	-3013.7871	6383.0473	0.0000	1.0000
K2PIG	58	-3014.8329	6391.3753	0.0000	1.0000
SYMG	61	-3013.1476	6406.7137	0.0000	1.0000
K2PI	57	-3028.1488	6411.7707	0.0000	1.0000
SYMIG	62	-3012.9148	6412.4846	0.0000	1.0000
SYMI	61	-3025.8274	6432.0733	0.0000	1.0000
HKY	59	-3042.3650	6452.6757	0.0000	1.0000
GTR	63	-3039.8059	6472.5030	0.0000	1.0000
F81G	59	-3069.0153	6505.9764	0.0000	1.0000
F81IG	60	-3068.5642	6511.3105	0.0000	1.0000
K2P	56	-3089.9682	6529.1731	0.0000	1.0000
F81I	59	-3081.1761	6530.2979	0.0000	1.0000
JC69G	56	-3097.1462	6543.5290	0.0000	1.0000
JC69IG	57	-3096.6398	6548.7527	0.0000	1.0000
SYM	60	-3089.4370	6553.0562	0.0000	1.0000
JC69I	56	-3109.3417	6567.9200	0.0000	1.0000
F81	58	-3143.5402	6648.7899	0.0000	1.0000
JC69	55	-3174.8851	6692.7705	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/7309_hg18.phylip = 1 - 511;
 Partition Dummy = 1:261-aggressive-phylip/7309_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/7309_hg18.phylip = 1 - 511;
 Partition Dummy = 1:261-aggressive-phylip/7309_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/7309_hg18.phylip = 1 - 511;
 Partition Dummy = 1:261-aggressive-phylip/7309_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

