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

Input data from file "591-tomodel-phylip/chr7_11511.phylip" (7 taxa, 620 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	-1862.9711	3758.8444	0.4414	0.4414
HKYI	16	-1863.1953	3759.2928	0.3528	0.7942
HKYIG	17	-1862.9763	3760.9693	0.1526	0.9467
GTRG	20	-1861.5788	3764.5600	0.0253	0.9720
GTRI	20	-1861.8604	3765.1232	0.0191	0.9912
GTRIG	21	-1861.5836	3766.7124	0.0086	0.9998
HKY	15	-1871.7558	3774.3063	0.0002	1.0000
GTR	19	-1870.4285	3780.1236	0.0000	1.0000
K2PG	13	-1889.6671	3805.9348	0.0000	1.0000
K2PI	13	-1889.8309	3806.2625	0.0000	1.0000
K2PIG	14	-1889.6803	3808.0548	0.0000	1.0000
SYMG	17	-1888.3528	3811.7221	0.0000	1.0000
SYMI	17	-1888.5250	3812.0665	0.0000	1.0000
SYMIG	18	-1888.3707	3813.8795	0.0000	1.0000
K2P	12	-1898.0832	3820.6805	0.0000	1.0000
SYM	16	-1896.7854	3826.4729	0.0000	1.0000
F81G	15	-1901.9078	3834.6104	0.0000	1.0000
F81I	15	-1902.0750	3834.9448	0.0000	1.0000
F81IG	16	-1901.9177	3836.7375	0.0000	1.0000
F81	14	-1910.3040	3849.3023	0.0000	1.0000
JC69G	12	-1923.8702	3872.2545	0.0000	1.0000
JC69I	12	-1924.0204	3872.5547	0.0000	1.0000
JC69IG	13	-1923.8739	3874.3485	0.0000	1.0000
JC69	11	-1932.2636	3886.9615	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1862.9711	3757.9423	0.4307	0.4307
HKYI	16	-1863.1953	3758.3906	0.3442	0.7749
HKYIG	17	-1862.9763	3759.9526	0.1576	0.9325
GTRG	20	-1861.5788	3763.1577	0.0317	0.9642
GTRI	20	-1861.8604	3763.7209	0.0240	0.9882
GTRIG	21	-1861.5836	3765.1672	0.0116	0.9998
HKY	15	-1871.7558	3773.5116	0.0002	1.0000
GTR	19	-1870.4285	3778.8570	0.0000	1.0000
K2PG	13	-1889.6671	3805.3342	0.0000	1.0000
K2PI	13	-1889.8309	3805.6618	0.0000	1.0000
K2PIG	14	-1889.6803	3807.3606	0.0000	1.0000
SYMG	17	-1888.3528	3810.7055	0.0000	1.0000
SYMI	17	-1888.5250	3811.0499	0.0000	1.0000
SYMIG	18	-1888.3707	3812.7414	0.0000	1.0000
K2P	12	-1898.0832	3820.1665	0.0000	1.0000
SYM	16	-1896.7854	3825.5708	0.0000	1.0000
F81G	15	-1901.9078	3833.8157	0.0000	1.0000
F81I	15	-1902.0750	3834.1501	0.0000	1.0000
F81IG	16	-1901.9177	3835.8353	0.0000	1.0000
F81	14	-1910.3040	3848.6080	0.0000	1.0000
JC69G	12	-1923.8702	3871.7405	0.0000	1.0000
JC69I	12	-1924.0204	3872.0407	0.0000	1.0000
JC69IG	13	-1923.8739	3873.7479	0.0000	1.0000
JC69	11	-1932.2636	3886.5273	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1862.9711	3828.8178	0.5426	0.5426
HKYI	16	-1863.1953	3829.2661	0.4336	0.9762
HKYIG	17	-1862.9763	3835.2579	0.0217	0.9979
HKY	15	-1871.7558	3839.9574	0.0021	1.0000
GTRG	20	-1861.5788	3851.7520	0.0000	1.0000
GTRI	20	-1861.8604	3852.3153	0.0000	1.0000
GTRIG	21	-1861.5836	3858.1913	0.0000	1.0000
K2PG	13	-1889.6671	3862.9205	0.0000	1.0000
GTR	19	-1870.4285	3863.0217	0.0000	1.0000
K2PI	13	-1889.8309	3863.2482	0.0000	1.0000
K2PIG	14	-1889.6803	3869.3766	0.0000	1.0000
K2P	12	-1898.0832	3873.3231	0.0000	1.0000
SYMG	17	-1888.3528	3886.0108	0.0000	1.0000
SYMI	17	-1888.5250	3886.3552	0.0000	1.0000
SYMIG	18	-1888.3707	3892.4764	0.0000	1.0000
SYM	16	-1896.7854	3896.4463	0.0000	1.0000
F81G	15	-1901.9078	3900.2615	0.0000	1.0000
F81I	15	-1902.0750	3900.5959	0.0000	1.0000
F81IG	16	-1901.9177	3906.7108	0.0000	1.0000
F81	14	-1910.3040	3910.6241	0.0000	1.0000
JC69G	12	-1923.8702	3924.8971	0.0000	1.0000
JC69I	12	-1924.0204	3925.1974	0.0000	1.0000
JC69IG	13	-1923.8739	3931.3342	0.0000	1.0000
JC69	11	-1932.2636	3935.2542	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_11511.phylip = 1 - 620;
 Partition Dummy = 1:591-tomodel-phylip/chr7_11511.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_11511.phylip = 1 - 620;
 Partition Dummy = 1:591-tomodel-phylip/chr7_11511.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_11511.phylip = 1 - 620;
 Partition Dummy = 1:591-tomodel-phylip/chr7_11511.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

