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

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

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
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
HKYI	16	-2410.1736	4853.1333	0.5616	0.5616
HKYIG	17	-2410.1858	4855.2573	0.1942	0.7558
HKYG	16	-2411.3784	4855.5428	0.1684	0.9242
GTRI	20	-2408.4144	4858.0496	0.0481	0.9722
GTRIG	21	-2408.4309	4860.2067	0.0163	0.9886
GTRG	20	-2409.8512	4860.9234	0.0114	1.0000
K2PI	13	-2438.1728	4902.8694	0.0000	1.0000
K2PG	13	-2439.0698	4904.6634	0.0000	1.0000
K2PIG	14	-2438.1806	4904.9664	0.0000	1.0000
SYMI	17	-2437.0409	4908.9675	0.0000	1.0000
SYMG	17	-2437.4858	4909.8573	0.0000	1.0000
SYMIG	18	-2437.0410	4911.0733	0.0000	1.0000
HKY	15	-2441.3701	4913.4329	0.0000	1.0000
GTR	19	-2440.4086	4919.9203	0.0000	1.0000
K2P	12	-2468.0505	4960.5493	0.0000	1.0000
SYM	16	-2466.9856	4966.7574	0.0000	1.0000
F81I	15	-2472.6068	4975.9063	0.0000	1.0000
F81IG	16	-2472.6278	4978.0416	0.0000	1.0000
F81G	15	-2474.0311	4978.7549	0.0000	1.0000
JC69I	12	-2495.7788	5016.0059	0.0000	1.0000
JC69IG	13	-2495.7935	5018.1107	0.0000	1.0000
JC69G	12	-2496.9192	5018.2867	0.0000	1.0000
F81	14	-2500.9410	5030.4871	0.0000	1.0000
JC69	11	-2522.5645	5067.5077	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-2410.1736	4852.3472	0.5455	0.5455
HKYIG	17	-2410.1858	4854.3716	0.1982	0.7437
HKYG	16	-2411.3784	4854.7567	0.1635	0.9072
GTRI	20	-2408.4144	4856.8287	0.0580	0.9652
GTRIG	21	-2408.4309	4858.8617	0.0210	0.9862
GTRG	20	-2409.8512	4859.7024	0.0138	1.0000
K2PI	13	-2438.1728	4902.3456	0.0000	1.0000
K2PG	13	-2439.0698	4904.1396	0.0000	1.0000
K2PIG	14	-2438.1806	4904.3612	0.0000	1.0000
SYMI	17	-2437.0409	4908.0818	0.0000	1.0000
SYMG	17	-2437.4858	4908.9716	0.0000	1.0000
SYMIG	18	-2437.0410	4910.0820	0.0000	1.0000
HKY	15	-2441.3701	4912.7402	0.0000	1.0000
GTR	19	-2440.4086	4918.8173	0.0000	1.0000
K2P	12	-2468.0505	4960.1010	0.0000	1.0000
SYM	16	-2466.9856	4965.9713	0.0000	1.0000
F81I	15	-2472.6068	4975.2137	0.0000	1.0000
F81IG	16	-2472.6278	4977.2555	0.0000	1.0000
F81G	15	-2474.0311	4978.0623	0.0000	1.0000
JC69I	12	-2495.7788	5015.5576	0.0000	1.0000
JC69IG	13	-2495.7935	5017.5869	0.0000	1.0000
JC69G	12	-2496.9192	5017.8384	0.0000	1.0000
F81	14	-2500.9410	5029.8819	0.0000	1.0000
JC69	11	-2522.5645	5067.1289	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2410.1736	4925.3689	0.7480	0.7480
HKYG	16	-2411.3784	4927.7784	0.2242	0.9722
HKYIG	17	-2410.1858	4931.9571	0.0278	1.0000
GTRI	20	-2408.4144	4948.1058	0.0000	1.0000
GTRG	20	-2409.8512	4950.9796	0.0000	1.0000
GTRIG	21	-2408.4309	4954.7027	0.0000	1.0000
K2PI	13	-2438.1728	4961.6758	0.0000	1.0000
K2PG	13	-2439.0698	4963.4697	0.0000	1.0000
K2PIG	14	-2438.1806	4968.2552	0.0000	1.0000
HKY	15	-2441.3701	4981.1981	0.0000	1.0000
SYMI	17	-2437.0409	4985.6673	0.0000	1.0000
SYMG	17	-2437.4858	4986.5572	0.0000	1.0000
SYMIG	18	-2437.0410	4992.2314	0.0000	1.0000
GTR	19	-2440.4086	5005.5306	0.0000	1.0000
K2P	12	-2468.0505	5014.8673	0.0000	1.0000
SYM	16	-2466.9856	5038.9930	0.0000	1.0000
F81I	15	-2472.6068	5043.6715	0.0000	1.0000
F81G	15	-2474.0311	5046.5201	0.0000	1.0000
F81IG	16	-2472.6278	5050.2772	0.0000	1.0000
JC69I	12	-2495.7788	5070.3239	0.0000	1.0000
JC69G	12	-2496.9192	5072.6047	0.0000	1.0000
JC69IG	13	-2495.7935	5076.9171	0.0000	1.0000
F81	14	-2500.9410	5093.7759	0.0000	1.0000
JC69	11	-2522.5645	5117.3313	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr2_14397.phylip = 1 - 709;
 Partition Dummy = 1:683-tomodel-phylip/chr2_14397.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;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr2_14397.phylip = 1 - 709;
 Partition Dummy = 1:683-tomodel-phylip/chr2_14397.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;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr2_14397.phylip = 1 - 709;
 Partition Dummy = 1:683-tomodel-phylip/chr2_14397.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

