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

Input data from file "183-phylip/chr8_4319.phylip" (29 taxa, 324 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	-1738.1229	3624.0786	0.6706	0.6706
HKYIG	61	-1737.6345	3626.1393	0.2393	0.9100
GTRG	64	-1734.2260	3628.5756	0.0708	0.9808
GTRIG	65	-1733.9636	3631.1831	0.0192	1.0000
HKYI	60	-1752.1245	3652.0818	0.0000	1.0000
K2PG	57	-1759.5860	3658.0291	0.0000	1.0000
K2PIG	58	-1759.0446	3659.9157	0.0000	1.0000
GTRI	64	-1751.1401	3662.4037	0.0000	1.0000
SYMG	61	-1757.3309	3665.5320	0.0000	1.0000
SYMIG	62	-1756.7419	3667.4148	0.0000	1.0000
K2PI	57	-1768.7674	3676.3919	0.0000	1.0000
SYMI	61	-1767.3244	3685.5190	0.0000	1.0000
F81G	59	-1803.7034	3752.2251	0.0000	1.0000
F81IG	60	-1803.1295	3754.0918	0.0000	1.0000
HKY	59	-1804.6869	3754.1920	0.0000	1.0000
GTR	63	-1802.4714	3761.9581	0.0000	1.0000
F81I	59	-1812.2636	3769.3454	0.0000	1.0000
JC69G	56	-1821.5518	3779.0137	0.0000	1.0000
JC69IG	57	-1821.0343	3780.9258	0.0000	1.0000
K2P	56	-1824.2998	3784.5097	0.0000	1.0000
SYM	60	-1820.5568	3788.9463	0.0000	1.0000
JC69I	56	-1829.4485	3794.8070	0.0000	1.0000
F81	58	-1866.9836	3875.7936	0.0000	1.0000
JC69	55	-1882.8226	3898.6302	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1738.1229	3596.2459	0.3409	0.3409
GTRG	64	-1734.2260	3596.4520	0.3076	0.6485
HKYIG	61	-1737.6345	3597.2691	0.2044	0.8529
GTRIG	65	-1733.9636	3597.9272	0.1471	1.0000
HKYI	60	-1752.1245	3624.2491	0.0000	1.0000
GTRI	64	-1751.1401	3630.2801	0.0000	1.0000
K2PG	57	-1759.5860	3633.1720	0.0000	1.0000
K2PIG	58	-1759.0446	3634.0893	0.0000	1.0000
SYMG	61	-1757.3309	3636.6618	0.0000	1.0000
SYMIG	62	-1756.7419	3637.4837	0.0000	1.0000
K2PI	57	-1768.7674	3651.5347	0.0000	1.0000
SYMI	61	-1767.3244	3656.6488	0.0000	1.0000
F81G	59	-1803.7034	3725.4069	0.0000	1.0000
F81IG	60	-1803.1295	3726.2591	0.0000	1.0000
HKY	59	-1804.6869	3727.3738	0.0000	1.0000
GTR	63	-1802.4714	3730.9427	0.0000	1.0000
F81I	59	-1812.2636	3742.5272	0.0000	1.0000
JC69G	56	-1821.5518	3755.1036	0.0000	1.0000
JC69IG	57	-1821.0343	3756.0686	0.0000	1.0000
K2P	56	-1824.2998	3760.5996	0.0000	1.0000
SYM	60	-1820.5568	3761.1136	0.0000	1.0000
JC69I	56	-1829.4485	3770.8969	0.0000	1.0000
F81	58	-1866.9836	3849.9672	0.0000	1.0000
JC69	55	-1882.8226	3875.6451	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1738.1229	3823.0905	0.9166	0.9166
HKYIG	61	-1737.6345	3827.8944	0.0830	0.9995
GTRG	64	-1734.2260	3838.4196	0.0004	1.0000
GTRIG	65	-1733.9636	3843.6756	0.0000	1.0000
K2PG	57	-1759.5860	3848.6743	0.0000	1.0000
HKYI	60	-1752.1245	3851.0937	0.0000	1.0000
K2PIG	58	-1759.0446	3853.3724	0.0000	1.0000
K2PI	57	-1768.7674	3867.0371	0.0000	1.0000
SYMG	61	-1757.3309	3867.2872	0.0000	1.0000
SYMIG	62	-1756.7419	3871.8898	0.0000	1.0000
GTRI	64	-1751.1401	3872.2477	0.0000	1.0000
SYMI	61	-1767.3244	3887.2741	0.0000	1.0000
F81G	59	-1803.7034	3948.4707	0.0000	1.0000
HKY	59	-1804.6869	3950.4377	0.0000	1.0000
F81IG	60	-1803.1295	3953.1037	0.0000	1.0000
F81I	59	-1812.2636	3965.5911	0.0000	1.0000
JC69G	56	-1821.5518	3966.8252	0.0000	1.0000
GTR	63	-1802.4714	3969.1296	0.0000	1.0000
JC69IG	57	-1821.0343	3971.5710	0.0000	1.0000
K2P	56	-1824.2998	3972.3212	0.0000	1.0000
JC69I	56	-1829.4485	3982.6186	0.0000	1.0000
SYM	60	-1820.5568	3987.9582	0.0000	1.0000
F81	58	-1866.9836	4069.2503	0.0000	1.0000
JC69	55	-1882.8226	4083.5860	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/chr8_4319.phylip = 1 - 324;
 Partition Dummy = 1:183-phylip/chr8_4319.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/chr8_4319.phylip = 1 - 324;
 Partition Dummy = 1:183-phylip/chr8_4319.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/chr8_4319.phylip = 1 - 324;
 Partition Dummy = 1:183-phylip/chr8_4319.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

