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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	60	-1971.1190	4082.5150	0.5812	0.5812
HKYIG	61	-1970.4612	4083.9335	0.2860	0.8671
HKYI	60	-1973.0702	4086.4173	0.0826	0.9497
K2PG	57	-1978.2776	4088.7200	0.0261	0.9758
K2PIG	58	-1977.8908	4090.6356	0.0100	0.9859
GTRG	64	-1969.9653	4091.2359	0.0074	0.9933
GTRIG	65	-1969.4553	4093.0117	0.0031	0.9963
K2PI	57	-1980.6847	4093.5342	0.0024	0.9987
GTRI	64	-1971.8958	4095.0970	0.0011	0.9998
SYMG	61	-1978.0045	4099.0202	0.0002	0.9999
SYMIG	62	-1977.5753	4100.9111	0.0001	1.0000
SYMI	61	-1980.3170	4103.6451	0.0000	1.0000
HKY	59	-1991.0012	4119.5604	0.0000	1.0000
GTR	63	-1989.3898	4127.3048	0.0000	1.0000
K2P	56	-1999.4229	4128.3362	0.0000	1.0000
SYM	60	-1999.1380	4138.5530	0.0000	1.0000
F81G	59	-2025.5562	4188.6704	0.0000	1.0000
F81IG	60	-2025.5609	4191.3989	0.0000	1.0000
F81I	59	-2027.2009	4191.9599	0.0000	1.0000
JC69G	56	-2031.9836	4193.4576	0.0000	1.0000
JC69IG	57	-2031.9832	4196.1312	0.0000	1.0000
JC69I	56	-2033.7067	4196.9037	0.0000	1.0000
F81	58	-2044.4769	4223.8078	0.0000	1.0000
JC69	55	-2051.6015	4230.0337	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1971.1190	4062.2380	0.5049	0.5049
HKYIG	61	-1970.4612	4062.9224	0.3586	0.8635
HKYI	60	-1973.0702	4066.1403	0.0718	0.9353
GTRG	64	-1969.9653	4067.9306	0.0293	0.9646
GTRIG	65	-1969.4553	4068.9106	0.0180	0.9826
K2PG	57	-1978.2776	4070.5551	0.0079	0.9904
K2PIG	58	-1977.8908	4071.7816	0.0043	0.9947
GTRI	64	-1971.8958	4071.7917	0.0043	0.9990
K2PI	57	-1980.6847	4075.3693	0.0007	0.9997
SYMG	61	-1978.0045	4078.0091	0.0002	0.9999
SYMIG	62	-1977.5753	4079.1506	0.0001	1.0000
SYMI	61	-1980.3170	4082.6339	0.0000	1.0000
HKY	59	-1991.0012	4100.0023	0.0000	1.0000
GTR	63	-1989.3898	4104.7796	0.0000	1.0000
K2P	56	-1999.4229	4110.8457	0.0000	1.0000
SYM	60	-1999.1380	4118.2760	0.0000	1.0000
F81G	59	-2025.5562	4169.1124	0.0000	1.0000
F81IG	60	-2025.5609	4171.1219	0.0000	1.0000
F81I	59	-2027.2009	4172.4019	0.0000	1.0000
JC69G	56	-2031.9836	4175.9671	0.0000	1.0000
JC69IG	57	-2031.9832	4177.9664	0.0000	1.0000
JC69I	56	-2033.7067	4179.4133	0.0000	1.0000
F81	58	-2044.4769	4204.9538	0.0000	1.0000
JC69	55	-2051.6015	4213.2031	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1978.2776	4301.1204	0.7435	0.7435
HKYG	60	-1971.1190	4304.9383	0.1102	0.8537
K2PI	57	-1980.6847	4305.9346	0.0670	0.9207
K2PIG	58	-1977.8908	4306.3919	0.0533	0.9740
HKYI	60	-1973.0702	4308.8406	0.0157	0.9896
HKYIG	61	-1970.4612	4309.6677	0.0104	1.0000
SYMG	61	-1978.0045	4324.7544	0.0000	1.0000
GTRG	64	-1969.9653	4326.8109	0.0000	1.0000
SYMI	61	-1980.3170	4329.3793	0.0000	1.0000
SYMIG	62	-1977.5753	4329.9410	0.0000	1.0000
GTRI	64	-1971.8958	4330.6720	0.0000	1.0000
GTRIG	65	-1969.4553	4331.8360	0.0000	1.0000
K2P	56	-1999.4229	4337.3660	0.0000	1.0000
HKY	59	-1991.0012	4338.6577	0.0000	1.0000
GTR	63	-1989.3898	4359.6150	0.0000	1.0000
SYM	60	-1999.1380	4360.9763	0.0000	1.0000
JC69G	56	-2031.9836	4402.4874	0.0000	1.0000
JC69I	56	-2033.7067	4405.9336	0.0000	1.0000
F81G	59	-2025.5562	4407.7677	0.0000	1.0000
JC69IG	57	-2031.9832	4408.5317	0.0000	1.0000
F81I	59	-2027.2009	4411.0572	0.0000	1.0000
F81IG	60	-2025.5609	4413.8222	0.0000	1.0000
JC69	55	-2051.6015	4435.6784	0.0000	1.0000
F81	58	-2044.4769	4439.5641	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_4777.phylip = 1 - 422;
 Partition Dummy = 1:183-phylip/chr7_4777.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_4777.phylip = 1 - 422;
 Partition Dummy = 1:183-phylip/chr7_4777.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 (K2PG)]

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

