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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-3096.5110	6225.0220	0.4562	0.4562
HKYIG	17	-3096.4609	6226.9219	0.1764	0.6326
HKYI	16	-3097.6946	6227.3893	0.1397	0.7723
GTRG	20	-3093.7315	6227.4631	0.1346	0.9069
GTRIG	21	-3093.6663	6229.3326	0.0529	0.9598
GTRI	20	-3094.9398	6229.8796	0.0402	1.0000
HKY	15	-3121.7935	6273.5869	0.0000	1.0000
GTR	19	-3118.4567	6274.9135	0.0000	1.0000
K2PG	13	-3131.1070	6288.2139	0.0000	1.0000
K2PIG	14	-3131.0547	6290.1094	0.0000	1.0000
K2PI	13	-3132.5494	6291.0988	0.0000	1.0000
SYMG	17	-3130.5556	6295.1112	0.0000	1.0000
SYMIG	18	-3130.5028	6297.0056	0.0000	1.0000
SYMI	17	-3132.1967	6298.3934	0.0000	1.0000
K2P	12	-3158.6927	6341.3855	0.0000	1.0000
SYM	16	-3158.0415	6348.0829	0.0000	1.0000
F81G	15	-3190.5287	6411.0573	0.0000	1.0000
F81IG	16	-3190.4911	6412.9822	0.0000	1.0000
F81I	15	-3192.2734	6414.5468	0.0000	1.0000
F81	14	-3215.0911	6458.1823	0.0000	1.0000
JC69G	12	-3219.5278	6463.0555	0.0000	1.0000
JC69IG	13	-3219.4812	6464.9624	0.0000	1.0000
JC69I	12	-3221.4681	6466.9362	0.0000	1.0000
JC69	11	-3245.3959	6512.7919	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3096.5110	6303.7211	0.7467	0.7467
HKYI	16	-3097.6946	6306.0884	0.2286	0.9753
HKYIG	17	-3096.4609	6310.5397	0.0247	1.0000
GTRG	20	-3093.7315	6325.8370	0.0000	1.0000
GTRI	20	-3094.9398	6328.2535	0.0000	1.0000
GTRIG	21	-3093.6663	6332.6252	0.0000	1.0000
HKY	15	-3121.7935	6347.3673	0.0000	1.0000
K2PG	13	-3131.1070	6352.1570	0.0000	1.0000
K2PI	13	-3132.5494	6355.0418	0.0000	1.0000
K2PIG	14	-3131.0547	6358.9711	0.0000	1.0000
GTR	19	-3118.4567	6368.3687	0.0000	1.0000
SYMG	17	-3130.5556	6378.7290	0.0000	1.0000
SYMI	17	-3132.1967	6382.0112	0.0000	1.0000
SYMIG	18	-3130.5028	6385.5421	0.0000	1.0000
K2P	12	-3158.6927	6400.4098	0.0000	1.0000
SYM	16	-3158.0415	6426.7820	0.0000	1.0000
F81G	15	-3190.5287	6484.8378	0.0000	1.0000
F81I	15	-3192.2734	6488.3272	0.0000	1.0000
F81IG	16	-3190.4911	6491.6813	0.0000	1.0000
JC69G	12	-3219.5278	6522.0798	0.0000	1.0000
JC69I	12	-3221.4681	6525.9605	0.0000	1.0000
F81	14	-3215.0911	6527.0440	0.0000	1.0000
JC69IG	13	-3219.4812	6528.9054	0.0000	1.0000
JC69	11	-3245.3959	6566.8975	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 683-tomodel-phylip/chr8_4405.phylip = 1 - 1011;
 Partition Dummy = 1:683-tomodel-phylip/chr8_4405.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 683-tomodel-phylip/chr8_4405.phylip = 1 - 1011;
 Partition Dummy = 1:683-tomodel-phylip/chr8_4405.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 (HKYG)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr8_4405.phylip = 1 - 1011;
 Partition Dummy = 1:683-tomodel-phylip/chr8_4405.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

