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

Input data from file "591-tomodel-phylip/chr1_5411.phylip" (7 taxa, 1058 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	-3501.3104	7042.6208	0.4696	0.4696
HKYG	16	-3505.9785	7043.9569	0.2407	0.7103
GTRIG	21	-3501.3105	7044.6210	0.1727	0.8830
HKYIG	17	-3505.9784	7045.9569	0.0886	0.9716
GTRI	20	-3504.7398	7049.4797	0.0152	0.9868
HKYI	16	-3508.8828	7049.7656	0.0132	1.0000
GTR	19	-3517.6740	7073.3480	0.0000	1.0000
HKY	15	-3522.5496	7075.0992	0.0000	1.0000
SYMG	17	-3583.2924	7200.5847	0.0000	1.0000
K2PG	13	-3587.4056	7200.8111	0.0000	1.0000
SYMIG	18	-3583.2924	7202.5847	0.0000	1.0000
K2PIG	14	-3587.4056	7202.8111	0.0000	1.0000
SYMI	17	-3585.4471	7204.8941	0.0000	1.0000
K2PI	13	-3589.4864	7204.9728	0.0000	1.0000
F81G	15	-3598.3308	7226.6617	0.0000	1.0000
F81IG	16	-3598.3309	7228.6617	0.0000	1.0000
F81I	15	-3599.8419	7229.6837	0.0000	1.0000
K2P	12	-3606.9045	7237.8089	0.0000	1.0000
SYM	16	-3602.9943	7237.9886	0.0000	1.0000
F81	14	-3612.7454	7253.4909	0.0000	1.0000
JC69G	12	-3664.3415	7352.6829	0.0000	1.0000
JC69IG	13	-3664.3415	7354.6830	0.0000	1.0000
JC69I	12	-3665.6501	7355.3002	0.0000	1.0000
JC69	11	-3680.4864	7382.9729	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3505.9785	7123.3831	0.9211	0.9211
HKYI	16	-3508.8828	7129.1917	0.0505	0.9716
HKYIG	17	-3505.9784	7130.3472	0.0283	0.9999
GTRG	20	-3501.3104	7141.9035	0.0001	1.0000
GTRI	20	-3504.7398	7148.7624	0.0000	1.0000
GTRIG	21	-3501.3105	7148.8678	0.0000	1.0000
HKY	15	-3522.5496	7149.5612	0.0000	1.0000
GTR	19	-3517.6740	7167.6665	0.0000	1.0000
K2PG	13	-3587.4056	7265.3449	0.0000	1.0000
K2PI	13	-3589.4864	7269.5066	0.0000	1.0000
K2PIG	14	-3587.4056	7272.3090	0.0000	1.0000
SYMG	17	-3583.2924	7284.9750	0.0000	1.0000
SYMI	17	-3585.4471	7289.2844	0.0000	1.0000
SYMIG	18	-3583.2924	7291.9392	0.0000	1.0000
K2P	12	-3606.9045	7297.3786	0.0000	1.0000
F81G	15	-3598.3308	7301.1237	0.0000	1.0000
F81I	15	-3599.8419	7304.1458	0.0000	1.0000
F81IG	16	-3598.3309	7308.0879	0.0000	1.0000
SYM	16	-3602.9943	7317.4148	0.0000	1.0000
F81	14	-3612.7454	7322.9888	0.0000	1.0000
JC69G	12	-3664.3415	7412.2526	0.0000	1.0000
JC69I	12	-3665.6501	7414.8698	0.0000	1.0000
JC69IG	13	-3664.3415	7419.2168	0.0000	1.0000
JC69	11	-3680.4864	7437.5784	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/chr1_5411.phylip = 1 - 1058;
 Partition Dummy = 1:591-tomodel-phylip/chr1_5411.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/chr1_5411.phylip = 1 - 1058;
 Partition Dummy = 1:591-tomodel-phylip/chr1_5411.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/chr1_5411.phylip = 1 - 1058;
 Partition Dummy = 1:591-tomodel-phylip/chr1_5411.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

