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

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

Minimum AIC  model: GTRG
Minimum AICc model: K2PG
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
K2PG	57	-797.5217	1745.3731	0.4689	0.4689
HKYG	60	-793.3813	1747.6565	0.1497	0.6186
K2PIG	58	-797.2354	1748.2829	0.1095	0.7281
K2PI	57	-799.3107	1748.9511	0.0784	0.8064
GTRG	64	-786.9728	1749.4885	0.0599	0.8663
SYMG	61	-792.5133	1749.5210	0.0589	0.9253
HKYIG	61	-793.0637	1750.6217	0.0340	0.9592
HKYI	60	-795.1895	1751.2729	0.0245	0.9838
GTRIG	65	-787.6747	1754.6598	0.0045	0.9883
SYMIG	62	-793.4105	1754.9566	0.0039	0.9922
GTRI	64	-789.7279	1754.9987	0.0038	0.9960
SYMI	61	-795.3101	1755.1147	0.0036	0.9996
K2P	56	-806.6465	1760.1783	0.0003	0.9999
HKY	59	-803.0191	1763.3716	0.0001	0.9999
SYM	60	-801.5612	1764.0163	0.0000	1.0000
GTR	63	-796.9325	1765.6833	0.0000	1.0000
JC69G	56	-816.9491	1780.7834	0.0000	1.0000
JC69IG	57	-816.6114	1783.5525	0.0000	1.0000
JC69I	56	-819.1577	1785.2006	0.0000	1.0000
F81G	59	-813.9585	1785.2503	0.0000	1.0000
F81IG	60	-813.6122	1788.1183	0.0000	1.0000
F81I	59	-816.1756	1789.6845	0.0000	1.0000
JC69	55	-826.3004	1796.0790	0.0000	1.0000
F81	58	-823.4516	1800.7154	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-786.9728	1701.9457	0.6496	0.6496
GTRIG	65	-787.6747	1705.3494	0.1185	0.7681
HKYG	60	-793.3813	1706.7626	0.0584	0.8265
SYMG	61	-792.5133	1707.0267	0.0512	0.8777
GTRI	64	-789.7279	1707.4558	0.0413	0.9190
HKYIG	61	-793.0637	1708.1273	0.0295	0.9485
K2PG	57	-797.5217	1709.0435	0.0187	0.9672
HKYI	60	-795.1895	1710.3790	0.0096	0.9768
K2PIG	58	-797.2354	1710.4708	0.0092	0.9860
SYMIG	62	-793.4105	1710.8210	0.0077	0.9936
SYMI	61	-795.3101	1712.6203	0.0031	0.9968
K2PI	57	-799.3107	1712.6215	0.0031	0.9999
GTR	63	-796.9325	1719.8651	0.0001	1.0000
SYM	60	-801.5612	1723.1224	0.0000	1.0000
HKY	59	-803.0191	1724.0382	0.0000	1.0000
K2P	56	-806.6465	1725.2930	0.0000	1.0000
JC69G	56	-816.9491	1745.8982	0.0000	1.0000
F81G	59	-813.9585	1745.9169	0.0000	1.0000
JC69IG	57	-816.6114	1747.2228	0.0000	1.0000
F81IG	60	-813.6122	1747.2245	0.0000	1.0000
JC69I	56	-819.1577	1750.3154	0.0000	1.0000
F81I	59	-816.1756	1750.3512	0.0000	1.0000
JC69	55	-826.3004	1762.6008	0.0000	1.0000
F81	58	-823.4516	1762.9032	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-797.5217	1907.4399	0.7819	0.7819
K2PI	57	-799.3107	1911.0179	0.1307	0.9125
K2PIG	58	-797.2354	1912.3478	0.0672	0.9798
HKYG	60	-793.3813	1915.6010	0.0132	0.9930
HKYI	60	-795.1895	1919.2173	0.0022	0.9951
SYMG	61	-792.5133	1919.3456	0.0020	0.9972
K2P	56	-806.6465	1920.2088	0.0013	0.9985
HKYIG	61	-793.0637	1920.4463	0.0012	0.9997
GTRG	64	-786.9728	1924.7066	0.0001	0.9998
SYMI	61	-795.3101	1924.9393	0.0001	0.9999
SYMIG	62	-793.4105	1926.6207	0.0001	1.0000
HKY	59	-803.0191	1929.3959	0.0000	1.0000
GTRI	64	-789.7279	1930.2167	0.0000	1.0000
GTRIG	65	-787.6747	1931.5910	0.0000	1.0000
SYM	60	-801.5612	1931.9608	0.0000	1.0000
GTR	63	-796.9325	1939.1454	0.0000	1.0000
JC69G	56	-816.9491	1940.8139	0.0000	1.0000
JC69I	56	-819.1577	1945.2312	0.0000	1.0000
JC69IG	57	-816.6114	1945.6192	0.0000	1.0000
F81G	59	-813.9585	1951.2746	0.0000	1.0000
JC69	55	-826.3004	1954.0359	0.0000	1.0000
F81I	59	-816.1756	1955.7089	0.0000	1.0000
F81IG	60	-813.6122	1956.0628	0.0000	1.0000
F81	58	-823.4516	1964.7803	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/chr6_9631.phylip = 1 - 240;
 Partition Dummy = 1:183-phylip/chr6_9631.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 (K2PG)]

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


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

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

