
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:16 08/01/2011

Input data from file "261-aggressive-phylip/13155_hg18.phylip" (29 taxa, 300 characters)

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: GTRI


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
GTRI	64	-782.0067	1727.4177	0.5287	0.5287
GTRG	64	-782.3095	1728.0233	0.3906	0.9193
GTRIG	65	-782.3093	1731.2853	0.0764	0.9957
HKYG	60	-793.9972	1738.6219	0.0020	0.9977
HKYI	60	-794.0200	1738.6677	0.0019	0.9996
HKYIG	61	-794.0038	1741.7891	0.0004	1.0000
GTR	63	-794.1926	1748.5548	0.0000	1.0000
K2PG	57	-806.4635	1754.2493	0.0000	1.0000
K2PI	57	-806.6648	1754.6520	0.0000	1.0000
HKY	59	-803.6982	1754.8964	0.0000	1.0000
SYMG	61	-800.6100	1755.0015	0.0000	1.0000
SYMI	61	-800.7817	1755.3449	0.0000	1.0000
K2PIG	58	-806.4182	1757.2348	0.0000	1.0000
SYMIG	62	-800.5982	1758.1584	0.0000	1.0000
K2P	56	-815.8531	1769.9779	0.0000	1.0000
SYM	60	-809.8555	1770.3386	0.0000	1.0000
F81G	59	-812.7644	1773.0287	0.0000	1.0000
F81I	59	-812.9309	1773.3618	0.0000	1.0000
F81IG	60	-812.7351	1776.0977	0.0000	1.0000
JC69G	56	-822.9429	1784.1573	0.0000	1.0000
JC69I	56	-823.1754	1784.6224	0.0000	1.0000
JC69IG	57	-822.8747	1787.0717	0.0000	1.0000
F81	58	-822.2047	1788.8077	0.0000	1.0000
JC69	55	-832.2863	1799.8185	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	64	-782.0067	1692.0134	0.4972	0.4972
GTRG	64	-782.3095	1692.6191	0.3673	0.8645
GTRIG	65	-782.3093	1694.6186	0.1351	0.9996
HKYG	60	-793.9972	1707.9943	0.0002	0.9998
HKYI	60	-794.0200	1708.0400	0.0002	0.9999
HKYIG	61	-794.0038	1710.0075	0.0001	1.0000
GTR	63	-794.1926	1714.3853	0.0000	1.0000
SYMG	61	-800.6100	1723.2200	0.0000	1.0000
SYMI	61	-800.7817	1723.5634	0.0000	1.0000
SYMIG	62	-800.5982	1725.1964	0.0000	1.0000
HKY	59	-803.6982	1725.3964	0.0000	1.0000
K2PG	57	-806.4635	1726.9269	0.0000	1.0000
K2PI	57	-806.6648	1727.3297	0.0000	1.0000
K2PIG	58	-806.4182	1728.8365	0.0000	1.0000
SYM	60	-809.8555	1739.7110	0.0000	1.0000
F81G	59	-812.7644	1743.5287	0.0000	1.0000
K2P	56	-815.8531	1743.7063	0.0000	1.0000
F81I	59	-812.9309	1743.8618	0.0000	1.0000
F81IG	60	-812.7351	1745.4701	0.0000	1.0000
JC69G	56	-822.9429	1757.8857	0.0000	1.0000
JC69I	56	-823.1754	1758.3508	0.0000	1.0000
JC69IG	57	-822.8747	1759.7494	0.0000	1.0000
F81	58	-822.2047	1760.4093	0.0000	1.0000
JC69	55	-832.2863	1774.5726	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRI	64	-782.0067	1929.0555	0.3402	0.3402
GTRG	64	-782.3095	1929.6611	0.2513	0.5916
HKYG	60	-793.9972	1930.2213	0.1899	0.7815
HKYI	60	-794.0200	1930.2670	0.1857	0.9672
GTRIG	65	-782.3093	1935.3645	0.0145	0.9817
HKYIG	61	-794.0038	1935.9383	0.0109	0.9926
K2PG	57	-806.4635	1938.0425	0.0038	0.9964
K2PI	57	-806.6648	1938.4453	0.0031	0.9995
K2PIG	58	-806.4182	1943.6558	0.0002	0.9997
HKY	59	-803.6982	1943.9195	0.0002	0.9999
GTR	63	-794.1926	1947.7236	0.0000	1.0000
SYMG	61	-800.6100	1949.1508	0.0000	1.0000
SYMI	61	-800.7817	1949.4941	0.0000	1.0000
K2P	56	-815.8531	1951.1181	0.0000	1.0000
SYMIG	62	-800.5982	1954.8309	0.0000	1.0000
SYM	60	-809.8555	1961.9379	0.0000	1.0000
F81G	59	-812.7644	1962.0519	0.0000	1.0000
F81I	59	-812.9309	1962.3850	0.0000	1.0000
JC69G	56	-822.9429	1965.2975	0.0000	1.0000
JC69I	56	-823.1754	1965.7626	0.0000	1.0000
F81IG	60	-812.7351	1967.6971	0.0000	1.0000
JC69IG	57	-822.8747	1970.8650	0.0000	1.0000
F81	58	-822.2047	1975.2287	0.0000	1.0000
JC69	55	-832.2863	1978.2806	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/13155_hg18.phylip = 1 - 300;
 Partition Dummy = 1:261-aggressive-phylip/13155_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/13155_hg18.phylip = 1 - 300;
 Partition Dummy = 1:261-aggressive-phylip/13155_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/13155_hg18.phylip = 1 - 300;
 Partition Dummy = 1:261-aggressive-phylip/13155_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

