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

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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
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
GTRIG	65	-2436.9336	5014.6869	0.7904	0.7904
GTRG	64	-2439.5366	5017.5518	0.1887	0.9790
GTRI	64	-2441.7570	5021.9926	0.0205	0.9995
HKYIG	61	-2449.1736	5029.8379	0.0004	0.9999
HKYG	60	-2451.9856	5033.1442	0.0001	1.0000
HKYI	60	-2453.8202	5036.8133	0.0000	1.0000
SYMIG	62	-2469.8467	5073.5074	0.0000	1.0000
K2PIG	58	-2474.5875	5073.7300	0.0000	1.0000
SYMG	61	-2472.2024	5075.8953	0.0000	1.0000
K2PG	57	-2476.9340	5076.1226	0.0000	1.0000
K2PI	57	-2478.0152	5078.2851	0.0000	1.0000
SYMI	61	-2473.4087	5078.3080	0.0000	1.0000
GTR	63	-2487.5477	5111.2388	0.0000	1.0000
HKY	59	-2500.8986	5128.6582	0.0000	1.0000
SYM	60	-2521.7746	5172.7222	0.0000	1.0000
F81IG	60	-2521.8525	5172.8779	0.0000	1.0000
F81G	59	-2523.9838	5174.8286	0.0000	1.0000
K2P	56	-2527.7038	5175.3677	0.0000	1.0000
F81I	59	-2525.2528	5177.3666	0.0000	1.0000
JC69IG	57	-2544.5043	5211.2632	0.0000	1.0000
JC69G	56	-2546.5941	5213.1484	0.0000	1.0000
JC69I	56	-2547.6665	5215.2932	0.0000	1.0000
F81	58	-2572.1356	5268.8262	0.0000	1.0000
JC69	55	-2595.9624	5309.5960	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-2436.9336	5003.8672	0.8174	0.8174
GTRG	64	-2439.5366	5007.0733	0.1645	0.9819
GTRI	64	-2441.7570	5011.5140	0.0179	0.9997
HKYIG	61	-2449.1736	5020.3473	0.0002	1.0000
HKYG	60	-2451.9856	5023.9712	0.0000	1.0000
HKYI	60	-2453.8202	5027.6403	0.0000	1.0000
SYMIG	62	-2469.8467	5063.6933	0.0000	1.0000
K2PIG	58	-2474.5875	5065.1750	0.0000	1.0000
SYMG	61	-2472.2024	5066.4048	0.0000	1.0000
K2PG	57	-2476.9340	5067.8679	0.0000	1.0000
SYMI	61	-2473.4087	5068.8174	0.0000	1.0000
K2PI	57	-2478.0152	5070.0304	0.0000	1.0000
GTR	63	-2487.5477	5101.0954	0.0000	1.0000
HKY	59	-2500.8986	5119.7972	0.0000	1.0000
SYM	60	-2521.7746	5163.5493	0.0000	1.0000
F81IG	60	-2521.8525	5163.7049	0.0000	1.0000
F81G	59	-2523.9838	5165.9675	0.0000	1.0000
K2P	56	-2527.7038	5167.4076	0.0000	1.0000
F81I	59	-2525.2528	5168.5056	0.0000	1.0000
JC69IG	57	-2544.5043	5203.0085	0.0000	1.0000
JC69G	56	-2546.5941	5205.1883	0.0000	1.0000
JC69I	56	-2547.6665	5207.3331	0.0000	1.0000
F81	58	-2572.1356	5260.2712	0.0000	1.0000
JC69	55	-2595.9624	5301.9248	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2451.9856	5309.3174	0.4406	0.4406
HKYIG	61	-2449.1736	5310.4492	0.2502	0.6907
GTRG	64	-2439.5366	5311.4425	0.1522	0.8430
HKYI	60	-2453.8202	5312.9865	0.0704	0.9133
GTRIG	65	-2436.9336	5312.9922	0.0701	0.9835
GTRI	64	-2441.7570	5315.8832	0.0165	1.0000
K2PG	57	-2476.9340	5338.9467	0.0000	1.0000
K2PIG	58	-2474.5875	5341.0096	0.0000	1.0000
K2PI	57	-2478.0152	5341.1092	0.0000	1.0000
SYMG	61	-2472.2024	5356.5067	0.0000	1.0000
SYMIG	62	-2469.8467	5358.5510	0.0000	1.0000
SYMI	61	-2473.4087	5358.9193	0.0000	1.0000
HKY	59	-2500.8986	5400.3875	0.0000	1.0000
GTR	63	-2487.5477	5400.7088	0.0000	1.0000
K2P	56	-2527.7038	5433.7307	0.0000	1.0000
F81G	59	-2523.9838	5446.5579	0.0000	1.0000
SYM	60	-2521.7746	5448.8954	0.0000	1.0000
F81IG	60	-2521.8525	5449.0511	0.0000	1.0000
F81I	59	-2525.2528	5449.0959	0.0000	1.0000
JC69G	56	-2546.5941	5471.5113	0.0000	1.0000
JC69I	56	-2547.6665	5473.6562	0.0000	1.0000
JC69IG	57	-2544.5043	5474.0874	0.0000	1.0000
F81	58	-2572.1356	5536.1058	0.0000	1.0000
JC69	55	-2595.9624	5563.4921	0.0000	1.0000



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

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

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


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

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


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

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

