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

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

Minimum AIC  model: GTR
Minimum AICc model: HKY
Minimum BIC  model: K2P


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
HKY	59	-336.0914	857.6114	0.4609	0.4609
K2P	56	-343.7284	858.5680	0.2857	0.7466
F81	58	-340.7269	862.0198	0.0509	0.7975
JC69	55	-347.9288	862.3713	0.0427	0.8402
HKYI	60	-336.0906	862.5658	0.0387	0.8789
HKYG	60	-336.0923	862.5691	0.0386	0.9175
K2PI	57	-343.7275	863.2494	0.0275	0.9450
K2PG	57	-343.7292	863.2529	0.0275	0.9725
F81I	59	-340.7263	866.8811	0.0045	0.9769
F81G	59	-340.7278	866.8842	0.0045	0.9814
JC69I	56	-347.9280	866.9672	0.0043	0.9857
JC69G	56	-347.9296	866.9704	0.0043	0.9900
HKYIG	61	-336.0996	867.6361	0.0031	0.9930
GTR	63	-331.0247	867.8910	0.0027	0.9957
K2PIG	58	-343.7359	868.0378	0.0025	0.9983
SYM	60	-340.4375	871.2596	0.0005	0.9988
JC69IG	57	-347.9367	871.6678	0.0004	0.9992
F81IG	60	-340.7355	871.8555	0.0004	0.9995
GTRI	64	-331.0261	873.2521	0.0002	0.9997
GTRG	64	-331.0277	873.2553	0.0002	0.9999
SYMI	61	-340.4369	876.3107	0.0000	0.9999
SYMG	61	-340.4386	876.3141	0.0000	1.0000
GTRIG	65	-331.0380	878.7426	0.0000	1.0000
SYMIG	62	-340.4468	881.4819	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTR	63	-331.0247	788.0495	0.3940	0.3940
GTRI	64	-331.0261	790.0521	0.1448	0.5388
GTRG	64	-331.0277	790.0553	0.1445	0.6833
HKY	59	-336.0914	790.1829	0.1356	0.8189
GTRIG	65	-331.0380	792.0760	0.0526	0.8716
HKYI	60	-336.0906	792.1812	0.0499	0.9215
HKYG	60	-336.0923	792.1845	0.0498	0.9713
HKYIG	61	-336.0996	794.1992	0.0182	0.9895
F81	58	-340.7269	797.4538	0.0036	0.9931
F81I	59	-340.7263	799.4525	0.0013	0.9944
F81G	59	-340.7278	799.4556	0.0013	0.9958
K2P	56	-343.7284	799.4569	0.0013	0.9971
SYM	60	-340.4375	800.8750	0.0006	0.9977
K2PI	57	-343.7275	801.4550	0.0005	0.9982
K2PG	57	-343.7292	801.4585	0.0005	0.9987
F81IG	60	-340.7355	801.4709	0.0005	0.9992
SYMI	61	-340.4369	802.8738	0.0002	0.9994
SYMG	61	-340.4386	802.8772	0.0002	0.9996
K2PIG	58	-343.7359	803.4718	0.0002	0.9998
SYMIG	62	-340.4468	804.8937	0.0001	0.9999
JC69	55	-347.9288	805.8576	0.0001	1.0000
JC69I	56	-347.9280	807.8561	0.0000	1.0000
JC69G	56	-347.9296	807.8592	0.0000	1.0000
JC69IG	57	-347.9367	809.8734	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2P	56	-343.7284	973.3898	0.3743	0.3743
HKY	59	-336.0914	973.4336	0.3662	0.7406
JC69	55	-347.9288	976.6846	0.0721	0.8127
F81	58	-340.7269	977.5986	0.0456	0.8583
K2PI	57	-343.7275	978.4939	0.0292	0.8875
K2PG	57	-343.7292	978.4974	0.0291	0.9166
HKYI	60	-336.0906	978.5379	0.0285	0.9451
HKYG	60	-336.0923	978.5412	0.0285	0.9736
JC69I	56	-347.9280	981.7890	0.0056	0.9792
JC69G	56	-347.9296	981.7922	0.0056	0.9848
F81I	59	-340.7263	982.7033	0.0036	0.9884
F81G	59	-340.7278	982.7064	0.0035	0.9919
K2PIG	58	-343.7359	983.6166	0.0023	0.9942
HKYIG	61	-336.0996	983.6619	0.0022	0.9964
GTR	63	-331.0247	983.7240	0.0021	0.9985
JC69IG	57	-347.9367	986.9123	0.0004	0.9990
SYM	60	-340.4375	987.2317	0.0004	0.9993
F81IG	60	-340.7355	987.8276	0.0003	0.9996
GTRI	64	-331.0261	988.8327	0.0002	0.9998
GTRG	64	-331.0277	988.8358	0.0002	0.9999
SYMI	61	-340.4369	992.3365	0.0000	1.0000
SYMG	61	-340.4386	992.3399	0.0000	1.0000
GTRIG	65	-331.0380	993.9624	0.0000	1.0000
SYMIG	62	-340.4468	997.4623	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (GTR)]

BEGIN MRBAYES;
 Charset 183-phylip/chrZ_4787.phylip = 1 - 165;
 Partition Dummy = 1:183-phylip/chrZ_4787.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=equal;
 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);
END;


[Mrbayes block for the best AICc model (HKY)]

BEGIN MRBAYES;
 Charset 183-phylip/chrZ_4787.phylip = 1 - 165;
 Partition Dummy = 1:183-phylip/chrZ_4787.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


[Mrbayes block for the best BIC model (K2P)]

BEGIN MRBAYES;
 Charset 183-phylip/chrZ_4787.phylip = 1 - 165;
 Partition Dummy = 1:183-phylip/chrZ_4787.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0);
END;


-------------------------------------------------------------
End of Output

