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

Input data from file "591-tomodel-phylip/chr3_17769.phylip" (7 taxa, 1076 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	-2704.9314	5449.8627	0.5288	0.5288
GTRIG	21	-2704.4977	5450.9953	0.3002	0.8289
GTRI	20	-2706.0686	5452.1373	0.1696	0.9985
HKYG	16	-2715.5259	5463.0518	0.0007	0.9992
HKYIG	17	-2714.8920	5463.7839	0.0005	0.9997
HKYI	16	-2716.5795	5465.1589	0.0003	1.0000
SYMG	17	-2742.4560	5518.9120	0.0000	1.0000
SYMIG	18	-2742.1688	5520.3376	0.0000	1.0000
SYMI	17	-2743.3029	5520.6057	0.0000	1.0000
GTR	19	-2742.0936	5522.1872	0.0000	1.0000
HKY	15	-2750.8793	5531.7586	0.0000	1.0000
K2PG	13	-2754.9170	5535.8340	0.0000	1.0000
K2PIG	14	-2754.3992	5536.7983	0.0000	1.0000
K2PI	13	-2755.8455	5537.6909	0.0000	1.0000
SYM	16	-2777.9505	5587.9011	0.0000	1.0000
K2P	12	-2788.8577	5601.7155	0.0000	1.0000
F81G	15	-2790.0975	5610.1951	0.0000	1.0000
F81IG	16	-2789.6653	5611.3306	0.0000	1.0000
F81I	15	-2791.0362	5612.0724	0.0000	1.0000
F81	14	-2820.7069	5669.4137	0.0000	1.0000
JC69G	12	-2825.6970	5675.3940	0.0000	1.0000
JC69IG	13	-2825.4124	5676.8248	0.0000	1.0000
JC69I	12	-2826.4767	5676.9535	0.0000	1.0000
JC69	11	-2855.3384	5732.6768	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2715.5259	5542.7479	0.6881	0.6881
HKYI	16	-2716.5795	5544.8550	0.2399	0.9280
HKYIG	17	-2714.8920	5548.4610	0.0395	0.9676
GTRG	20	-2704.9314	5549.4828	0.0237	0.9913
GTRI	20	-2706.0686	5551.7574	0.0076	0.9989
GTRIG	21	-2704.4977	5555.5965	0.0011	1.0000
K2PG	13	-2754.9170	5600.5871	0.0000	1.0000
K2PI	13	-2755.8455	5602.4440	0.0000	1.0000
SYMG	17	-2742.4560	5603.5891	0.0000	1.0000
SYMI	17	-2743.3029	5605.2828	0.0000	1.0000
HKY	15	-2750.8793	5606.4737	0.0000	1.0000
K2PIG	14	-2754.3992	5606.5324	0.0000	1.0000
SYMIG	18	-2742.1688	5609.9957	0.0000	1.0000
GTR	19	-2742.0936	5616.8263	0.0000	1.0000
K2P	12	-2788.8577	5661.4875	0.0000	1.0000
SYM	16	-2777.9505	5667.5972	0.0000	1.0000
F81G	15	-2790.0975	5684.9102	0.0000	1.0000
F81I	15	-2791.0362	5686.7875	0.0000	1.0000
F81IG	16	-2789.6653	5691.0267	0.0000	1.0000
JC69G	12	-2825.6970	5735.1661	0.0000	1.0000
JC69I	12	-2826.4767	5736.7255	0.0000	1.0000
F81	14	-2820.7069	5739.1478	0.0000	1.0000
JC69IG	13	-2825.4124	5741.5779	0.0000	1.0000
JC69	11	-2855.3384	5787.4678	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/chr3_17769.phylip = 1 - 1076;
 Partition Dummy = 1:591-tomodel-phylip/chr3_17769.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/chr3_17769.phylip = 1 - 1076;
 Partition Dummy = 1:591-tomodel-phylip/chr3_17769.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/chr3_17769.phylip = 1 - 1076;
 Partition Dummy = 1:591-tomodel-phylip/chr3_17769.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

