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

Input data from file "683-tomodel-phylip/chr12_2162.phylip" (7 taxa, 1021 characters)

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-3576.1526	7184.3053	0.2845	0.2845
HKYIG	17	-3575.5274	7185.0547	0.1956	0.4801
GTRI	20	-3572.6761	7185.3521	0.1686	0.6487
HKYG	16	-3577.0152	7186.0304	0.1201	0.7687
GTRIG	21	-3572.0311	7186.0623	0.1182	0.8869
GTRG	20	-3573.0752	7186.1505	0.1131	1.0000
HKY	15	-3590.5035	7211.0071	0.0000	1.0000
GTR	19	-3586.7414	7211.4829	0.0000	1.0000
K2PI	13	-3598.3551	7222.7101	0.0000	1.0000
K2PIG	14	-3597.4860	7222.9720	0.0000	1.0000
K2PG	13	-3598.5269	7223.0537	0.0000	1.0000
SYMG	17	-3595.9139	7225.8277	0.0000	1.0000
SYMI	17	-3596.1293	7226.2586	0.0000	1.0000
SYMIG	18	-3595.2520	7226.5039	0.0000	1.0000
K2P	12	-3611.7063	7247.4127	0.0000	1.0000
SYM	16	-3609.0194	7250.0388	0.0000	1.0000
F81G	15	-3663.4908	7356.9815	0.0000	1.0000
F81I	15	-3663.5042	7357.0084	0.0000	1.0000
F81IG	16	-3664.2862	7360.5724	0.0000	1.0000
F81	14	-3674.4952	7376.9903	0.0000	1.0000
JC69G	12	-3680.2773	7384.5545	0.0000	1.0000
JC69I	12	-3680.4328	7384.8655	0.0000	1.0000
JC69IG	13	-3681.3096	7388.6192	0.0000	1.0000
JC69	11	-3691.2458	7404.4916	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-3576.1526	7263.1619	0.6754	0.6754
HKYG	16	-3577.0152	7264.8870	0.2851	0.9604
HKYIG	17	-3575.5274	7268.8399	0.0395	0.9999
GTRI	20	-3572.6761	7283.9229	0.0000	1.0000
GTRG	20	-3573.0752	7284.7213	0.0000	1.0000
HKY	15	-3590.5035	7284.9351	0.0000	1.0000
K2PI	13	-3598.3551	7286.7811	0.0000	1.0000
K2PG	13	-3598.5269	7287.1247	0.0000	1.0000
GTRIG	21	-3572.0311	7289.5616	0.0000	1.0000
K2PIG	14	-3597.4860	7291.9715	0.0000	1.0000
GTR	19	-3586.7414	7305.1251	0.0000	1.0000
K2P	12	-3611.7063	7306.5551	0.0000	1.0000
SYMG	17	-3595.9139	7309.6129	0.0000	1.0000
SYMI	17	-3596.1293	7310.0437	0.0000	1.0000
SYMIG	18	-3595.2520	7315.2176	0.0000	1.0000
SYM	16	-3609.0194	7328.8954	0.0000	1.0000
F81G	15	-3663.4908	7430.9096	0.0000	1.0000
F81I	15	-3663.5042	7430.9365	0.0000	1.0000
F81IG	16	-3664.2862	7439.4290	0.0000	1.0000
JC69G	12	-3680.2773	7443.6970	0.0000	1.0000
JC69I	12	-3680.4328	7444.0080	0.0000	1.0000
F81	14	-3674.4952	7445.9898	0.0000	1.0000
JC69IG	13	-3681.3096	7452.6902	0.0000	1.0000
JC69	11	-3691.2458	7458.7056	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr12_2162.phylip = 1 - 1021;
 Partition Dummy = 1:683-tomodel-phylip/chr12_2162.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr12_2162.phylip = 1 - 1021;
 Partition Dummy = 1:683-tomodel-phylip/chr12_2162.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best BIC model (HKYI)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr12_2162.phylip = 1 - 1021;
 Partition Dummy = 1:683-tomodel-phylip/chr12_2162.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


-------------------------------------------------------------
End of Output

