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

Input data from file "683-tomodel-phylip/chr5_10353.phylip" (7 taxa, 1182 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	-3539.6545	7119.3090	0.4847	0.4847
GTRIG	21	-3539.4764	7120.9528	0.2131	0.6978
HKYG	16	-3544.6187	7121.2374	0.1848	0.8827
HKYIG	17	-3544.4366	7122.8733	0.0816	0.9642
GTRI	20	-3542.6100	7125.2201	0.0252	0.9895
HKYI	16	-3547.4826	7126.9651	0.0105	1.0000
GTR	19	-3577.1241	7192.2482	0.0000	1.0000
HKY	15	-3581.8764	7193.7527	0.0000	1.0000
K2PG	13	-3602.0276	7230.0552	0.0000	1.0000
SYMG	17	-3598.6407	7231.2814	0.0000	1.0000
K2PIG	14	-3601.8126	7231.6253	0.0000	1.0000
SYMIG	18	-3598.4303	7232.8606	0.0000	1.0000
K2PI	13	-3604.8210	7235.6421	0.0000	1.0000
SYMI	17	-3601.4694	7236.9388	0.0000	1.0000
K2P	12	-3641.5010	7307.0021	0.0000	1.0000
SYM	16	-3639.1034	7310.2068	0.0000	1.0000
F81G	15	-3643.8524	7317.7048	0.0000	1.0000
F81IG	16	-3643.7653	7319.5305	0.0000	1.0000
F81I	15	-3646.0453	7322.0905	0.0000	1.0000
F81	14	-3678.2458	7384.4916	0.0000	1.0000
JC69G	12	-3691.4268	7406.8536	0.0000	1.0000
JC69IG	13	-3691.3205	7408.6411	0.0000	1.0000
JC69I	12	-3693.7101	7411.4201	0.0000	1.0000
JC69	11	-3727.3091	7476.6181	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3544.6187	7202.4368	0.9157	0.9157
HKYI	16	-3547.4826	7208.1645	0.0522	0.9679
HKYIG	17	-3544.4366	7209.1476	0.0320	0.9999
GTRG	20	-3539.6545	7220.8083	0.0001	1.0000
GTRI	20	-3542.6100	7226.7194	0.0000	1.0000
GTRIG	21	-3539.4764	7227.5270	0.0000	1.0000
HKY	15	-3581.8764	7269.8772	0.0000	1.0000
GTR	19	-3577.1241	7288.6725	0.0000	1.0000
K2PG	13	-3602.0276	7296.0298	0.0000	1.0000
K2PI	13	-3604.8210	7301.6166	0.0000	1.0000
K2PIG	14	-3601.8126	7302.6747	0.0000	1.0000
SYMG	17	-3598.6407	7317.5558	0.0000	1.0000
SYMI	17	-3601.4694	7323.2132	0.0000	1.0000
SYMIG	18	-3598.4303	7324.2099	0.0000	1.0000
K2P	12	-3641.5010	7367.9017	0.0000	1.0000
SYM	16	-3639.1034	7391.4062	0.0000	1.0000
F81G	15	-3643.8524	7393.8292	0.0000	1.0000
F81I	15	-3646.0453	7398.2149	0.0000	1.0000
F81IG	16	-3643.7653	7400.7299	0.0000	1.0000
F81	14	-3678.2458	7455.5411	0.0000	1.0000
JC69G	12	-3691.4268	7467.7531	0.0000	1.0000
JC69I	12	-3693.7101	7472.3197	0.0000	1.0000
JC69IG	13	-3691.3205	7474.6156	0.0000	1.0000
JC69	11	-3727.3091	7532.4427	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/chr5_10353.phylip = 1 - 1182;
 Partition Dummy = 1:683-tomodel-phylip/chr5_10353.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/chr5_10353.phylip = 1 - 1182;
 Partition Dummy = 1:683-tomodel-phylip/chr5_10353.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/chr5_10353.phylip = 1 - 1182;
 Partition Dummy = 1:683-tomodel-phylip/chr5_10353.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

