
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
21:40 04/04/2011

Input data from file "591-tomodel-phylip/chr4_4926.phylip" (7 taxa, 726 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	20	-2164.3645	4369.9206	0.4747	0.4747
GTRI	20	-2164.6725	4370.5365	0.3489	0.8236
GTRIG	21	-2164.3237	4371.9600	0.1712	0.9948
HKYG	16	-2173.9170	4380.6012	0.0023	0.9971
HKYI	16	-2174.0161	4380.7995	0.0021	0.9991
HKYIG	17	-2173.8349	4382.5341	0.0009	1.0000
GTR	19	-2184.5760	4408.2285	0.0000	1.0000
HKY	15	-2192.6766	4416.0293	0.0000	1.0000
F81I	15	-2199.2874	4429.2508	0.0000	1.0000
F81G	15	-2199.5093	4429.6948	0.0000	1.0000
F81IG	16	-2199.5162	4431.7997	0.0000	1.0000
K2PI	13	-2205.9530	4438.4172	0.0000	1.0000
K2PG	13	-2206.1123	4438.7359	0.0000	1.0000
SYMI	17	-2202.3153	4439.4949	0.0000	1.0000
SYMG	17	-2202.5150	4439.8945	0.0000	1.0000
K2PIG	14	-2205.9061	4440.4029	0.0000	1.0000
SYMIG	18	-2202.2670	4441.5016	0.0000	1.0000
F81	14	-2217.5912	4463.7730	0.0000	1.0000
K2P	12	-2225.0081	4474.4537	0.0000	1.0000
SYM	16	-2221.9213	4476.6099	0.0000	1.0000
JC69I	12	-2226.8685	4478.1746	0.0000	1.0000
JC69G	12	-2227.2453	4478.9283	0.0000	1.0000
JC69IG	13	-2226.8713	4480.2538	0.0000	1.0000
JC69	11	-2245.4904	4513.3504	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2164.3645	4368.7291	0.4701	0.4701
GTRI	20	-2164.6725	4369.3450	0.3455	0.8156
GTRIG	21	-2164.3237	4370.6475	0.1802	0.9958
HKYG	16	-2173.9170	4379.8339	0.0018	0.9976
HKYI	16	-2174.0161	4380.0322	0.0017	0.9993
HKYIG	17	-2173.8349	4381.6697	0.0007	1.0000
GTR	19	-2184.5760	4407.1520	0.0000	1.0000
HKY	15	-2192.6766	4415.3533	0.0000	1.0000
F81I	15	-2199.2874	4428.5748	0.0000	1.0000
F81G	15	-2199.5093	4429.0187	0.0000	1.0000
F81IG	16	-2199.5162	4431.0325	0.0000	1.0000
K2PI	13	-2205.9530	4437.9060	0.0000	1.0000
K2PG	13	-2206.1123	4438.2246	0.0000	1.0000
SYMI	17	-2202.3153	4438.6305	0.0000	1.0000
SYMG	17	-2202.5150	4439.0301	0.0000	1.0000
K2PIG	14	-2205.9061	4439.8122	0.0000	1.0000
SYMIG	18	-2202.2670	4440.5341	0.0000	1.0000
F81	14	-2217.5912	4463.1823	0.0000	1.0000
K2P	12	-2225.0081	4474.0162	0.0000	1.0000
SYM	16	-2221.9213	4475.8426	0.0000	1.0000
JC69I	12	-2226.8685	4477.7370	0.0000	1.0000
JC69G	12	-2227.2453	4478.4907	0.0000	1.0000
JC69IG	13	-2226.8713	4479.7426	0.0000	1.0000
JC69	11	-2245.4904	4512.9807	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2173.9170	4453.2347	0.5017	0.5017
HKYI	16	-2174.0161	4453.4330	0.4543	0.9560
HKYIG	17	-2173.8349	4459.6581	0.0202	0.9762
GTRG	20	-2164.3645	4460.4801	0.0134	0.9896
GTRI	20	-2164.6725	4461.0960	0.0098	0.9995
GTRIG	21	-2164.3237	4466.9860	0.0005	1.0000
HKY	15	-2192.6766	4484.1665	0.0000	1.0000
GTR	19	-2184.5760	4494.3155	0.0000	1.0000
F81I	15	-2199.2874	4497.3880	0.0000	1.0000
K2PI	13	-2205.9530	4497.5441	0.0000	1.0000
F81G	15	-2199.5093	4497.8320	0.0000	1.0000
K2PG	13	-2206.1123	4497.8628	0.0000	1.0000
K2PIG	14	-2205.9061	4504.0379	0.0000	1.0000
F81IG	16	-2199.5162	4504.4333	0.0000	1.0000
SYMI	17	-2202.3153	4516.6189	0.0000	1.0000
SYMG	17	-2202.5150	4517.0184	0.0000	1.0000
SYMIG	18	-2202.2670	4523.1100	0.0000	1.0000
F81	14	-2217.5912	4527.4080	0.0000	1.0000
K2P	12	-2225.0081	4529.0668	0.0000	1.0000
JC69I	12	-2226.8685	4532.7876	0.0000	1.0000
JC69G	12	-2227.2453	4533.5413	0.0000	1.0000
JC69IG	13	-2226.8713	4539.3807	0.0000	1.0000
SYM	16	-2221.9213	4549.2434	0.0000	1.0000
JC69	11	-2245.4904	4563.4438	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 591-tomodel-phylip/chr4_4926.phylip = 1 - 726;
 Partition Dummy = 1:591-tomodel-phylip/chr4_4926.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 591-tomodel-phylip/chr4_4926.phylip = 1 - 726;
 Partition Dummy = 1:591-tomodel-phylip/chr4_4926.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 591-tomodel-phylip/chr4_4926.phylip = 1 - 726;
 Partition Dummy = 1:591-tomodel-phylip/chr4_4926.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

