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

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

Minimum AIC  model: GTRG
Minimum AICc model: HKYG
Minimum BIC  model: K2PG


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	60	-860.0438	1881.9161	0.3378	0.3378
HKYI	60	-860.7439	1883.3164	0.1677	0.5056
GTRG	64	-853.4975	1883.6499	0.1420	0.6475
GTRI	64	-853.9851	1884.6252	0.0872	0.7347
HKYIG	61	-859.9497	1885.3706	0.0601	0.7948
HKY	59	-863.6036	1885.4345	0.0582	0.8529
K2PG	57	-867.2490	1885.6441	0.0524	0.9053
GTR	63	-856.6804	1886.2445	0.0388	0.9441
K2PI	57	-868.1250	1887.3962	0.0218	0.9659
GTRIG	65	-853.5046	1887.4798	0.0209	0.9869
K2PIG	58	-867.2478	1889.1622	0.0090	0.9959
K2P	56	-872.0803	1891.8253	0.0024	0.9983
SYMG	61	-864.0573	1893.5859	0.0010	0.9992
SYMI	61	-864.8025	1895.0762	0.0005	0.9997
SYMIG	62	-864.0617	1897.2794	0.0002	0.9999
SYM	60	-867.9001	1897.6288	0.0001	1.0000
F81G	59	-878.6501	1915.5275	0.0000	1.0000
F81I	59	-879.4292	1917.0857	0.0000	1.0000
JC69G	56	-885.5667	1918.7981	0.0000	1.0000
F81IG	60	-878.6526	1919.1338	0.0000	1.0000
F81	58	-882.2907	1919.2480	0.0000	1.0000
JC69I	56	-886.4712	1920.6073	0.0000	1.0000
JC69IG	57	-885.5716	1922.2893	0.0000	1.0000
JC69	55	-890.1034	1924.4290	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-853.4975	1834.9949	0.4444	0.4444
GTRI	64	-853.9851	1835.9702	0.2729	0.7173
GTRIG	65	-853.5046	1837.0092	0.1623	0.8796
GTR	63	-856.6804	1839.3608	0.0501	0.9297
HKYG	60	-860.0438	1840.0876	0.0348	0.9646
HKYI	60	-860.7439	1841.4878	0.0173	0.9819
HKYIG	61	-859.9497	1841.8993	0.0141	0.9959
HKY	59	-863.6036	1845.2072	0.0027	0.9986
K2PG	57	-867.2490	1848.4980	0.0005	0.9991
SYMG	61	-864.0573	1850.1146	0.0002	0.9994
K2PI	57	-868.1250	1850.2501	0.0002	0.9996
K2PIG	58	-867.2478	1850.4956	0.0002	0.9998
SYMI	61	-864.8025	1851.6050	0.0001	0.9999
SYMIG	62	-864.0617	1852.1234	0.0001	1.0000
SYM	60	-867.9001	1855.8002	0.0000	1.0000
K2P	56	-872.0803	1856.1605	0.0000	1.0000
F81G	59	-878.6501	1875.3003	0.0000	1.0000
F81I	59	-879.4292	1876.8584	0.0000	1.0000
F81IG	60	-878.6526	1877.3052	0.0000	1.0000
F81	58	-882.2907	1880.5814	0.0000	1.0000
JC69G	56	-885.5667	1883.1333	0.0000	1.0000
JC69I	56	-886.4712	1884.9425	0.0000	1.0000
JC69IG	57	-885.5716	1885.1433	0.0000	1.0000
JC69	55	-890.1034	1890.2068	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-867.2490	2045.9365	0.4237	0.4237
K2PI	57	-868.1250	2047.6885	0.1764	0.6001
HKYG	60	-860.0438	2047.9175	0.1573	0.7574
HKYI	60	-860.7439	2049.3177	0.0781	0.8356
HKY	59	-863.6036	2049.5733	0.0688	0.9043
K2P	56	-872.0803	2050.1351	0.0519	0.9562
K2PIG	58	-867.2478	2051.3978	0.0276	0.9838
HKYIG	61	-859.9497	2053.1930	0.0113	0.9951
GTRG	64	-853.4975	2056.6802	0.0020	0.9971
GTR	63	-856.6804	2057.5822	0.0013	0.9983
GTRI	64	-853.9851	2057.6555	0.0012	0.9995
SYMG	61	-864.0573	2061.4084	0.0002	0.9997
GTRIG	65	-853.5046	2062.1582	0.0001	0.9998
SYMI	61	-864.8025	2062.8987	0.0001	0.9999
SYM	60	-867.9001	2063.6301	0.0001	1.0000
SYMIG	62	-864.0617	2066.8810	0.0000	1.0000
JC69G	56	-885.5667	2077.1079	0.0000	1.0000
JC69I	56	-886.4712	2078.9170	0.0000	1.0000
F81G	59	-878.6501	2079.6663	0.0000	1.0000
JC69	55	-890.1034	2080.7175	0.0000	1.0000
F81I	59	-879.4292	2081.2245	0.0000	1.0000
F81	58	-882.2907	2081.4836	0.0000	1.0000
JC69IG	57	-885.5716	2082.5817	0.0000	1.0000
F81IG	60	-878.6526	2085.1351	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chrZ_11569.phylip = 1 - 236;
 Partition Dummy = 1:183-phylip/chrZ_11569.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

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

BEGIN MRBAYES;
 Charset 183-phylip/chrZ_11569.phylip = 1 - 236;
 Partition Dummy = 1:183-phylip/chrZ_11569.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

