
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:11 04/04/2011

Input data from file "683-tomodel-phylip/chr1_18908.phylip" (7 taxa, 747 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	-2755.3730	5543.4911	0.3983	0.3983
HKYI	16	-2755.3831	5543.5114	0.3943	0.7926
HKYIG	17	-2755.1585	5545.1565	0.1732	0.9658
GTRI	20	-2754.4824	5550.1219	0.0145	0.9803
GTRG	20	-2754.5579	5550.2728	0.0134	0.9937
GTRIG	21	-2754.2910	5551.8566	0.0061	0.9998
HKY	15	-2763.8759	5558.4085	0.0002	1.0000
GTR	19	-2762.8733	5564.7920	0.0000	1.0000
K2PG	13	-2807.9954	5642.4874	0.0000	1.0000
K2PI	13	-2808.2843	5643.0652	0.0000	1.0000
K2PIG	14	-2807.9361	5644.4459	0.0000	1.0000
SYMG	17	-2805.6231	5646.0856	0.0000	1.0000
SYMI	17	-2806.1005	5647.0405	0.0000	1.0000
SYMIG	18	-2805.6059	5648.1514	0.0000	1.0000
K2P	12	-2818.2231	5660.8713	0.0000	1.0000
SYM	16	-2816.6899	5666.1250	0.0000	1.0000
F81I	15	-2844.3644	5719.3854	0.0000	1.0000
F81G	15	-2844.5515	5719.7596	0.0000	1.0000
F81IG	16	-2844.3070	5721.3593	0.0000	1.0000
F81	14	-2851.6194	5731.8125	0.0000	1.0000
JC69I	12	-2885.8823	5796.1897	0.0000	1.0000
JC69G	12	-2885.9300	5796.2850	0.0000	1.0000
JC69IG	13	-2885.7899	5798.0765	0.0000	1.0000
JC69	11	-2893.9137	5810.1866	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2755.3730	5542.7459	0.3918	0.3918
HKYI	16	-2755.3831	5542.7662	0.3879	0.7797
HKYIG	17	-2755.1585	5544.3170	0.1786	0.9583
GTRI	20	-2754.4824	5548.9649	0.0175	0.9758
GTRG	20	-2754.5579	5549.1158	0.0162	0.9920
GTRIG	21	-2754.2910	5550.5821	0.0078	0.9998
HKY	15	-2763.8759	5557.7518	0.0002	1.0000
GTR	19	-2762.8733	5563.7466	0.0000	1.0000
K2PG	13	-2807.9954	5641.9908	0.0000	1.0000
K2PI	13	-2808.2843	5642.5686	0.0000	1.0000
K2PIG	14	-2807.9361	5643.8722	0.0000	1.0000
SYMG	17	-2805.6231	5645.2461	0.0000	1.0000
SYMI	17	-2806.1005	5646.2010	0.0000	1.0000
SYMIG	18	-2805.6059	5647.2118	0.0000	1.0000
K2P	12	-2818.2231	5660.4462	0.0000	1.0000
SYM	16	-2816.6899	5665.3798	0.0000	1.0000
F81I	15	-2844.3644	5718.7287	0.0000	1.0000
F81G	15	-2844.5515	5719.1030	0.0000	1.0000
F81IG	16	-2844.3070	5720.6141	0.0000	1.0000
F81	14	-2851.6194	5731.2387	0.0000	1.0000
JC69I	12	-2885.8823	5795.7647	0.0000	1.0000
JC69G	12	-2885.9300	5795.8600	0.0000	1.0000
JC69IG	13	-2885.7899	5797.5799	0.0000	1.0000
JC69	11	-2893.9137	5809.8274	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2755.3730	5616.6030	0.4900	0.4900
HKYI	16	-2755.3831	5616.6232	0.4851	0.9751
HKYIG	17	-2755.1585	5622.7901	0.0222	0.9973
HKY	15	-2763.8759	5626.9928	0.0027	1.0000
GTRI	20	-2754.4824	5641.2862	0.0000	1.0000
GTRG	20	-2754.5579	5641.4371	0.0000	1.0000
GTRIG	21	-2754.2910	5647.5194	0.0000	1.0000
GTR	19	-2762.8733	5651.4519	0.0000	1.0000
K2PG	13	-2807.9954	5701.9996	0.0000	1.0000
K2PI	13	-2808.2843	5702.5775	0.0000	1.0000
K2PIG	14	-2807.9361	5708.4971	0.0000	1.0000
K2P	12	-2818.2231	5715.8390	0.0000	1.0000
SYMG	17	-2805.6231	5723.7192	0.0000	1.0000
SYMI	17	-2806.1005	5724.6741	0.0000	1.0000
SYMIG	18	-2805.6059	5730.3010	0.0000	1.0000
SYM	16	-2816.6899	5739.2368	0.0000	1.0000
F81I	15	-2844.3644	5787.9697	0.0000	1.0000
F81G	15	-2844.5515	5788.3440	0.0000	1.0000
F81IG	16	-2844.3070	5794.4711	0.0000	1.0000
F81	14	-2851.6194	5795.8637	0.0000	1.0000
JC69I	12	-2885.8823	5851.1574	0.0000	1.0000
JC69G	12	-2885.9300	5851.2528	0.0000	1.0000
JC69IG	13	-2885.7899	5857.5887	0.0000	1.0000
JC69	11	-2893.9137	5860.6041	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 683-tomodel-phylip/chr1_18908.phylip = 1 - 747;
 Partition Dummy = 1:683-tomodel-phylip/chr1_18908.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 683-tomodel-phylip/chr1_18908.phylip = 1 - 747;
 Partition Dummy = 1:683-tomodel-phylip/chr1_18908.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 683-tomodel-phylip/chr1_18908.phylip = 1 - 747;
 Partition Dummy = 1:683-tomodel-phylip/chr1_18908.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

