
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
21:30 04/04/2011

Input data from file "591-tomodel-phylip/chr20_2151.phylip" (7 taxa, 369 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	16	-864.8372	1763.2199	0.3739	0.3739
HKYI	16	-865.0122	1763.5698	0.3139	0.6877
HKYIG	17	-864.7339	1765.2115	0.1381	0.8258
GTRG	20	-861.9771	1766.3680	0.0775	0.9033
GTRI	20	-862.1889	1766.7916	0.0627	0.9660
GTRIG	21	-861.7786	1768.2199	0.0307	0.9967
HKY	15	-872.0914	1775.5425	0.0008	0.9975
K2PG	13	-874.6702	1776.3658	0.0005	0.9980
K2PI	13	-874.8759	1776.7772	0.0004	0.9984
SYMG	17	-870.7367	1777.2170	0.0003	0.9987
SYMI	17	-871.0025	1777.7486	0.0003	0.9990
F81G	15	-873.2316	1777.8229	0.0003	0.9993
F81I	15	-873.3622	1778.0843	0.0002	0.9995
K2PIG	14	-874.5753	1778.3371	0.0002	0.9997
SYMIG	18	-870.5414	1779.0371	0.0001	0.9998
GTR	19	-869.8455	1779.8686	0.0001	0.9999
F81IG	16	-873.1809	1779.9072	0.0001	1.0000
K2P	12	-881.3151	1787.5065	0.0000	1.0000
F81	14	-879.8943	1788.9751	0.0000	1.0000
JC69G	12	-882.1535	1789.1834	0.0000	1.0000
SYM	16	-877.9587	1789.4629	0.0000	1.0000
JC69I	12	-882.3163	1789.5091	0.0000	1.0000
JC69IG	13	-882.0952	1791.2157	0.0000	1.0000
JC69	11	-888.5088	1799.7570	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-864.8372	1761.6744	0.3358	0.3358
HKYI	16	-865.0122	1762.0243	0.2819	0.6177
HKYIG	17	-864.7339	1763.4679	0.1370	0.7547
GTRG	20	-861.9771	1763.9542	0.1074	0.8621
GTRI	20	-862.1889	1764.3778	0.0869	0.9490
GTRIG	21	-861.7786	1765.5571	0.0482	0.9972
HKY	15	-872.0914	1774.1827	0.0006	0.9979
K2PG	13	-874.6702	1775.3404	0.0004	0.9982
SYMG	17	-870.7367	1775.4734	0.0003	0.9986
K2PI	13	-874.8759	1775.7519	0.0003	0.9989
SYMI	17	-871.0025	1776.0051	0.0003	0.9991
F81G	15	-873.2316	1776.4632	0.0002	0.9993
F81I	15	-873.3622	1776.7245	0.0002	0.9995
SYMIG	18	-870.5414	1777.0829	0.0002	0.9997
K2PIG	14	-874.5753	1777.1507	0.0001	0.9998
GTR	19	-869.8455	1777.6910	0.0001	0.9999
F81IG	16	-873.1809	1778.3618	0.0001	1.0000
K2P	12	-881.3151	1786.6301	0.0000	1.0000
F81	14	-879.8943	1787.7886	0.0000	1.0000
SYM	16	-877.9587	1787.9175	0.0000	1.0000
JC69G	12	-882.1535	1788.3070	0.0000	1.0000
JC69I	12	-882.3163	1788.6327	0.0000	1.0000
JC69IG	13	-882.0952	1790.1904	0.0000	1.0000
JC69	11	-888.5088	1799.0175	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-864.8372	1824.2472	0.3775	0.3775
HKYI	16	-865.0122	1824.5971	0.3169	0.6944
K2PG	13	-874.6702	1826.1808	0.1436	0.8380
K2PI	13	-874.8759	1826.5922	0.1169	0.9549
HKYIG	17	-864.7339	1829.9514	0.0218	0.9767
K2PIG	14	-874.5753	1831.9018	0.0082	0.9849
HKY	15	-872.0914	1832.8447	0.0051	0.9900
K2P	12	-881.3151	1833.5597	0.0036	0.9936
F81G	15	-873.2316	1835.1251	0.0016	0.9952
JC69G	12	-882.1535	1835.2365	0.0016	0.9968
F81I	15	-873.3622	1835.3864	0.0014	0.9982
JC69I	12	-882.3163	1835.5622	0.0013	0.9995
F81IG	16	-873.1809	1840.9345	0.0001	0.9996
JC69IG	13	-882.0952	1841.0307	0.0001	0.9997
SYMG	17	-870.7367	1841.9570	0.0001	0.9998
JC69	11	-888.5088	1842.0363	0.0001	0.9998
GTRG	20	-861.9771	1842.1702	0.0000	0.9999
SYMI	17	-871.0025	1842.4886	0.0000	0.9999
F81	14	-879.8943	1842.5398	0.0000	1.0000
GTRI	20	-862.1889	1842.5937	0.0000	1.0000
SYMIG	18	-870.5414	1847.4772	0.0000	1.0000
GTRIG	21	-861.7786	1847.6838	0.0000	1.0000
SYM	16	-877.9587	1850.4902	0.0000	1.0000
GTR	19	-869.8455	1851.9961	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 591-tomodel-phylip/chr20_2151.phylip = 1 - 369;
 Partition Dummy = 1:591-tomodel-phylip/chr20_2151.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 591-tomodel-phylip/chr20_2151.phylip = 1 - 369;
 Partition Dummy = 1:591-tomodel-phylip/chr20_2151.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 591-tomodel-phylip/chr20_2151.phylip = 1 - 369;
 Partition Dummy = 1:591-tomodel-phylip/chr20_2151.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

