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

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

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


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	-832.1290	1822.2941	0.2677	0.2677
GTRG	64	-832.1805	1822.3971	0.2542	0.5219
HKYI	60	-838.3783	1822.8063	0.2072	0.7291
HKYG	60	-838.6237	1823.2973	0.1621	0.8912
GTRIG	65	-832.0839	1825.2547	0.0609	0.9521
HKYIG	61	-838.3792	1825.7727	0.0470	0.9991
K2PI	57	-849.1316	1835.5448	0.0004	0.9995
K2PG	57	-849.3730	1836.0276	0.0003	0.9998
K2PIG	58	-849.1325	1838.4487	0.0001	0.9998
SYMI	61	-844.8217	1838.6578	0.0001	0.9999
SYMG	61	-844.9180	1838.8503	0.0001	1.0000
SYMIG	62	-844.8190	1841.6380	0.0000	1.0000
F81I	59	-852.6901	1848.4865	0.0000	1.0000
F81G	59	-852.9727	1849.0519	0.0000	1.0000
F81IG	60	-852.6908	1851.4315	0.0000	1.0000
GTR	63	-849.1631	1853.3335	0.0000	1.0000
HKY	59	-856.4136	1855.9336	0.0000	1.0000
JC69I	56	-862.4039	1859.2078	0.0000	1.0000
JC69G	56	-862.7319	1859.8638	0.0000	1.0000
JC69IG	57	-862.4048	1862.0912	0.0000	1.0000
K2P	56	-866.9908	1868.3817	0.0000	1.0000
SYM	60	-861.6182	1869.2861	0.0000	1.0000
F81	58	-871.5722	1883.3282	0.0000	1.0000
JC69	55	-880.8658	1893.2700	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	64	-832.1290	1792.2580	0.3904	0.3904
GTRG	64	-832.1805	1792.3610	0.3708	0.7611
GTRIG	65	-832.0839	1794.1678	0.1502	0.9114
HKYI	60	-838.3783	1796.7565	0.0412	0.9526
HKYG	60	-838.6237	1797.2475	0.0322	0.9848
HKYIG	61	-838.3792	1798.7584	0.0151	0.9999
SYMI	61	-844.8217	1811.6435	0.0000	0.9999
SYMG	61	-844.9180	1811.8360	0.0000	1.0000
K2PI	57	-849.1316	1812.2631	0.0000	1.0000
K2PG	57	-849.3730	1812.7459	0.0000	1.0000
SYMIG	62	-844.8190	1813.6380	0.0000	1.0000
K2PIG	58	-849.1325	1814.2650	0.0000	1.0000
F81I	59	-852.6901	1823.3802	0.0000	1.0000
F81G	59	-852.9727	1823.9455	0.0000	1.0000
GTR	63	-849.1631	1824.3263	0.0000	1.0000
F81IG	60	-852.6908	1825.3817	0.0000	1.0000
HKY	59	-856.4136	1830.8272	0.0000	1.0000
JC69I	56	-862.4039	1836.8078	0.0000	1.0000
JC69G	56	-862.7319	1837.4638	0.0000	1.0000
JC69IG	57	-862.4048	1838.8095	0.0000	1.0000
SYM	60	-861.6182	1843.2363	0.0000	1.0000
K2P	56	-866.9908	1845.9817	0.0000	1.0000
F81	58	-871.5722	1859.1445	0.0000	1.0000
JC69	55	-880.8658	1871.7315	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	60	-838.3783	2026.8452	0.4775	0.4775
HKYG	60	-838.6237	2027.3361	0.3736	0.8511
K2PI	57	-849.1316	2030.8473	0.0646	0.9157
K2PG	57	-849.3730	2031.3301	0.0507	0.9664
HKYIG	61	-838.3792	2032.6819	0.0258	0.9922
K2PIG	58	-849.1325	2036.6840	0.0035	0.9957
GTRI	64	-832.1290	2037.6859	0.0021	0.9978
GTRG	64	-832.1805	2037.7889	0.0020	0.9998
GTRIG	65	-832.0839	2043.4305	0.0001	0.9999
SYMI	61	-844.8217	2045.5670	0.0000	0.9999
SYMG	61	-844.9180	2045.7595	0.0000	1.0000
F81I	59	-852.6901	2049.6340	0.0000	1.0000
F81G	59	-852.9727	2050.1993	0.0000	1.0000
SYMIG	62	-844.8190	2051.3962	0.0000	1.0000
JC69I	56	-862.4039	2051.5572	0.0000	1.0000
JC69G	56	-862.7319	2052.2132	0.0000	1.0000
F81IG	60	-852.6908	2055.4703	0.0000	1.0000
HKY	59	-856.4136	2057.0811	0.0000	1.0000
JC69IG	57	-862.4048	2057.3937	0.0000	1.0000
K2P	56	-866.9908	2060.7311	0.0000	1.0000
GTR	63	-849.1631	2065.9193	0.0000	1.0000
SYM	60	-861.6182	2073.3249	0.0000	1.0000
F81	58	-871.5722	2081.5635	0.0000	1.0000
JC69	55	-880.8658	2082.6461	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 183-phylip/chr4_7577.phylip = 1 - 342;
 Partition Dummy = 1:183-phylip/chr4_7577.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 183-phylip/chr4_7577.phylip = 1 - 342;
 Partition Dummy = 1:183-phylip/chr4_7577.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 (HKYI)]

BEGIN MRBAYES;
 Charset 183-phylip/chr4_7577.phylip = 1 - 342;
 Partition Dummy = 1:183-phylip/chr4_7577.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

