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

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

Minimum AIC  model: GTRG
Minimum AICc model: K2P
Minimum BIC  model: K2P


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
K2P	56	-821.3751	1787.6574	0.3969	0.3969
K2PG	57	-820.3277	1788.9145	0.2117	0.6086
K2PI	57	-820.3513	1788.9617	0.2068	0.8153
K2PIG	58	-820.3301	1792.3060	0.0388	0.8542
SYM	60	-816.9416	1792.4096	0.0369	0.8910
SYMG	61	-815.6467	1793.3145	0.0235	0.9145
HKY	59	-819.1700	1793.4080	0.0224	0.9369
SYMI	61	-815.9224	1793.8660	0.0178	0.9547
HKYG	60	-817.9814	1794.4891	0.0130	0.9677
HKYI	60	-818.0487	1794.6238	0.0122	0.9799
GTR	63	-813.5871	1796.2973	0.0053	0.9852
GTRG	64	-811.9205	1796.5722	0.0046	0.9898
SYMIG	62	-815.6461	1796.8454	0.0040	0.9938
GTRI	64	-812.2846	1797.3003	0.0032	0.9970
HKYIG	61	-817.9909	1798.0029	0.0022	0.9993
GTRIG	65	-811.9201	1800.2185	0.0007	1.0000
JC69	55	-839.5043	1820.5983	0.0000	1.0000
JC69G	56	-838.2462	1821.3997	0.0000	1.0000
JC69I	56	-838.5216	1821.9504	0.0000	1.0000
JC69IG	57	-838.2457	1824.7505	0.0000	1.0000
F81	58	-837.6585	1826.9629	0.0000	1.0000
F81G	59	-836.3459	1827.7598	0.0000	1.0000
F81I	59	-836.6371	1828.3422	0.0000	1.0000
F81IG	60	-836.3454	1831.2171	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-811.9205	1751.8411	0.1934	0.1934
GTRI	64	-812.2846	1752.5691	0.1344	0.3278
GTR	63	-813.5871	1753.1743	0.0993	0.4270
SYMG	61	-815.6467	1753.2933	0.0936	0.5206
GTRIG	65	-811.9201	1753.8401	0.0712	0.5918
SYMI	61	-815.9224	1753.8448	0.0710	0.6628
SYM	60	-816.9416	1753.8833	0.0697	0.7324
K2PG	57	-820.3277	1754.6555	0.0473	0.7798
K2PI	57	-820.3513	1754.7026	0.0462	0.8260
K2P	56	-821.3751	1754.7502	0.0452	0.8712
SYMIG	62	-815.6461	1755.2922	0.0344	0.9056
HKYG	60	-817.9814	1755.9628	0.0246	0.9302
HKYI	60	-818.0487	1756.0974	0.0230	0.9533
HKY	59	-819.1700	1756.3400	0.0204	0.9736
K2PIG	58	-820.3301	1756.6602	0.0174	0.9910
HKYIG	61	-817.9909	1757.9818	0.0090	1.0000
JC69G	56	-838.2462	1788.4924	0.0000	1.0000
JC69	55	-839.5043	1789.0085	0.0000	1.0000
JC69I	56	-838.5216	1789.0432	0.0000	1.0000
JC69IG	57	-838.2457	1790.4914	0.0000	1.0000
F81G	59	-836.3459	1790.6917	0.0000	1.0000
F81I	59	-836.6371	1791.2741	0.0000	1.0000
F81	58	-837.6585	1791.3171	0.0000	1.0000
F81IG	60	-836.3454	1792.6908	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2P	56	-821.3751	1952.1756	0.7288	0.7288
K2PG	57	-820.3277	1955.6063	0.1311	0.8599
K2PI	57	-820.3513	1955.6534	0.1281	0.9880
K2PIG	58	-820.3301	1961.1364	0.0083	0.9963
HKY	59	-819.1700	1964.3417	0.0017	0.9979
SYM	60	-816.9416	1965.4105	0.0010	0.9989
HKYG	60	-817.9814	1967.4900	0.0003	0.9992
HKYI	60	-818.0487	1967.6246	0.0003	0.9996
SYMG	61	-815.6467	1968.3459	0.0002	0.9998
SYMI	61	-815.9224	1968.8974	0.0002	1.0000
HKYIG	61	-817.9909	1973.0344	0.0000	1.0000
SYMIG	62	-815.6461	1973.8703	0.0000	1.0000
GTR	63	-813.5871	1975.2778	0.0000	1.0000
GTRG	64	-811.9205	1977.4700	0.0000	1.0000
GTRI	64	-812.2846	1978.1981	0.0000	1.0000
JC69	55	-839.5043	1982.9084	0.0000	1.0000
GTRIG	65	-811.9201	1982.9946	0.0000	1.0000
JC69G	56	-838.2462	1985.9178	0.0000	1.0000
JC69I	56	-838.5216	1986.4685	0.0000	1.0000
JC69IG	57	-838.2457	1991.4423	0.0000	1.0000
F81	58	-837.6585	1995.7934	0.0000	1.0000
F81G	59	-836.3459	1998.6935	0.0000	1.0000
F81I	59	-836.6371	1999.2758	0.0000	1.0000
F81IG	60	-836.3454	2004.2180	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/chr4_13431.phylip = 1 - 251;
 Partition Dummy = 1:183-phylip/chr4_13431.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 (K2P)]

BEGIN MRBAYES;
 Charset 183-phylip/chr4_13431.phylip = 1 - 251;
 Partition Dummy = 1:183-phylip/chr4_13431.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr4_13431.phylip = 1 - 251;
 Partition Dummy = 1:183-phylip/chr4_13431.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0);
END;


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

