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

Input data from file "591-tomodel-phylip/chr2_22050.phylip" (7 taxa, 316 characters)

Minimum AIC  model: GTRI
Minimum AICc model: HKYI
Minimum BIC  model: HKY


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
HKYI	16	-889.4505	1812.7205	0.1972	0.1972
HKYG	16	-889.5461	1812.9116	0.1792	0.3764
GTRI	20	-885.0410	1812.9294	0.1776	0.5540
GTRG	20	-885.0431	1812.9338	0.1772	0.7313
HKY	15	-891.3487	1814.2973	0.0896	0.8209
HKYIG	17	-889.4523	1814.9583	0.0644	0.8853
GTR	19	-887.2925	1815.1526	0.0584	0.9438
GTRIG	21	-885.0437	1815.2303	0.0562	1.0000
K2PI	13	-905.8138	1838.8328	0.0000	1.0000
K2PG	13	-905.8935	1838.9924	0.0000	1.0000
SYMG	17	-901.5867	1839.2271	0.0000	1.0000
SYMI	17	-901.6408	1839.3353	0.0000	1.0000
K2P	12	-907.4469	1839.9236	0.0000	1.0000
SYM	16	-903.5326	1840.8847	0.0000	1.0000
K2PIG	14	-905.8197	1841.0347	0.0000	1.0000
SYMIG	18	-901.5901	1841.4831	0.0000	1.0000
F81	14	-919.8218	1869.0390	0.0000	1.0000
F81I	15	-918.7271	1869.0542	0.0000	1.0000
F81G	15	-918.8246	1869.2492	0.0000	1.0000
F81IG	16	-918.7293	1871.2780	0.0000	1.0000
JC69	11	-933.3316	1889.5317	0.0000	1.0000
JC69I	12	-932.2841	1889.5980	0.0000	1.0000
JC69G	12	-932.3597	1889.7491	0.0000	1.0000
JC69IG	13	-932.2863	1891.7778	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-885.0410	1810.0819	0.2256	0.2256
GTRG	20	-885.0431	1810.0863	0.2251	0.4507
HKYI	16	-889.4505	1810.9011	0.1498	0.6005
HKYG	16	-889.5461	1811.0922	0.1361	0.7367
GTRIG	21	-885.0437	1812.0875	0.0828	0.8194
GTR	19	-887.2925	1812.5851	0.0645	0.8840
HKY	15	-891.3487	1812.6973	0.0610	0.9450
HKYIG	17	-889.4523	1812.9046	0.0550	1.0000
SYMG	17	-901.5867	1837.1734	0.0000	1.0000
SYMI	17	-901.6408	1837.2816	0.0000	1.0000
K2PI	13	-905.8138	1837.6275	0.0000	1.0000
K2PG	13	-905.8935	1837.7871	0.0000	1.0000
K2P	12	-907.4469	1838.8939	0.0000	1.0000
SYM	16	-903.5326	1839.0653	0.0000	1.0000
SYMIG	18	-901.5901	1839.1801	0.0000	1.0000
K2PIG	14	-905.8197	1839.6394	0.0000	1.0000
F81I	15	-918.7271	1867.4542	0.0000	1.0000
F81	14	-919.8218	1867.6437	0.0000	1.0000
F81G	15	-918.8246	1867.6492	0.0000	1.0000
F81IG	16	-918.7293	1869.4586	0.0000	1.0000
JC69I	12	-932.2841	1888.5683	0.0000	1.0000
JC69	11	-933.3316	1888.6633	0.0000	1.0000
JC69G	12	-932.3597	1888.7194	0.0000	1.0000
JC69IG	13	-932.2863	1890.5725	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKY	15	-891.3487	1869.0334	0.5748	0.5748
HKYI	16	-889.4505	1870.9929	0.2158	0.7906
HKYG	16	-889.5461	1871.1841	0.1961	0.9867
HKYIG	17	-889.4523	1876.7522	0.0121	0.9988
GTR	19	-887.2925	1883.9442	0.0003	0.9991
K2P	12	-907.4469	1883.9628	0.0003	0.9994
GTRI	20	-885.0410	1885.1968	0.0002	0.9996
GTRG	20	-885.0431	1885.2011	0.0002	0.9998
K2PI	13	-905.8138	1886.4522	0.0001	0.9999
K2PG	13	-905.8935	1886.6117	0.0001	1.0000
GTRIG	21	-885.0437	1890.9581	0.0000	1.0000
K2PIG	14	-905.8197	1892.2198	0.0000	1.0000
SYM	16	-903.5326	1899.1572	0.0000	1.0000
SYMG	17	-901.5867	1901.0210	0.0000	1.0000
SYMI	17	-901.6408	1901.1292	0.0000	1.0000
SYMIG	18	-901.5901	1906.7835	0.0000	1.0000
F81	14	-919.8218	1920.2241	0.0000	1.0000
F81I	15	-918.7271	1923.7903	0.0000	1.0000
F81G	15	-918.8246	1923.9854	0.0000	1.0000
F81IG	16	-918.7293	1929.5504	0.0000	1.0000
JC69	11	-933.3316	1929.9764	0.0000	1.0000
JC69I	12	-932.2841	1933.6372	0.0000	1.0000
JC69G	12	-932.3597	1933.7883	0.0000	1.0000
JC69IG	13	-932.2863	1939.3972	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 591-tomodel-phylip/chr2_22050.phylip = 1 - 316;
 Partition Dummy = 1:591-tomodel-phylip/chr2_22050.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 (HKYI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_22050.phylip = 1 - 316;
 Partition Dummy = 1:591-tomodel-phylip/chr2_22050.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;


[Mrbayes block for the best BIC model (HKY)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_22050.phylip = 1 - 316;
 Partition Dummy = 1:591-tomodel-phylip/chr2_22050.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


-------------------------------------------------------------
End of Output

