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

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

Minimum AIC  model: GTRI
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	-866.0037	1882.1385	0.3911	0.3911
K2PI	57	-866.2337	1882.5986	0.3107	0.7018
K2PIG	58	-866.0148	1885.6340	0.0681	0.7699
HKYI	60	-862.4856	1885.6378	0.0680	0.8379
HKYG	60	-862.5589	1885.7845	0.0632	0.9010
K2P	56	-869.7893	1886.2743	0.0494	0.9505
HKYIG	61	-862.4610	1889.1789	0.0116	0.9620
HKY	59	-866.2251	1889.5662	0.0095	0.9716
SYMG	61	-862.8743	1890.0056	0.0077	0.9792
SYMI	61	-863.0318	1890.3206	0.0065	0.9858
GTRI	64	-857.5844	1890.4415	0.0062	0.9919
GTRG	64	-857.8294	1890.9316	0.0048	0.9967
SYMIG	62	-862.8864	1893.6604	0.0012	0.9980
GTRIG	65	-857.5861	1894.2007	0.0009	0.9989
SYM	60	-866.9902	1894.6471	0.0008	0.9997
GTR	63	-862.3748	1896.3088	0.0003	1.0000
JC69G	56	-884.1909	1915.0775	0.0000	1.0000
JC69I	56	-884.6313	1915.9582	0.0000	1.0000
JC69IG	57	-884.2030	1918.5372	0.0000	1.0000
JC69	55	-887.8615	1919.0203	0.0000	1.0000
F81G	59	-882.0253	1921.1666	0.0000	1.0000
F81I	59	-882.4653	1922.0467	0.0000	1.0000
F81	58	-885.4718	1924.5481	0.0000	1.0000
F81IG	60	-882.0362	1924.7391	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	64	-857.5844	1843.1687	0.2563	0.2563
GTRG	64	-857.8294	1843.6589	0.2006	0.4569
HKYI	60	-862.4856	1844.9711	0.1041	0.5610
HKYG	60	-862.5589	1845.1179	0.0967	0.6577
GTRIG	65	-857.5861	1845.1721	0.0941	0.7519
K2PG	57	-866.0037	1846.0073	0.0620	0.8139
K2PI	57	-866.2337	1846.4674	0.0493	0.8631
HKYIG	61	-862.4610	1846.9219	0.0392	0.9024
SYMG	61	-862.8743	1847.7486	0.0260	0.9283
K2PIG	58	-866.0148	1848.0296	0.0226	0.9509
SYMI	61	-863.0318	1848.0637	0.0222	0.9731
SYMIG	62	-862.8864	1849.7728	0.0094	0.9825
HKY	59	-866.2251	1850.4502	0.0067	0.9892
GTR	63	-862.3748	1850.7495	0.0058	0.9950
K2P	56	-869.7893	1851.5787	0.0038	0.9988
SYM	60	-866.9902	1853.9804	0.0012	1.0000
JC69G	56	-884.1909	1880.3818	0.0000	1.0000
JC69I	56	-884.6313	1881.2625	0.0000	1.0000
F81G	59	-882.0253	1882.0505	0.0000	1.0000
JC69IG	57	-884.2030	1882.4060	0.0000	1.0000
F81I	59	-882.4653	1882.9307	0.0000	1.0000
F81IG	60	-882.0362	1884.0724	0.0000	1.0000
JC69	55	-887.8615	1885.7230	0.0000	1.0000
F81	58	-885.4718	1886.9437	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-866.0037	2044.6408	0.4479	0.4479
K2PI	57	-866.2337	2045.1008	0.3558	0.8037
K2P	56	-869.7893	2046.7273	0.1578	0.9615
K2PIG	58	-866.0148	2050.1478	0.0285	0.9901
HKYI	60	-862.4856	2054.0590	0.0040	0.9941
HKYG	60	-862.5589	2054.2057	0.0038	0.9978
HKY	59	-866.2251	2056.0532	0.0015	0.9993
HKYIG	61	-862.4610	2059.4945	0.0003	0.9996
SYMG	61	-862.8743	2060.3212	0.0002	0.9998
SYMI	61	-863.0318	2060.6363	0.0002	0.9999
SYM	60	-866.9902	2063.0682	0.0000	1.0000
SYMIG	62	-862.8864	2065.8302	0.0000	1.0000
GTRI	64	-857.5844	2066.1957	0.0000	1.0000
GTRG	64	-857.8294	2066.6859	0.0000	1.0000
GTR	63	-862.3748	2070.2917	0.0000	1.0000
GTRIG	65	-857.5861	2071.6839	0.0000	1.0000
JC69G	56	-884.1909	2075.5305	0.0000	1.0000
JC69I	56	-884.6313	2076.4112	0.0000	1.0000
JC69	55	-887.8615	2077.3869	0.0000	1.0000
JC69IG	57	-884.2030	2081.0394	0.0000	1.0000
F81G	59	-882.0253	2087.6536	0.0000	1.0000
F81I	59	-882.4653	2088.5337	0.0000	1.0000
F81	58	-885.4718	2089.0619	0.0000	1.0000
F81IG	60	-882.0362	2093.1602	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/7244_hg18.phylip = 1 - 241;
 Partition Dummy = 1:261-aggressive-phylip/7244_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 (K2PG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/7244_hg18.phylip = 1 - 241;
 Partition Dummy = 1:261-aggressive-phylip/7244_hg18.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 261-aggressive-phylip/7244_hg18.phylip = 1 - 241;
 Partition Dummy = 1:261-aggressive-phylip/7244_hg18.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

