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

Input data from file "591-tomodel-phylip/chr7_2669.phylip" (7 taxa, 1122 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	-2660.9833	5361.9666	0.3046	0.3046
GTRI	20	-2661.3701	5362.7402	0.2069	0.5115
HKYG	16	-2665.6035	5363.2070	0.1638	0.6753
HKYI	16	-2665.7952	5363.5903	0.1353	0.8106
GTRIG	21	-2660.8582	5363.7163	0.1270	0.9376
HKYIG	17	-2665.5689	5365.1379	0.0624	1.0000
GTR	19	-2684.5821	5407.1641	0.0000	1.0000
HKY	15	-2688.8916	5407.7831	0.0000	1.0000
K2PG	13	-2715.3592	5456.7183	0.0000	1.0000
K2PI	13	-2715.4419	5456.8838	0.0000	1.0000
K2PIG	14	-2715.3692	5458.7383	0.0000	1.0000
SYMG	17	-2712.5177	5459.0353	0.0000	1.0000
SYMI	17	-2712.5986	5459.1971	0.0000	1.0000
SYMIG	18	-2712.5210	5461.0420	0.0000	1.0000
F81G	15	-2720.5325	5471.0650	0.0000	1.0000
F81I	15	-2720.6911	5471.3823	0.0000	1.0000
F81IG	16	-2720.5118	5473.0235	0.0000	1.0000
K2P	12	-2739.5293	5503.0586	0.0000	1.0000
SYM	16	-2736.7107	5505.4215	0.0000	1.0000
F81	14	-2742.6745	5513.3490	0.0000	1.0000
JC69G	12	-2765.4580	5554.9159	0.0000	1.0000
JC69I	12	-2765.5811	5555.1623	0.0000	1.0000
JC69IG	13	-2765.4502	5556.9005	0.0000	1.0000
JC69	11	-2788.1792	5598.3584	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2665.6035	5443.5728	0.5386	0.5386
HKYI	16	-2665.7952	5443.9562	0.4447	0.9833
HKYIG	17	-2665.5689	5450.5266	0.0166	0.9999
GTRG	20	-2660.9833	5462.4240	0.0000	1.0000
GTRI	20	-2661.3701	5463.1976	0.0000	1.0000
GTRIG	21	-2660.8582	5469.1965	0.0000	1.0000
HKY	15	-2688.8916	5483.1261	0.0000	1.0000
GTR	19	-2684.5821	5502.5986	0.0000	1.0000
K2PG	13	-2715.3592	5522.0156	0.0000	1.0000
K2PI	13	-2715.4419	5522.1810	0.0000	1.0000
K2PIG	14	-2715.3692	5529.0585	0.0000	1.0000
SYMG	17	-2712.5177	5544.4241	0.0000	1.0000
SYMI	17	-2712.5986	5544.5859	0.0000	1.0000
F81G	15	-2720.5325	5546.4080	0.0000	1.0000
F81I	15	-2720.6911	5546.7253	0.0000	1.0000
SYMIG	18	-2712.5210	5551.4536	0.0000	1.0000
F81IG	16	-2720.5118	5553.3894	0.0000	1.0000
K2P	12	-2739.5293	5563.3330	0.0000	1.0000
F81	14	-2742.6745	5583.6691	0.0000	1.0000
SYM	16	-2736.7107	5585.7873	0.0000	1.0000
JC69G	12	-2765.4580	5615.1904	0.0000	1.0000
JC69I	12	-2765.5811	5615.4367	0.0000	1.0000
JC69IG	13	-2765.4502	5622.1977	0.0000	1.0000
JC69	11	-2788.1792	5653.6099	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/chr7_2669.phylip = 1 - 1122;
 Partition Dummy = 1:591-tomodel-phylip/chr7_2669.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/chr7_2669.phylip = 1 - 1122;
 Partition Dummy = 1:591-tomodel-phylip/chr7_2669.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 591-tomodel-phylip/chr7_2669.phylip = 1 - 1122;
 Partition Dummy = 1:591-tomodel-phylip/chr7_2669.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

