
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
21:35 04/04/2011

Input data from file "591-tomodel-phylip/chr2_26586.phylip" (7 taxa, 1064 characters)

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: K2PG


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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-3437.8710	6907.7420	0.4294	0.4294
K2PG	13	-3441.5233	6909.0466	0.2236	0.6530
HKYIG	17	-3437.8701	6909.7401	0.1581	0.8111
K2PIG	14	-3441.5350	6911.0700	0.0813	0.8924
GTRG	20	-3436.4060	6912.8119	0.0340	0.9265
SYMG	17	-3439.7953	6913.5907	0.0231	0.9495
HKYI	16	-3441.0334	6914.0669	0.0182	0.9677
GTRIG	21	-3436.4078	6914.8155	0.0125	0.9802
K2PI	13	-3444.6705	6915.3409	0.0096	0.9898
SYMIG	18	-3439.7935	6915.5869	0.0085	0.9983
GTRI	20	-3439.9371	6919.8742	0.0010	0.9993
SYMI	17	-3443.3023	6920.6046	0.0007	1.0000
HKY	15	-3469.2658	6968.5316	0.0000	1.0000
GTR	19	-3465.7040	6969.4080	0.0000	1.0000
K2P	12	-3473.0650	6970.1301	0.0000	1.0000
SYM	16	-3471.0165	6974.0330	0.0000	1.0000
F81G	15	-3505.3974	7040.7948	0.0000	1.0000
JC69G	12	-3508.6186	7041.2371	0.0000	1.0000
F81IG	16	-3505.3973	7042.7946	0.0000	1.0000
JC69IG	13	-3508.6185	7043.2370	0.0000	1.0000
F81I	15	-3507.9065	7045.8130	0.0000	1.0000
JC69I	12	-3511.1126	7046.2252	0.0000	1.0000
F81	14	-3533.7907	7095.5813	0.0000	1.0000
JC69	11	-3537.0361	7096.0721	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-3441.5233	6973.6538	0.9307	0.9307
K2PI	13	-3444.6705	6979.9482	0.0400	0.9707
K2PIG	14	-3441.5350	6980.6471	0.0282	0.9989
HKYG	16	-3437.8710	6987.2586	0.0010	0.9999
HKYI	16	-3441.0334	6993.5835	0.0000	1.0000
HKYIG	17	-3437.8701	6994.2266	0.0000	1.0000
SYMG	17	-3439.7953	6998.0771	0.0000	1.0000
SYMIG	18	-3439.7935	7005.0432	0.0000	1.0000
SYMI	17	-3443.3023	7005.0911	0.0000	1.0000
GTRG	20	-3436.4060	7012.2077	0.0000	1.0000
GTRIG	21	-3436.4078	7019.1811	0.0000	1.0000
GTRI	20	-3439.9371	7019.2700	0.0000	1.0000
K2P	12	-3473.0650	7029.7675	0.0000	1.0000
HKY	15	-3469.2658	7043.0785	0.0000	1.0000
SYM	16	-3471.0165	7053.5496	0.0000	1.0000
GTR	19	-3465.7040	7063.8341	0.0000	1.0000
JC69G	12	-3508.6186	7100.8746	0.0000	1.0000
JC69I	12	-3511.1126	7105.8627	0.0000	1.0000
JC69IG	13	-3508.6185	7107.8442	0.0000	1.0000
F81G	15	-3505.3974	7115.3417	0.0000	1.0000
F81I	15	-3507.9065	7120.3599	0.0000	1.0000
F81IG	16	-3505.3973	7122.3113	0.0000	1.0000
JC69	11	-3537.0361	7150.7398	0.0000	1.0000
F81	14	-3533.7907	7165.1584	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 591-tomodel-phylip/chr2_26586.phylip = 1 - 1064;
 Partition Dummy = 1:591-tomodel-phylip/chr2_26586.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 591-tomodel-phylip/chr2_26586.phylip = 1 - 1064;
 Partition Dummy = 1:591-tomodel-phylip/chr2_26586.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 591-tomodel-phylip/chr2_26586.phylip = 1 - 1064;
 Partition Dummy = 1:591-tomodel-phylip/chr2_26586.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

