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

Input data from file "683-tomodel-phylip/chr19_2459.phylip" (7 taxa, 1103 characters)

Minimum AIC  model: HKYG
Minimum AICc model: K2PG
Minimum BIC  model: K2PG


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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-3684.6909	7401.3817	0.2064	0.2064
K2PG	13	-3687.7583	7401.5166	0.1930	0.3994
K2PI	13	-3687.8406	7401.6812	0.1777	0.5771
HKYI	16	-3684.9917	7401.9834	0.1528	0.7299
K2PIG	14	-3687.4548	7402.9096	0.0962	0.8261
HKYIG	17	-3684.7655	7403.5310	0.0705	0.8966
GTRG	20	-3682.2756	7404.5512	0.0423	0.9389
GTRIG	21	-3682.2238	7406.4477	0.0164	0.9553
GTRI	20	-3683.2465	7406.4930	0.0160	0.9713
SYMG	17	-3686.3602	7406.7204	0.0143	0.9856
SYMI	17	-3686.8801	7407.7601	0.0085	0.9941
SYMIG	18	-3686.2511	7408.5021	0.0059	1.0000
GTR	19	-3741.5868	7521.1736	0.0000	1.0000
HKY	15	-3745.6916	7521.3833	0.0000	1.0000
K2P	12	-3749.9211	7523.8422	0.0000	1.0000
SYM	16	-3746.6582	7525.3164	0.0000	1.0000
F81I	15	-3810.7383	7651.4765	0.0000	1.0000
JC69I	12	-3813.7729	7651.5458	0.0000	1.0000
F81IG	16	-3810.7373	7653.4747	0.0000	1.0000
JC69IG	13	-3813.7723	7653.5447	0.0000	1.0000
F81G	15	-3811.8814	7653.7628	0.0000	1.0000
JC69G	12	-3814.9068	7653.8135	0.0000	1.0000
F81	14	-3865.5356	7759.0712	0.0000	1.0000
JC69	11	-3868.7753	7759.5506	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-3687.7583	7466.5919	0.5095	0.5095
K2PI	13	-3687.8406	7466.7564	0.4692	0.9787
K2PIG	14	-3687.4548	7472.9906	0.0208	0.9995
HKYG	16	-3684.6909	7481.4744	0.0003	0.9998
HKYI	16	-3684.9917	7482.0760	0.0002	1.0000
HKYIG	17	-3684.7655	7488.6294	0.0000	1.0000
SYMG	17	-3686.3602	7491.8188	0.0000	1.0000
SYMI	17	-3686.8801	7492.8586	0.0000	1.0000
SYMIG	18	-3686.2511	7498.6063	0.0000	1.0000
GTRG	20	-3682.2756	7504.6670	0.0000	1.0000
GTRI	20	-3683.2465	7506.6088	0.0000	1.0000
GTRIG	21	-3682.2238	7511.5692	0.0000	1.0000
K2P	12	-3749.9211	7583.9117	0.0000	1.0000
HKY	15	-3745.6916	7596.4701	0.0000	1.0000
SYM	16	-3746.6582	7605.4090	0.0000	1.0000
GTR	19	-3741.5868	7616.2836	0.0000	1.0000
JC69I	12	-3813.7729	7711.6152	0.0000	1.0000
JC69G	12	-3814.9068	7713.8830	0.0000	1.0000
JC69IG	13	-3813.7723	7718.6199	0.0000	1.0000
F81I	15	-3810.7383	7726.5634	0.0000	1.0000
F81G	15	-3811.8814	7728.8496	0.0000	1.0000
F81IG	16	-3810.7373	7733.5673	0.0000	1.0000
JC69	11	-3868.7753	7814.6142	0.0000	1.0000
F81	14	-3865.5356	7829.1522	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/chr19_2459.phylip = 1 - 1103;
 Partition Dummy = 1:683-tomodel-phylip/chr19_2459.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 (K2PG)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr19_2459.phylip = 1 - 1103;
 Partition Dummy = 1:683-tomodel-phylip/chr19_2459.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr19_2459.phylip = 1 - 1103;
 Partition Dummy = 1:683-tomodel-phylip/chr19_2459.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

