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

Input data from file "591-tomodel-phylip/chr7_10865.phylip" (7 taxa, 626 characters)

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	16	-2023.3977	4079.6887	0.5881	0.5881
HKYIG	17	-2023.3977	4081.8020	0.2044	0.7926
HKYI	16	-2025.1144	4083.1220	0.1057	0.8982
GTRG	20	-2021.2797	4083.9478	0.0699	0.9682
GTRIG	21	-2021.2797	4086.0891	0.0240	0.9921
GTRI	20	-2023.4651	4088.3185	0.0079	1.0000
HKY	15	-2036.2514	4103.2897	0.0000	1.0000
GTR	19	-2033.4511	4106.1562	0.0000	1.0000
K2PG	13	-2056.6797	4139.9542	0.0000	1.0000
K2PIG	14	-2056.6807	4142.0489	0.0000	1.0000
K2PI	13	-2057.8702	4142.3352	0.0000	1.0000
SYMG	17	-2055.4396	4145.8857	0.0000	1.0000
SYMIG	18	-2055.4397	4148.0062	0.0000	1.0000
SYMI	17	-2056.7200	4148.4466	0.0000	1.0000
K2P	12	-2069.8223	4164.1535	0.0000	1.0000
SYM	16	-2068.4415	4169.7762	0.0000	1.0000
F81G	15	-2071.6011	4173.9890	0.0000	1.0000
F81I	15	-2073.2394	4177.2657	0.0000	1.0000
F81IG	16	-2072.4659	4177.8251	0.0000	1.0000
F81	14	-2083.6283	4195.9441	0.0000	1.0000
JC69G	12	-2099.5188	4223.5466	0.0000	1.0000
JC69I	12	-2100.4166	4225.3422	0.0000	1.0000
JC69IG	13	-2099.5217	4225.6382	0.0000	1.0000
JC69	11	-2110.9717	4244.3734	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2023.3977	4078.7954	0.5640	0.5640
HKYIG	17	-2023.3977	4080.7954	0.2075	0.7715
HKYI	16	-2025.1144	4082.2288	0.1013	0.8728
GTRG	20	-2021.2797	4082.5593	0.0859	0.9587
GTRIG	21	-2021.2797	4084.5593	0.0316	0.9903
GTRI	20	-2023.4651	4086.9301	0.0097	1.0000
HKY	15	-2036.2514	4102.5029	0.0000	1.0000
GTR	19	-2033.4511	4104.9021	0.0000	1.0000
K2PG	13	-2056.6797	4139.3594	0.0000	1.0000
K2PIG	14	-2056.6807	4141.3615	0.0000	1.0000
K2PI	13	-2057.8702	4141.7404	0.0000	1.0000
SYMG	17	-2055.4396	4144.8792	0.0000	1.0000
SYMIG	18	-2055.4397	4146.8794	0.0000	1.0000
SYMI	17	-2056.7200	4147.4400	0.0000	1.0000
K2P	12	-2069.8223	4163.6445	0.0000	1.0000
SYM	16	-2068.4415	4168.8829	0.0000	1.0000
F81G	15	-2071.6011	4173.2021	0.0000	1.0000
F81I	15	-2073.2394	4176.4788	0.0000	1.0000
F81IG	16	-2072.4659	4176.9318	0.0000	1.0000
F81	14	-2083.6283	4195.2567	0.0000	1.0000
JC69G	12	-2099.5188	4223.0376	0.0000	1.0000
JC69I	12	-2100.4166	4224.8332	0.0000	1.0000
JC69IG	13	-2099.5217	4225.0435	0.0000	1.0000
JC69	11	-2110.9717	4243.9435	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2023.3977	4149.8250	0.8199	0.8199
HKYI	16	-2025.1144	4153.2584	0.1473	0.9672
HKYIG	17	-2023.3977	4156.2643	0.0328	0.9999
HKY	15	-2036.2514	4169.0931	0.0001	1.0000
GTRG	20	-2021.2797	4171.3463	0.0000	1.0000
GTRI	20	-2023.4651	4175.7171	0.0000	1.0000
GTRIG	21	-2021.2797	4177.7857	0.0000	1.0000
GTR	19	-2033.4511	4189.2498	0.0000	1.0000
K2PG	13	-2056.6797	4197.0710	0.0000	1.0000
K2PI	13	-2057.8702	4199.4520	0.0000	1.0000
K2PIG	14	-2056.6807	4203.5124	0.0000	1.0000
K2P	12	-2069.8223	4216.9167	0.0000	1.0000
SYMG	17	-2055.4396	4220.3481	0.0000	1.0000
SYMI	17	-2056.7200	4222.9089	0.0000	1.0000
SYMIG	18	-2055.4397	4226.7877	0.0000	1.0000
F81G	15	-2071.6011	4239.7924	0.0000	1.0000
SYM	16	-2068.4415	4239.9125	0.0000	1.0000
F81I	15	-2073.2394	4243.0691	0.0000	1.0000
F81IG	16	-2072.4659	4247.9614	0.0000	1.0000
F81	14	-2083.6283	4257.4076	0.0000	1.0000
JC69G	12	-2099.5188	4276.3098	0.0000	1.0000
JC69I	12	-2100.4166	4278.1054	0.0000	1.0000
JC69IG	13	-2099.5217	4282.7550	0.0000	1.0000
JC69	11	-2110.9717	4292.7763	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 591-tomodel-phylip/chr7_10865.phylip = 1 - 626;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10865.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 591-tomodel-phylip/chr7_10865.phylip = 1 - 626;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10865.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 591-tomodel-phylip/chr7_10865.phylip = 1 - 626;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10865.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

