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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	16	-2356.2147	4745.3184	0.4645	0.4645
HKYIG	17	-2356.2153	4747.4322	0.1614	0.6259
HKYI	16	-2357.3735	4747.6360	0.1458	0.7717
GTRG	20	-2353.1600	4747.7015	0.1411	0.9127
GTRIG	21	-2353.1600	4749.8422	0.0484	0.9611
GTRI	20	-2354.4655	4750.3126	0.0382	0.9994
HKY	15	-2364.0682	4758.9194	0.0005	0.9999
GTR	19	-2361.2250	4761.6979	0.0001	1.0000
F81G	15	-2388.0873	4806.9577	0.0000	1.0000
F81I	15	-2389.1462	4809.0753	0.0000	1.0000
F81IG	16	-2388.0977	4809.0843	0.0000	1.0000
F81	14	-2395.6181	4819.9202	0.0000	1.0000
K2PG	13	-2415.4399	4857.4717	0.0000	1.0000
K2PIG	14	-2415.4399	4859.5638	0.0000	1.0000
K2PI	13	-2416.5920	4859.7759	0.0000	1.0000
SYMG	17	-2412.8870	4860.7756	0.0000	1.0000
SYMIG	18	-2412.8870	4862.8953	0.0000	1.0000
SYMI	17	-2414.1947	4863.3910	0.0000	1.0000
K2P	12	-2424.0611	4872.6287	0.0000	1.0000
SYM	16	-2422.0786	4877.0461	0.0000	1.0000
JC69G	12	-2437.7589	4900.0244	0.0000	1.0000
JC69I	12	-2438.7949	4902.0964	0.0000	1.0000
JC69IG	13	-2437.7589	4902.1097	0.0000	1.0000
JC69	11	-2445.7557	4913.9393	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2356.2147	4744.4295	0.4311	0.4311
GTRG	20	-2353.1600	4746.3199	0.1675	0.5986
HKYIG	17	-2356.2153	4746.4306	0.1585	0.7571
HKYI	16	-2357.3735	4746.7471	0.1353	0.8924
GTRIG	21	-2353.1600	4748.3199	0.0616	0.9540
GTRI	20	-2354.4655	4748.9310	0.0454	0.9994
HKY	15	-2364.0682	4758.1364	0.0005	0.9999
GTR	19	-2361.2250	4760.4500	0.0001	1.0000
F81G	15	-2388.0873	4806.1746	0.0000	1.0000
F81IG	16	-2388.0977	4808.1954	0.0000	1.0000
F81I	15	-2389.1462	4808.2923	0.0000	1.0000
F81	14	-2395.6181	4819.2361	0.0000	1.0000
K2PG	13	-2415.4399	4856.8798	0.0000	1.0000
K2PIG	14	-2415.4399	4858.8798	0.0000	1.0000
K2PI	13	-2416.5920	4859.1841	0.0000	1.0000
SYMG	17	-2412.8870	4859.7740	0.0000	1.0000
SYMIG	18	-2412.8870	4861.7740	0.0000	1.0000
SYMI	17	-2414.1947	4862.3894	0.0000	1.0000
K2P	12	-2424.0611	4872.1222	0.0000	1.0000
SYM	16	-2422.0786	4876.1572	0.0000	1.0000
JC69G	12	-2437.7589	4899.5179	0.0000	1.0000
JC69IG	13	-2437.7589	4901.5179	0.0000	1.0000
JC69I	12	-2438.7949	4901.5899	0.0000	1.0000
JC69	11	-2445.7557	4913.5114	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2356.2147	4815.5356	0.7334	0.7334
HKYI	16	-2357.3735	4817.8532	0.2302	0.9636
HKYIG	17	-2356.2153	4821.9808	0.0292	0.9928
HKY	15	-2364.0682	4824.7984	0.0071	0.9999
GTRG	20	-2353.1600	4835.2025	0.0000	1.0000
GTRI	20	-2354.4655	4837.8136	0.0000	1.0000
GTRIG	21	-2353.1600	4841.6467	0.0000	1.0000
GTR	19	-2361.2250	4844.8885	0.0000	1.0000
F81G	15	-2388.0873	4872.8366	0.0000	1.0000
F81I	15	-2389.1462	4874.9543	0.0000	1.0000
F81IG	16	-2388.0977	4879.3015	0.0000	1.0000
F81	14	-2395.6181	4881.4540	0.0000	1.0000
K2PG	13	-2415.4399	4914.6535	0.0000	1.0000
K2PI	13	-2416.5920	4916.9578	0.0000	1.0000
K2PIG	14	-2415.4399	4921.0976	0.0000	1.0000
K2P	12	-2424.0611	4925.4518	0.0000	1.0000
SYMG	17	-2412.8870	4935.3242	0.0000	1.0000
SYMI	17	-2414.1947	4937.9396	0.0000	1.0000
SYMIG	18	-2412.8870	4941.7684	0.0000	1.0000
SYM	16	-2422.0786	4947.2633	0.0000	1.0000
JC69G	12	-2437.7589	4952.8475	0.0000	1.0000
JC69I	12	-2438.7949	4954.9195	0.0000	1.0000
JC69IG	13	-2437.7589	4959.2916	0.0000	1.0000
JC69	11	-2445.7557	4962.3969	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr13_714.phylip = 1 - 629;
 Partition Dummy = 1:591-tomodel-phylip/chr13_714.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;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr13_714.phylip = 1 - 629;
 Partition Dummy = 1:591-tomodel-phylip/chr13_714.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;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr13_714.phylip = 1 - 629;
 Partition Dummy = 1:591-tomodel-phylip/chr13_714.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

