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

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

Minimum AIC  model: GTRG
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	-854.2150	1865.6456	0.6577	0.6577
K2PI	57	-855.4679	1868.1515	0.1879	0.8456
K2PIG	58	-854.2081	1869.4425	0.0985	0.9441
HKYG	60	-851.4617	1871.7235	0.0315	0.9756
HKYI	60	-852.6581	1874.1162	0.0095	0.9851
SYMG	61	-851.1823	1875.1297	0.0057	0.9909
HKYIG	61	-851.4557	1875.6765	0.0044	0.9952
SYMI	61	-852.3023	1877.3696	0.0019	0.9971
GTRG	64	-846.5405	1878.0673	0.0013	0.9984
SYMIG	62	-851.1811	1879.1460	0.0008	0.9992
GTRI	64	-847.7076	1880.4016	0.0004	0.9996
K2P	56	-864.2396	1881.9338	0.0002	0.9998
GTRIG	65	-846.5409	1882.2542	0.0002	1.0000
HKY	59	-860.6420	1886.1714	0.0000	1.0000
SYM	60	-860.5397	1889.8793	0.0000	1.0000
GTR	63	-855.6452	1892.1476	0.0000	1.0000
JC69G	56	-869.6320	1892.7185	0.0000	1.0000
JC69I	56	-870.7660	1894.9866	0.0000	1.0000
JC69IG	57	-869.6263	1896.4682	0.0000	1.0000
F81G	59	-867.8116	1900.5106	0.0000	1.0000
F81I	59	-868.9198	1902.7269	0.0000	1.0000
F81IG	60	-867.8081	1904.4161	0.0000	1.0000
JC69	55	-879.2759	1908.2938	0.0000	1.0000
F81	58	-877.2147	1915.4558	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-846.5405	1821.0810	0.2843	0.2843
K2PG	57	-854.2150	1822.4299	0.1448	0.4292
HKYG	60	-851.4617	1822.9235	0.1132	0.5424
GTRIG	65	-846.5409	1823.0818	0.1046	0.6469
GTRI	64	-847.7076	1823.4153	0.0885	0.7354
SYMG	61	-851.1823	1824.3646	0.0551	0.7905
K2PIG	58	-854.2081	1824.4161	0.0537	0.8441
HKYIG	61	-851.4557	1824.9114	0.0419	0.8860
K2PI	57	-855.4679	1824.9358	0.0414	0.9274
HKYI	60	-852.6581	1825.3162	0.0342	0.9616
SYMIG	62	-851.1811	1826.3622	0.0203	0.9819
SYMI	61	-852.3023	1826.6045	0.0180	0.9999
GTR	63	-855.6452	1837.2905	0.0001	0.9999
HKY	59	-860.6420	1839.2840	0.0000	1.0000
K2P	56	-864.2396	1840.4793	0.0000	1.0000
SYM	60	-860.5397	1841.0793	0.0000	1.0000
JC69G	56	-869.6320	1851.2639	0.0000	1.0000
JC69IG	57	-869.6263	1853.2525	0.0000	1.0000
JC69I	56	-870.7660	1853.5320	0.0000	1.0000
F81G	59	-867.8116	1853.6232	0.0000	1.0000
F81IG	60	-867.8081	1855.6161	0.0000	1.0000
F81I	59	-868.9198	1855.8395	0.0000	1.0000
JC69	55	-879.2759	1868.5518	0.0000	1.0000
F81	58	-877.2147	1870.4294	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-854.2150	2013.4858	0.7335	0.7335
K2PI	57	-855.4679	2015.9917	0.2095	0.9430
K2PIG	58	-854.2081	2018.8239	0.0508	0.9939
HKYG	60	-851.4617	2024.0349	0.0038	0.9976
HKYI	60	-852.6581	2026.4277	0.0011	0.9988
K2P	56	-864.2396	2028.1833	0.0005	0.9992
SYMG	61	-851.1823	2028.8280	0.0003	0.9996
HKYIG	61	-851.4557	2029.3747	0.0003	0.9998
SYMI	61	-852.3023	2031.0679	0.0001	1.0000
SYMIG	62	-851.1811	2034.1774	0.0000	1.0000
GTRG	64	-846.5405	2035.5999	0.0000	1.0000
HKY	59	-860.6420	2037.0436	0.0000	1.0000
GTRI	64	-847.7076	2037.9342	0.0000	1.0000
JC69G	56	-869.6320	2038.9680	0.0000	1.0000
GTRIG	65	-846.5409	2040.9526	0.0000	1.0000
JC69I	56	-870.7660	2041.2361	0.0000	1.0000
SYM	60	-860.5397	2042.1908	0.0000	1.0000
JC69IG	57	-869.6263	2044.3084	0.0000	1.0000
GTR	63	-855.6452	2048.4575	0.0000	1.0000
F81G	59	-867.8116	2051.3828	0.0000	1.0000
JC69	55	-879.2759	2052.9040	0.0000	1.0000
F81I	59	-868.9198	2053.5991	0.0000	1.0000
F81IG	60	-867.8081	2056.7276	0.0000	1.0000
F81	58	-877.2147	2064.8372	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr8_1574.phylip = 1 - 211;
 Partition Dummy = 1:183-phylip/chr8_1574.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr8_1574.phylip = 1 - 211;
 Partition Dummy = 1:183-phylip/chr8_1574.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 183-phylip/chr8_1574.phylip = 1 - 211;
 Partition Dummy = 1:183-phylip/chr8_1574.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

