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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKYG


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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2416.8548	4865.7097	0.4022	0.4022
HKYI	16	-2417.1626	4866.3252	0.2956	0.6978
HKYIG	17	-2416.8617	4867.7234	0.1469	0.8448
GTRG	20	-2414.5426	4869.0852	0.0744	0.9192
GTRI	20	-2414.8810	4869.7620	0.0530	0.9722
GTRIG	21	-2414.5474	4871.0949	0.0272	0.9994
HKY	15	-2424.5484	4879.0968	0.0005	0.9999
GTR	19	-2422.2403	4882.4806	0.0001	1.0000
K2PG	13	-2449.9657	4925.9314	0.0000	1.0000
K2PI	13	-2450.1961	4926.3921	0.0000	1.0000
K2PIG	14	-2449.9770	4927.9541	0.0000	1.0000
SYMG	17	-2448.1935	4930.3870	0.0000	1.0000
SYMI	17	-2448.4070	4930.8141	0.0000	1.0000
SYMIG	18	-2448.2045	4932.4091	0.0000	1.0000
K2P	12	-2457.5822	4939.1643	0.0000	1.0000
SYM	16	-2455.7683	4943.5366	0.0000	1.0000
F81G	15	-2471.6657	4973.3314	0.0000	1.0000
F81I	15	-2471.7653	4973.5305	0.0000	1.0000
F81IG	16	-2471.6711	4975.3423	0.0000	1.0000
F81	14	-2478.0355	4984.0709	0.0000	1.0000
JC69G	12	-2500.0577	5024.1155	0.0000	1.0000
JC69I	12	-2500.1254	5024.2509	0.0000	1.0000
JC69IG	13	-2500.0616	5026.1232	0.0000	1.0000
JC69	11	-2506.5008	5035.0017	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2416.8548	4943.2777	0.5613	0.5613
HKYI	16	-2417.1626	4943.8933	0.4126	0.9740
HKYIG	17	-2416.8617	4950.1395	0.0182	0.9921
HKY	15	-2424.5484	4951.8169	0.0079	1.0000
GTRG	20	-2414.5426	4966.0453	0.0000	1.0000
GTRI	20	-2414.8810	4966.7221	0.0000	1.0000
GTRIG	21	-2414.5474	4972.9030	0.0000	1.0000
GTR	19	-2422.2403	4974.5927	0.0000	1.0000
K2PG	13	-2449.9657	4988.9554	0.0000	1.0000
K2PI	13	-2450.1961	4989.4162	0.0000	1.0000
K2PIG	14	-2449.9770	4995.8262	0.0000	1.0000
K2P	12	-2457.5822	4997.3404	0.0000	1.0000
SYMG	17	-2448.1935	5012.8031	0.0000	1.0000
SYMI	17	-2448.4070	5013.2302	0.0000	1.0000
SYMIG	18	-2448.2045	5019.6732	0.0000	1.0000
SYM	16	-2455.7683	5021.1047	0.0000	1.0000
F81G	15	-2471.6657	5046.0515	0.0000	1.0000
F81I	15	-2471.7653	5046.2506	0.0000	1.0000
F81	14	-2478.0355	5051.9430	0.0000	1.0000
F81IG	16	-2471.6711	5052.9103	0.0000	1.0000
JC69G	12	-2500.0577	5082.2916	0.0000	1.0000
JC69I	12	-2500.1254	5082.4269	0.0000	1.0000
JC69	11	-2506.5008	5088.3298	0.0000	1.0000
JC69IG	13	-2500.0616	5089.1473	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr8_4403.phylip = 1 - 942;
 Partition Dummy = 1:683-tomodel-phylip/chr8_4403.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;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr8_4403.phylip = 1 - 942;
 Partition Dummy = 1:683-tomodel-phylip/chr8_4403.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;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr8_4403.phylip = 1 - 942;
 Partition Dummy = 1:683-tomodel-phylip/chr8_4403.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

