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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: HKYI


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2062.9952	4165.9905	0.4232	0.4232
GTRI	20	-2063.0020	4166.0039	0.4204	0.8436
GTRIG	21	-2062.9936	4167.9872	0.1559	0.9995
HKYI	16	-2074.6778	4181.3557	0.0002	0.9997
HKYG	16	-2074.6953	4181.3907	0.0002	0.9999
HKYIG	17	-2074.6946	4183.3893	0.0001	1.0000
SYMG	17	-2078.1715	4190.3429	0.0000	1.0000
SYMI	17	-2078.1823	4190.3646	0.0000	1.0000
SYMIG	18	-2078.1841	4192.3682	0.0000	1.0000
K2PI	13	-2086.0407	4198.0813	0.0000	1.0000
K2PG	13	-2086.0473	4198.0946	0.0000	1.0000
K2PIG	14	-2086.0509	4200.1018	0.0000	1.0000
GTR	19	-2091.7130	4221.4259	0.0000	1.0000
HKY	15	-2102.1971	4234.3943	0.0000	1.0000
F81I	15	-2104.8484	4239.6969	0.0000	1.0000
F81G	15	-2105.0318	4240.0636	0.0000	1.0000
F81IG	16	-2104.8802	4241.7603	0.0000	1.0000
SYM	16	-2104.9935	4241.9869	0.0000	1.0000
K2P	12	-2112.4743	4248.9487	0.0000	1.0000
JC69I	12	-2114.5975	4253.1950	0.0000	1.0000
JC69G	12	-2114.5999	4253.1999	0.0000	1.0000
JC69IG	13	-2114.5999	4255.1999	0.0000	1.0000
F81	14	-2130.5228	4289.0455	0.0000	1.0000
JC69	11	-2139.7826	4301.5652	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2074.6778	4257.8707	0.3397	0.3397
HKYG	16	-2074.6953	4257.9057	0.3338	0.6736
K2PI	13	-2086.0407	4260.2498	0.1034	0.7770
K2PG	13	-2086.0473	4260.2631	0.1027	0.8797
GTRG	20	-2062.9952	4261.6343	0.0517	0.9315
GTRI	20	-2063.0020	4261.6478	0.0514	0.9829
HKYIG	17	-2074.6946	4264.6865	0.0112	0.9941
K2PIG	14	-2086.0509	4267.0525	0.0034	0.9976
GTRIG	21	-2062.9936	4268.4132	0.0017	0.9993
SYMG	17	-2078.1715	4271.6402	0.0003	0.9996
SYMI	17	-2078.1823	4271.6619	0.0003	1.0000
SYMIG	18	-2078.1841	4278.4476	0.0000	1.0000
HKY	15	-2102.1971	4306.1272	0.0000	1.0000
K2P	12	-2112.4743	4306.3350	0.0000	1.0000
JC69I	12	-2114.5975	4310.5813	0.0000	1.0000
JC69G	12	-2114.5999	4310.5862	0.0000	1.0000
F81I	15	-2104.8484	4311.4297	0.0000	1.0000
F81G	15	-2105.0318	4311.7965	0.0000	1.0000
GTR	19	-2091.7130	4312.2876	0.0000	1.0000
JC69IG	13	-2114.5999	4317.3684	0.0000	1.0000
F81IG	16	-2104.8802	4318.2754	0.0000	1.0000
SYM	16	-2104.9935	4318.5020	0.0000	1.0000
JC69	11	-2139.7826	4354.1693	0.0000	1.0000
F81	14	-2130.5228	4355.9962	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/chr8_4340.phylip = 1 - 882;
 Partition Dummy = 1:683-tomodel-phylip/chr8_4340.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/chr8_4340.phylip = 1 - 882;
 Partition Dummy = 1:683-tomodel-phylip/chr8_4340.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 (HKYI)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr8_4340.phylip = 1 - 882;
 Partition Dummy = 1:683-tomodel-phylip/chr8_4340.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

