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

Input data from file "591-tomodel-phylip/chr8_675.phylip" (7 taxa, 806 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	-2303.0860	4647.2420	0.5291	0.5291
GTRI	20	-2303.7483	4648.5666	0.2728	0.8019
GTRIG	21	-2303.1001	4649.3788	0.1818	0.9836
HKYG	16	-2311.4244	4655.5384	0.0084	0.9920
HKYI	16	-2311.9252	4656.5400	0.0051	0.9971
HKYIG	17	-2311.4273	4657.6312	0.0029	1.0000
GTR	19	-2317.9703	4674.9076	0.0000	1.0000
HKY	15	-2326.2827	4683.1730	0.0000	1.0000
F81G	15	-2344.4391	4719.4857	0.0000	1.0000
F81I	15	-2344.8464	4720.3003	0.0000	1.0000
F81IG	16	-2344.4465	4721.5824	0.0000	1.0000
F81	14	-2356.8780	4742.2870	0.0000	1.0000
SYMG	17	-2372.3659	4779.5085	0.0000	1.0000
SYMI	17	-2372.9037	4780.5839	0.0000	1.0000
SYMIG	18	-2372.3689	4781.6069	0.0000	1.0000
K2PG	13	-2378.6468	4783.7532	0.0000	1.0000
K2PI	13	-2379.0590	4784.5777	0.0000	1.0000
K2PIG	14	-2378.6567	4785.8444	0.0000	1.0000
SYM	16	-2383.9062	4800.5020	0.0000	1.0000
K2P	12	-2390.7097	4805.8128	0.0000	1.0000
JC69G	12	-2401.3376	4827.0686	0.0000	1.0000
JC69I	12	-2401.6943	4827.7819	0.0000	1.0000
JC69IG	13	-2401.3466	4829.1529	0.0000	1.0000
JC69	11	-2412.7183	4847.7691	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2303.0860	4646.1720	0.5252	0.5252
GTRI	20	-2303.7483	4647.4966	0.2708	0.7960
GTRIG	21	-2303.1001	4648.2003	0.1905	0.9865
HKYG	16	-2311.4244	4654.8489	0.0069	0.9933
HKYI	16	-2311.9252	4655.8505	0.0042	0.9975
HKYIG	17	-2311.4273	4656.8546	0.0025	1.0000
GTR	19	-2317.9703	4673.9406	0.0000	1.0000
HKY	15	-2326.2827	4682.5654	0.0000	1.0000
F81G	15	-2344.4391	4718.8781	0.0000	1.0000
F81I	15	-2344.8464	4719.6927	0.0000	1.0000
F81IG	16	-2344.4465	4720.8929	0.0000	1.0000
F81	14	-2356.8780	4741.7560	0.0000	1.0000
SYMG	17	-2372.3659	4778.7319	0.0000	1.0000
SYMI	17	-2372.9037	4779.8073	0.0000	1.0000
SYMIG	18	-2372.3689	4780.7378	0.0000	1.0000
K2PG	13	-2378.6468	4783.2936	0.0000	1.0000
K2PI	13	-2379.0590	4784.1181	0.0000	1.0000
K2PIG	14	-2378.6567	4785.3134	0.0000	1.0000
SYM	16	-2383.9062	4799.8125	0.0000	1.0000
K2P	12	-2390.7097	4805.4194	0.0000	1.0000
JC69G	12	-2401.3376	4826.6751	0.0000	1.0000
JC69I	12	-2401.6943	4827.3885	0.0000	1.0000
JC69IG	13	-2401.3466	4828.6933	0.0000	1.0000
JC69	11	-2412.7183	4847.4366	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2311.4244	4729.9222	0.6056	0.6056
HKYI	16	-2311.9252	4730.9238	0.3670	0.9727
HKYIG	17	-2311.4273	4736.6200	0.0213	0.9939
GTRG	20	-2303.0860	4740.0137	0.0039	0.9978
GTRI	20	-2303.7483	4741.3382	0.0020	0.9999
GTRIG	21	-2303.1001	4746.7340	0.0001	1.0000
HKY	15	-2326.2827	4752.9467	0.0000	1.0000
GTR	19	-2317.9703	4763.0902	0.0000	1.0000
F81G	15	-2344.4391	4789.2594	0.0000	1.0000
F81I	15	-2344.8464	4790.0740	0.0000	1.0000
F81IG	16	-2344.4465	4795.9663	0.0000	1.0000
F81	14	-2356.8780	4807.4452	0.0000	1.0000
K2PG	13	-2378.6468	4844.2907	0.0000	1.0000
K2PI	13	-2379.0590	4845.1152	0.0000	1.0000
K2PIG	14	-2378.6567	4851.0026	0.0000	1.0000
SYMG	17	-2372.3659	4858.4973	0.0000	1.0000
SYMI	17	-2372.9037	4859.5727	0.0000	1.0000
K2P	12	-2390.7097	4861.7244	0.0000	1.0000
SYMIG	18	-2372.3689	4865.1953	0.0000	1.0000
SYM	16	-2383.9062	4874.8858	0.0000	1.0000
JC69G	12	-2401.3376	4882.9801	0.0000	1.0000
JC69I	12	-2401.6943	4883.6935	0.0000	1.0000
JC69IG	13	-2401.3466	4889.6904	0.0000	1.0000
JC69	11	-2412.7183	4899.0495	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/chr8_675.phylip = 1 - 806;
 Partition Dummy = 1:591-tomodel-phylip/chr8_675.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/chr8_675.phylip = 1 - 806;
 Partition Dummy = 1:591-tomodel-phylip/chr8_675.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/chr8_675.phylip = 1 - 806;
 Partition Dummy = 1:591-tomodel-phylip/chr8_675.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

