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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKYG


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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2755.9326	5543.8652	0.4418	0.4418
HKYI	16	-2756.1952	5544.3904	0.3398	0.7815
HKYIG	17	-2755.9024	5545.8047	0.1675	0.9490
GTRG	20	-2754.8538	5549.7076	0.0238	0.9728
GTRI	20	-2755.1260	5550.2519	0.0181	0.9910
GTRIG	21	-2754.8235	5551.6470	0.0090	1.0000
HKY	15	-2786.3110	5602.6221	0.0000	1.0000
GTR	19	-2785.2421	5608.4842	0.0000	1.0000
F81G	15	-2812.6600	5655.3201	0.0000	1.0000
F81I	15	-2812.7149	5655.4298	0.0000	1.0000
F81IG	16	-2812.6657	5657.3315	0.0000	1.0000
SYMG	17	-2812.8117	5659.6233	0.0000	1.0000
K2PG	13	-2817.1917	5660.3835	0.0000	1.0000
SYMI	17	-2813.3496	5660.6992	0.0000	1.0000
SYMIG	18	-2812.6457	5661.2915	0.0000	1.0000
K2PIG	14	-2817.1181	5662.2363	0.0000	1.0000
K2PI	13	-2818.4551	5662.9103	0.0000	1.0000
F81	14	-2841.9051	5711.8103	0.0000	1.0000
K2P	12	-2848.8051	5721.6103	0.0000	1.0000
SYM	16	-2844.8411	5721.6821	0.0000	1.0000
JC69G	12	-2867.3094	5758.6189	0.0000	1.0000
JC69I	12	-2867.4570	5758.9139	0.0000	1.0000
JC69IG	13	-2866.7943	5759.5886	0.0000	1.0000
JC69	11	-2895.5859	5813.1719	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2755.9326	5623.6058	0.5554	0.5554
HKYI	16	-2756.1952	5624.1310	0.4271	0.9826
HKYIG	17	-2755.9024	5630.5292	0.0174	1.0000
GTRG	20	-2754.8538	5649.3834	0.0000	1.0000
GTRI	20	-2755.1260	5649.9277	0.0000	1.0000
GTRIG	21	-2754.8235	5656.3066	0.0000	1.0000
HKY	15	-2786.3110	5677.3789	0.0000	1.0000
GTR	19	-2785.2421	5703.1762	0.0000	1.0000
K2PG	13	-2817.1917	5725.1727	0.0000	1.0000
K2PI	13	-2818.4551	5727.6996	0.0000	1.0000
F81G	15	-2812.6600	5730.0769	0.0000	1.0000
F81I	15	-2812.7149	5730.1866	0.0000	1.0000
K2PIG	14	-2817.1181	5732.0093	0.0000	1.0000
F81IG	16	-2812.6657	5737.0721	0.0000	1.0000
SYMG	17	-2812.8117	5744.3478	0.0000	1.0000
SYMI	17	-2813.3496	5745.4236	0.0000	1.0000
SYMIG	18	-2812.6457	5750.9997	0.0000	1.0000
K2P	12	-2848.8051	5781.4158	0.0000	1.0000
F81	14	-2841.9051	5781.5833	0.0000	1.0000
SYM	16	-2844.8411	5801.4228	0.0000	1.0000
JC69G	12	-2867.3094	5818.4244	0.0000	1.0000
JC69I	12	-2867.4570	5818.7194	0.0000	1.0000
JC69IG	13	-2866.7943	5824.3778	0.0000	1.0000
JC69	11	-2895.5859	5867.9936	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (HKYG)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr7_4885.phylip = 1 - 1079;
 Partition Dummy = 1:683-tomodel-phylip/chr7_4885.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;


[Mrbayes block for the best AICc model (HKYG)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr7_4885.phylip = 1 - 1079;
 Partition Dummy = 1:683-tomodel-phylip/chr7_4885.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;


[Mrbayes block for the best BIC model (HKYG)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr7_4885.phylip = 1 - 1079;
 Partition Dummy = 1:683-tomodel-phylip/chr7_4885.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

