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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	20	-2247.1741	4535.4196	0.3718	0.3718
GTRI	20	-2247.8821	4536.8355	0.1831	0.5549
HKYG	16	-2252.0825	4536.8553	0.1813	0.7363
HKYI	16	-2252.5003	4537.6909	0.1194	0.8557
GTRIG	21	-2247.4643	4538.1087	0.0969	0.9526
HKYIG	17	-2252.3801	4539.5378	0.0474	1.0000
K2PG	13	-2281.7573	4589.9747	0.0000	1.0000
K2PI	13	-2282.1745	4590.8091	0.0000	1.0000
K2PIG	14	-2282.6941	4593.9199	0.0000	1.0000
SYMG	17	-2281.1908	4597.1593	0.0000	1.0000
SYMI	17	-2281.6188	4598.0153	0.0000	1.0000
SYMIG	18	-2281.6123	4600.0948	0.0000	1.0000
HKY	15	-2288.8301	4608.2687	0.0000	1.0000
GTR	19	-2284.7073	4608.3827	0.0000	1.0000
K2P	12	-2318.5242	4661.4424	0.0000	1.0000
F81G	15	-2318.5512	4667.7109	0.0000	1.0000
SYM	16	-2317.5880	4667.8663	0.0000	1.0000
F81I	15	-2318.8982	4668.4048	0.0000	1.0000
F81IG	16	-2318.5405	4669.7714	0.0000	1.0000
JC69G	12	-2342.6918	4709.7776	0.0000	1.0000
JC69I	12	-2342.9269	4710.2478	0.0000	1.0000
JC69IG	13	-2342.7150	4711.8902	0.0000	1.0000
F81	14	-2350.6502	4729.8320	0.0000	1.0000
JC69	11	-2374.9633	4772.2595	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2247.1741	4534.3482	0.3927	0.3927
GTRI	20	-2247.8821	4535.7641	0.1934	0.5861
HKYG	16	-2252.0825	4536.1649	0.1583	0.7444
GTRIG	21	-2247.4643	4536.9286	0.1081	0.8525
HKYI	16	-2252.5003	4537.0005	0.1043	0.9568
HKYIG	17	-2252.3801	4538.7601	0.0432	1.0000
K2PG	13	-2281.7573	4589.5145	0.0000	1.0000
K2PI	13	-2282.1745	4590.3490	0.0000	1.0000
K2PIG	14	-2282.6941	4593.3882	0.0000	1.0000
SYMG	17	-2281.1908	4596.3816	0.0000	1.0000
SYMI	17	-2281.6188	4597.2377	0.0000	1.0000
SYMIG	18	-2281.6123	4599.2246	0.0000	1.0000
GTR	19	-2284.7073	4607.4146	0.0000	1.0000
HKY	15	-2288.8301	4607.6603	0.0000	1.0000
K2P	12	-2318.5242	4661.0484	0.0000	1.0000
F81G	15	-2318.5512	4667.1025	0.0000	1.0000
SYM	16	-2317.5880	4667.1759	0.0000	1.0000
F81I	15	-2318.8982	4667.7965	0.0000	1.0000
F81IG	16	-2318.5405	4669.0811	0.0000	1.0000
JC69G	12	-2342.6918	4709.3836	0.0000	1.0000
JC69I	12	-2342.9269	4709.8539	0.0000	1.0000
JC69IG	13	-2342.7150	4711.4301	0.0000	1.0000
F81	14	-2350.6502	4729.3003	0.0000	1.0000
JC69	11	-2374.9633	4771.9266	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2252.0825	4611.2184	0.5935	0.5935
HKYI	16	-2252.5003	4612.0540	0.3908	0.9843
HKYIG	17	-2252.3801	4618.5045	0.0155	0.9998
GTRG	20	-2247.1741	4628.1650	0.0001	0.9999
GTRI	20	-2247.8821	4629.5810	0.0001	1.0000
GTRIG	21	-2247.4643	4635.4363	0.0000	1.0000
K2PG	13	-2281.7573	4650.4954	0.0000	1.0000
K2PI	13	-2282.1745	4651.3299	0.0000	1.0000
K2PIG	14	-2282.6941	4659.0600	0.0000	1.0000
SYMG	17	-2281.1908	4676.1260	0.0000	1.0000
SYMI	17	-2281.6188	4676.9820	0.0000	1.0000
HKY	15	-2288.8301	4678.0229	0.0000	1.0000
SYMIG	18	-2281.6123	4683.6598	0.0000	1.0000
GTR	19	-2284.7073	4696.5406	0.0000	1.0000
K2P	12	-2318.5242	4717.3385	0.0000	1.0000
F81G	15	-2318.5512	4737.4651	0.0000	1.0000
F81I	15	-2318.8982	4738.1591	0.0000	1.0000
SYM	16	-2317.5880	4742.2294	0.0000	1.0000
F81IG	16	-2318.5405	4744.1345	0.0000	1.0000
JC69G	12	-2342.6918	4765.6737	0.0000	1.0000
JC69I	12	-2342.9269	4766.1440	0.0000	1.0000
JC69IG	13	-2342.7150	4772.4110	0.0000	1.0000
F81	14	-2350.6502	4794.9721	0.0000	1.0000
JC69	11	-2374.9633	4823.5258	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (GTRG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chrZ_4782.phylip = 1 - 805;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_4782.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


[Mrbayes block for the best AICc model (GTRG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chrZ_4782.phylip = 1 - 805;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_4782.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,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/chrZ_4782.phylip = 1 - 805;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_4782.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

