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

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

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


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	-691.9662	1570.4779	0.5797	0.5797
HKYI	60	-692.8666	1572.2786	0.2356	0.8152
HKYIG	61	-691.9689	1575.3323	0.0512	0.8664
GTRG	64	-684.6737	1575.8380	0.0397	0.9061
HKY	59	-697.0972	1575.9782	0.0370	0.9432
K2PG	57	-702.2161	1576.9455	0.0228	0.9660
K2PI	57	-702.5971	1577.7074	0.0156	0.9816
GTRI	64	-686.1263	1578.7431	0.0093	0.9909
GTRIG	65	-684.6730	1581.0604	0.0029	0.9939
GTR	63	-689.9374	1581.2393	0.0027	0.9965
K2PIG	58	-702.1929	1581.4929	0.0024	0.9989
K2P	56	-708.0382	1584.0763	0.0006	0.9995
SYMG	61	-697.2009	1585.7962	0.0003	0.9998
SYMI	61	-697.7437	1586.8819	0.0002	1.0000
SYMIG	62	-697.1958	1590.7249	0.0000	1.0000
SYM	60	-702.1473	1590.8400	0.0000	1.0000
F81G	59	-719.9049	1621.5935	0.0000	1.0000
F81I	59	-720.3514	1622.4867	0.0000	1.0000
F81IG	60	-719.8963	1626.3380	0.0000	1.0000
JC69G	56	-730.3110	1628.6219	0.0000	1.0000
F81	58	-725.7957	1628.6986	0.0000	1.0000
JC69I	56	-730.5737	1629.1474	0.0000	1.0000
JC69IG	57	-730.2822	1633.0777	0.0000	1.0000
JC69	55	-735.8456	1635.2564	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-684.6737	1497.3474	0.5942	0.5942
GTRIG	65	-684.6730	1499.3461	0.2187	0.8130
GTRI	64	-686.1263	1500.2525	0.1390	0.9520
HKYG	60	-691.9662	1503.9324	0.0221	0.9741
HKYI	60	-692.8666	1505.7332	0.0090	0.9831
GTR	63	-689.9374	1505.8748	0.0084	0.9914
HKYIG	61	-691.9689	1505.9378	0.0081	0.9995
HKY	59	-697.0972	1512.1944	0.0004	0.9999
SYMG	61	-697.2009	1516.4017	0.0000	0.9999
SYMI	61	-697.7437	1517.4874	0.0000	1.0000
SYMIG	62	-697.1958	1518.3916	0.0000	1.0000
K2PG	57	-702.2161	1518.4323	0.0000	1.0000
K2PI	57	-702.5971	1519.1942	0.0000	1.0000
K2PIG	58	-702.1929	1520.3857	0.0000	1.0000
SYM	60	-702.1473	1524.2945	0.0000	1.0000
K2P	56	-708.0382	1528.0763	0.0000	1.0000
F81G	59	-719.9049	1557.8097	0.0000	1.0000
F81I	59	-720.3514	1558.7029	0.0000	1.0000
F81IG	60	-719.8963	1559.7926	0.0000	1.0000
F81	58	-725.7957	1567.5915	0.0000	1.0000
JC69G	56	-730.3110	1572.6219	0.0000	1.0000
JC69I	56	-730.5737	1573.1474	0.0000	1.0000
JC69IG	57	-730.2822	1574.5644	0.0000	1.0000
JC69	55	-735.8456	1581.6911	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-691.9662	1692.4322	0.5642	0.5642
HKYI	60	-692.8666	1694.2330	0.2293	0.7936
K2PG	57	-702.2161	1697.5071	0.0446	0.8382
HKY	59	-697.0972	1697.5526	0.0436	0.8818
HKYIG	61	-691.9689	1697.5793	0.0430	0.9248
K2PI	57	-702.5971	1698.2690	0.0305	0.9553
GTRG	64	-684.6737	1698.4139	0.0284	0.9837
GTRI	64	-686.1263	1701.3190	0.0066	0.9903
K2PIG	58	-702.1929	1702.6022	0.0035	0.9938
GTRIG	65	-684.6730	1703.5542	0.0022	0.9960
GTR	63	-689.9374	1703.7996	0.0019	0.9979
K2P	56	-708.0382	1704.0095	0.0017	0.9996
SYMG	61	-697.2009	1708.0432	0.0002	0.9998
SYMI	61	-697.7437	1709.1288	0.0001	1.0000
SYM	60	-702.1473	1712.7944	0.0000	1.0000
SYMIG	62	-697.1958	1713.1747	0.0000	1.0000
F81G	59	-719.9049	1743.1679	0.0000	1.0000
F81I	59	-720.3514	1744.0610	0.0000	1.0000
F81IG	60	-719.8963	1748.2924	0.0000	1.0000
JC69G	56	-730.3110	1748.5551	0.0000	1.0000
JC69I	56	-730.5737	1749.0806	0.0000	1.0000
F81	58	-725.7957	1749.8080	0.0000	1.0000
JC69IG	57	-730.2822	1753.6392	0.0000	1.0000
JC69	55	-735.8456	1754.4826	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/chr3_3107.phylip = 1 - 171;
 Partition Dummy = 1:183-phylip/chr3_3107.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/chr3_3107.phylip = 1 - 171;
 Partition Dummy = 1:183-phylip/chr3_3107.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 (HKYG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr3_3107.phylip = 1 - 171;
 Partition Dummy = 1:183-phylip/chr3_3107.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;


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

