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

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

Minimum AIC  model: GTRG
Minimum AICc model: K2P
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
K2P	56	-299.7842	850.3509	0.8335	0.8335
JC69	55	-307.2625	855.5887	0.0607	0.8943
K2PG	57	-297.4355	855.8044	0.0545	0.9488
K2PI	57	-297.6811	856.2956	0.0427	0.9915
JC69G	56	-304.9843	860.7511	0.0046	0.9961
JC69I	56	-305.2191	861.2208	0.0036	0.9997
K2PIG	58	-297.4331	866.4117	0.0003	1.0000
JC69IG	57	-304.9840	870.9014	0.0000	1.0000
HKY	59	-296.4605	875.5721	0.0000	1.0000
F81	58	-304.0729	879.6913	0.0000	1.0000
HKYG	60	-294.0352	882.3561	0.0000	1.0000
HKYI	60	-294.2026	882.6909	0.0000	1.0000
SYM	60	-294.9917	884.2692	0.0000	1.0000
F81G	59	-301.8824	886.4160	0.0000	1.0000
F81I	59	-302.0725	886.7961	0.0000	1.0000
SYMG	61	-292.8901	892.2679	0.0000	1.0000
SYMI	61	-293.0337	892.5553	0.0000	1.0000
HKYIG	61	-294.0395	894.5668	0.0000	1.0000
F81IG	60	-301.8861	898.0578	0.0000	1.0000
SYMIG	62	-292.8934	905.0869	0.0000	1.0000
GTR	63	-290.8004	914.3701	0.0000	1.0000
GTRG	64	-288.6944	924.3361	0.0000	1.0000
GTRI	64	-288.7752	924.4977	0.0000	1.0000
GTRIG	65	-288.6984	939.2887	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-288.6944	705.3888	0.2242	0.2242
GTRI	64	-288.7752	705.5504	0.2068	0.4309
GTRIG	65	-288.6984	707.3968	0.0821	0.5130
GTR	63	-290.8004	707.6008	0.0742	0.5872
SYMG	61	-292.8901	707.7801	0.0678	0.6550
SYMI	61	-293.0337	708.0675	0.0587	0.7137
HKYG	60	-294.0352	708.0703	0.0586	0.7724
HKYI	60	-294.2026	708.4052	0.0496	0.8220
K2PG	57	-297.4355	708.8711	0.0393	0.8613
K2PI	57	-297.6811	709.3623	0.0307	0.8920
SYMIG	62	-292.8934	709.7869	0.0249	0.9169
SYM	60	-294.9917	709.9835	0.0225	0.9394
HKYIG	61	-294.0395	710.0790	0.0215	0.9609
K2PIG	58	-297.4331	710.8663	0.0145	0.9754
HKY	59	-296.4605	710.9210	0.0141	0.9895
K2P	56	-299.7842	711.5683	0.0102	0.9997
F81G	59	-301.8824	721.7648	0.0001	0.9998
JC69G	56	-304.9843	721.9685	0.0001	0.9998
F81I	59	-302.0725	722.1449	0.0001	0.9999
JC69I	56	-305.2191	722.4382	0.0000	0.9999
F81IG	60	-301.8861	723.7721	0.0000	0.9999
JC69IG	57	-304.9840	723.9680	0.0000	1.0000
F81	58	-304.0729	724.1459	0.0000	1.0000
JC69	55	-307.2625	724.5249	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-297.4355	859.0506	0.3358	0.3358
K2P	56	-299.7842	859.1132	0.3255	0.6613
K2PI	57	-297.6811	859.5418	0.2627	0.9239
K2PIG	58	-297.4331	863.6806	0.0332	0.9571
HKYG	60	-294.0352	866.1541	0.0096	0.9667
HKY	59	-296.4605	866.3700	0.0086	0.9754
HKYI	60	-294.2026	866.4889	0.0081	0.9835
SYM	60	-294.9917	868.0672	0.0037	0.9872
SYMG	61	-292.8901	868.4986	0.0030	0.9902
SYMI	61	-293.0337	868.7859	0.0026	0.9928
JC69	55	-307.2625	869.4350	0.0019	0.9946
JC69G	56	-304.9843	869.5133	0.0018	0.9964
JC69I	56	-305.2191	869.9830	0.0014	0.9979
HKYIG	61	-294.0395	870.7975	0.0009	0.9988
SYMIG	62	-292.8934	873.1401	0.0003	0.9991
GTR	63	-290.8004	873.5887	0.0002	0.9993
GTRG	64	-288.6944	874.0114	0.0002	0.9995
JC69IG	57	-304.9840	874.1476	0.0002	0.9997
GTRI	64	-288.7752	874.1730	0.0002	0.9999
F81	58	-304.0729	876.9602	0.0000	0.9999
F81G	59	-301.8824	877.2139	0.0000	0.9999
F81I	59	-302.0725	877.5939	0.0000	1.0000
GTRIG	65	-288.6984	878.6542	0.0000	1.0000
F81IG	60	-301.8861	881.8558	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/chr2_25984.phylip = 1 - 103;
 Partition Dummy = 1:183-phylip/chr2_25984.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 (K2P)]

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


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

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

