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

Input data from file "591-tomodel-phylip/chr13_2246.phylip" (7 taxa, 554 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	-2185.2224	4403.4578	0.3760	0.3760
GTRG	20	-2181.3860	4404.3480	0.2409	0.6169
HKYIG	17	-2185.1080	4405.3577	0.1454	0.7623
GTRIG	21	-2181.2432	4406.2233	0.0943	0.8566
HKYI	16	-2186.6345	4406.2821	0.0916	0.9482
GTRI	20	-2182.9249	4407.4257	0.0517	0.9999
K2PG	13	-2197.7669	4422.2079	0.0000	1.0000
K2PIG	14	-2197.7038	4424.1867	0.0000	1.0000
K2PI	13	-2199.2092	4425.0925	0.0000	1.0000
SYMG	17	-2195.1539	4425.4497	0.0000	1.0000
SYMIG	18	-2195.1656	4427.6097	0.0000	1.0000
SYMI	17	-2197.1566	4429.4550	0.0000	1.0000
HKY	15	-2218.8574	4468.6069	0.0000	1.0000
GTR	19	-2215.3255	4470.0743	0.0000	1.0000
K2P	12	-2231.9132	4488.4031	0.0000	1.0000
SYM	16	-2229.9478	4492.9086	0.0000	1.0000
F81G	15	-2270.4498	4571.7917	0.0000	1.0000
F81I	15	-2270.7197	4572.3316	0.0000	1.0000
F81IG	16	-2270.2126	4573.4381	0.0000	1.0000
JC69G	12	-2280.1133	4584.8033	0.0000	1.0000
JC69I	12	-2280.6055	4585.7878	0.0000	1.0000
JC69IG	13	-2279.9348	4586.5437	0.0000	1.0000
F81	14	-2300.7930	4630.3652	0.0000	1.0000
JC69	11	-2310.4137	4643.3145	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2185.2224	4402.4447	0.3281	0.3281
GTRG	20	-2181.3860	4402.7720	0.2786	0.6067
HKYIG	17	-2185.1080	4404.2159	0.1353	0.7420
GTRIG	21	-2181.2432	4404.4865	0.1182	0.8602
HKYI	16	-2186.6345	4405.2691	0.0799	0.9402
GTRI	20	-2182.9249	4405.8497	0.0598	1.0000
K2PG	13	-2197.7669	4421.5338	0.0000	1.0000
K2PIG	14	-2197.7038	4423.4075	0.0000	1.0000
SYMG	17	-2195.1539	4424.3079	0.0000	1.0000
K2PI	13	-2199.2092	4424.4184	0.0000	1.0000
SYMIG	18	-2195.1656	4426.3312	0.0000	1.0000
SYMI	17	-2197.1566	4428.3132	0.0000	1.0000
HKY	15	-2218.8574	4467.7147	0.0000	1.0000
GTR	19	-2215.3255	4468.6511	0.0000	1.0000
K2P	12	-2231.9132	4487.8264	0.0000	1.0000
SYM	16	-2229.9478	4491.8956	0.0000	1.0000
F81G	15	-2270.4498	4570.8995	0.0000	1.0000
F81I	15	-2270.7197	4571.4394	0.0000	1.0000
F81IG	16	-2270.2126	4572.4251	0.0000	1.0000
JC69G	12	-2280.1133	4584.2266	0.0000	1.0000
JC69I	12	-2280.6055	4585.2111	0.0000	1.0000
JC69IG	13	-2279.9348	4585.8696	0.0000	1.0000
F81	14	-2300.7930	4629.5860	0.0000	1.0000
JC69	11	-2310.4137	4642.8274	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2185.2224	4471.5194	0.7402	0.7402
HKYI	16	-2186.6345	4474.3437	0.1803	0.9205
HKYIG	17	-2185.1080	4477.6077	0.0353	0.9558
K2PG	13	-2197.7669	4477.6569	0.0344	0.9902
K2PI	13	-2199.2092	4480.5416	0.0081	0.9983
K2PIG	14	-2197.7038	4483.8478	0.0016	0.9999
GTRG	20	-2181.3860	4489.1153	0.0001	1.0000
GTRI	20	-2182.9249	4492.1930	0.0000	1.0000
GTRIG	21	-2181.2432	4495.1470	0.0000	1.0000
SYMG	17	-2195.1539	4497.6997	0.0000	1.0000
SYMI	17	-2197.1566	4501.7050	0.0000	1.0000
SYMIG	18	-2195.1656	4504.0402	0.0000	1.0000
HKY	15	-2218.8574	4532.4722	0.0000	1.0000
K2P	12	-2231.9132	4539.6324	0.0000	1.0000
GTR	19	-2215.3255	4550.6772	0.0000	1.0000
SYM	16	-2229.9478	4560.9702	0.0000	1.0000
F81G	15	-2270.4498	4635.6570	0.0000	1.0000
JC69G	12	-2280.1133	4636.0326	0.0000	1.0000
F81I	15	-2270.7197	4636.1969	0.0000	1.0000
JC69I	12	-2280.6055	4637.0171	0.0000	1.0000
F81IG	16	-2270.2126	4641.4997	0.0000	1.0000
JC69IG	13	-2279.9348	4641.9927	0.0000	1.0000
F81	14	-2300.7930	4690.0263	0.0000	1.0000
JC69	11	-2310.4137	4690.3162	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/chr13_2246.phylip = 1 - 554;
 Partition Dummy = 1:591-tomodel-phylip/chr13_2246.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/chr13_2246.phylip = 1 - 554;
 Partition Dummy = 1:591-tomodel-phylip/chr13_2246.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/chr13_2246.phylip = 1 - 554;
 Partition Dummy = 1:591-tomodel-phylip/chr13_2246.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

