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

Input data from file "683-tomodel-phylip/chr1_32375.phylip" (7 taxa, 794 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	-2698.0368	5437.1602	0.5673	0.5673
GTRIG	21	-2698.0404	5439.2777	0.1968	0.7641
GTRI	20	-2699.1664	5439.4195	0.1833	0.9474
HKYG	16	-2705.2958	5443.2917	0.0264	0.9738
HKYI	16	-2705.7718	5444.2437	0.0164	0.9902
HKYIG	17	-2705.2482	5445.2851	0.0098	1.0000
GTR	19	-2724.7950	5488.5719	0.0000	1.0000
HKY	15	-2730.9650	5492.5470	0.0000	1.0000
K2PI	13	-2744.4331	5515.3329	0.0000	1.0000
K2PG	13	-2744.5913	5515.6493	0.0000	1.0000
SYMI	17	-2741.0769	5516.9425	0.0000	1.0000
K2PIG	14	-2744.2408	5517.0207	0.0000	1.0000
SYMG	17	-2741.2592	5517.3070	0.0000	1.0000
SYMIG	18	-2740.9147	5518.7120	0.0000	1.0000
F81G	15	-2751.9347	5534.4863	0.0000	1.0000
F81I	15	-2751.9557	5534.5284	0.0000	1.0000
F81IG	16	-2751.8062	5536.3126	0.0000	1.0000
K2P	12	-2770.1161	5564.6318	0.0000	1.0000
SYM	16	-2767.0926	5566.8854	0.0000	1.0000
F81	14	-2774.8306	5578.2004	0.0000	1.0000
JC69I	12	-2785.1641	5594.7277	0.0000	1.0000
JC69G	12	-2785.4909	5595.3812	0.0000	1.0000
JC69IG	13	-2785.0850	5596.6367	0.0000	1.0000
JC69	11	-2808.8204	5639.9783	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2698.0368	5436.0735	0.5660	0.5660
GTRIG	21	-2698.0404	5438.0808	0.2075	0.7735
GTRI	20	-2699.1664	5438.3328	0.1829	0.9563
HKYG	16	-2705.2958	5442.5916	0.0217	0.9781
HKYI	16	-2705.7718	5443.5436	0.0135	0.9916
HKYIG	17	-2705.2482	5444.4965	0.0084	1.0000
GTR	19	-2724.7950	5487.5900	0.0000	1.0000
HKY	15	-2730.9650	5491.9300	0.0000	1.0000
K2PI	13	-2744.4331	5514.8662	0.0000	1.0000
K2PG	13	-2744.5913	5515.1826	0.0000	1.0000
SYMI	17	-2741.0769	5516.1538	0.0000	1.0000
K2PIG	14	-2744.2408	5516.4816	0.0000	1.0000
SYMG	17	-2741.2592	5516.5183	0.0000	1.0000
SYMIG	18	-2740.9147	5517.8294	0.0000	1.0000
F81G	15	-2751.9347	5533.8694	0.0000	1.0000
F81I	15	-2751.9557	5533.9114	0.0000	1.0000
F81IG	16	-2751.8062	5535.6125	0.0000	1.0000
K2P	12	-2770.1161	5564.2323	0.0000	1.0000
SYM	16	-2767.0926	5566.1853	0.0000	1.0000
F81	14	-2774.8306	5577.6612	0.0000	1.0000
JC69I	12	-2785.1641	5594.3282	0.0000	1.0000
JC69G	12	-2785.4909	5594.9817	0.0000	1.0000
JC69IG	13	-2785.0850	5596.1701	0.0000	1.0000
JC69	11	-2808.8204	5639.6407	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2705.2958	5517.4249	0.6019	0.6019
HKYI	16	-2705.7718	5518.3769	0.3739	0.9758
HKYIG	17	-2705.2482	5524.0069	0.0224	0.9982
GTRG	20	-2698.0368	5529.6152	0.0014	0.9995
GTRI	20	-2699.1664	5531.8745	0.0004	1.0000
GTRIG	21	-2698.0404	5536.2995	0.0000	1.0000
HKY	15	-2730.9650	5562.0863	0.0000	1.0000
K2PI	13	-2744.4331	5575.6683	0.0000	1.0000
K2PG	13	-2744.5913	5575.9847	0.0000	1.0000
GTR	19	-2724.7950	5576.4545	0.0000	1.0000
K2PIG	14	-2744.2408	5581.9607	0.0000	1.0000
SYMI	17	-2741.0769	5595.6642	0.0000	1.0000
SYMG	17	-2741.2592	5596.0287	0.0000	1.0000
SYMIG	18	-2740.9147	5602.0169	0.0000	1.0000
F81G	15	-2751.9347	5604.0256	0.0000	1.0000
F81I	15	-2751.9557	5604.0677	0.0000	1.0000
F81IG	16	-2751.8062	5610.4458	0.0000	1.0000
K2P	12	-2770.1161	5620.3573	0.0000	1.0000
SYM	16	-2767.0926	5641.0186	0.0000	1.0000
F81	14	-2774.8306	5643.1404	0.0000	1.0000
JC69I	12	-2785.1641	5650.4532	0.0000	1.0000
JC69G	12	-2785.4909	5651.1067	0.0000	1.0000
JC69IG	13	-2785.0850	5656.9721	0.0000	1.0000
JC69	11	-2808.8204	5691.0887	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 683-tomodel-phylip/chr1_32375.phylip = 1 - 794;
 Partition Dummy = 1:683-tomodel-phylip/chr1_32375.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 683-tomodel-phylip/chr1_32375.phylip = 1 - 794;
 Partition Dummy = 1:683-tomodel-phylip/chr1_32375.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 683-tomodel-phylip/chr1_32375.phylip = 1 - 794;
 Partition Dummy = 1:683-tomodel-phylip/chr1_32375.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

