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

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

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


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
HKYI	16	-2407.6731	4848.0245	0.4623	0.4623
HKYG	16	-2408.7350	4850.1483	0.1599	0.6222
HKYIG	17	-2407.6986	4850.1613	0.1588	0.7810
GTRI	20	-2404.7530	4850.5587	0.1302	0.9112
GTRG	20	-2405.8273	4852.7072	0.0445	0.9557
GTRIG	21	-2404.7782	4852.7157	0.0443	1.0000
HKY	15	-2421.5283	4873.6544	0.0000	1.0000
GTR	19	-2419.1218	4877.1948	0.0000	1.0000
K2PI	13	-2425.5858	4877.6237	0.0000	1.0000
K2PIG	14	-2425.6155	4879.7534	0.0000	1.0000
K2PG	13	-2426.7729	4879.9980	0.0000	1.0000
SYMI	17	-2424.4678	4883.6996	0.0000	1.0000
SYMIG	18	-2424.4986	4885.8521	0.0000	1.0000
SYMG	17	-2425.6989	4886.1618	0.0000	1.0000
K2P	12	-2438.7145	4901.8160	0.0000	1.0000
SYM	16	-2437.5445	4907.7674	0.0000	1.0000
F81I	15	-2461.1284	4952.8545	0.0000	1.0000
F81G	15	-2462.1582	4954.9142	0.0000	1.0000
F81IG	16	-2461.1541	4954.9865	0.0000	1.0000
JC69I	12	-2475.4048	4975.1967	0.0000	1.0000
F81	14	-2473.9101	4976.3427	0.0000	1.0000
JC69IG	13	-2475.4319	4977.3160	0.0000	1.0000
JC69G	12	-2476.4769	4977.3409	0.0000	1.0000
JC69	11	-2487.9423	4998.2118	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-2407.6731	4847.3462	0.4382	0.4382
HKYIG	17	-2407.6986	4849.3973	0.1572	0.5954
HKYG	16	-2408.7350	4849.4700	0.1515	0.7469
GTRI	20	-2404.7530	4849.5061	0.1488	0.8958
GTRIG	21	-2404.7782	4851.5564	0.0534	0.9492
GTRG	20	-2405.8273	4851.6546	0.0508	1.0000
HKY	15	-2421.5283	4873.0567	0.0000	1.0000
GTR	19	-2419.1218	4876.2437	0.0000	1.0000
K2PI	13	-2425.5858	4877.1716	0.0000	1.0000
K2PIG	14	-2425.6155	4879.2310	0.0000	1.0000
K2PG	13	-2426.7729	4879.5458	0.0000	1.0000
SYMI	17	-2424.4678	4882.9355	0.0000	1.0000
SYMIG	18	-2424.4986	4884.9971	0.0000	1.0000
SYMG	17	-2425.6989	4885.3978	0.0000	1.0000
K2P	12	-2438.7145	4901.4289	0.0000	1.0000
SYM	16	-2437.5445	4907.0891	0.0000	1.0000
F81I	15	-2461.1284	4952.2568	0.0000	1.0000
F81IG	16	-2461.1541	4954.3082	0.0000	1.0000
F81G	15	-2462.1582	4954.3165	0.0000	1.0000
JC69I	12	-2475.4048	4974.8096	0.0000	1.0000
F81	14	-2473.9101	4975.8203	0.0000	1.0000
JC69IG	13	-2475.4319	4976.8639	0.0000	1.0000
JC69G	12	-2476.4769	4976.9538	0.0000	1.0000
JC69	11	-2487.9423	4997.8847	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2407.6731	4922.6756	0.7244	0.7244
HKYG	16	-2408.7350	4924.7993	0.2505	0.9749
HKYIG	17	-2407.6986	4929.4347	0.0247	0.9996
K2PI	13	-2425.5858	4938.3767	0.0003	0.9999
K2PG	13	-2426.7729	4940.7509	0.0001	0.9999
GTRI	20	-2404.7530	4943.6677	0.0000	1.0000
HKY	15	-2421.5283	4943.6779	0.0000	1.0000
K2PIG	14	-2425.6155	4945.1442	0.0000	1.0000
GTRG	20	-2405.8273	4945.8163	0.0000	1.0000
GTRIG	21	-2404.7782	4950.4261	0.0000	1.0000
K2P	12	-2438.7145	4957.9259	0.0000	1.0000
SYMI	17	-2424.4678	4962.9729	0.0000	1.0000
SYMG	17	-2425.6989	4965.4352	0.0000	1.0000
GTR	19	-2419.1218	4965.6973	0.0000	1.0000
SYMIG	18	-2424.4986	4969.7427	0.0000	1.0000
SYM	16	-2437.5445	4982.4184	0.0000	1.0000
F81I	15	-2461.1284	5022.8780	0.0000	1.0000
F81G	15	-2462.1582	5024.9377	0.0000	1.0000
F81IG	16	-2461.1541	5029.6376	0.0000	1.0000
JC69I	12	-2475.4048	5031.3066	0.0000	1.0000
JC69G	12	-2476.4769	5033.4508	0.0000	1.0000
JC69IG	13	-2475.4319	5038.0690	0.0000	1.0000
F81	14	-2473.9101	5041.7335	0.0000	1.0000
JC69	11	-2487.9423	5049.6736	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/chr2_22620.phylip = 1 - 819;
 Partition Dummy = 1:683-tomodel-phylip/chr2_22620.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/chr2_22620.phylip = 1 - 819;
 Partition Dummy = 1:683-tomodel-phylip/chr2_22620.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/chr2_22620.phylip = 1 - 819;
 Partition Dummy = 1:683-tomodel-phylip/chr2_22620.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

