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

Input data from file "683-tomodel-phylip/chr5_3332.phylip" (7 taxa, 1165 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	-4190.0449	8420.0898	0.6831	0.6831
GTRIG	21	-4190.0448	8422.0896	0.2513	0.9344
GTRI	20	-4192.6779	8425.3557	0.0491	0.9835
HKYG	16	-4198.1444	8428.2889	0.0113	0.9948
HKYIG	17	-4198.1445	8430.2891	0.0042	0.9990
HKYI	16	-4200.5867	8433.1734	0.0010	1.0000
GTR	19	-4210.5204	8459.0408	0.0000	1.0000
HKY	15	-4218.6667	8467.3335	0.0000	1.0000
SYMG	17	-4223.2912	8480.5824	0.0000	1.0000
SYMIG	18	-4223.2943	8482.5887	0.0000	1.0000
SYMI	17	-4225.0079	8484.0158	0.0000	1.0000
K2PG	13	-4230.5979	8487.1958	0.0000	1.0000
K2PIG	14	-4230.6028	8489.2057	0.0000	1.0000
K2PI	13	-4232.3441	8490.6882	0.0000	1.0000
SYM	16	-4241.1689	8514.3378	0.0000	1.0000
K2P	12	-4249.7511	8523.5022	0.0000	1.0000
F81G	15	-4307.1142	8644.2285	0.0000	1.0000
F81IG	16	-4307.1236	8646.2471	0.0000	1.0000
F81I	15	-4308.4192	8646.8384	0.0000	1.0000
F81	14	-4324.6450	8677.2899	0.0000	1.0000
JC69G	12	-4332.6738	8689.3475	0.0000	1.0000
JC69I	12	-4333.6835	8691.3670	0.0000	1.0000
JC69IG	13	-4332.6841	8691.3682	0.0000	1.0000
JC69	11	-4349.2818	8720.5635	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-4198.1444	8509.2565	0.8937	0.8937
HKYI	16	-4200.5867	8514.1410	0.0777	0.9714
HKYIG	17	-4198.1445	8516.3172	0.0262	0.9976
GTRG	20	-4190.0449	8521.2993	0.0022	0.9998
GTRI	20	-4192.6779	8526.5652	0.0002	0.9999
GTRIG	21	-4190.0448	8528.3596	0.0001	1.0000
HKY	15	-4218.6667	8543.2406	0.0000	1.0000
K2PG	13	-4230.5979	8552.9820	0.0000	1.0000
GTR	19	-4210.5204	8555.1899	0.0000	1.0000
K2PI	13	-4232.3441	8556.4744	0.0000	1.0000
K2PIG	14	-4230.6028	8560.0523	0.0000	1.0000
SYMG	17	-4223.2912	8566.6105	0.0000	1.0000
SYMI	17	-4225.0079	8570.0439	0.0000	1.0000
SYMIG	18	-4223.2943	8573.6772	0.0000	1.0000
K2P	12	-4249.7511	8584.2279	0.0000	1.0000
SYM	16	-4241.1689	8595.3054	0.0000	1.0000
F81G	15	-4307.1142	8720.1356	0.0000	1.0000
F81I	15	-4308.4192	8722.7456	0.0000	1.0000
F81IG	16	-4307.1236	8727.2147	0.0000	1.0000
F81	14	-4324.6450	8748.1366	0.0000	1.0000
JC69G	12	-4332.6738	8750.0732	0.0000	1.0000
JC69I	12	-4333.6835	8752.0927	0.0000	1.0000
JC69IG	13	-4332.6841	8757.1544	0.0000	1.0000
JC69	11	-4349.2818	8776.2288	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 683-tomodel-phylip/chr5_3332.phylip = 1 - 1165;
 Partition Dummy = 1:683-tomodel-phylip/chr5_3332.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 683-tomodel-phylip/chr5_3332.phylip = 1 - 1165;
 Partition Dummy = 1:683-tomodel-phylip/chr5_3332.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 683-tomodel-phylip/chr5_3332.phylip = 1 - 1165;
 Partition Dummy = 1:683-tomodel-phylip/chr5_3332.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

