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

Input data from file "591-tomodel-phylip/chr5_10794.phylip" (7 taxa, 853 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	-2274.2511	4580.5023	0.5038	0.5038
HKYI	16	-2274.7895	4581.5791	0.2940	0.7978
HKYIG	17	-2274.2501	4582.5002	0.1855	0.9833
GTRG	20	-2274.2736	4588.5471	0.0090	0.9923
GTRI	20	-2275.0085	4590.0169	0.0043	0.9967
GTRIG	21	-2274.2666	4590.5332	0.0033	1.0000
HKY	15	-2288.5767	4607.1534	0.0000	1.0000
GTR	19	-2286.5772	4611.1544	0.0000	1.0000
K2PG	13	-2300.7296	4627.4593	0.0000	1.0000
K2PI	13	-2301.2681	4628.5362	0.0000	1.0000
K2PIG	14	-2300.7277	4629.4554	0.0000	1.0000
SYMIG	18	-2298.3776	4632.7551	0.0000	1.0000
SYMG	17	-2300.4271	4634.8542	0.0000	1.0000
SYMI	17	-2300.9962	4635.9924	0.0000	1.0000
K2P	12	-2315.0866	4654.1731	0.0000	1.0000
SYM	16	-2312.5772	4657.1543	0.0000	1.0000
F81G	15	-2331.2030	4692.4059	0.0000	1.0000
F81IG	16	-2331.2005	4694.4010	0.0000	1.0000
F81I	15	-2333.8933	4697.7867	0.0000	1.0000
F81	14	-2341.6884	4711.3768	0.0000	1.0000
JC69G	12	-2353.8873	4731.7747	0.0000	1.0000
JC69I	12	-2354.4114	4732.8229	0.0000	1.0000
JC69IG	13	-2353.8849	4733.7698	0.0000	1.0000
JC69	11	-2364.3609	4750.7218	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2274.2511	4656.4824	0.6181	0.6181
HKYI	16	-2274.7895	4657.5592	0.3607	0.9788
HKYIG	17	-2274.2501	4663.2291	0.0212	1.0000
HKY	15	-2288.5767	4678.3848	0.0000	1.0000
GTRG	20	-2274.2736	4683.5223	0.0000	1.0000
GTRI	20	-2275.0085	4684.9921	0.0000	1.0000
K2PG	13	-2300.7296	4689.1931	0.0000	1.0000
GTRIG	21	-2274.2666	4690.2572	0.0000	1.0000
K2PI	13	-2301.2681	4690.2700	0.0000	1.0000
K2PIG	14	-2300.7277	4695.9381	0.0000	1.0000
GTR	19	-2286.5772	4701.3808	0.0000	1.0000
K2P	12	-2315.0866	4711.1582	0.0000	1.0000
SYMG	17	-2300.4271	4715.5831	0.0000	1.0000
SYMI	17	-2300.9962	4716.7213	0.0000	1.0000
SYMIG	18	-2298.3776	4718.2328	0.0000	1.0000
SYM	16	-2312.5772	4733.1345	0.0000	1.0000
F81G	15	-2331.2030	4763.6373	0.0000	1.0000
F81I	15	-2333.8933	4769.0181	0.0000	1.0000
F81IG	16	-2331.2005	4770.3811	0.0000	1.0000
F81	14	-2341.6884	4777.8594	0.0000	1.0000
JC69G	12	-2353.8873	4788.7598	0.0000	1.0000
JC69I	12	-2354.4114	4789.8080	0.0000	1.0000
JC69IG	13	-2353.8849	4795.5036	0.0000	1.0000
JC69	11	-2364.3609	4802.9582	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_10794.phylip = 1 - 853;
 Partition Dummy = 1:591-tomodel-phylip/chr5_10794.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_10794.phylip = 1 - 853;
 Partition Dummy = 1:591-tomodel-phylip/chr5_10794.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_10794.phylip = 1 - 853;
 Partition Dummy = 1:591-tomodel-phylip/chr5_10794.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

