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

Input data from file "261-aggressive-phylip/7207_hg18.phylip" (29 taxa, 729 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	-2139.9964	4420.5230	0.5835	0.5835
GTRIG	65	-2139.8725	4422.6862	0.1978	0.7814
GTRI	64	-2141.5268	4423.5838	0.1263	0.9077
HKYG	60	-2147.1075	4425.1732	0.0571	0.9647
HKYIG	61	-2146.9560	4427.2523	0.0202	0.9849
HKYI	60	-2148.4358	4427.8296	0.0151	1.0000
GTR	63	-2158.6788	4455.4840	0.0000	1.0000
HKY	59	-2165.8961	4460.3752	0.0000	1.0000
K2PG	57	-2177.8663	4479.5866	0.0000	1.0000
K2PIG	58	-2177.6840	4481.5830	0.0000	1.0000
K2PI	57	-2179.0791	4482.0121	0.0000	1.0000
SYMG	61	-2174.5766	4482.4934	0.0000	1.0000
SYMIG	62	-2174.3899	4484.5095	0.0000	1.0000
SYMI	61	-2175.8391	4485.0186	0.0000	1.0000
K2P	56	-2197.9005	4517.3010	0.0000	1.0000
SYM	60	-2194.1406	4519.2393	0.0000	1.0000
F81G	59	-2197.2325	4523.0480	0.0000	1.0000
F81IG	60	-2197.0614	4525.0809	0.0000	1.0000
F81I	59	-2198.4871	4525.5572	0.0000	1.0000
F81	58	-2216.1988	4558.6124	0.0000	1.0000
JC69G	56	-2225.1816	4571.8633	0.0000	1.0000
JC69IG	57	-2224.9945	4573.8429	0.0000	1.0000
JC69I	56	-2226.4096	4574.3193	0.0000	1.0000
JC69	55	-2244.7279	4608.6089	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2139.9964	4407.9928	0.5854	0.5854
GTRIG	65	-2139.8725	4409.7450	0.2438	0.8291
GTRI	64	-2141.5268	4411.0537	0.1267	0.9558
HKYG	60	-2147.1075	4414.2151	0.0261	0.9819
HKYIG	61	-2146.9560	4415.9120	0.0112	0.9931
HKYI	60	-2148.4358	4416.8715	0.0069	1.0000
GTR	63	-2158.6788	4443.3577	0.0000	1.0000
HKY	59	-2165.8961	4449.7923	0.0000	1.0000
K2PG	57	-2177.8663	4469.7326	0.0000	1.0000
SYMG	61	-2174.5766	4471.1531	0.0000	1.0000
K2PIG	58	-2177.6840	4471.3681	0.0000	1.0000
K2PI	57	-2179.0791	4472.1582	0.0000	1.0000
SYMIG	62	-2174.3899	4472.7797	0.0000	1.0000
SYMI	61	-2175.8391	4473.6783	0.0000	1.0000
K2P	56	-2197.9005	4507.8010	0.0000	1.0000
SYM	60	-2194.1406	4508.2812	0.0000	1.0000
F81G	59	-2197.2325	4512.4650	0.0000	1.0000
F81IG	60	-2197.0614	4514.1228	0.0000	1.0000
F81I	59	-2198.4871	4514.9742	0.0000	1.0000
F81	58	-2216.1988	4548.3975	0.0000	1.0000
JC69G	56	-2225.1816	4562.3633	0.0000	1.0000
JC69IG	57	-2224.9945	4563.9890	0.0000	1.0000
JC69I	56	-2226.4096	4564.8193	0.0000	1.0000
JC69	55	-2244.7279	4599.4558	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2147.1075	4689.7155	0.7628	0.7628
HKYI	60	-2148.4358	4692.3719	0.2021	0.9649
HKYIG	61	-2146.9560	4696.0041	0.0329	0.9978
GTRG	64	-2139.9964	4701.8600	0.0018	0.9995
GTRI	64	-2141.5268	4704.9208	0.0004	0.9999
GTRIG	65	-2139.8725	4708.2038	0.0001	1.0000
HKY	59	-2165.8961	4720.7010	0.0000	1.0000
K2PG	57	-2177.8663	4731.4580	0.0000	1.0000
GTR	63	-2158.6788	4732.6331	0.0000	1.0000
K2PI	57	-2179.0791	4733.8836	0.0000	1.0000
K2PIG	58	-2177.6840	4737.6852	0.0000	1.0000
SYMG	61	-2174.5766	4751.2452	0.0000	1.0000
SYMI	61	-2175.8391	4753.7704	0.0000	1.0000
SYMIG	62	-2174.3899	4757.4635	0.0000	1.0000
K2P	56	-2197.9005	4764.9347	0.0000	1.0000
F81G	59	-2197.2325	4783.3738	0.0000	1.0000
SYM	60	-2194.1406	4783.7816	0.0000	1.0000
F81I	59	-2198.4871	4785.8830	0.0000	1.0000
F81IG	60	-2197.0614	4789.6233	0.0000	1.0000
F81	58	-2216.1988	4814.7146	0.0000	1.0000
JC69G	56	-2225.1816	4819.4970	0.0000	1.0000
JC69I	56	-2226.4096	4821.9530	0.0000	1.0000
JC69IG	57	-2224.9945	4825.7144	0.0000	1.0000
JC69	55	-2244.7279	4851.9979	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/7207_hg18.phylip = 1 - 729;
 Partition Dummy = 1:261-aggressive-phylip/7207_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/7207_hg18.phylip = 1 - 729;
 Partition Dummy = 1:261-aggressive-phylip/7207_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/7207_hg18.phylip = 1 - 729;
 Partition Dummy = 1:261-aggressive-phylip/7207_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

