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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: HKYI


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
GTRI	20	-2132.4772	4306.0641	0.4694	0.4694
GTRG	20	-2133.0299	4307.1694	0.2701	0.7395
GTRIG	21	-2132.4846	4308.1915	0.1620	0.9015
HKYI	16	-2138.9641	4310.6430	0.0476	0.9491
HKYG	16	-2139.2931	4311.3011	0.0342	0.9833
HKYIG	17	-2138.9669	4312.7391	0.0167	1.0000
GTR	19	-2145.4803	4329.9633	0.0000	1.0000
HKY	15	-2151.0408	4332.7115	0.0000	1.0000
F81I	15	-2184.1222	4398.8743	0.0000	1.0000
F81G	15	-2184.4457	4399.5214	0.0000	1.0000
F81IG	16	-2184.1251	4400.9650	0.0000	1.0000
K2PI	13	-2192.4487	4411.3738	0.0000	1.0000
K2PG	13	-2192.7779	4412.0323	0.0000	1.0000
K2PIG	14	-2192.4512	4413.4529	0.0000	1.0000
SYMI	17	-2191.7287	4418.2626	0.0000	1.0000
F81	14	-2195.0974	4418.7453	0.0000	1.0000
SYMG	17	-2192.2007	4419.2066	0.0000	1.0000
SYMIG	18	-2191.7348	4420.3709	0.0000	1.0000
K2P	12	-2204.4466	4433.3011	0.0000	1.0000
SYM	16	-2203.6285	4439.9718	0.0000	1.0000
JC69I	12	-2229.2398	4482.8874	0.0000	1.0000
JC69G	12	-2229.5784	4483.5647	0.0000	1.0000
JC69IG	13	-2229.2419	4484.9602	0.0000	1.0000
JC69	11	-2240.3639	4503.0723	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-2132.4772	4304.9545	0.4730	0.4730
GTRG	20	-2133.0299	4306.0598	0.2722	0.7452
GTRIG	21	-2132.4846	4306.9692	0.1727	0.9179
HKYI	16	-2138.9641	4309.9282	0.0393	0.9573
HKYG	16	-2139.2931	4310.5863	0.0283	0.9856
HKYIG	17	-2138.9669	4311.9338	0.0144	1.0000
GTR	19	-2145.4803	4328.9606	0.0000	1.0000
HKY	15	-2151.0408	4332.0816	0.0000	1.0000
F81I	15	-2184.1222	4398.2444	0.0000	1.0000
F81G	15	-2184.4457	4398.8915	0.0000	1.0000
F81IG	16	-2184.1251	4400.2501	0.0000	1.0000
K2PI	13	-2192.4487	4410.8974	0.0000	1.0000
K2PG	13	-2192.7779	4411.5558	0.0000	1.0000
K2PIG	14	-2192.4512	4412.9024	0.0000	1.0000
SYMI	17	-2191.7287	4417.4573	0.0000	1.0000
F81	14	-2195.0974	4418.1948	0.0000	1.0000
SYMG	17	-2192.2007	4418.4013	0.0000	1.0000
SYMIG	18	-2191.7348	4419.4697	0.0000	1.0000
K2P	12	-2204.4466	4432.8933	0.0000	1.0000
SYM	16	-2203.6285	4439.2569	0.0000	1.0000
JC69I	12	-2229.2398	4482.4796	0.0000	1.0000
JC69G	12	-2229.5784	4483.1568	0.0000	1.0000
JC69IG	13	-2229.2419	4484.4837	0.0000	1.0000
JC69	11	-2240.3639	4502.7277	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2138.9641	4384.4358	0.5691	0.5691
HKYG	16	-2139.2931	4385.0939	0.4095	0.9786
HKYIG	17	-2138.9669	4391.0982	0.0203	0.9989
GTRI	20	-2132.4772	4398.0890	0.0006	0.9995
GTRG	20	-2133.0299	4399.1943	0.0004	0.9999
HKY	15	-2151.0408	4401.9325	0.0001	1.0000
GTRIG	21	-2132.4846	4404.7605	0.0000	1.0000
GTR	19	-2145.4803	4417.4384	0.0000	1.0000
F81I	15	-2184.1222	4468.0953	0.0000	1.0000
F81G	15	-2184.4457	4468.7424	0.0000	1.0000
K2PI	13	-2192.4487	4471.4348	0.0000	1.0000
K2PG	13	-2192.7779	4472.0933	0.0000	1.0000
F81IG	16	-2184.1251	4474.7578	0.0000	1.0000
K2PIG	14	-2192.4512	4478.0966	0.0000	1.0000
F81	14	-2195.0974	4483.3890	0.0000	1.0000
K2P	12	-2204.4466	4488.7740	0.0000	1.0000
SYMI	17	-2191.7287	4496.6217	0.0000	1.0000
SYMG	17	-2192.2007	4497.5657	0.0000	1.0000
SYMIG	18	-2191.7348	4503.2908	0.0000	1.0000
SYM	16	-2203.6285	4513.7645	0.0000	1.0000
JC69I	12	-2229.2398	4538.3603	0.0000	1.0000
JC69G	12	-2229.5784	4539.0376	0.0000	1.0000
JC69IG	13	-2229.2419	4545.0212	0.0000	1.0000
JC69	11	-2240.3639	4553.9517	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (GTRI)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr2_8747.phylip = 1 - 778;
 Partition Dummy = 1:683-tomodel-phylip/chr2_8747.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (GTRI)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr2_8747.phylip = 1 - 778;
 Partition Dummy = 1:683-tomodel-phylip/chr2_8747.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best BIC model (HKYI)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr2_8747.phylip = 1 - 778;
 Partition Dummy = 1:683-tomodel-phylip/chr2_8747.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


-------------------------------------------------------------
End of Output

