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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	64	-1816.6717	3778.9334	0.4113	0.4113
GTRIG	65	-1815.8013	3779.7806	0.2693	0.6806
HKYG	60	-1822.6889	3780.7238	0.1680	0.8486
HKYIG	61	-1821.6751	3781.2410	0.1297	0.9784
GTRI	64	-1820.1996	3785.9891	0.0121	0.9904
HKYI	60	-1825.5564	3786.4587	0.0096	1.0000
F81G	59	-1847.7102	3828.2321	0.0000	1.0000
F81IG	60	-1846.5967	3828.5393	0.0000	1.0000
F81I	59	-1850.2975	3833.4067	0.0000	1.0000
SYMG	61	-1861.2241	3860.3390	0.0000	1.0000
SYMIG	62	-1860.2171	3860.8806	0.0000	1.0000
SYMI	61	-1865.4634	3868.8176	0.0000	1.0000
HKY	59	-1871.7665	3876.3447	0.0000	1.0000
K2PG	57	-1876.0231	3879.8212	0.0000	1.0000
GTR	63	-1868.4981	3880.0089	0.0000	1.0000
K2PIG	58	-1875.2846	3880.8574	0.0000	1.0000
K2PI	57	-1879.0826	3885.9402	0.0000	1.0000
F81	58	-1893.5288	3917.3457	0.0000	1.0000
JC69G	56	-1905.7478	3936.7680	0.0000	1.0000
JC69IG	57	-1905.0962	3937.9674	0.0000	1.0000
JC69I	56	-1908.0106	3941.2936	0.0000	1.0000
SYM	60	-1908.5076	3952.3612	0.0000	1.0000
K2P	56	-1920.6487	3966.5697	0.0000	1.0000
JC69	55	-1947.7092	4018.1985	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1816.6717	3761.3435	0.4580	0.4580
GTRIG	65	-1815.8013	3761.6027	0.4024	0.8604
HKYIG	61	-1821.6751	3765.3502	0.0618	0.9222
HKYG	60	-1822.6889	3765.3779	0.0609	0.9831
GTRI	64	-1820.1996	3768.3993	0.0135	0.9965
HKYI	60	-1825.5564	3771.1127	0.0035	1.0000
F81IG	60	-1846.5967	3813.1933	0.0000	1.0000
F81G	59	-1847.7102	3813.4203	0.0000	1.0000
F81I	59	-1850.2975	3818.5950	0.0000	1.0000
SYMIG	62	-1860.2171	3844.4343	0.0000	1.0000
SYMG	61	-1861.2241	3844.4483	0.0000	1.0000
SYMI	61	-1865.4634	3852.9269	0.0000	1.0000
HKY	59	-1871.7665	3861.5330	0.0000	1.0000
GTR	63	-1868.4981	3862.9962	0.0000	1.0000
K2PG	57	-1876.0231	3866.0462	0.0000	1.0000
K2PIG	58	-1875.2846	3866.5693	0.0000	1.0000
K2PI	57	-1879.0826	3872.1652	0.0000	1.0000
F81	58	-1893.5288	3903.0576	0.0000	1.0000
JC69G	56	-1905.7478	3923.4956	0.0000	1.0000
JC69IG	57	-1905.0962	3924.1924	0.0000	1.0000
JC69I	56	-1908.0106	3928.0213	0.0000	1.0000
SYM	60	-1908.5076	3937.0153	0.0000	1.0000
K2P	56	-1920.6487	3953.2973	0.0000	1.0000
JC69	55	-1947.7092	4005.4184	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1822.6889	4022.6494	0.8494	0.8494
HKYIG	61	-1821.6751	4026.9096	0.1009	0.9504
HKYI	60	-1825.5564	4028.3843	0.0483	0.9986
GTRG	64	-1816.6717	4035.7664	0.0012	0.9998
GTRIG	65	-1815.8013	4040.3135	0.0001	1.0000
GTRI	64	-1820.1996	4042.8222	0.0000	1.0000
F81G	59	-1847.7102	4066.4040	0.0000	1.0000
F81IG	60	-1846.5967	4070.4649	0.0000	1.0000
F81I	59	-1850.2975	4071.5786	0.0000	1.0000
SYMG	61	-1861.2241	4106.0077	0.0000	1.0000
SYMIG	62	-1860.2171	4110.2815	0.0000	1.0000
K2PG	57	-1876.0231	4110.4541	0.0000	1.0000
SYMI	61	-1865.4634	4114.4862	0.0000	1.0000
HKY	59	-1871.7665	4114.5167	0.0000	1.0000
K2PIG	58	-1875.2846	4115.2651	0.0000	1.0000
K2PI	57	-1879.0826	4116.5731	0.0000	1.0000
GTR	63	-1868.4981	4133.1313	0.0000	1.0000
F81	58	-1893.5288	4151.7534	0.0000	1.0000
JC69G	56	-1905.7478	4163.6157	0.0000	1.0000
JC69I	56	-1908.0106	4168.1414	0.0000	1.0000
JC69IG	57	-1905.0962	4168.6004	0.0000	1.0000
K2P	56	-1920.6487	4193.4174	0.0000	1.0000
SYM	60	-1908.5076	4194.2868	0.0000	1.0000
JC69	55	-1947.7092	4241.2506	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/chr5_2303.phylip = 1 - 538;
 Partition Dummy = 1:183-phylip/chr5_2303.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 (GTRG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr5_2303.phylip = 1 - 538;
 Partition Dummy = 1:183-phylip/chr5_2303.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 BIC model (HKYG)]

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

