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

Input data from file "261-aggressive-phylip/11261_hg18.phylip" (29 taxa, 296 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	-864.0780	1879.3049	0.4913	0.4913
HKYI	60	-864.6766	1880.5022	0.2700	0.7612
HKYIG	61	-863.5732	1881.4711	0.1663	0.9275
F81G	59	-868.7929	1885.5858	0.0213	0.9488
GTRG	64	-861.0136	1886.0446	0.0169	0.9657
F81I	59	-869.3862	1886.7725	0.0117	0.9774
GTRI	64	-861.6436	1887.3046	0.0090	0.9864
F81IG	60	-868.3084	1887.7657	0.0071	0.9936
GTRIG	65	-860.5689	1888.4422	0.0051	0.9987
K2PG	57	-875.6945	1893.1704	0.0005	0.9991
K2PI	57	-876.2922	1894.3659	0.0003	0.9994
SYMI	61	-870.2011	1894.7271	0.0002	0.9996
K2PIG	58	-875.2803	1895.4383	0.0002	0.9998
SYMIG	62	-869.1977	1895.9233	0.0001	0.9999
JC69G	56	-879.6306	1897.9725	0.0000	0.9999
JC69I	56	-880.2273	1899.1658	0.0000	1.0000
JC69IG	57	-879.2177	1900.2168	0.0000	1.0000
SYMG	61	-872.9994	1900.3237	0.0000	1.0000
HKY	59	-877.7871	1903.5743	0.0000	1.0000
GTR	63	-873.7289	1908.2164	0.0000	1.0000
F81	58	-882.3241	1909.5257	0.0000	1.0000
SYM	60	-881.4468	1914.0425	0.0000	1.0000
K2P	56	-888.4669	1915.6452	0.0000	1.0000
JC69	55	-892.3533	1920.3732	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-864.0780	1848.1559	0.3290	0.3290
HKYIG	61	-863.5732	1849.1463	0.2005	0.5294
HKYI	60	-864.6766	1849.3533	0.1808	0.7102
GTRG	64	-861.0136	1850.0273	0.1291	0.8393
GTRIG	65	-860.5689	1851.1378	0.0741	0.9134
GTRI	64	-861.6436	1851.2873	0.0687	0.9821
F81G	59	-868.7929	1855.5858	0.0080	0.9901
F81IG	60	-868.3084	1856.6167	0.0048	0.9949
F81I	59	-869.3862	1856.7725	0.0044	0.9993
SYMIG	62	-869.1977	1862.3954	0.0003	0.9996
SYMI	61	-870.2011	1862.4023	0.0003	0.9998
K2PG	57	-875.6945	1865.3889	0.0001	0.9999
K2PIG	58	-875.2803	1866.5606	0.0000	0.9999
K2PI	57	-876.2922	1866.5844	0.0000	1.0000
SYMG	61	-872.9994	1867.9989	0.0000	1.0000
JC69G	56	-879.6306	1871.2612	0.0000	1.0000
JC69IG	57	-879.2177	1872.4353	0.0000	1.0000
JC69I	56	-880.2273	1872.4545	0.0000	1.0000
GTR	63	-873.7289	1873.4578	0.0000	1.0000
HKY	59	-877.7871	1873.5743	0.0000	1.0000
F81	58	-882.3241	1880.6481	0.0000	1.0000
SYM	60	-881.4468	1882.8935	0.0000	1.0000
K2P	56	-888.4669	1888.9339	0.0000	1.0000
JC69	55	-892.3533	1894.7066	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-864.0780	2069.5775	0.4998	0.4998
HKYI	60	-864.6766	2070.7748	0.2747	0.7745
F81G	59	-868.7929	2073.3170	0.0771	0.8515
HKYIG	61	-863.5732	2074.2583	0.0481	0.8997
F81I	59	-869.3862	2074.5037	0.0426	0.9422
K2PG	57	-875.6945	2075.7394	0.0229	0.9652
K2PI	57	-876.2922	2076.9348	0.0126	0.9778
JC69G	56	-879.6306	2077.9214	0.0077	0.9855
F81IG	60	-868.3084	2078.0383	0.0073	0.9928
JC69I	56	-880.2273	2079.1146	0.0042	0.9970
K2PIG	58	-875.2803	2080.6015	0.0020	0.9990
JC69IG	57	-879.2177	2082.7858	0.0007	0.9997
GTRG	64	-861.0136	2086.2103	0.0001	0.9998
GTRI	64	-861.6436	2087.4703	0.0001	0.9999
SYMI	61	-870.2011	2087.5142	0.0001	1.0000
GTRIG	65	-860.5689	2091.0112	0.0000	1.0000
SYMIG	62	-869.1977	2091.1977	0.0000	1.0000
HKY	59	-877.7871	2091.3055	0.0000	1.0000
SYMG	61	-872.9994	2093.1108	0.0000	1.0000
F81	58	-882.3241	2094.6889	0.0000	1.0000
K2P	56	-888.4669	2095.5940	0.0000	1.0000
JC69	55	-892.3533	2097.6763	0.0000	1.0000
SYM	60	-881.4468	2104.3151	0.0000	1.0000
GTR	63	-873.7289	2105.9504	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 261-aggressive-phylip/11261_hg18.phylip = 1 - 296;
 Partition Dummy = 1:261-aggressive-phylip/11261_hg18.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 261-aggressive-phylip/11261_hg18.phylip = 1 - 296;
 Partition Dummy = 1:261-aggressive-phylip/11261_hg18.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 261-aggressive-phylip/11261_hg18.phylip = 1 - 296;
 Partition Dummy = 1:261-aggressive-phylip/11261_hg18.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

