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

Input data from file "591-tomodel-phylip/chr12_1169.phylip" (7 taxa, 801 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	-2678.8480	5390.3899	0.4135	0.4135
HKYI	16	-2678.9522	5390.5983	0.3726	0.7861
HKYIG	17	-2678.8862	5392.5541	0.1401	0.9262
GTRG	20	-2677.1539	5395.3847	0.0340	0.9603
GTRI	20	-2677.3378	5395.7524	0.0283	0.9886
GTRIG	21	-2677.1921	5397.5704	0.0114	1.0000
SYMG	17	-2703.5620	5441.9056	0.0000	1.0000
SYMIG	18	-2703.5834	5444.0415	0.0000	1.0000
K2PG	13	-2708.8493	5444.1611	0.0000	1.0000
SYMI	17	-2704.8362	5444.4540	0.0000	1.0000
K2PI	13	-2709.4648	5445.3922	0.0000	1.0000
K2PIG	14	-2708.8978	5446.3299	0.0000	1.0000
HKY	15	-2709.4956	5449.6026	0.0000	1.0000
GTR	19	-2706.7339	5452.4410	0.0000	1.0000
F81G	15	-2726.6796	5483.9706	0.0000	1.0000
F81I	15	-2727.1177	5484.8470	0.0000	1.0000
F81IG	16	-2726.7481	5486.1901	0.0000	1.0000
SYM	16	-2736.1046	5504.9030	0.0000	1.0000
K2P	12	-2740.5806	5505.5571	0.0000	1.0000
JC69G	12	-2752.9737	5530.3434	0.0000	1.0000
JC69I	12	-2753.7547	5531.9054	0.0000	1.0000
JC69IG	13	-2753.0238	5532.5101	0.0000	1.0000
F81	14	-2756.1569	5540.8482	0.0000	1.0000
JC69	11	-2784.7397	5591.8139	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2678.8480	5389.6960	0.4044	0.4044
HKYI	16	-2678.9522	5389.9044	0.3643	0.7687
HKYIG	17	-2678.8862	5391.7725	0.1432	0.9119
GTRG	20	-2677.1539	5394.3078	0.0403	0.9522
GTRI	20	-2677.3378	5394.6755	0.0335	0.9857
GTRIG	21	-2677.1921	5396.3842	0.0143	1.0000
SYMG	17	-2703.5620	5441.1240	0.0000	1.0000
SYMIG	18	-2703.5834	5443.1669	0.0000	1.0000
SYMI	17	-2704.8362	5443.6724	0.0000	1.0000
K2PG	13	-2708.8493	5443.6986	0.0000	1.0000
K2PI	13	-2709.4648	5444.9296	0.0000	1.0000
K2PIG	14	-2708.8978	5445.7955	0.0000	1.0000
HKY	15	-2709.4956	5448.9911	0.0000	1.0000
GTR	19	-2706.7339	5451.4679	0.0000	1.0000
F81G	15	-2726.6796	5483.3592	0.0000	1.0000
F81I	15	-2727.1177	5484.2355	0.0000	1.0000
F81IG	16	-2726.7481	5485.4962	0.0000	1.0000
SYM	16	-2736.1046	5504.2091	0.0000	1.0000
K2P	12	-2740.5806	5505.1611	0.0000	1.0000
JC69G	12	-2752.9737	5529.9474	0.0000	1.0000
JC69I	12	-2753.7547	5531.5095	0.0000	1.0000
JC69IG	13	-2753.0238	5532.0476	0.0000	1.0000
F81	14	-2756.1569	5540.3139	0.0000	1.0000
JC69	11	-2784.7397	5591.4793	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2678.8480	5464.6698	0.5168	0.5168
HKYI	16	-2678.9522	5464.8782	0.4656	0.9824
HKYIG	17	-2678.8862	5471.4321	0.0176	1.0000
GTRG	20	-2677.1539	5488.0250	0.0000	1.0000
GTRI	20	-2677.3378	5488.3927	0.0000	1.0000
GTRIG	21	-2677.1921	5494.7873	0.0000	1.0000
K2PG	13	-2708.8493	5504.6148	0.0000	1.0000
K2PI	13	-2709.4648	5505.8458	0.0000	1.0000
K2PIG	14	-2708.8978	5511.3976	0.0000	1.0000
HKY	15	-2709.4956	5519.2790	0.0000	1.0000
SYMG	17	-2703.5620	5520.7837	0.0000	1.0000
SYMI	17	-2704.8362	5523.3321	0.0000	1.0000
SYMIG	18	-2703.5834	5527.5124	0.0000	1.0000
GTR	19	-2706.7339	5540.4992	0.0000	1.0000
F81G	15	-2726.6796	5553.6471	0.0000	1.0000
F81I	15	-2727.1177	5554.5234	0.0000	1.0000
F81IG	16	-2726.7481	5560.4700	0.0000	1.0000
K2P	12	-2740.5806	5561.3915	0.0000	1.0000
SYM	16	-2736.1046	5579.1829	0.0000	1.0000
JC69G	12	-2752.9737	5586.1778	0.0000	1.0000
JC69I	12	-2753.7547	5587.7398	0.0000	1.0000
JC69IG	13	-2753.0238	5592.9638	0.0000	1.0000
F81	14	-2756.1569	5605.9159	0.0000	1.0000
JC69	11	-2784.7397	5643.0238	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/chr12_1169.phylip = 1 - 801;
 Partition Dummy = 1:591-tomodel-phylip/chr12_1169.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/chr12_1169.phylip = 1 - 801;
 Partition Dummy = 1:591-tomodel-phylip/chr12_1169.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/chr12_1169.phylip = 1 - 801;
 Partition Dummy = 1:591-tomodel-phylip/chr12_1169.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

