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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: JC69G


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	-771.8411	1708.0458	0.2850	0.2850
K2PG	57	-777.4403	1708.2377	0.2589	0.5439
JC69G	56	-779.8545	1709.4841	0.1388	0.6828
K2PI	57	-778.5524	1710.4618	0.0852	0.7679
F81G	59	-775.3463	1711.3433	0.0548	0.8227
JC69I	56	-780.9519	1711.6790	0.0463	0.8691
K2PIG	58	-777.3844	1711.7509	0.0447	0.9138
JC69IG	57	-779.8083	1712.9737	0.0243	0.9380
HKYI	60	-774.4842	1713.3321	0.0203	0.9583
F81I	59	-776.3704	1713.3914	0.0197	0.9780
HKYIG	61	-773.2657	1714.6533	0.0105	0.9884
F81IG	60	-775.2988	1714.9612	0.0090	0.9974
GTRG	64	-770.5085	1720.6939	0.0005	0.9979
GTRI	64	-770.7596	1721.1963	0.0004	0.9983
SYMG	61	-776.5450	1721.2120	0.0004	0.9987
K2P	56	-785.8411	1721.4573	0.0003	0.9991
HKY	59	-780.5188	1721.6883	0.0003	0.9994
JC69	55	-788.0698	1722.3749	0.0002	0.9996
SYMI	61	-777.5547	1723.2314	0.0001	0.9997
GTRIG	65	-770.0033	1723.6315	0.0001	0.9999
F81	58	-783.7633	1724.5086	0.0001	0.9999
SYMIG	62	-776.4949	1724.9161	0.0001	1.0000
SYM	60	-785.1040	1734.5717	0.0000	1.0000
GTR	63	-779.5751	1734.9281	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-771.8411	1663.6822	0.6123	0.6123
HKYIG	61	-773.2657	1668.5313	0.0542	0.6665
F81G	59	-775.3463	1668.6927	0.0500	0.7165
K2PG	57	-777.4403	1668.8806	0.0455	0.7620
HKYI	60	-774.4842	1668.9685	0.0436	0.8056
GTRG	64	-770.5085	1669.0169	0.0425	0.8481
GTRI	64	-770.7596	1669.5193	0.0331	0.8812
GTRIG	65	-770.0033	1670.0065	0.0259	0.9071
F81IG	60	-775.2988	1670.5976	0.0193	0.9264
F81I	59	-776.3704	1670.7408	0.0180	0.9444
K2PIG	58	-777.3844	1670.7688	0.0177	0.9621
K2PI	57	-778.5524	1671.1047	0.0150	0.9770
JC69G	56	-779.8545	1671.7090	0.0111	0.9881
JC69IG	57	-779.8083	1673.6166	0.0043	0.9924
JC69I	56	-780.9519	1673.9039	0.0037	0.9961
SYMG	61	-776.5450	1675.0901	0.0020	0.9981
SYMIG	62	-776.4949	1676.9897	0.0008	0.9989
SYMI	61	-777.5547	1677.1095	0.0007	0.9996
HKY	59	-780.5188	1679.0377	0.0003	0.9999
F81	58	-783.7633	1683.5266	0.0000	0.9999
K2P	56	-785.8411	1683.6822	0.0000	1.0000
GTR	63	-779.5751	1685.1503	0.0000	1.0000
JC69	55	-788.0698	1686.1396	0.0000	1.0000
SYM	60	-785.1040	1690.2080	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
JC69G	56	-779.8545	1863.2590	0.3916	0.3916
K2PG	57	-777.4403	1863.8511	0.2913	0.6829
JC69I	56	-780.9519	1865.4538	0.1307	0.8136
K2PI	57	-778.5524	1866.0752	0.0958	0.9094
JC69IG	57	-779.8083	1868.5871	0.0273	0.9367
HKYG	60	-771.8411	1868.9143	0.0232	0.9599
K2PIG	58	-777.3844	1869.1598	0.0205	0.9803
F81G	59	-775.3463	1870.5042	0.0105	0.9908
F81I	59	-776.3704	1872.5524	0.0038	0.9946
HKYI	60	-774.4842	1874.2006	0.0016	0.9962
JC69	55	-788.0698	1874.2691	0.0016	0.9978
K2P	56	-785.8411	1875.2322	0.0010	0.9988
F81IG	60	-775.2988	1875.8297	0.0007	0.9995
HKYIG	61	-773.2657	1877.1839	0.0004	0.9999
HKY	59	-780.5188	1880.8492	0.0001	0.9999
F81	58	-783.7633	1881.9176	0.0000	1.0000
SYMG	61	-776.5450	1883.7427	0.0000	1.0000
SYMI	61	-777.5547	1885.7621	0.0000	1.0000
GTRG	64	-770.5085	1887.9311	0.0000	1.0000
GTRI	64	-770.7596	1888.4335	0.0000	1.0000
SYMIG	62	-776.4949	1889.0629	0.0000	1.0000
GTRIG	65	-770.0033	1892.3413	0.0000	1.0000
SYM	60	-785.1040	1895.4401	0.0000	1.0000
GTR	63	-779.5751	1900.6440	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 183-phylip/chrZ_11557.phylip = 1 - 226;
 Partition Dummy = 1:183-phylip/chrZ_11557.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 183-phylip/chrZ_11557.phylip = 1 - 226;
 Partition Dummy = 1:183-phylip/chrZ_11557.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 (JC69G)]

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


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

