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

Input data from file "591-tomodel-phylip/chr3_2614.phylip" (7 taxa, 686 characters)

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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
HKYI	16	-2262.7820	4558.3771	0.5390	0.5390
HKYG	16	-2263.7924	4560.3980	0.1962	0.7352
HKYIG	17	-2262.7685	4560.4532	0.1909	0.9261
GTRI	20	-2261.0662	4563.3956	0.0438	0.9699
GTRIG	21	-2261.0168	4565.4251	0.0159	0.9858
GTRG	20	-2262.1947	4565.6525	0.0142	1.0000
K2PI	13	-2301.5221	4629.5859	0.0000	1.0000
HKY	15	-2300.0823	4630.8809	0.0000	1.0000
K2PIG	14	-2301.5358	4631.6976	0.0000	1.0000
K2PG	13	-2303.0260	4632.5936	0.0000	1.0000
SYMI	17	-2299.2102	4633.3366	0.0000	1.0000
SYMIG	18	-2299.2243	4635.4740	0.0000	1.0000
SYMG	17	-2300.7085	4636.3332	0.0000	1.0000
GTR	19	-2298.8271	4636.7954	0.0000	1.0000
F81I	15	-2319.8524	4670.4213	0.0000	1.0000
F81IG	16	-2319.8729	4672.5590	0.0000	1.0000
F81G	15	-2321.5809	4673.8781	0.0000	1.0000
K2P	12	-2339.4793	4703.4221	0.0000	1.0000
SYM	16	-2337.5401	4707.8935	0.0000	1.0000
JC69I	12	-2351.5475	4727.5586	0.0000	1.0000
JC69IG	13	-2351.5756	4729.6929	0.0000	1.0000
JC69G	12	-2353.4862	4731.4360	0.0000	1.0000
F81	14	-2354.1378	4736.9015	0.0000	1.0000
JC69	11	-2386.2471	4794.8858	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-2262.7820	4557.5640	0.5233	0.5233
HKYIG	17	-2262.7685	4559.5371	0.1951	0.7184
HKYG	16	-2263.7924	4559.5848	0.1905	0.9089
GTRI	20	-2261.0662	4562.1325	0.0533	0.9622
GTRIG	21	-2261.0168	4564.0335	0.0206	0.9828
GTRG	20	-2262.1947	4564.3893	0.0172	1.0000
K2PI	13	-2301.5221	4629.0442	0.0000	1.0000
HKY	15	-2300.0823	4630.1645	0.0000	1.0000
K2PIG	14	-2301.5358	4631.0717	0.0000	1.0000
K2PG	13	-2303.0260	4632.0519	0.0000	1.0000
SYMI	17	-2299.2102	4632.4205	0.0000	1.0000
SYMIG	18	-2299.2243	4634.4485	0.0000	1.0000
SYMG	17	-2300.7085	4635.4170	0.0000	1.0000
GTR	19	-2298.8271	4635.6542	0.0000	1.0000
F81I	15	-2319.8524	4669.7049	0.0000	1.0000
F81IG	16	-2319.8729	4671.7459	0.0000	1.0000
F81G	15	-2321.5809	4673.1617	0.0000	1.0000
K2P	12	-2339.4793	4702.9585	0.0000	1.0000
SYM	16	-2337.5401	4707.0803	0.0000	1.0000
JC69I	12	-2351.5475	4727.0950	0.0000	1.0000
JC69IG	13	-2351.5756	4729.1512	0.0000	1.0000
JC69G	12	-2353.4862	4730.9724	0.0000	1.0000
F81	14	-2354.1378	4736.2755	0.0000	1.0000
JC69	11	-2386.2471	4794.4941	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2262.7820	4630.0580	0.7129	0.7129
HKYG	16	-2263.7924	4632.0788	0.2595	0.9724
HKYIG	17	-2262.7685	4636.5620	0.0276	1.0000
GTRI	20	-2261.0662	4652.7500	0.0000	1.0000
GTRG	20	-2262.1947	4655.0069	0.0000	1.0000
GTRIG	21	-2261.0168	4659.1820	0.0000	1.0000
K2PI	13	-2301.5221	4687.9456	0.0000	1.0000
K2PG	13	-2303.0260	4690.9533	0.0000	1.0000
K2PIG	14	-2301.5358	4694.5039	0.0000	1.0000
HKY	15	-2300.0823	4698.1277	0.0000	1.0000
SYMI	17	-2299.2102	4709.4454	0.0000	1.0000
SYMG	17	-2300.7085	4712.4420	0.0000	1.0000
SYMIG	18	-2299.2243	4716.0043	0.0000	1.0000
GTR	19	-2298.8271	4721.7409	0.0000	1.0000
F81I	15	-2319.8524	4737.6680	0.0000	1.0000
F81G	15	-2321.5809	4741.1249	0.0000	1.0000
F81IG	16	-2319.8729	4744.2399	0.0000	1.0000
K2P	12	-2339.4793	4757.3291	0.0000	1.0000
SYM	16	-2337.5401	4779.5743	0.0000	1.0000
JC69I	12	-2351.5475	4781.4656	0.0000	1.0000
JC69G	12	-2353.4862	4785.3429	0.0000	1.0000
JC69IG	13	-2351.5756	4788.0526	0.0000	1.0000
F81	14	-2354.1378	4799.7078	0.0000	1.0000
JC69	11	-2386.2471	4844.3338	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_2614.phylip = 1 - 686;
 Partition Dummy = 1:591-tomodel-phylip/chr3_2614.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_2614.phylip = 1 - 686;
 Partition Dummy = 1:591-tomodel-phylip/chr3_2614.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_2614.phylip = 1 - 686;
 Partition Dummy = 1:591-tomodel-phylip/chr3_2614.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


-------------------------------------------------------------
End of Output

