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

Input data from file "591-tomodel-phylip/chr5_8793.phylip" (7 taxa, 1148 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	-3729.9836	7491.9671	0.6567	0.6567
HKYIG	17	-3730.0061	7494.0123	0.2362	0.8929
HKYI	16	-3732.0498	7496.0995	0.0832	0.9761
GTRG	20	-3729.6967	7499.3934	0.0160	0.9922
GTRIG	21	-3729.7083	7501.4167	0.0058	0.9980
GTRI	20	-3731.7692	7503.5384	0.0020	1.0000
HKY	15	-3763.0753	7556.1507	0.0000	1.0000
GTR	19	-3762.6352	7563.2703	0.0000	1.0000
F81G	15	-3793.5948	7617.1896	0.0000	1.0000
F81IG	16	-3793.6056	7619.2112	0.0000	1.0000
F81I	15	-3795.5274	7621.0548	0.0000	1.0000
SYMG	17	-3814.9455	7663.8910	0.0000	1.0000
SYMIG	18	-3814.9284	7665.8567	0.0000	1.0000
SYMI	17	-3817.7093	7669.4186	0.0000	1.0000
K2PG	13	-3824.0672	7674.1345	0.0000	1.0000
K2PIG	14	-3824.0594	7676.1188	0.0000	1.0000
F81	14	-3824.3298	7676.6595	0.0000	1.0000
K2PI	13	-3826.3906	7678.7812	0.0000	1.0000
SYM	16	-3849.6302	7731.2604	0.0000	1.0000
K2P	12	-3857.5392	7739.0784	0.0000	1.0000
JC69G	12	-3874.8165	7773.6329	0.0000	1.0000
JC69IG	13	-3874.8155	7775.6310	0.0000	1.0000
JC69I	12	-3876.9345	7777.8690	0.0000	1.0000
JC69	11	-3905.8998	7833.7995	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3729.9836	7572.6995	0.8654	0.8654
HKYI	16	-3732.0498	7576.8320	0.1096	0.9750
HKYIG	17	-3730.0061	7579.7905	0.0250	1.0000
GTRG	20	-3729.6967	7600.3089	0.0000	1.0000
GTRI	20	-3731.7692	7604.4539	0.0000	1.0000
GTRIG	21	-3729.7083	7607.3780	0.0000	1.0000
HKY	15	-3763.0753	7631.8373	0.0000	1.0000
GTR	19	-3762.6352	7659.1401	0.0000	1.0000
F81G	15	-3793.5948	7692.8762	0.0000	1.0000
F81I	15	-3795.5274	7696.7415	0.0000	1.0000
F81IG	16	-3793.6056	7699.9436	0.0000	1.0000
K2PG	13	-3824.0672	7739.7296	0.0000	1.0000
K2PI	13	-3826.3906	7744.3763	0.0000	1.0000
K2PIG	14	-3824.0594	7746.7597	0.0000	1.0000
F81	14	-3824.3298	7747.3004	0.0000	1.0000
SYMG	17	-3814.9455	7749.6692	0.0000	1.0000
SYMI	17	-3817.7093	7755.1968	0.0000	1.0000
SYMIG	18	-3814.9284	7756.6807	0.0000	1.0000
K2P	12	-3857.5392	7799.6277	0.0000	1.0000
SYM	16	-3849.6302	7811.9928	0.0000	1.0000
JC69G	12	-3874.8165	7834.1823	0.0000	1.0000
JC69I	12	-3876.9345	7838.4184	0.0000	1.0000
JC69IG	13	-3874.8155	7841.2261	0.0000	1.0000
JC69	11	-3905.8998	7889.3030	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 591-tomodel-phylip/chr5_8793.phylip = 1 - 1148;
 Partition Dummy = 1:591-tomodel-phylip/chr5_8793.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 591-tomodel-phylip/chr5_8793.phylip = 1 - 1148;
 Partition Dummy = 1:591-tomodel-phylip/chr5_8793.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 591-tomodel-phylip/chr5_8793.phylip = 1 - 1148;
 Partition Dummy = 1:591-tomodel-phylip/chr5_8793.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

