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

Input data from file "591-tomodel-phylip/chrZ_6769.phylip" (7 taxa, 663 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	-2028.7378	4090.3178	0.4287	0.4287
HKYI	16	-2029.0131	4090.8683	0.3255	0.7542
HKYIG	17	-2028.6717	4092.2922	0.1597	0.9140
GTRG	20	-2026.9873	4095.2831	0.0358	0.9498
GTRI	20	-2027.0055	4095.3195	0.0352	0.9850
GTRIG	21	-2026.7888	4097.0190	0.0150	1.0000
K2PG	13	-2043.4874	4113.5357	0.0000	1.0000
K2PI	13	-2043.7994	4114.1596	0.0000	1.0000
K2PIG	14	-2043.4112	4115.4705	0.0000	1.0000
SYMG	17	-2041.2695	4117.4879	0.0000	1.0000
SYMI	17	-2041.5655	4118.0799	0.0000	1.0000
SYMIG	18	-2041.1484	4119.3589	0.0000	1.0000
HKY	15	-2051.5894	4133.9207	0.0000	1.0000
GTR	19	-2049.3322	4137.8463	0.0000	1.0000
F81G	15	-2060.7667	4152.2753	0.0000	1.0000
F81I	15	-2061.0075	4152.7569	0.0000	1.0000
F81IG	16	-2060.6933	4154.2287	0.0000	1.0000
K2P	12	-2066.1201	4156.7202	0.0000	1.0000
SYM	16	-2063.8247	4160.4914	0.0000	1.0000
JC69G	12	-2073.8099	4172.0998	0.0000	1.0000
JC69I	12	-2074.1088	4172.6977	0.0000	1.0000
JC69IG	13	-2073.7344	4174.0297	0.0000	1.0000
F81	14	-2082.4907	4193.6295	0.0000	1.0000
JC69	11	-2095.5768	4213.5590	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2028.7378	4089.4757	0.4152	0.4152
HKYI	16	-2029.0131	4090.0262	0.3153	0.7304
HKYIG	17	-2028.6717	4091.3434	0.1632	0.8936
GTRG	20	-2026.9873	4093.9746	0.0438	0.9374
GTRI	20	-2027.0055	4094.0111	0.0430	0.9803
GTRIG	21	-2026.7888	4095.5775	0.0196	1.0000
K2PG	13	-2043.4874	4112.9748	0.0000	1.0000
K2PI	13	-2043.7994	4113.5987	0.0000	1.0000
K2PIG	14	-2043.4112	4114.8224	0.0000	1.0000
SYMG	17	-2041.2695	4116.5390	0.0000	1.0000
SYMI	17	-2041.5655	4117.1310	0.0000	1.0000
SYMIG	18	-2041.1484	4118.2968	0.0000	1.0000
HKY	15	-2051.5894	4133.1788	0.0000	1.0000
GTR	19	-2049.3322	4136.6644	0.0000	1.0000
F81G	15	-2060.7667	4151.5334	0.0000	1.0000
F81I	15	-2061.0075	4152.0150	0.0000	1.0000
F81IG	16	-2060.6933	4153.3866	0.0000	1.0000
K2P	12	-2066.1201	4156.2402	0.0000	1.0000
SYM	16	-2063.8247	4159.6493	0.0000	1.0000
JC69G	12	-2073.8099	4171.6198	0.0000	1.0000
JC69I	12	-2074.1088	4172.2177	0.0000	1.0000
JC69IG	13	-2073.7344	4173.4689	0.0000	1.0000
F81	14	-2082.4907	4192.9814	0.0000	1.0000
JC69	11	-2095.5768	4213.1535	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2028.7378	4161.4241	0.5516	0.5516
HKYI	16	-2029.0131	4161.9746	0.4189	0.9705
HKYIG	17	-2028.6717	4167.7886	0.0229	0.9934
K2PG	13	-2043.4874	4171.4329	0.0037	0.9971
K2PI	13	-2043.7994	4172.0568	0.0027	0.9998
K2PIG	14	-2043.4112	4177.7772	0.0002	1.0000
GTRG	20	-2026.9873	4183.9101	0.0000	1.0000
GTRI	20	-2027.0055	4183.9466	0.0000	1.0000
GTRIG	21	-2026.7888	4190.0098	0.0000	1.0000
SYMG	17	-2041.2695	4192.9842	0.0000	1.0000
SYMI	17	-2041.5655	4193.5762	0.0000	1.0000
SYMIG	18	-2041.1484	4199.2387	0.0000	1.0000
HKY	15	-2051.5894	4200.6304	0.0000	1.0000
K2P	12	-2066.1201	4210.2015	0.0000	1.0000
F81G	15	-2060.7667	4218.9851	0.0000	1.0000
F81I	15	-2061.0075	4219.4667	0.0000	1.0000
GTR	19	-2049.3322	4222.1031	0.0000	1.0000
F81IG	16	-2060.6933	4225.3350	0.0000	1.0000
JC69G	12	-2073.8099	4225.5811	0.0000	1.0000
JC69I	12	-2074.1088	4226.1790	0.0000	1.0000
SYM	16	-2063.8247	4231.5977	0.0000	1.0000
JC69IG	13	-2073.7344	4231.9269	0.0000	1.0000
F81	14	-2082.4907	4255.9362	0.0000	1.0000
JC69	11	-2095.5768	4262.6180	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/chrZ_6769.phylip = 1 - 663;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_6769.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/chrZ_6769.phylip = 1 - 663;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_6769.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/chrZ_6769.phylip = 1 - 663;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_6769.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

