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

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

Minimum AIC  model: HKYG
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	-850.1396	1863.5929	0.5918	0.5918
HKYI	60	-850.8239	1864.9613	0.2985	0.8903
HKYIG	61	-850.0473	1867.1183	0.1015	0.9918
HKY	59	-856.4774	1872.6019	0.0065	0.9984
GTRG	64	-848.9780	1876.3802	0.0010	0.9994
GTRI	64	-849.7294	1877.8830	0.0005	0.9998
GTRIG	65	-848.9012	1880.1194	0.0002	1.0000
GTR	63	-855.1280	1884.8343	0.0000	1.0000
K2PG	57	-874.2650	1900.9718	0.0000	1.0000
K2PI	57	-875.0416	1902.5250	0.0000	1.0000
K2PIG	58	-874.1654	1904.3541	0.0000	1.0000
F81G	59	-873.0034	1905.6538	0.0000	1.0000
F81I	59	-873.6459	1906.9389	0.0000	1.0000
K2P	56	-879.3030	1907.5076	0.0000	1.0000
SYMG	61	-870.5823	1908.1883	0.0000	1.0000
F81IG	60	-872.9086	1909.1307	0.0000	1.0000
SYMI	61	-872.6301	1912.2841	0.0000	1.0000
SYMIG	62	-871.4668	1913.7120	0.0000	1.0000
F81	58	-879.2359	1914.4951	0.0000	1.0000
SYM	60	-876.8121	1916.9379	0.0000	1.0000
JC69G	56	-891.2136	1931.3290	0.0000	1.0000
JC69I	56	-893.8019	1936.5056	0.0000	1.0000
JC69IG	57	-892.9695	1938.3809	0.0000	1.0000
JC69	55	-896.9913	1939.3848	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-850.1396	1820.2793	0.4944	0.4944
HKYI	60	-850.8239	1821.6477	0.2494	0.7439
HKYIG	61	-850.0473	1822.0945	0.1995	0.9434
GTRG	64	-848.9780	1825.9559	0.0289	0.9723
GTRI	64	-849.7294	1827.4588	0.0136	0.9860
GTRIG	65	-848.9012	1827.8023	0.0115	0.9975
HKY	59	-856.4774	1830.9549	0.0024	0.9998
GTR	63	-855.1280	1836.2559	0.0002	1.0000
K2PG	57	-874.2650	1862.5300	0.0000	1.0000
SYMG	61	-870.5823	1863.1645	0.0000	1.0000
F81G	59	-873.0034	1864.0068	0.0000	1.0000
K2PI	57	-875.0416	1864.0832	0.0000	1.0000
K2PIG	58	-874.1654	1864.3307	0.0000	1.0000
F81I	59	-873.6459	1865.2919	0.0000	1.0000
F81IG	60	-872.9086	1865.8171	0.0000	1.0000
SYMIG	62	-871.4668	1866.9336	0.0000	1.0000
SYMI	61	-872.6301	1867.2603	0.0000	1.0000
K2P	56	-879.3030	1870.6059	0.0000	1.0000
SYM	60	-876.8121	1873.6243	0.0000	1.0000
F81	58	-879.2359	1874.4717	0.0000	1.0000
JC69G	56	-891.2136	1894.4273	0.0000	1.0000
JC69I	56	-893.8019	1899.6039	0.0000	1.0000
JC69IG	57	-892.9695	1899.9390	0.0000	1.0000
JC69	55	-896.9913	1903.9825	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-850.1396	2026.5641	0.6235	0.6235
HKYI	60	-850.8239	2027.9325	0.3146	0.9381
HKYIG	61	-850.0473	2031.8174	0.0451	0.9832
HKY	59	-856.4774	2033.8015	0.0167	0.9999
GTRG	64	-848.9780	2045.9930	0.0000	1.0000
GTRI	64	-849.7294	2047.4958	0.0000	1.0000
GTRIG	65	-848.9012	2051.2775	0.0000	1.0000
GTR	63	-855.1280	2052.8549	0.0000	1.0000
K2PG	57	-874.2650	2058.5005	0.0000	1.0000
K2PI	57	-875.0416	2060.0537	0.0000	1.0000
K2P	56	-879.3030	2063.1383	0.0000	1.0000
K2PIG	58	-874.1654	2063.7393	0.0000	1.0000
F81G	59	-873.0034	2066.8535	0.0000	1.0000
F81I	59	-873.6459	2068.1385	0.0000	1.0000
F81IG	60	-872.9086	2072.1019	0.0000	1.0000
SYMG	61	-870.5823	2072.8874	0.0000	1.0000
F81	58	-879.2359	2073.8803	0.0000	1.0000
SYMI	61	-872.6301	2076.9831	0.0000	1.0000
SYM	60	-876.8121	2079.9090	0.0000	1.0000
SYMIG	62	-871.4668	2080.0945	0.0000	1.0000
JC69G	56	-891.2136	2086.9597	0.0000	1.0000
JC69I	56	-893.8019	2092.1363	0.0000	1.0000
JC69	55	-896.9913	2093.0769	0.0000	1.0000
JC69IG	57	-892.9695	2095.9095	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/chr3_3177.phylip = 1 - 230;
 Partition Dummy = 1:183-phylip/chr3_3177.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/chr3_3177.phylip = 1 - 230;
 Partition Dummy = 1:183-phylip/chr3_3177.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_3177.phylip = 1 - 230;
 Partition Dummy = 1:183-phylip/chr3_3177.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

