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

Input data from file "683-tomodel-phylip/chr7_10735.phylip" (7 taxa, 1139 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	-3710.0787	7460.1574	0.3806	0.3806
GTRI	20	-3710.1446	7460.2893	0.3563	0.7369
GTRIG	21	-3709.9387	7461.8774	0.1611	0.8980
HKYI	16	-3716.1756	7464.3513	0.0468	0.9447
HKYG	16	-3716.4108	7464.8216	0.0370	0.9817
HKYIG	17	-3716.1129	7466.2257	0.0183	1.0000
GTR	19	-3727.5088	7493.0176	0.0000	1.0000
HKY	15	-3733.7933	7497.5866	0.0000	1.0000
K2PI	13	-3772.3834	7570.7668	0.0000	1.0000
K2PG	13	-3772.5017	7571.0034	0.0000	1.0000
K2PIG	14	-3772.2838	7572.5675	0.0000	1.0000
SYMG	17	-3770.4871	7574.9742	0.0000	1.0000
SYMI	17	-3770.5355	7575.0710	0.0000	1.0000
SYMIG	18	-3770.3528	7576.7055	0.0000	1.0000
K2P	12	-3791.6588	7607.3177	0.0000	1.0000
SYM	16	-3789.6186	7611.2372	0.0000	1.0000
F81I	15	-3814.0538	7658.1076	0.0000	1.0000
F81G	15	-3814.7265	7659.4529	0.0000	1.0000
F81IG	16	-3814.0565	7660.1131	0.0000	1.0000
F81	14	-3830.0016	7688.0031	0.0000	1.0000
JC69I	12	-3858.7138	7741.4276	0.0000	1.0000
JC69G	12	-3859.1783	7742.3566	0.0000	1.0000
JC69IG	13	-3858.7025	7743.4050	0.0000	1.0000
JC69	11	-3875.7230	7773.4461	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-3716.1756	7544.9578	0.5486	0.5486
HKYG	16	-3716.4108	7545.4281	0.4337	0.9823
HKYIG	17	-3716.1129	7551.8701	0.0173	0.9996
GTRG	20	-3710.0787	7560.9156	0.0002	0.9998
GTRI	20	-3710.1446	7561.0474	0.0002	1.0000
GTRIG	21	-3709.9387	7567.6734	0.0000	1.0000
HKY	15	-3733.7933	7573.1552	0.0000	1.0000
GTR	19	-3727.5088	7588.7378	0.0000	1.0000
K2PI	13	-3772.3834	7636.2596	0.0000	1.0000
K2PG	13	-3772.5017	7636.4962	0.0000	1.0000
K2PIG	14	-3772.2838	7643.0982	0.0000	1.0000
SYMG	17	-3770.4871	7660.6186	0.0000	1.0000
SYMI	17	-3770.5355	7660.7154	0.0000	1.0000
SYMIG	18	-3770.3528	7667.3878	0.0000	1.0000
K2P	12	-3791.6588	7667.7726	0.0000	1.0000
SYM	16	-3789.6186	7691.8437	0.0000	1.0000
F81I	15	-3814.0538	7733.6761	0.0000	1.0000
F81G	15	-3814.7265	7735.0215	0.0000	1.0000
F81IG	16	-3814.0565	7740.7196	0.0000	1.0000
F81	14	-3830.0016	7758.5338	0.0000	1.0000
JC69I	12	-3858.7138	7801.8825	0.0000	1.0000
JC69G	12	-3859.1783	7802.8115	0.0000	1.0000
JC69IG	13	-3858.7025	7808.8978	0.0000	1.0000
JC69	11	-3875.7230	7828.8630	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_10735.phylip = 1 - 1139;
 Partition Dummy = 1:683-tomodel-phylip/chr7_10735.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_10735.phylip = 1 - 1139;
 Partition Dummy = 1:683-tomodel-phylip/chr7_10735.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/chr7_10735.phylip = 1 - 1139;
 Partition Dummy = 1:683-tomodel-phylip/chr7_10735.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

