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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
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
GTRI	20	-2096.4878	4234.0794	0.2714	0.2714
GTRG	20	-2096.6098	4234.3234	0.2403	0.5117
HKYI	16	-2101.1789	4235.0689	0.1655	0.6772
HKYG	16	-2101.1845	4235.0800	0.1646	0.8418
GTRIG	21	-2096.4561	4236.1280	0.0975	0.9393
HKYIG	17	-2101.1373	4237.0757	0.0607	1.0000
GTR	19	-2106.6573	4252.3121	0.0000	1.0000
HKY	15	-2111.4529	4253.5325	0.0000	1.0000
K2PI	13	-2149.6927	4325.8594	0.0000	1.0000
K2PG	13	-2149.7934	4326.0607	0.0000	1.0000
SYMI	17	-2145.7169	4326.2348	0.0000	1.0000
SYMG	17	-2146.0386	4326.8782	0.0000	1.0000
K2PIG	14	-2149.6405	4327.8285	0.0000	1.0000
SYMIG	18	-2145.7141	4328.3247	0.0000	1.0000
F81I	15	-2155.0075	4340.6416	0.0000	1.0000
F81G	15	-2155.1616	4340.9499	0.0000	1.0000
F81IG	16	-2154.9858	4342.6827	0.0000	1.0000
K2P	12	-2160.9439	4346.2934	0.0000	1.0000
SYM	16	-2157.1203	4346.9517	0.0000	1.0000
F81	14	-2164.4784	4357.5044	0.0000	1.0000
JC69I	12	-2196.0465	4416.4987	0.0000	1.0000
JC69G	12	-2196.2055	4416.8168	0.0000	1.0000
JC69IG	13	-2196.0197	4418.5133	0.0000	1.0000
JC69	11	-2205.7568	4433.8564	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-2096.4878	4232.9756	0.2893	0.2893
GTRG	20	-2096.6098	4233.2196	0.2561	0.5454
HKYI	16	-2101.1789	4234.3577	0.1450	0.6904
HKYG	16	-2101.1845	4234.3689	0.1441	0.8345
GTRIG	21	-2096.4561	4234.9122	0.1099	0.9444
HKYIG	17	-2101.1373	4236.2747	0.0556	1.0000
GTR	19	-2106.6573	4251.3147	0.0000	1.0000
HKY	15	-2111.4529	4252.9059	0.0000	1.0000
K2PI	13	-2149.6927	4325.3855	0.0000	1.0000
SYMI	17	-2145.7169	4325.4338	0.0000	1.0000
K2PG	13	-2149.7934	4325.5867	0.0000	1.0000
SYMG	17	-2146.0386	4326.0771	0.0000	1.0000
K2PIG	14	-2149.6405	4327.2809	0.0000	1.0000
SYMIG	18	-2145.7141	4327.4283	0.0000	1.0000
F81I	15	-2155.0075	4340.0150	0.0000	1.0000
F81G	15	-2155.1616	4340.3232	0.0000	1.0000
F81IG	16	-2154.9858	4341.9716	0.0000	1.0000
K2P	12	-2160.9439	4345.8877	0.0000	1.0000
SYM	16	-2157.1203	4346.2406	0.0000	1.0000
F81	14	-2164.4784	4356.9568	0.0000	1.0000
JC69I	12	-2196.0465	4416.0930	0.0000	1.0000
JC69G	12	-2196.2055	4416.4111	0.0000	1.0000
JC69IG	13	-2196.0197	4418.0394	0.0000	1.0000
JC69	11	-2205.7568	4433.5135	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2101.1789	4308.9474	0.4919	0.4919
HKYG	16	-2101.1845	4308.9586	0.4891	0.9810
HKYIG	17	-2101.1373	4315.5262	0.0183	0.9994
HKY	15	-2111.4529	4322.8337	0.0005	0.9998
GTRI	20	-2096.4878	4326.2127	0.0001	0.9999
GTRG	20	-2096.6098	4326.4567	0.0001	1.0000
GTRIG	21	-2096.4561	4332.8112	0.0000	1.0000
GTR	19	-2106.6573	4339.8899	0.0000	1.0000
K2PI	13	-2149.6927	4385.9896	0.0000	1.0000
K2PG	13	-2149.7934	4386.1909	0.0000	1.0000
K2PIG	14	-2149.6405	4392.5469	0.0000	1.0000
K2P	12	-2160.9439	4401.8300	0.0000	1.0000
SYMI	17	-2145.7169	4404.6853	0.0000	1.0000
SYMG	17	-2146.0386	4405.3287	0.0000	1.0000
F81I	15	-2155.0075	4409.9428	0.0000	1.0000
F81G	15	-2155.1616	4410.2510	0.0000	1.0000
SYMIG	18	-2145.7141	4411.3416	0.0000	1.0000
F81IG	16	-2154.9858	4416.5612	0.0000	1.0000
SYM	16	-2157.1203	4420.8303	0.0000	1.0000
F81	14	-2164.4784	4422.2228	0.0000	1.0000
JC69I	12	-2196.0465	4472.0352	0.0000	1.0000
JC69G	12	-2196.2055	4472.3533	0.0000	1.0000
JC69IG	13	-2196.0197	4478.6435	0.0000	1.0000
JC69	11	-2205.7568	4484.7939	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr1_6051.phylip = 1 - 782;
 Partition Dummy = 1:683-tomodel-phylip/chr1_6051.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr1_6051.phylip = 1 - 782;
 Partition Dummy = 1:683-tomodel-phylip/chr1_6051.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,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/chr1_6051.phylip = 1 - 782;
 Partition Dummy = 1:683-tomodel-phylip/chr1_6051.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

