
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
00:52 02/17/2011

Input data from file "183-phylip/chr5_12015.phylip" (29 taxa, 300 characters)

Minimum AIC  model: HKYIG
Minimum AICc model: HKYI
Minimum BIC  model: K2PI


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	60	-854.5259	1859.6794	0.3744	0.3744
HKYIG	61	-853.3179	1860.4172	0.2589	0.6332
HKYG	60	-854.9429	1860.5134	0.2467	0.8799
K2PI	57	-861.0471	1863.4166	0.0578	0.9377
K2PIG	58	-860.2820	1864.9623	0.0267	0.9644
K2PG	57	-861.9023	1865.1268	0.0246	0.9890
GTRI	64	-852.5552	1868.5146	0.0045	0.9935
GTRG	64	-853.1001	1869.6044	0.0026	0.9961
GTRIG	65	-851.6356	1869.9379	0.0022	0.9983
SYMI	61	-858.9675	1871.7166	0.0009	0.9992
SYMIG	62	-858.1451	1873.2523	0.0004	0.9996
SYMG	61	-859.8874	1873.5563	0.0004	1.0000
F81I	59	-873.6332	1894.7663	0.0000	1.0000
F81G	59	-874.4245	1896.3490	0.0000	1.0000
F81IG	60	-873.5273	1897.6821	0.0000	1.0000
JC69I	56	-880.8838	1900.0392	0.0000	1.0000
JC69G	56	-881.7106	1901.6929	0.0000	1.0000
JC69IG	57	-880.7192	1902.7608	0.0000	1.0000
HKY	59	-888.1426	1923.7852	0.0000	1.0000
K2P	56	-895.5581	1929.3877	0.0000	1.0000
GTR	63	-886.9080	1933.9854	0.0000	1.0000
SYM	60	-894.0040	1938.6355	0.0000	1.0000
F81	58	-906.9646	1958.3276	0.0000	1.0000
JC69	55	-914.1601	1963.5661	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYIG	61	-853.3179	1828.6357	0.3737	0.3737
HKYI	60	-854.5259	1829.0518	0.3035	0.6773
HKYG	60	-854.9429	1829.8858	0.2000	0.8773
GTRI	64	-852.5552	1833.1103	0.0399	0.9172
GTRIG	65	-851.6356	1833.2713	0.0368	0.9540
GTRG	64	-853.1001	1834.2001	0.0231	0.9772
K2PI	57	-861.0471	1836.0943	0.0090	0.9861
K2PIG	58	-860.2820	1836.5640	0.0071	0.9932
K2PG	57	-861.9023	1837.8045	0.0038	0.9971
SYMI	61	-858.9675	1839.9351	0.0013	0.9984
SYMIG	62	-858.1451	1840.2902	0.0011	0.9995
SYMG	61	-859.8874	1841.7748	0.0005	1.0000
F81I	59	-873.6332	1865.2663	0.0000	1.0000
F81G	59	-874.4245	1866.8490	0.0000	1.0000
F81IG	60	-873.5273	1867.0545	0.0000	1.0000
JC69I	56	-880.8838	1873.7676	0.0000	1.0000
JC69G	56	-881.7106	1875.4213	0.0000	1.0000
JC69IG	57	-880.7192	1875.4385	0.0000	1.0000
HKY	59	-888.1426	1894.2852	0.0000	1.0000
GTR	63	-886.9080	1899.8159	0.0000	1.0000
K2P	56	-895.5581	1903.1161	0.0000	1.0000
SYM	60	-894.0040	1908.0079	0.0000	1.0000
F81	58	-906.9646	1929.9292	0.0000	1.0000
JC69	55	-914.1601	1938.3202	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	57	-861.0471	2047.2099	0.5581	0.5581
K2PG	57	-861.9023	2048.9201	0.2373	0.7955
HKYI	60	-854.5259	2051.2787	0.0730	0.8685
K2PIG	58	-860.2820	2051.3834	0.0693	0.9377
HKYG	60	-854.9429	2052.1128	0.0481	0.9858
HKYIG	61	-853.3179	2054.5664	0.0141	0.9999
SYMI	61	-858.9675	2065.8658	0.0000	1.0000
SYMG	61	-859.8874	2067.7055	0.0000	1.0000
SYMIG	62	-858.1451	2069.9248	0.0000	1.0000
GTRI	64	-852.5552	2070.1524	0.0000	1.0000
GTRG	64	-853.1001	2071.2422	0.0000	1.0000
GTRIG	65	-851.6356	2074.0171	0.0000	1.0000
JC69I	56	-880.8838	2081.1795	0.0000	1.0000
JC69G	56	-881.7106	2082.8331	0.0000	1.0000
F81I	59	-873.6332	2083.7895	0.0000	1.0000
F81G	59	-874.4245	2085.3721	0.0000	1.0000
JC69IG	57	-880.7192	2086.5541	0.0000	1.0000
F81IG	60	-873.5273	2089.2815	0.0000	1.0000
K2P	56	-895.5581	2110.5280	0.0000	1.0000
HKY	59	-888.1426	2112.8083	0.0000	1.0000
SYM	60	-894.0040	2130.2349	0.0000	1.0000
GTR	63	-886.9080	2133.1542	0.0000	1.0000
JC69	55	-914.1601	2142.0282	0.0000	1.0000
F81	58	-906.9646	2144.7486	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr5_12015.phylip = 1 - 300;
 Partition Dummy = 1:183-phylip/chr5_12015.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr5_12015.phylip = 1 - 300;
 Partition Dummy = 1:183-phylip/chr5_12015.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 (K2PI)]

BEGIN MRBAYES;
 Charset 183-phylip/chr5_12015.phylip = 1 - 300;
 Partition Dummy = 1:183-phylip/chr5_12015.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

