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

Input data from file "683-tomodel-phylip/chrZ_7476.phylip" (7 taxa, 737 characters)

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKYG


Note: number of parameters compared to sample size
is low (Nchar/Nparams < 40). Consider the use of AICc.


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

Model	df	lnL		AICc		wAICc	cwAICc
HKYG	16	-2709.2364	5451.2284	0.4445	0.4445
HKYI	16	-2709.5140	5451.7836	0.3368	0.7813
HKYIG	17	-2709.0066	5452.8644	0.1962	0.9774
GTRG	20	-2708.7930	5458.7592	0.0103	0.9877
GTRI	20	-2709.0637	5459.3005	0.0079	0.9956
GTRIG	21	-2708.5766	5460.4455	0.0044	1.0000
K2PG	13	-2725.6818	5477.8670	0.0000	1.0000
K2PI	13	-2726.2767	5479.0569	0.0000	1.0000
K2PIG	14	-2725.6049	5479.7915	0.0000	1.0000
SYMG	17	-2724.3122	5483.4756	0.0000	1.0000
SYMI	17	-2724.8186	5484.4883	0.0000	1.0000
SYMIG	18	-2724.2090	5485.3707	0.0000	1.0000
HKY	15	-2731.5246	5493.7150	0.0000	1.0000
GTR	19	-2731.0162	5501.0923	0.0000	1.0000
K2P	12	-2748.3162	5521.0634	0.0000	1.0000
SYM	16	-2746.4693	5525.6941	0.0000	1.0000
F81G	15	-2799.7803	5630.2264	0.0000	1.0000
F81IG	16	-2799.7804	5632.3163	0.0000	1.0000
F81I	15	-2801.6411	5633.9480	0.0000	1.0000
JC69G	12	-2815.1754	5654.7817	0.0000	1.0000
JC69IG	13	-2815.1754	5656.8543	0.0000	1.0000
JC69I	12	-2816.9395	5658.3100	0.0000	1.0000
F81	14	-2821.1737	5670.9290	0.0000	1.0000
JC69	11	-2836.1606	5694.6854	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2709.2364	5450.4728	0.4379	0.4379
HKYI	16	-2709.5140	5451.0281	0.3317	0.7696
HKYIG	17	-2709.0066	5452.0133	0.2027	0.9723
GTRG	20	-2708.7930	5457.5860	0.0125	0.9848
GTRI	20	-2709.0637	5458.1273	0.0095	0.9943
GTRIG	21	-2708.5766	5459.1532	0.0057	1.0000
K2PG	13	-2725.6818	5477.3635	0.0000	1.0000
K2PI	13	-2726.2767	5478.5534	0.0000	1.0000
K2PIG	14	-2725.6049	5479.2098	0.0000	1.0000
SYMG	17	-2724.3122	5482.6245	0.0000	1.0000
SYMI	17	-2724.8186	5483.6372	0.0000	1.0000
SYMIG	18	-2724.2090	5484.4181	0.0000	1.0000
HKY	15	-2731.5246	5493.0492	0.0000	1.0000
GTR	19	-2731.0162	5500.0323	0.0000	1.0000
K2P	12	-2748.3162	5520.6324	0.0000	1.0000
SYM	16	-2746.4693	5524.9385	0.0000	1.0000
F81G	15	-2799.7803	5629.5606	0.0000	1.0000
F81IG	16	-2799.7804	5631.5607	0.0000	1.0000
F81I	15	-2801.6411	5633.2823	0.0000	1.0000
JC69G	12	-2815.1754	5654.3508	0.0000	1.0000
JC69IG	13	-2815.1754	5656.3509	0.0000	1.0000
JC69I	12	-2816.9395	5657.8791	0.0000	1.0000
F81	14	-2821.1737	5670.3473	0.0000	1.0000
JC69	11	-2836.1606	5694.3213	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2709.2364	5524.1142	0.5536	0.5536
HKYI	16	-2709.5140	5524.6695	0.4194	0.9731
HKYIG	17	-2709.0066	5530.2573	0.0257	0.9987
K2PG	13	-2725.6818	5537.1971	0.0008	0.9995
K2PI	13	-2726.2767	5538.3870	0.0004	1.0000
K2PIG	14	-2725.6049	5543.6461	0.0000	1.0000
GTRG	20	-2708.7930	5549.6378	0.0000	1.0000
GTRI	20	-2709.0637	5550.1791	0.0000	1.0000
GTRIG	21	-2708.5766	5555.8076	0.0000	1.0000
SYMG	17	-2724.3122	5560.8685	0.0000	1.0000
SYMI	17	-2724.8186	5561.8812	0.0000	1.0000
HKY	15	-2731.5246	5562.0880	0.0000	1.0000
SYMIG	18	-2724.2090	5567.2646	0.0000	1.0000
K2P	12	-2748.3162	5575.8635	0.0000	1.0000
GTR	19	-2731.0162	5587.4815	0.0000	1.0000
SYM	16	-2746.4693	5598.5799	0.0000	1.0000
F81G	15	-2799.7803	5698.5995	0.0000	1.0000
F81I	15	-2801.6411	5702.3211	0.0000	1.0000
F81IG	16	-2799.7804	5705.2021	0.0000	1.0000
JC69G	12	-2815.1754	5709.5818	0.0000	1.0000
JC69I	12	-2816.9395	5713.1101	0.0000	1.0000
JC69IG	13	-2815.1754	5716.1845	0.0000	1.0000
F81	14	-2821.1737	5734.7836	0.0000	1.0000
JC69	11	-2836.1606	5744.9498	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 683-tomodel-phylip/chrZ_7476.phylip = 1 - 737;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_7476.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 683-tomodel-phylip/chrZ_7476.phylip = 1 - 737;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_7476.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 683-tomodel-phylip/chrZ_7476.phylip = 1 - 737;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_7476.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

