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

Input data from file "183-phylip/chr7_4670.phylip" (29 taxa, 606 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	-1844.4932	3822.4176	0.5596	0.5596
HKYI	60	-1845.4883	3824.4077	0.2069	0.7665
HKYIG	61	-1844.4927	3824.8898	0.1626	0.9291
HKY	59	-1848.7448	3828.4567	0.0273	0.9564
GTRG	64	-1842.6145	3828.6079	0.0253	0.9817
GTRI	64	-1843.6760	3830.7309	0.0088	0.9905
GTRIG	65	-1842.6141	3831.1170	0.0072	0.9977
GTR	63	-1846.6572	3834.1926	0.0016	0.9993
K2PG	57	-1855.3611	3836.7878	0.0004	0.9997
K2PI	57	-1856.4525	3838.9706	0.0001	0.9998
K2PIG	58	-1855.3605	3839.2329	0.0001	0.9999
SYMG	61	-1853.4344	3842.7733	0.0000	1.0000
K2P	56	-1859.6428	3842.9139	0.0000	1.0000
SYMI	61	-1854.6361	3845.1767	0.0000	1.0000
SYMIG	62	-1853.4340	3845.2546	0.0000	1.0000
SYM	60	-1857.5433	3848.5179	0.0000	1.0000
F81G	59	-1885.2416	3901.4502	0.0000	1.0000
F81I	59	-1886.2490	3903.4650	0.0000	1.0000
F81IG	60	-1885.2412	3903.9136	0.0000	1.0000
F81	58	-1889.1355	3906.7829	0.0000	1.0000
JC69G	56	-1894.5599	3912.7481	0.0000	1.0000
JC69I	56	-1895.6161	3914.8606	0.0000	1.0000
JC69IG	57	-1894.5593	3915.1842	0.0000	1.0000
JC69	55	-1898.4834	3918.1668	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1844.4932	3808.9864	0.5029	0.5029
HKYI	60	-1845.4883	3810.9765	0.1859	0.6888
HKYIG	61	-1844.4927	3810.9854	0.1851	0.8739
GTRG	64	-1842.6145	3813.2289	0.0603	0.9342
GTRIG	65	-1842.6141	3815.2282	0.0222	0.9564
GTRI	64	-1843.6760	3815.3520	0.0209	0.9773
HKY	59	-1848.7448	3815.4896	0.0195	0.9967
GTR	63	-1846.6572	3819.3144	0.0029	0.9996
K2PG	57	-1855.3611	3824.7221	0.0002	0.9998
K2PIG	58	-1855.3605	3826.7210	0.0001	0.9999
K2PI	57	-1856.4525	3826.9049	0.0001	1.0000
SYMG	61	-1853.4344	3828.8689	0.0000	1.0000
SYMIG	62	-1853.4340	3830.8679	0.0000	1.0000
SYMI	61	-1854.6361	3831.2723	0.0000	1.0000
K2P	56	-1859.6428	3831.2855	0.0000	1.0000
SYM	60	-1857.5433	3835.0867	0.0000	1.0000
F81G	59	-1885.2416	3888.4832	0.0000	1.0000
F81IG	60	-1885.2412	3890.4824	0.0000	1.0000
F81I	59	-1886.2490	3890.4980	0.0000	1.0000
F81	58	-1889.1355	3894.2710	0.0000	1.0000
JC69G	56	-1894.5599	3901.1197	0.0000	1.0000
JC69IG	57	-1894.5593	3903.1185	0.0000	1.0000
JC69I	56	-1895.6161	3903.2322	0.0000	1.0000
JC69	55	-1898.4834	3906.9668	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1844.4932	4073.3992	0.4446	0.4446
HKYI	60	-1845.4883	4075.3893	0.1644	0.6090
HKY	59	-1848.7448	4075.4956	0.1559	0.7649
K2PG	57	-1855.3611	4075.9143	0.1264	0.8913
K2P	56	-1859.6428	4078.0708	0.0430	0.9343
K2PI	57	-1856.4525	4078.0971	0.0424	0.9768
HKYIG	61	-1844.4927	4079.8051	0.0181	0.9948
K2PIG	58	-1855.3605	4082.3201	0.0051	1.0000
GTRG	64	-1842.6145	4095.2693	0.0000	1.0000
GTR	63	-1846.6572	4096.9479	0.0000	1.0000
GTRI	64	-1843.6760	4097.3923	0.0000	1.0000
SYMG	61	-1853.4344	4097.6885	0.0000	1.0000
SYM	60	-1857.5433	4099.4995	0.0000	1.0000
SYMI	61	-1854.6361	4100.0920	0.0000	1.0000
GTRIG	65	-1842.6141	4101.6754	0.0000	1.0000
SYMIG	62	-1853.4340	4104.0945	0.0000	1.0000
JC69G	56	-1894.5599	4147.9050	0.0000	1.0000
F81G	59	-1885.2416	4148.4891	0.0000	1.0000
JC69	55	-1898.4834	4149.3452	0.0000	1.0000
F81	58	-1889.1355	4149.8700	0.0000	1.0000
JC69I	56	-1895.6161	4150.0175	0.0000	1.0000
F81I	59	-1886.2490	4150.5039	0.0000	1.0000
JC69IG	57	-1894.5593	4154.3107	0.0000	1.0000
F81IG	60	-1885.2412	4154.8952	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 183-phylip/chr7_4670.phylip = 1 - 606;
 Partition Dummy = 1:183-phylip/chr7_4670.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 183-phylip/chr7_4670.phylip = 1 - 606;
 Partition Dummy = 1:183-phylip/chr7_4670.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 183-phylip/chr7_4670.phylip = 1 - 606;
 Partition Dummy = 1:183-phylip/chr7_4670.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

