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

Input data from file "683-tomodel-phylip/chr7_3352.phylip" (7 taxa, 1067 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	-4648.0392	9336.0784	0.4336	0.4336
GTRI	20	-4648.2928	9336.5856	0.3365	0.7700
GTRIG	21	-4647.7157	9337.4314	0.2204	0.9905
HKYG	16	-4656.2740	9344.5480	0.0063	0.9968
HKYIG	17	-4656.2813	9346.5626	0.0023	0.9991
HKYI	16	-4658.1681	9348.3362	0.0009	1.0000
SYMG	17	-4668.6562	9371.3125	0.0000	1.0000
SYMI	17	-4668.7216	9371.4431	0.0000	1.0000
SYMIG	18	-4668.3278	9372.6556	0.0000	1.0000
K2PG	13	-4676.4781	9378.9562	0.0000	1.0000
K2PIG	14	-4676.4797	9380.9594	0.0000	1.0000
K2PI	13	-4678.1441	9382.2883	0.0000	1.0000
GTR	19	-4672.4322	9382.8643	0.0000	1.0000
HKY	15	-4681.0345	9392.0689	0.0000	1.0000
SYM	16	-4689.7305	9411.4610	0.0000	1.0000
K2P	12	-4699.1635	9422.3271	0.0000	1.0000
F81I	15	-4777.6178	9585.2356	0.0000	1.0000
F81G	15	-4778.2967	9586.5935	0.0000	1.0000
F81IG	16	-4777.5457	9587.0913	0.0000	1.0000
JC69I	12	-4790.7312	9605.4624	0.0000	1.0000
JC69G	12	-4791.4601	9606.9203	0.0000	1.0000
JC69IG	13	-4790.6710	9607.3420	0.0000	1.0000
F81	14	-4796.7847	9621.5694	0.0000	1.0000
JC69	11	-4809.5204	9641.0408	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-4656.2740	9424.1097	0.8425	0.8425
HKYI	16	-4658.1681	9427.8979	0.1268	0.9693
HKYIG	17	-4656.2813	9431.0969	0.0256	0.9949
GTRG	20	-4648.0392	9435.5305	0.0028	0.9977
GTRI	20	-4648.2928	9436.0377	0.0022	0.9998
GTRIG	21	-4647.7157	9441.8561	0.0001	0.9999
K2PG	13	-4676.4781	9443.6001	0.0000	1.0000
K2PI	13	-4678.1441	9446.9322	0.0000	1.0000
K2PIG	14	-4676.4797	9450.5759	0.0000	1.0000
SYMG	17	-4668.6562	9455.8468	0.0000	1.0000
SYMI	17	-4668.7216	9455.9774	0.0000	1.0000
SYMIG	18	-4668.3278	9462.1626	0.0000	1.0000
HKY	15	-4681.0345	9466.6580	0.0000	1.0000
GTR	19	-4672.4322	9477.3438	0.0000	1.0000
K2P	12	-4699.1635	9481.9983	0.0000	1.0000
SYM	16	-4689.7305	9491.0227	0.0000	1.0000
F81I	15	-4777.6178	9659.8247	0.0000	1.0000
F81G	15	-4778.2967	9661.1826	0.0000	1.0000
JC69I	12	-4790.7312	9665.1336	0.0000	1.0000
JC69G	12	-4791.4601	9666.5916	0.0000	1.0000
F81IG	16	-4777.5457	9666.6530	0.0000	1.0000
JC69IG	13	-4790.6710	9671.9858	0.0000	1.0000
F81	14	-4796.7847	9691.1859	0.0000	1.0000
JC69	11	-4809.5204	9695.7394	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/chr7_3352.phylip = 1 - 1067;
 Partition Dummy = 1:683-tomodel-phylip/chr7_3352.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/chr7_3352.phylip = 1 - 1067;
 Partition Dummy = 1:683-tomodel-phylip/chr7_3352.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/chr7_3352.phylip = 1 - 1067;
 Partition Dummy = 1:683-tomodel-phylip/chr7_3352.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

