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

Input data from file "591-tomodel-phylip/chr1_32290.phylip" (7 taxa, 694 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	-1854.3961	3741.5956	0.3918	0.3918
HKYI	16	-1854.5164	3741.8363	0.3474	0.7392
HKYIG	17	-1854.3994	3743.7042	0.1365	0.8757
GTRG	20	-1852.0867	3745.4215	0.0579	0.9336
GTRI	20	-1852.3056	3745.8594	0.0465	0.9800
GTRIG	21	-1852.0875	3747.5499	0.0200	1.0000
K2PG	13	-1883.3148	3793.1649	0.0000	1.0000
K2PI	13	-1883.4684	3793.4721	0.0000	1.0000
K2PIG	14	-1883.3148	3795.2481	0.0000	1.0000
HKY	15	-1882.7169	3796.1418	0.0000	1.0000
GTR	19	-1880.2823	3799.6922	0.0000	1.0000
SYMG	17	-1882.9682	3800.8417	0.0000	1.0000
SYMI	17	-1883.1511	3801.2074	0.0000	1.0000
SYMIG	18	-1882.9695	3802.9524	0.0000	1.0000
F81G	15	-1890.3352	3811.3784	0.0000	1.0000
F81I	15	-1890.6512	3812.0103	0.0000	1.0000
F81IG	16	-1890.3312	3813.4660	0.0000	1.0000
K2P	12	-1911.2007	3846.8595	0.0000	1.0000
SYM	16	-1910.0903	3852.9842	0.0000	1.0000
JC69G	12	-1914.8860	3854.2301	0.0000	1.0000
JC69I	12	-1915.2178	3854.8938	0.0000	1.0000
JC69IG	13	-1914.8611	3856.2574	0.0000	1.0000
F81	14	-1917.2415	3863.1015	0.0000	1.0000
JC69	11	-1941.1806	3904.7482	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1854.3961	3740.7921	0.3768	0.3768
HKYI	16	-1854.5164	3741.0327	0.3341	0.7110
HKYIG	17	-1854.3994	3742.7989	0.1382	0.8491
GTRG	20	-1852.0867	3744.1733	0.0695	0.9186
GTRI	20	-1852.3056	3744.6112	0.0558	0.9745
GTRIG	21	-1852.0875	3746.1749	0.0255	1.0000
K2PG	13	-1883.3148	3792.6296	0.0000	1.0000
K2PI	13	-1883.4684	3792.9368	0.0000	1.0000
K2PIG	14	-1883.3148	3794.6296	0.0000	1.0000
HKY	15	-1882.7169	3795.4338	0.0000	1.0000
GTR	19	-1880.2823	3798.5646	0.0000	1.0000
SYMG	17	-1882.9682	3799.9364	0.0000	1.0000
SYMI	17	-1883.1511	3800.3021	0.0000	1.0000
SYMIG	18	-1882.9695	3801.9390	0.0000	1.0000
F81G	15	-1890.3352	3810.6705	0.0000	1.0000
F81I	15	-1890.6512	3811.3024	0.0000	1.0000
F81IG	16	-1890.3312	3812.6625	0.0000	1.0000
K2P	12	-1911.2007	3846.4013	0.0000	1.0000
SYM	16	-1910.0903	3852.1807	0.0000	1.0000
JC69G	12	-1914.8860	3853.7720	0.0000	1.0000
JC69I	12	-1915.2178	3854.4356	0.0000	1.0000
JC69IG	13	-1914.8611	3855.7221	0.0000	1.0000
F81	14	-1917.2415	3862.4829	0.0000	1.0000
JC69	11	-1941.1806	3904.3611	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1854.3961	3813.4717	0.5196	0.5196
HKYI	16	-1854.5164	3813.7123	0.4607	0.9803
HKYIG	17	-1854.3994	3820.0209	0.0197	1.0000
GTRG	20	-1852.0867	3835.0228	0.0000	1.0000
GTRI	20	-1852.3056	3835.4607	0.0000	1.0000
GTRIG	21	-1852.0875	3841.5669	0.0000	1.0000
K2PG	13	-1883.3148	3851.6817	0.0000	1.0000
K2PI	13	-1883.4684	3851.9889	0.0000	1.0000
K2PIG	14	-1883.3148	3858.2242	0.0000	1.0000
HKY	15	-1882.7169	3863.5709	0.0000	1.0000
SYMG	17	-1882.9682	3877.1584	0.0000	1.0000
SYMI	17	-1883.1511	3877.5241	0.0000	1.0000
F81G	15	-1890.3352	3878.8075	0.0000	1.0000
F81I	15	-1890.6512	3879.4394	0.0000	1.0000
SYMIG	18	-1882.9695	3883.7035	0.0000	1.0000
GTR	19	-1880.2823	3884.8716	0.0000	1.0000
F81IG	16	-1890.3312	3885.3421	0.0000	1.0000
K2P	12	-1911.2007	3900.9110	0.0000	1.0000
JC69G	12	-1914.8860	3908.2816	0.0000	1.0000
JC69I	12	-1915.2178	3908.9453	0.0000	1.0000
JC69IG	13	-1914.8611	3914.7743	0.0000	1.0000
SYM	16	-1910.0903	3924.8602	0.0000	1.0000
F81	14	-1917.2415	3926.0775	0.0000	1.0000
JC69	11	-1941.1806	3954.3283	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/chr1_32290.phylip = 1 - 694;
 Partition Dummy = 1:591-tomodel-phylip/chr1_32290.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/chr1_32290.phylip = 1 - 694;
 Partition Dummy = 1:591-tomodel-phylip/chr1_32290.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/chr1_32290.phylip = 1 - 694;
 Partition Dummy = 1:591-tomodel-phylip/chr1_32290.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

