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

Input data from file "591-tomodel-phylip/chrZ_11558.phylip" (7 taxa, 618 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	-2469.9456	4972.7964	0.4639	0.4639
HKYIG	17	-2469.7739	4974.5678	0.1913	0.6552
HKYI	16	-2470.9333	4974.7718	0.1728	0.8280
GTRG	20	-2467.1912	4975.7894	0.1039	0.9319
GTRIG	21	-2467.0428	4977.6359	0.0413	0.9731
GTRI	20	-2468.5435	4978.4940	0.0269	1.0000
HKY	15	-2497.6695	5026.1363	0.0000	1.0000
GTR	19	-2494.3819	5028.0347	0.0000	1.0000
F81I	15	-2510.8239	5052.4452	0.0000	1.0000
F81G	15	-2511.0674	5052.9321	0.0000	1.0000
F81IG	16	-2510.5020	5053.9091	0.0000	1.0000
F81	14	-2536.6749	5102.0464	0.0000	1.0000
SYMG	17	-2554.9833	5144.9865	0.0000	1.0000
SYMIG	18	-2554.9753	5147.0924	0.0000	1.0000
SYMI	17	-2556.4318	5147.8836	0.0000	1.0000
K2PG	13	-2563.3148	5153.2322	0.0000	1.0000
K2PI	13	-2564.1472	5154.8971	0.0000	1.0000
K2PIG	14	-2563.2087	5155.1138	0.0000	1.0000
JC69I	12	-2588.5456	5201.6069	0.0000	1.0000
JC69G	12	-2588.5905	5201.6967	0.0000	1.0000
JC69IG	13	-2588.1233	5202.8493	0.0000	1.0000
SYM	16	-2586.2141	5205.3334	0.0000	1.0000
K2P	12	-2592.6754	5209.8665	0.0000	1.0000
JC69	11	-2616.0500	5254.5356	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2469.9456	4971.8912	0.4359	0.4359
HKYIG	17	-2469.7739	4973.5478	0.1904	0.6263
HKYI	16	-2470.9333	4973.8667	0.1623	0.7886
GTRG	20	-2467.1912	4974.3823	0.1254	0.9140
GTRIG	21	-2467.0428	4976.0856	0.0535	0.9676
GTRI	20	-2468.5435	4977.0870	0.0324	1.0000
HKY	15	-2497.6695	5025.3389	0.0000	1.0000
GTR	19	-2494.3819	5026.7638	0.0000	1.0000
F81I	15	-2510.8239	5051.6478	0.0000	1.0000
F81G	15	-2511.0674	5052.1348	0.0000	1.0000
F81IG	16	-2510.5020	5053.0040	0.0000	1.0000
F81	14	-2536.6749	5101.3498	0.0000	1.0000
SYMG	17	-2554.9833	5143.9665	0.0000	1.0000
SYMIG	18	-2554.9753	5145.9505	0.0000	1.0000
SYMI	17	-2556.4318	5146.8636	0.0000	1.0000
K2PG	13	-2563.3148	5152.6295	0.0000	1.0000
K2PI	13	-2564.1472	5154.2945	0.0000	1.0000
K2PIG	14	-2563.2087	5154.4173	0.0000	1.0000
JC69I	12	-2588.5456	5201.0912	0.0000	1.0000
JC69G	12	-2588.5905	5201.1810	0.0000	1.0000
JC69IG	13	-2588.1233	5202.2466	0.0000	1.0000
SYM	16	-2586.2141	5204.4282	0.0000	1.0000
K2P	12	-2592.6754	5209.3508	0.0000	1.0000
JC69	11	-2616.0500	5254.0999	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2469.9456	5042.7151	0.7041	0.7041
HKYI	16	-2470.9333	5044.6905	0.2622	0.9663
HKYIG	17	-2469.7739	5048.7981	0.0336	1.0000
GTRG	20	-2467.1912	5062.9121	0.0000	1.0000
GTRI	20	-2468.5435	5065.6167	0.0000	1.0000
GTRIG	21	-2467.0428	5069.0418	0.0000	1.0000
HKY	15	-2497.6695	5091.7362	0.0000	1.0000
GTR	19	-2494.3819	5110.8671	0.0000	1.0000
F81I	15	-2510.8239	5118.0452	0.0000	1.0000
F81G	15	-2511.0674	5118.5321	0.0000	1.0000
F81IG	16	-2510.5020	5123.8278	0.0000	1.0000
F81	14	-2536.6749	5163.3207	0.0000	1.0000
K2PG	13	-2563.3148	5210.1739	0.0000	1.0000
K2PI	13	-2564.1472	5211.8388	0.0000	1.0000
K2PIG	14	-2563.2087	5216.3881	0.0000	1.0000
SYMG	17	-2554.9833	5219.2168	0.0000	1.0000
SYMI	17	-2556.4318	5222.1139	0.0000	1.0000
SYMIG	18	-2554.9753	5225.6273	0.0000	1.0000
JC69I	12	-2588.5456	5254.2090	0.0000	1.0000
JC69G	12	-2588.5905	5254.2989	0.0000	1.0000
JC69IG	13	-2588.1233	5259.7910	0.0000	1.0000
K2P	12	-2592.6754	5262.4687	0.0000	1.0000
SYM	16	-2586.2141	5275.2520	0.0000	1.0000
JC69	11	-2616.0500	5302.7913	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_11558.phylip = 1 - 618;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_11558.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_11558.phylip = 1 - 618;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_11558.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_11558.phylip = 1 - 618;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_11558.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

