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

Input data from file "683-tomodel-phylip/chr2_8735.phylip" (7 taxa, 780 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	-1857.9605	3748.6340	0.4094	0.4094
HKYI	16	-1858.5148	3749.7425	0.2352	0.6445
HKYIG	17	-1857.7880	3750.3792	0.1710	0.8156
GTRG	20	-1855.2264	3751.5595	0.0948	0.9104
GTRI	20	-1855.8730	3752.8527	0.0497	0.9600
GTRIG	21	-1855.0344	3753.2878	0.0400	1.0000
HKY	15	-1871.3701	3773.3684	0.0000	1.0000
GTR	19	-1868.6746	3776.3492	0.0000	1.0000
K2PG	13	-1886.7103	3799.8957	0.0000	1.0000
K2PI	13	-1887.0176	3800.5104	0.0000	1.0000
K2PIG	14	-1886.5819	3801.7128	0.0000	1.0000
SYMG	17	-1885.5521	3805.9074	0.0000	1.0000
SYMI	17	-1885.8167	3806.4366	0.0000	1.0000
SYMIG	18	-1885.4367	3807.7722	0.0000	1.0000
K2P	12	-1901.0710	3826.5489	0.0000	1.0000
SYM	16	-1899.9672	3832.6474	0.0000	1.0000
F81G	15	-1905.9511	3842.5305	0.0000	1.0000
F81I	15	-1906.3008	3843.2299	0.0000	1.0000
F81IG	16	-1905.8513	3844.4157	0.0000	1.0000
F81	14	-1919.0045	3866.5581	0.0000	1.0000
JC69G	12	-1931.4618	3887.3305	0.0000	1.0000
JC69I	12	-1931.7276	3887.8620	0.0000	1.0000
JC69IG	13	-1931.3731	3889.2214	0.0000	1.0000
JC69	11	-1944.9930	3912.3297	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1857.9605	3747.9210	0.3896	0.3896
HKYI	16	-1858.5148	3749.0295	0.2238	0.6134
HKYIG	17	-1857.7880	3749.5761	0.1703	0.7836
GTRG	20	-1855.2264	3750.4528	0.1098	0.8935
GTRI	20	-1855.8730	3751.7460	0.0575	0.9510
GTRIG	21	-1855.0344	3752.0688	0.0490	1.0000
HKY	15	-1871.3701	3772.7402	0.0000	1.0000
GTR	19	-1868.6746	3775.3492	0.0000	1.0000
K2PG	13	-1886.7103	3799.4205	0.0000	1.0000
K2PI	13	-1887.0176	3800.0352	0.0000	1.0000
K2PIG	14	-1886.5819	3801.1638	0.0000	1.0000
SYMG	17	-1885.5521	3805.1043	0.0000	1.0000
SYMI	17	-1885.8167	3805.6335	0.0000	1.0000
SYMIG	18	-1885.4367	3806.8734	0.0000	1.0000
K2P	12	-1901.0710	3826.1421	0.0000	1.0000
SYM	16	-1899.9672	3831.9344	0.0000	1.0000
F81G	15	-1905.9511	3841.9022	0.0000	1.0000
F81I	15	-1906.3008	3842.6017	0.0000	1.0000
F81IG	16	-1905.8513	3843.7027	0.0000	1.0000
F81	14	-1919.0045	3866.0091	0.0000	1.0000
JC69G	12	-1931.4618	3886.9237	0.0000	1.0000
JC69I	12	-1931.7276	3887.4553	0.0000	1.0000
JC69IG	13	-1931.3731	3888.7462	0.0000	1.0000
JC69	11	-1944.9930	3911.9860	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1857.9605	3822.4697	0.6184	0.6184
HKYI	16	-1858.5148	3823.5782	0.3553	0.9736
HKYIG	17	-1857.7880	3828.7841	0.0263	0.9999
HKY	15	-1871.3701	3842.6296	0.0000	1.0000
GTRG	20	-1855.2264	3843.6387	0.0000	1.0000
GTRI	20	-1855.8730	3844.9318	0.0000	1.0000
GTRIG	21	-1855.0344	3849.9140	0.0000	1.0000
K2PG	13	-1886.7103	3859.9913	0.0000	1.0000
K2PI	13	-1887.0176	3860.6061	0.0000	1.0000
GTR	19	-1868.6746	3863.8758	0.0000	1.0000
K2PIG	14	-1886.5819	3866.3939	0.0000	1.0000
K2P	12	-1901.0710	3882.0536	0.0000	1.0000
SYMG	17	-1885.5521	3884.3123	0.0000	1.0000
SYMI	17	-1885.8167	3884.8415	0.0000	1.0000
SYMIG	18	-1885.4367	3890.7407	0.0000	1.0000
SYM	16	-1899.9672	3906.4831	0.0000	1.0000
F81G	15	-1905.9511	3911.7916	0.0000	1.0000
F81I	15	-1906.3008	3912.4911	0.0000	1.0000
F81IG	16	-1905.8513	3918.2514	0.0000	1.0000
F81	14	-1919.0045	3931.2392	0.0000	1.0000
JC69G	12	-1931.4618	3942.8352	0.0000	1.0000
JC69I	12	-1931.7276	3943.3668	0.0000	1.0000
JC69IG	13	-1931.3731	3949.3170	0.0000	1.0000
JC69	11	-1944.9930	3963.2382	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/chr2_8735.phylip = 1 - 780;
 Partition Dummy = 1:683-tomodel-phylip/chr2_8735.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/chr2_8735.phylip = 1 - 780;
 Partition Dummy = 1:683-tomodel-phylip/chr2_8735.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/chr2_8735.phylip = 1 - 780;
 Partition Dummy = 1:683-tomodel-phylip/chr2_8735.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

