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

Input data from file "591-tomodel-phylip/chr8_4221.phylip" (7 taxa, 744 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	-2378.4583	4789.6649	0.3666	0.3666
GTRI	20	-2374.4822	4790.1263	0.2910	0.6576
HKYIG	17	-2378.4792	4791.8013	0.1260	0.7836
GTRIG	21	-2374.4952	4792.2702	0.0996	0.8832
HKYG	16	-2380.1543	4793.0569	0.0672	0.9505
GTRG	20	-2376.2529	4793.6676	0.0495	1.0000
GTR	19	-2410.3649	4859.7795	0.0000	1.0000
HKY	15	-2416.3818	4863.4230	0.0000	1.0000
K2PI	13	-2434.1993	4894.8973	0.0000	1.0000
K2PG	13	-2435.3967	4897.2921	0.0000	1.0000
K2PIG	14	-2436.0096	4900.5953	0.0000	1.0000
SYMI	17	-2433.8713	4902.5857	0.0000	1.0000
SYMG	17	-2434.9932	4904.8294	0.0000	1.0000
SYMIG	18	-2435.1667	4907.2768	0.0000	1.0000
F81I	15	-2447.5702	4925.7997	0.0000	1.0000
F81G	15	-2448.1873	4927.0340	0.0000	1.0000
F81IG	16	-2447.4743	4927.6968	0.0000	1.0000
K2P	12	-2470.2856	4964.9981	0.0000	1.0000
SYM	16	-2470.1454	4973.0391	0.0000	1.0000
F81	14	-2481.0243	4990.6247	0.0000	1.0000
JC69I	12	-2492.7028	5009.8323	0.0000	1.0000
JC69G	12	-2493.4197	5011.2663	0.0000	1.0000
JC69IG	13	-2493.3260	5013.1506	0.0000	1.0000
JC69	11	-2527.2360	5076.8326	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-2378.4583	4788.9166	0.3289	0.3289
GTRI	20	-2374.4822	4788.9645	0.3211	0.6499
HKYIG	17	-2378.4792	4790.9583	0.1185	0.7684
GTRIG	21	-2374.4952	4790.9904	0.1166	0.8850
HKYG	16	-2380.1543	4792.3086	0.0603	0.9453
GTRG	20	-2376.2529	4792.5058	0.0547	1.0000
GTR	19	-2410.3649	4858.7298	0.0000	1.0000
HKY	15	-2416.3818	4862.7636	0.0000	1.0000
K2PI	13	-2434.1993	4894.3987	0.0000	1.0000
K2PG	13	-2435.3967	4896.7934	0.0000	1.0000
K2PIG	14	-2436.0096	4900.0192	0.0000	1.0000
SYMI	17	-2433.8713	4901.7427	0.0000	1.0000
SYMG	17	-2434.9932	4903.9864	0.0000	1.0000
SYMIG	18	-2435.1667	4906.3333	0.0000	1.0000
F81I	15	-2447.5702	4925.1404	0.0000	1.0000
F81G	15	-2448.1873	4926.3747	0.0000	1.0000
F81IG	16	-2447.4743	4926.9485	0.0000	1.0000
K2P	12	-2470.2856	4964.5713	0.0000	1.0000
SYM	16	-2470.1454	4972.2908	0.0000	1.0000
F81	14	-2481.0243	4990.0486	0.0000	1.0000
JC69I	12	-2492.7028	5009.4055	0.0000	1.0000
JC69G	12	-2493.4197	5010.8395	0.0000	1.0000
JC69IG	13	-2493.3260	5012.6520	0.0000	1.0000
JC69	11	-2527.2360	5076.4719	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2378.4583	4862.7093	0.8200	0.8200
HKYG	16	-2380.1543	4866.1012	0.1504	0.9705
HKYIG	17	-2378.4792	4869.3630	0.0294	0.9999
GTRI	20	-2374.4822	4881.2053	0.0001	1.0000
GTRG	20	-2376.2529	4884.7466	0.0000	1.0000
GTRIG	21	-2374.4952	4887.8433	0.0000	1.0000
HKY	15	-2416.3818	4931.9442	0.0000	1.0000
GTR	19	-2410.3649	4946.3586	0.0000	1.0000
K2PI	13	-2434.1993	4954.3552	0.0000	1.0000
K2PG	13	-2435.3967	4956.7500	0.0000	1.0000
K2PIG	14	-2436.0096	4964.5878	0.0000	1.0000
SYMI	17	-2433.8713	4980.1474	0.0000	1.0000
SYMG	17	-2434.9932	4982.3911	0.0000	1.0000
SYMIG	18	-2435.1667	4989.3501	0.0000	1.0000
F81I	15	-2447.5702	4994.3210	0.0000	1.0000
F81G	15	-2448.1873	4995.5553	0.0000	1.0000
F81IG	16	-2447.4743	5000.7412	0.0000	1.0000
K2P	12	-2470.2856	5019.9158	0.0000	1.0000
SYM	16	-2470.1454	5046.0835	0.0000	1.0000
F81	14	-2481.0243	5054.6171	0.0000	1.0000
JC69I	12	-2492.7028	5064.7500	0.0000	1.0000
JC69G	12	-2493.4197	5066.1840	0.0000	1.0000
JC69IG	13	-2493.3260	5072.6085	0.0000	1.0000
JC69	11	-2527.2360	5127.2044	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 591-tomodel-phylip/chr8_4221.phylip = 1 - 744;
 Partition Dummy = 1:591-tomodel-phylip/chr8_4221.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 591-tomodel-phylip/chr8_4221.phylip = 1 - 744;
 Partition Dummy = 1:591-tomodel-phylip/chr8_4221.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 591-tomodel-phylip/chr8_4221.phylip = 1 - 744;
 Partition Dummy = 1:591-tomodel-phylip/chr8_4221.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

