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

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

Minimum AIC  model: F81I
Minimum AICc model: F81I
Minimum BIC  model: F81I


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
F81I	59	-839.4775	1814.6108	0.5427	0.5427
F81IG	60	-839.2003	1816.7006	0.1909	0.7336
HKYI	60	-839.4236	1817.1471	0.1527	0.8862
HKYIG	61	-839.1274	1819.2122	0.0544	0.9406
F81G	59	-841.9899	1819.6356	0.0440	0.9846
HKYG	60	-841.9065	1822.1130	0.0127	0.9974
GTRI	64	-838.5396	1826.0893	0.0017	0.9991
GTRIG	65	-838.1530	1828.0275	0.0007	0.9998
GTRG	64	-840.6748	1830.3596	0.0002	1.0000
F81	58	-851.4024	1835.8297	0.0000	1.0000
HKY	59	-851.3014	1838.2586	0.0000	1.0000
JC69I	56	-859.0704	1845.9428	0.0000	1.0000
GTR	63	-849.9074	1846.1271	0.0000	1.0000
JC69IG	57	-858.8364	1848.0797	0.0000	1.0000
K2PI	57	-859.0477	1848.5024	0.0000	1.0000
K2PIG	58	-858.8030	1850.6308	0.0000	1.0000
JC69G	56	-861.7626	1851.3272	0.0000	1.0000
K2PG	57	-861.7361	1853.8791	0.0000	1.0000
SYMI	61	-858.0527	1857.0627	0.0000	1.0000
SYMIG	62	-857.7329	1859.0940	0.0000	1.0000
SYMG	61	-860.5366	1862.0305	0.0000	1.0000
JC69	55	-871.2604	1867.7306	0.0000	1.0000
K2P	56	-871.2318	1870.2655	0.0000	1.0000
SYM	60	-869.9236	1878.1472	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
F81I	59	-839.4775	1796.9550	0.4529	0.4529
F81IG	60	-839.2003	1798.4006	0.2199	0.6728
HKYI	60	-839.4236	1798.8471	0.1759	0.8486
HKYIG	61	-839.1274	1800.2548	0.0870	0.9356
F81G	59	-841.9899	1801.9798	0.0367	0.9724
HKYG	60	-841.9065	1803.8130	0.0147	0.9870
GTRI	64	-838.5396	1805.0792	0.0078	0.9948
GTRIG	65	-838.1530	1806.3060	0.0042	0.9991
GTRG	64	-840.6748	1809.3495	0.0009	1.0000
F81	58	-851.4024	1818.8048	0.0000	1.0000
HKY	59	-851.3014	1820.6027	0.0000	1.0000
GTR	63	-849.9074	1825.8147	0.0000	1.0000
JC69I	56	-859.0704	1830.1408	0.0000	1.0000
JC69IG	57	-858.8364	1831.6727	0.0000	1.0000
K2PI	57	-859.0477	1832.0954	0.0000	1.0000
K2PIG	58	-858.8030	1833.6059	0.0000	1.0000
JC69G	56	-861.7626	1835.5252	0.0000	1.0000
K2PG	57	-861.7361	1837.4721	0.0000	1.0000
SYMI	61	-858.0527	1838.1053	0.0000	1.0000
SYMIG	62	-857.7329	1839.4658	0.0000	1.0000
SYMG	61	-860.5366	1843.0731	0.0000	1.0000
JC69	55	-871.2604	1852.5208	0.0000	1.0000
K2P	56	-871.2318	1854.4635	0.0000	1.0000
SYM	60	-869.9236	1859.8472	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
F81I	59	-839.4775	2040.8255	0.8340	0.8340
F81G	59	-841.9899	2045.8503	0.0676	0.9016
F81IG	60	-839.2003	2046.4044	0.0513	0.9529
HKYI	60	-839.4236	2046.8510	0.0410	0.9939
HKYG	60	-841.9065	2051.8169	0.0034	0.9973
HKYIG	61	-839.1274	2052.3921	0.0026	0.9998
F81	58	-851.4024	2058.5419	0.0001	1.0000
JC69I	56	-859.0704	2061.6111	0.0000	1.0000
HKY	59	-851.3014	2064.4732	0.0000	1.0000
JC69G	56	-861.7626	2066.9955	0.0000	1.0000
JC69IG	57	-858.8364	2067.2764	0.0000	1.0000
K2PI	57	-859.0477	2067.6991	0.0000	1.0000
GTRI	64	-838.5396	2069.6167	0.0000	1.0000
K2PG	57	-861.7361	2073.0758	0.0000	1.0000
K2PIG	58	-858.8030	2073.3430	0.0000	1.0000
GTRG	64	-840.6748	2073.8870	0.0000	1.0000
GTRIG	65	-838.1530	2074.9768	0.0000	1.0000
JC69	55	-871.2604	2079.8577	0.0000	1.0000
K2P	56	-871.2318	2085.9338	0.0000	1.0000
GTR	63	-849.9074	2086.2188	0.0000	1.0000
SYMI	61	-858.0527	2090.2426	0.0000	1.0000
SYMG	61	-860.5366	2095.2104	0.0000	1.0000
SYMIG	62	-857.7329	2095.7365	0.0000	1.0000
SYM	60	-869.9236	2107.8510	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/7680_hg18.phylip = 1 - 461;
 Partition Dummy = 1:261-aggressive-phylip/7680_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=1 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/7680_hg18.phylip = 1 - 461;
 Partition Dummy = 1:261-aggressive-phylip/7680_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=1 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/7680_hg18.phylip = 1 - 461;
 Partition Dummy = 1:261-aggressive-phylip/7680_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=1 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

