
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:12 04/04/2011

Input data from file "683-tomodel-phylip/chr2_8600.phylip" (7 taxa, 1008 characters)

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: HKYG


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2980.7576	6001.5152	0.7095	0.7095
GTRIG	21	-2980.7531	6003.5062	0.2622	0.9717
GTRI	20	-2984.1710	6008.3419	0.0234	0.9950
HKYG	16	-2990.0660	6012.1320	0.0035	0.9986
HKYIG	17	-2990.0654	6014.1308	0.0013	0.9999
HKYI	16	-2993.2267	6018.4534	0.0001	1.0000
GTR	19	-3004.9376	6047.8752	0.0000	1.0000
HKY	15	-3013.1990	6056.3979	0.0000	1.0000
SYMG	17	-3017.7380	6069.4760	0.0000	1.0000
K2PG	13	-3022.1020	6070.2040	0.0000	1.0000
SYMIG	18	-3017.7329	6071.4658	0.0000	1.0000
K2PIG	14	-3022.1002	6072.2005	0.0000	1.0000
SYMI	17	-3021.0690	6076.1381	0.0000	1.0000
K2PI	13	-3025.1990	6076.3979	0.0000	1.0000
SYM	16	-3040.9323	6113.8646	0.0000	1.0000
K2P	12	-3045.6224	6115.2448	0.0000	1.0000
F81G	15	-3052.7848	6135.5696	0.0000	1.0000
F81IG	16	-3052.7848	6137.5695	0.0000	1.0000
F81I	15	-3055.2375	6140.4750	0.0000	1.0000
F81	14	-3074.5184	6177.0368	0.0000	1.0000
JC69G	12	-3081.0687	6186.1375	0.0000	1.0000
JC69IG	13	-3081.0687	6188.1374	0.0000	1.0000
JC69I	12	-3083.6554	6191.3107	0.0000	1.0000
JC69	11	-3102.8281	6227.6562	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2990.0660	6090.7835	0.9213	0.9213
HKYI	16	-2993.2267	6097.1050	0.0391	0.9603
HKYIG	17	-2990.0654	6097.6981	0.0290	0.9894
GTRG	20	-2980.7576	6099.8297	0.0100	0.9994
GTRI	20	-2984.1710	6106.6564	0.0003	0.9997
GTRIG	21	-2980.7531	6106.7364	0.0003	1.0000
HKY	15	-3013.1990	6130.1338	0.0000	1.0000
K2PG	13	-3022.1020	6134.1084	0.0000	1.0000
K2PI	13	-3025.1990	6140.3023	0.0000	1.0000
K2PIG	14	-3022.1002	6141.0206	0.0000	1.0000
GTR	19	-3004.9376	6141.2739	0.0000	1.0000
SYMG	17	-3017.7380	6153.0433	0.0000	1.0000
SYMI	17	-3021.0690	6159.7054	0.0000	1.0000
SYMIG	18	-3017.7329	6159.9488	0.0000	1.0000
K2P	12	-3045.6224	6174.2335	0.0000	1.0000
SYM	16	-3040.9323	6192.5162	0.0000	1.0000
F81G	15	-3052.7848	6209.3055	0.0000	1.0000
F81I	15	-3055.2375	6214.2108	0.0000	1.0000
F81IG	16	-3052.7848	6216.2211	0.0000	1.0000
JC69G	12	-3081.0687	6245.1261	0.0000	1.0000
F81	14	-3074.5184	6245.8570	0.0000	1.0000
JC69I	12	-3083.6554	6250.2994	0.0000	1.0000
JC69IG	13	-3081.0687	6252.0418	0.0000	1.0000
JC69	11	-3102.8281	6281.7292	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 683-tomodel-phylip/chr2_8600.phylip = 1 - 1008;
 Partition Dummy = 1:683-tomodel-phylip/chr2_8600.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 683-tomodel-phylip/chr2_8600.phylip = 1 - 1008;
 Partition Dummy = 1:683-tomodel-phylip/chr2_8600.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 683-tomodel-phylip/chr2_8600.phylip = 1 - 1008;
 Partition Dummy = 1:683-tomodel-phylip/chr2_8600.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

