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

Input data from file "591-tomodel-phylip/chr3_5441.phylip" (7 taxa, 744 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	-2242.6673	4526.4965	0.5150	0.5150
GTRIG	21	-2241.7366	4526.7531	0.4530	0.9680
GTRI	20	-2245.4565	4532.0748	0.0317	0.9997
HKYG	16	-2254.9180	4542.5843	0.0002	0.9999
HKYIG	17	-2254.1155	4543.0740	0.0001	1.0000
HKYI	16	-2257.3265	4547.4013	0.0000	1.0000
SYMG	17	-2272.4403	4579.7235	0.0000	1.0000
SYMIG	18	-2271.4957	4579.9348	0.0000	1.0000
K2PG	13	-2278.1012	4582.7010	0.0000	1.0000
K2PIG	14	-2277.2307	4583.0375	0.0000	1.0000
SYMI	17	-2275.2024	4585.2477	0.0000	1.0000
K2PI	13	-2280.7350	4587.9686	0.0000	1.0000
GTR	19	-2290.7678	4620.5853	0.0000	1.0000
F81G	15	-2298.7800	4628.2194	0.0000	1.0000
F81IG	16	-2298.2808	4629.3099	0.0000	1.0000
F81I	15	-2300.7113	4632.0820	0.0000	1.0000
HKY	15	-2301.1107	4632.8808	0.0000	1.0000
JC69G	12	-2317.2441	4658.9151	0.0000	1.0000
JC69IG	13	-2316.7078	4659.9142	0.0000	1.0000
JC69I	12	-2319.3286	4663.0839	0.0000	1.0000
SYM	16	-2317.7185	4668.1853	0.0000	1.0000
K2P	12	-2323.4329	4671.2927	0.0000	1.0000
F81	14	-2341.1761	4710.9283	0.0000	1.0000
JC69	11	-2359.4532	4741.2670	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2242.6673	4525.3347	0.5013	0.5013
GTRIG	21	-2241.7366	4525.4733	0.4677	0.9689
GTRI	20	-2245.4565	4530.9130	0.0308	0.9997
HKYG	16	-2254.9180	4541.8360	0.0001	0.9999
HKYIG	17	-2254.1155	4542.2310	0.0001	1.0000
HKYI	16	-2257.3265	4546.6530	0.0000	1.0000
SYMG	17	-2272.4403	4578.8805	0.0000	1.0000
SYMIG	18	-2271.4957	4578.9914	0.0000	1.0000
K2PG	13	-2278.1012	4582.2024	0.0000	1.0000
K2PIG	14	-2277.2307	4582.4614	0.0000	1.0000
SYMI	17	-2275.2024	4584.4048	0.0000	1.0000
K2PI	13	-2280.7350	4587.4700	0.0000	1.0000
GTR	19	-2290.7678	4619.5356	0.0000	1.0000
F81G	15	-2298.7800	4627.5600	0.0000	1.0000
F81IG	16	-2298.2808	4628.5616	0.0000	1.0000
F81I	15	-2300.7113	4631.4227	0.0000	1.0000
HKY	15	-2301.1107	4632.2215	0.0000	1.0000
JC69G	12	-2317.2441	4658.4883	0.0000	1.0000
JC69IG	13	-2316.7078	4659.4156	0.0000	1.0000
JC69I	12	-2319.3286	4662.6571	0.0000	1.0000
SYM	16	-2317.7185	4667.4370	0.0000	1.0000
K2P	12	-2323.4329	4670.8658	0.0000	1.0000
F81	14	-2341.1761	4710.3521	0.0000	1.0000
JC69	11	-2359.4532	4740.9063	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2254.9180	4615.6287	0.6219	0.6219
GTRG	20	-2242.6673	4617.5755	0.2350	0.8569
HKYI	16	-2257.3265	4620.4457	0.0559	0.9128
HKYIG	17	-2254.1155	4620.6357	0.0509	0.9637
GTRIG	21	-2241.7366	4622.3262	0.0218	0.9856
GTRI	20	-2245.4565	4623.1538	0.0144	1.0000
K2PG	13	-2278.1012	4642.1589	0.0000	1.0000
K2PIG	14	-2277.2307	4647.0299	0.0000	1.0000
K2PI	13	-2280.7350	4647.4265	0.0000	1.0000
SYMG	17	-2272.4403	4657.2852	0.0000	1.0000
SYMIG	18	-2271.4957	4662.0081	0.0000	1.0000
SYMI	17	-2275.2024	4662.8095	0.0000	1.0000
F81G	15	-2298.7800	4696.7407	0.0000	1.0000
F81I	15	-2300.7113	4700.6033	0.0000	1.0000
HKY	15	-2301.1107	4701.4021	0.0000	1.0000
F81IG	16	-2298.2808	4702.3543	0.0000	1.0000
GTR	19	-2290.7678	4707.1644	0.0000	1.0000
JC69G	12	-2317.2441	4713.8328	0.0000	1.0000
JC69I	12	-2319.3286	4718.0016	0.0000	1.0000
JC69IG	13	-2316.7078	4719.3721	0.0000	1.0000
K2P	12	-2323.4329	4726.2103	0.0000	1.0000
SYM	16	-2317.7185	4741.2297	0.0000	1.0000
F81	14	-2341.1761	4774.9207	0.0000	1.0000
JC69	11	-2359.4532	4791.6388	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/chr3_5441.phylip = 1 - 744;
 Partition Dummy = 1:591-tomodel-phylip/chr3_5441.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/chr3_5441.phylip = 1 - 744;
 Partition Dummy = 1:591-tomodel-phylip/chr3_5441.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/chr3_5441.phylip = 1 - 744;
 Partition Dummy = 1:591-tomodel-phylip/chr3_5441.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

