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

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

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


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	-2051.3859	4135.5159	0.3885	0.3885
HKYI	16	-2051.9374	4136.6189	0.2238	0.6123
HKYIG	17	-2051.0670	4136.9723	0.1876	0.7998
GTRG	20	-2048.5846	4138.3247	0.0954	0.8952
GTRIG	21	-2047.9392	4139.1511	0.0631	0.9583
GTRI	20	-2049.6088	4140.3730	0.0343	0.9926
K2PG	13	-2059.3887	4145.2733	0.0030	0.9955
K2PI	13	-2059.7148	4145.9255	0.0021	0.9976
K2PIG	14	-2059.1579	4146.8888	0.0013	0.9990
SYMG	17	-2057.0060	4148.8504	0.0005	0.9995
SYMIG	18	-2056.3987	4149.7356	0.0003	0.9998
SYMI	17	-2057.7986	4150.4355	0.0002	1.0000
GTR	19	-2088.2381	4215.5202	0.0000	1.0000
HKY	15	-2096.5437	4223.7431	0.0000	1.0000
SYM	16	-2099.1242	4230.9926	0.0000	1.0000
F81G	15	-2102.8633	4236.3824	0.0000	1.0000
K2P	12	-2106.3827	4237.1898	0.0000	1.0000
F81IG	16	-2102.3465	4237.4372	0.0000	1.0000
F81I	15	-2103.9084	4238.4725	0.0000	1.0000
JC69G	12	-2112.1331	4248.6907	0.0000	1.0000
JC69IG	13	-2111.5826	4249.6612	0.0000	1.0000
JC69I	12	-2113.1611	4250.7466	0.0000	1.0000
F81	14	-2143.5617	4315.6964	0.0000	1.0000
JC69	11	-2153.0368	4328.4323	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2051.3859	4134.7718	0.3675	0.3675
HKYI	16	-2051.9374	4135.8747	0.2117	0.5792
HKYIG	17	-2051.0670	4136.1339	0.1860	0.7652
GTRG	20	-2048.5846	4137.1692	0.1108	0.8760
GTRIG	21	-2047.9392	4137.8783	0.0777	0.9538
GTRI	20	-2049.6088	4139.2176	0.0398	0.9936
K2PG	13	-2059.3887	4144.7774	0.0025	0.9960
K2PI	13	-2059.7148	4145.4295	0.0018	0.9978
K2PIG	14	-2059.1579	4146.3158	0.0011	0.9990
SYMG	17	-2057.0060	4148.0120	0.0005	0.9994
SYMIG	18	-2056.3987	4148.7973	0.0003	0.9998
SYMI	17	-2057.7986	4149.5972	0.0002	1.0000
GTR	19	-2088.2381	4214.4763	0.0000	1.0000
HKY	15	-2096.5437	4223.0874	0.0000	1.0000
SYM	16	-2099.1242	4230.2484	0.0000	1.0000
F81G	15	-2102.8633	4235.7266	0.0000	1.0000
F81IG	16	-2102.3465	4236.6930	0.0000	1.0000
K2P	12	-2106.3827	4236.7653	0.0000	1.0000
F81I	15	-2103.9084	4237.8168	0.0000	1.0000
JC69G	12	-2112.1331	4248.2662	0.0000	1.0000
JC69IG	13	-2111.5826	4249.1653	0.0000	1.0000
JC69I	12	-2113.1611	4250.3222	0.0000	1.0000
F81	14	-2143.5617	4315.1234	0.0000	1.0000
JC69	11	-2153.0368	4328.0736	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-2059.3887	4204.8037	0.4986	0.4986
K2PI	13	-2059.7148	4205.4558	0.3599	0.8585
HKYG	16	-2051.3859	4208.6502	0.0729	0.9314
HKYI	16	-2051.9374	4209.7532	0.0420	0.9734
K2PIG	14	-2059.1579	4210.9595	0.0230	0.9963
HKYIG	17	-2051.0670	4214.6298	0.0037	1.0000
SYMG	17	-2057.0060	4226.5079	0.0000	1.0000
SYMI	17	-2057.7986	4228.0930	0.0000	1.0000
GTRG	20	-2048.5846	4229.5173	0.0000	1.0000
GTRI	20	-2049.6088	4231.5656	0.0000	1.0000
SYMIG	18	-2056.3987	4231.9106	0.0000	1.0000
GTRIG	21	-2047.9392	4234.8438	0.0000	1.0000
K2P	12	-2106.3827	4292.1742	0.0000	1.0000
HKY	15	-2096.5437	4292.3484	0.0000	1.0000
GTR	19	-2088.2381	4302.2069	0.0000	1.0000
JC69G	12	-2112.1331	4303.6751	0.0000	1.0000
SYM	16	-2099.1242	4304.1268	0.0000	1.0000
F81G	15	-2102.8633	4304.9877	0.0000	1.0000
JC69I	12	-2113.1611	4305.7310	0.0000	1.0000
F81I	15	-2103.9084	4307.0778	0.0000	1.0000
JC69IG	13	-2111.5826	4309.1915	0.0000	1.0000
F81IG	16	-2102.3465	4310.5714	0.0000	1.0000
JC69	11	-2153.0368	4378.8651	0.0000	1.0000
F81	14	-2143.5617	4379.7670	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/chr8_9165.phylip = 1 - 748;
 Partition Dummy = 1:591-tomodel-phylip/chr8_9165.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/chr8_9165.phylip = 1 - 748;
 Partition Dummy = 1:591-tomodel-phylip/chr8_9165.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 (K2PG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr8_9165.phylip = 1 - 748;
 Partition Dummy = 1:591-tomodel-phylip/chr8_9165.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

