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

Input data from file "591-tomodel-phylip/chr10_4885.phylip" (7 taxa, 737 characters)

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


Note: number of parameters compared to sample size
is low (Nchar/Nparams < 40). Consider the use of AICc.


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

Model	df	lnL		AICc		wAICc	cwAICc
GTRG	20	-2624.8336	5290.8403	0.6904	0.6904
GTRIG	21	-2624.8209	5292.9341	0.2423	0.9327
HKYG	16	-2631.9985	5296.7525	0.0359	0.9686
GTRI	20	-2628.4809	5298.1349	0.0180	0.9866
HKYIG	17	-2631.9775	5298.8063	0.0129	0.9995
HKYI	16	-2636.2323	5305.2202	0.0005	1.0000
SYMG	17	-2649.3011	5333.4534	0.0000	1.0000
SYMIG	18	-2649.2873	5335.5272	0.0000	1.0000
SYMI	17	-2651.7522	5338.3556	0.0000	1.0000
K2PG	13	-2657.1072	5340.7178	0.0000	1.0000
K2PIG	14	-2657.0752	5342.7322	0.0000	1.0000
GTR	19	-2653.4555	5345.9711	0.0000	1.0000
K2PI	13	-2660.5210	5347.5455	0.0000	1.0000
HKY	15	-2661.9393	5354.5443	0.0000	1.0000
SYM	16	-2680.2566	5393.2688	0.0000	1.0000
K2P	12	-2688.0653	5400.5615	0.0000	1.0000
F81G	15	-2703.4302	5437.5260	0.0000	1.0000
F81IG	16	-2703.4280	5439.6116	0.0000	1.0000
F81I	15	-2706.1325	5442.9307	0.0000	1.0000
JC69G	12	-2724.0183	5472.4675	0.0000	1.0000
JC69IG	13	-2724.0170	5474.5375	0.0000	1.0000
JC69I	12	-2726.4061	5477.2432	0.0000	1.0000
F81	14	-2730.4904	5489.5625	0.0000	1.0000
JC69	11	-2751.9432	5526.2505	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2624.8336	5289.6671	0.6862	0.6862
GTRIG	21	-2624.8209	5291.6418	0.2557	0.9418
HKYG	16	-2631.9985	5295.9970	0.0290	0.9708
GTRI	20	-2628.4809	5296.9618	0.0179	0.9887
HKYIG	17	-2631.9775	5297.9551	0.0109	0.9996
HKYI	16	-2636.2323	5304.4646	0.0004	1.0000
SYMG	17	-2649.3011	5332.6023	0.0000	1.0000
SYMIG	18	-2649.2873	5334.5746	0.0000	1.0000
SYMI	17	-2651.7522	5337.5045	0.0000	1.0000
K2PG	13	-2657.1072	5340.2144	0.0000	1.0000
K2PIG	14	-2657.0752	5342.1504	0.0000	1.0000
GTR	19	-2653.4555	5344.9111	0.0000	1.0000
K2PI	13	-2660.5210	5347.0420	0.0000	1.0000
HKY	15	-2661.9393	5353.8785	0.0000	1.0000
SYM	16	-2680.2566	5392.5133	0.0000	1.0000
K2P	12	-2688.0653	5400.1305	0.0000	1.0000
F81G	15	-2703.4302	5436.8603	0.0000	1.0000
F81IG	16	-2703.4280	5438.8560	0.0000	1.0000
F81I	15	-2706.1325	5442.2650	0.0000	1.0000
JC69G	12	-2724.0183	5472.0366	0.0000	1.0000
JC69IG	13	-2724.0170	5474.0341	0.0000	1.0000
JC69I	12	-2726.4061	5476.8122	0.0000	1.0000
F81	14	-2730.4904	5488.9808	0.0000	1.0000
JC69	11	-2751.9432	5525.8864	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2631.9985	5369.6384	0.9482	0.9482
HKYIG	17	-2631.9775	5376.1991	0.0357	0.9839
HKYI	16	-2636.2323	5378.1060	0.0137	0.9976
GTRG	20	-2624.8336	5381.7189	0.0023	0.9999
GTRIG	21	-2624.8209	5388.2961	0.0001	0.9999
GTRI	20	-2628.4809	5389.0135	0.0001	1.0000
K2PG	13	-2657.1072	5400.0480	0.0000	1.0000
K2PIG	14	-2657.0752	5406.5867	0.0000	1.0000
K2PI	13	-2660.5210	5406.8757	0.0000	1.0000
SYMG	17	-2649.3011	5410.8463	0.0000	1.0000
SYMI	17	-2651.7522	5415.7485	0.0000	1.0000
SYMIG	18	-2649.2873	5417.4212	0.0000	1.0000
HKY	15	-2661.9393	5422.9173	0.0000	1.0000
GTR	19	-2653.4555	5432.3602	0.0000	1.0000
K2P	12	-2688.0653	5455.3616	0.0000	1.0000
SYM	16	-2680.2566	5466.1547	0.0000	1.0000
F81G	15	-2703.4302	5505.8991	0.0000	1.0000
F81I	15	-2706.1325	5511.3038	0.0000	1.0000
F81IG	16	-2703.4280	5512.4974	0.0000	1.0000
JC69G	12	-2724.0183	5527.2676	0.0000	1.0000
JC69I	12	-2726.4061	5532.0433	0.0000	1.0000
JC69IG	13	-2724.0170	5533.8677	0.0000	1.0000
F81	14	-2730.4904	5553.4170	0.0000	1.0000
JC69	11	-2751.9432	5576.5148	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/chr10_4885.phylip = 1 - 737;
 Partition Dummy = 1:591-tomodel-phylip/chr10_4885.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/chr10_4885.phylip = 1 - 737;
 Partition Dummy = 1:591-tomodel-phylip/chr10_4885.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/chr10_4885.phylip = 1 - 737;
 Partition Dummy = 1:591-tomodel-phylip/chr10_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

