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

Input data from file "183-phylip/chr9_5246.phylip" (29 taxa, 738 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	-1950.9952	4042.3530	0.4457	0.4457
GTRI	64	-1951.3439	4043.0503	0.3145	0.7601
GTRIG	65	-1950.7271	4044.2221	0.1750	0.9352
HKYG	60	-1958.4019	4047.6162	0.0321	0.9672
HKYI	60	-1959.0677	4048.9478	0.0165	0.9837
HKYIG	61	-1957.8942	4048.9777	0.0162	0.9999
K2PG	57	-1969.1292	4061.9819	0.0000	1.0000
K2PI	57	-1969.6413	4063.0061	0.0000	1.0000
K2PIG	58	-1968.7924	4063.6643	0.0000	1.0000
SYMG	61	-1966.5525	4066.2944	0.0000	1.0000
SYMI	61	-1967.0105	4067.2103	0.0000	1.0000
SYMIG	62	-1966.2388	4068.0509	0.0000	1.0000
F81G	59	-1976.9116	4082.2658	0.0000	1.0000
F81I	59	-1977.3457	4083.1338	0.0000	1.0000
F81IG	60	-1976.6382	4084.0889	0.0000	1.0000
GTR	63	-1973.2066	4084.3776	0.0000	1.0000
HKY	59	-1979.2585	4086.9594	0.0000	1.0000
JC69G	56	-1987.3674	4096.1093	0.0000	1.0000
JC69I	56	-1987.8591	4097.0925	0.0000	1.0000
JC69IG	57	-1987.0648	4097.8531	0.0000	1.0000
K2P	56	-1989.1236	4099.6216	0.0000	1.0000
SYM	60	-1987.0277	4104.8679	0.0000	1.0000
F81	58	-1996.6504	4119.3803	0.0000	1.0000
JC69	55	-2006.8857	4132.8037	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1950.9952	4029.9904	0.4430	0.4430
GTRI	64	-1951.3439	4030.6877	0.3126	0.7557
GTRIG	65	-1950.7271	4031.4542	0.2131	0.9688
HKYG	60	-1958.4019	4036.8038	0.0147	0.9835
HKYIG	61	-1957.8942	4037.7883	0.0090	0.9924
HKYI	60	-1959.0677	4038.1354	0.0075	1.0000
K2PG	57	-1969.1292	4052.2584	0.0000	1.0000
K2PI	57	-1969.6413	4053.2826	0.0000	1.0000
K2PIG	58	-1968.7924	4053.5848	0.0000	1.0000
SYMG	61	-1966.5525	4055.1051	0.0000	1.0000
SYMI	61	-1967.0105	4056.0210	0.0000	1.0000
SYMIG	62	-1966.2388	4056.4776	0.0000	1.0000
F81G	59	-1976.9116	4071.8233	0.0000	1.0000
GTR	63	-1973.2066	4072.4132	0.0000	1.0000
F81I	59	-1977.3457	4072.6913	0.0000	1.0000
F81IG	60	-1976.6382	4073.2765	0.0000	1.0000
HKY	59	-1979.2585	4076.5170	0.0000	1.0000
JC69G	56	-1987.3674	4086.7349	0.0000	1.0000
JC69I	56	-1987.8591	4087.7181	0.0000	1.0000
JC69IG	57	-1987.0648	4088.1295	0.0000	1.0000
K2P	56	-1989.1236	4090.2471	0.0000	1.0000
SYM	60	-1987.0277	4094.0555	0.0000	1.0000
F81	58	-1996.6504	4109.3008	0.0000	1.0000
JC69	55	-2006.8857	4123.7715	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1958.4019	4313.0404	0.4336	0.4336
HKYI	60	-1959.0677	4314.3720	0.2228	0.6564
K2PG	57	-1969.1292	4314.6832	0.1907	0.8471
K2PI	57	-1969.6413	4315.7074	0.1143	0.9613
HKYIG	61	-1957.8942	4318.6289	0.0265	0.9879
K2PIG	58	-1968.7924	4320.6135	0.0098	0.9977
GTRG	64	-1950.9952	4324.6428	0.0013	0.9990
GTRI	64	-1951.3439	4325.3401	0.0009	0.9999
GTRIG	65	-1950.7271	4330.7106	0.0001	1.0000
SYMG	61	-1966.5525	4335.9457	0.0000	1.0000
SYMI	61	-1967.0105	4336.8616	0.0000	1.0000
SYMIG	62	-1966.2388	4341.9221	0.0000	1.0000
F81G	59	-1976.9116	4343.4560	0.0000	1.0000
F81I	59	-1977.3457	4344.3240	0.0000	1.0000
JC69G	56	-1987.3674	4344.5557	0.0000	1.0000
JC69I	56	-1987.8591	4345.5390	0.0000	1.0000
K2P	56	-1989.1236	4348.0680	0.0000	1.0000
HKY	59	-1979.2585	4348.1496	0.0000	1.0000
F81IG	60	-1976.6382	4349.5131	0.0000	1.0000
JC69IG	57	-1987.0648	4350.5543	0.0000	1.0000
GTR	63	-1973.2066	4362.4617	0.0000	1.0000
SYM	60	-1987.0277	4370.2921	0.0000	1.0000
F81	58	-1996.6504	4376.3295	0.0000	1.0000
JC69	55	-2006.8857	4376.9884	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/chr9_5246.phylip = 1 - 738;
 Partition Dummy = 1:183-phylip/chr9_5246.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/chr9_5246.phylip = 1 - 738;
 Partition Dummy = 1:183-phylip/chr9_5246.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/chr9_5246.phylip = 1 - 738;
 Partition Dummy = 1:183-phylip/chr9_5246.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

