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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: HKYG


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2557.0482	5154.0965	0.4518	0.4518
GTRIG	21	-2556.0617	5154.1233	0.4458	0.8977
GTRI	20	-2558.6229	5157.2459	0.0936	0.9912
HKYIG	17	-2564.7950	5163.5901	0.0039	0.9951
HKYG	16	-2565.8069	5163.6138	0.0039	0.9990
HKYI	16	-2567.1745	5166.3490	0.0010	1.0000
GTR	19	-2581.5830	5201.1660	0.0000	1.0000
HKY	15	-2590.5083	5211.0166	0.0000	1.0000
F81IG	16	-2601.5281	5235.0561	0.0000	1.0000
F81G	15	-2602.6945	5235.3890	0.0000	1.0000
F81I	15	-2604.1637	5238.3273	0.0000	1.0000
K2PIG	14	-2623.2477	5274.4955	0.0000	1.0000
K2PG	13	-2624.3580	5274.7161	0.0000	1.0000
SYMIG	18	-2619.6602	5275.3204	0.0000	1.0000
SYMG	17	-2620.7478	5275.4956	0.0000	1.0000
K2PI	13	-2625.5408	5277.0816	0.0000	1.0000
SYMI	17	-2622.0011	5278.0022	0.0000	1.0000
F81	14	-2627.8240	5283.6479	0.0000	1.0000
SYM	16	-2645.9698	5323.9396	0.0000	1.0000
K2P	12	-2649.9892	5323.9783	0.0000	1.0000
JC69IG	13	-2654.5965	5335.1929	0.0000	1.0000
JC69G	12	-2655.8227	5335.6454	0.0000	1.0000
JC69I	12	-2657.0648	5338.1296	0.0000	1.0000
JC69	11	-2681.5704	5385.1408	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2565.8069	5243.5900	0.7437	0.7437
HKYI	16	-2567.1745	5246.3252	0.1894	0.9331
HKYIG	17	-2564.7950	5248.5647	0.0618	0.9949
GTRG	20	-2557.0482	5254.0667	0.0039	0.9989
GTRI	20	-2558.6229	5257.2161	0.0008	0.9997
GTRIG	21	-2556.0617	5259.0920	0.0003	1.0000
HKY	15	-2590.5083	5285.9942	0.0000	1.0000
GTR	19	-2581.5830	5296.1377	0.0000	1.0000
F81G	15	-2602.6945	5310.3666	0.0000	1.0000
F81I	15	-2604.1637	5313.3050	0.0000	1.0000
F81IG	16	-2601.5281	5315.0323	0.0000	1.0000
K2PG	13	-2624.3580	5339.6967	0.0000	1.0000
K2PI	13	-2625.5408	5342.0622	0.0000	1.0000
K2PIG	14	-2623.2477	5344.4746	0.0000	1.0000
F81	14	-2627.8240	5353.6271	0.0000	1.0000
SYMG	17	-2620.7478	5360.4702	0.0000	1.0000
SYMI	17	-2622.0011	5362.9769	0.0000	1.0000
SYMIG	18	-2619.6602	5365.2936	0.0000	1.0000
K2P	12	-2649.9892	5383.9605	0.0000	1.0000
JC69G	12	-2655.8227	5395.6275	0.0000	1.0000
JC69I	12	-2657.0648	5398.1117	0.0000	1.0000
JC69IG	13	-2654.5965	5400.1735	0.0000	1.0000
SYM	16	-2645.9698	5403.9158	0.0000	1.0000
JC69	11	-2681.5704	5440.1244	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/chr5_3418.phylip = 1 - 1095;
 Partition Dummy = 1:591-tomodel-phylip/chr5_3418.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/chr5_3418.phylip = 1 - 1095;
 Partition Dummy = 1:591-tomodel-phylip/chr5_3418.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/chr5_3418.phylip = 1 - 1095;
 Partition Dummy = 1:591-tomodel-phylip/chr5_3418.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

