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

Input data from file "591-tomodel-phylip/chr2_8146.phylip" (7 taxa, 920 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	-3245.2504	6530.5008	0.4219	0.4219
GTRI	20	-3245.3585	6530.7170	0.3787	0.8006
GTRIG	21	-3245.0158	6532.0316	0.1962	0.9968
HKYG	16	-3255.0349	6542.0697	0.0013	0.9981
HKYI	16	-3255.0487	6542.0973	0.0013	0.9994
HKYIG	17	-3254.7942	6543.5884	0.0006	1.0000
GTR	19	-3257.6604	6553.3209	0.0000	1.0000
HKY	15	-3267.1059	6564.2119	0.0000	1.0000
SYMG	17	-3333.4468	6700.8936	0.0000	1.0000
K2PG	13	-3337.7798	6701.5595	0.0000	1.0000
SYMI	17	-3333.8051	6701.6103	0.0000	1.0000
K2PI	13	-3337.9020	6701.8041	0.0000	1.0000
SYMIG	18	-3333.3971	6702.7943	0.0000	1.0000
K2PIG	14	-3337.6378	6703.2755	0.0000	1.0000
K2P	12	-3351.2070	6726.4140	0.0000	1.0000
SYM	16	-3347.2545	6726.5091	0.0000	1.0000
F81G	15	-3364.8512	6759.7023	0.0000	1.0000
F81I	15	-3365.3449	6760.6898	0.0000	1.0000
F81IG	16	-3364.8419	6761.6839	0.0000	1.0000
F81	14	-3375.0327	6778.0653	0.0000	1.0000
JC69G	12	-3426.2332	6876.4664	0.0000	1.0000
JC69I	12	-3426.9001	6877.8002	0.0000	1.0000
JC69IG	13	-3426.2358	6878.4715	0.0000	1.0000
JC69	11	-3437.4945	6896.9890	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3255.0349	6619.2597	0.4833	0.4833
HKYI	16	-3255.0487	6619.2873	0.4767	0.9600
HKYIG	17	-3254.7942	6625.6028	0.0203	0.9803
GTRG	20	-3245.2504	6626.9883	0.0101	0.9904
GTRI	20	-3245.3585	6627.2045	0.0091	0.9995
GTRIG	21	-3245.0158	6633.3434	0.0004	0.9999
HKY	15	-3267.1059	6636.5775	0.0001	1.0000
GTR	19	-3257.6604	6644.9840	0.0000	1.0000
K2PG	13	-3337.7798	6764.2764	0.0000	1.0000
K2PI	13	-3337.9020	6764.5209	0.0000	1.0000
K2PIG	14	-3337.6378	6770.8168	0.0000	1.0000
SYMG	17	-3333.4468	6782.9080	0.0000	1.0000
SYMI	17	-3333.8051	6783.6246	0.0000	1.0000
K2P	12	-3351.2070	6784.3065	0.0000	1.0000
SYMIG	18	-3333.3971	6789.6330	0.0000	1.0000
SYM	16	-3347.2545	6803.6991	0.0000	1.0000
F81G	15	-3364.8512	6832.0679	0.0000	1.0000
F81I	15	-3365.3449	6833.0554	0.0000	1.0000
F81IG	16	-3364.8419	6838.8739	0.0000	1.0000
F81	14	-3375.0327	6845.6065	0.0000	1.0000
JC69G	12	-3426.2332	6934.3589	0.0000	1.0000
JC69I	12	-3426.9001	6935.6927	0.0000	1.0000
JC69IG	13	-3426.2358	6941.1884	0.0000	1.0000
JC69	11	-3437.4945	6950.0571	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/chr2_8146.phylip = 1 - 920;
 Partition Dummy = 1:591-tomodel-phylip/chr2_8146.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/chr2_8146.phylip = 1 - 920;
 Partition Dummy = 1:591-tomodel-phylip/chr2_8146.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/chr2_8146.phylip = 1 - 920;
 Partition Dummy = 1:591-tomodel-phylip/chr2_8146.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

