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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: K2PI


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	20	-3153.9012	6348.8306	0.4706	0.4706
GTRI	20	-3154.2910	6349.6101	0.3187	0.7893
GTRIG	21	-3153.6549	6350.4421	0.2102	0.9996
SYMG	17	-3165.4974	6365.7412	0.0001	0.9997
SYMI	17	-3165.5003	6365.7469	0.0001	0.9998
HKYI	16	-3166.8430	6366.3485	0.0001	0.9999
HKYG	16	-3167.1368	6366.9363	0.0001	0.9999
SYMIG	18	-3165.1040	6367.0431	0.0001	1.0000
HKYIG	17	-3166.5690	6367.8843	0.0000	1.0000
K2PI	13	-3174.4444	6375.3306	0.0000	1.0000
K2PG	13	-3174.8317	6376.1052	0.0000	1.0000
K2PIG	14	-3174.2068	6376.9239	0.0000	1.0000
GTR	19	-3178.3544	6395.6379	0.0000	1.0000
SYM	16	-3187.6699	6408.0025	0.0000	1.0000
HKY	15	-3189.9890	6410.5620	0.0000	1.0000
K2P	12	-3197.3231	6419.0245	0.0000	1.0000
F81G	15	-3238.3684	6507.3206	0.0000	1.0000
F81I	15	-3238.3798	6507.3436	0.0000	1.0000
F81IG	16	-3237.9943	6508.6511	0.0000	1.0000
JC69I	12	-3244.2038	6512.7859	0.0000	1.0000
JC69G	12	-3244.2115	6512.8012	0.0000	1.0000
JC69IG	13	-3243.8293	6514.1003	0.0000	1.0000
F81	14	-3259.4832	6547.4766	0.0000	1.0000
JC69	11	-3264.9382	6552.1961	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-3153.9012	6347.8024	0.4654	0.4654
GTRI	20	-3154.2910	6348.5819	0.3152	0.7806
GTRIG	21	-3153.6549	6349.3097	0.2190	0.9996
SYMG	17	-3165.4974	6364.9949	0.0001	0.9997
SYMI	17	-3165.5003	6365.0006	0.0001	0.9998
HKYI	16	-3166.8430	6365.6859	0.0001	0.9999
SYMIG	18	-3165.1040	6366.2079	0.0000	0.9999
HKYG	16	-3167.1368	6366.2737	0.0000	1.0000
HKYIG	17	-3166.5690	6367.1380	0.0000	1.0000
K2PI	13	-3174.4444	6374.8888	0.0000	1.0000
K2PG	13	-3174.8317	6375.6635	0.0000	1.0000
K2PIG	14	-3174.2068	6376.4135	0.0000	1.0000
GTR	19	-3178.3544	6394.7088	0.0000	1.0000
SYM	16	-3187.6699	6407.3399	0.0000	1.0000
HKY	15	-3189.9890	6409.9780	0.0000	1.0000
K2P	12	-3197.3231	6418.6463	0.0000	1.0000
F81G	15	-3238.3684	6506.7367	0.0000	1.0000
F81I	15	-3238.3798	6506.7597	0.0000	1.0000
F81IG	16	-3237.9943	6507.9885	0.0000	1.0000
JC69I	12	-3244.2038	6512.4077	0.0000	1.0000
JC69G	12	-3244.2115	6512.4230	0.0000	1.0000
JC69IG	13	-3243.8293	6513.6586	0.0000	1.0000
F81	14	-3259.4832	6546.9663	0.0000	1.0000
JC69	11	-3264.9382	6551.8764	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	13	-3174.4444	6436.3921	0.5058	0.5058
K2PG	13	-3174.8317	6437.1667	0.3433	0.8491
HKYI	16	-3166.8430	6441.3822	0.0417	0.8908
HKYG	16	-3167.1368	6441.9700	0.0311	0.9219
GTRG	20	-3153.9012	6442.4228	0.0248	0.9467
K2PIG	14	-3174.2068	6442.6478	0.0222	0.9689
GTRI	20	-3154.2910	6443.2023	0.0168	0.9857
SYMG	17	-3165.4974	6445.4222	0.0055	0.9912
SYMI	17	-3165.5003	6445.4279	0.0055	0.9967
HKYIG	17	-3166.5690	6447.5653	0.0019	0.9986
GTRIG	21	-3153.6549	6448.6611	0.0011	0.9997
SYMIG	18	-3165.1040	6451.3662	0.0003	1.0000
K2P	12	-3197.3231	6475.4185	0.0000	1.0000
HKY	15	-3189.9890	6480.9433	0.0000	1.0000
SYM	16	-3187.6699	6483.0361	0.0000	1.0000
GTR	19	-3178.3544	6484.5981	0.0000	1.0000
JC69I	12	-3244.2038	6569.1799	0.0000	1.0000
JC69G	12	-3244.2115	6569.1952	0.0000	1.0000
JC69IG	13	-3243.8293	6575.1618	0.0000	1.0000
F81G	15	-3238.3684	6577.7020	0.0000	1.0000
F81I	15	-3238.3798	6577.7249	0.0000	1.0000
F81IG	16	-3237.9943	6583.6848	0.0000	1.0000
JC69	11	-3264.9382	6603.9176	0.0000	1.0000
F81	14	-3259.4832	6613.2006	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 591-tomodel-phylip/chr9_6388.phylip = 1 - 838;
 Partition Dummy = 1:591-tomodel-phylip/chr9_6388.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 591-tomodel-phylip/chr9_6388.phylip = 1 - 838;
 Partition Dummy = 1:591-tomodel-phylip/chr9_6388.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 (K2PI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr9_6388.phylip = 1 - 838;
 Partition Dummy = 1:591-tomodel-phylip/chr9_6388.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

