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

Input data from file "591-tomodel-phylip/chr1_8428.phylip" (7 taxa, 773 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	-2359.4326	4751.5848	0.2621	0.2621
HKYG	16	-2359.6305	4751.9807	0.2150	0.4771
GTRG	20	-2355.5878	4752.2926	0.1840	0.6611
GTRI	20	-2355.6695	4752.4560	0.1695	0.8306
HKYIG	17	-2359.3914	4753.5934	0.0960	0.9266
GTRIG	21	-2355.4503	4754.1310	0.0734	1.0000
K2PI	13	-2373.1923	4772.8642	0.0000	1.0000
K2PG	13	-2373.4244	4773.3284	0.0000	1.0000
K2PIG	14	-2373.1660	4774.8862	0.0000	1.0000
SYMI	17	-2372.0977	4779.0059	0.0000	1.0000
SYMG	17	-2372.3004	4779.4114	0.0000	1.0000
SYMIG	18	-2372.0606	4781.0284	0.0000	1.0000
HKY	15	-2380.5117	4791.6575	0.0000	1.0000
GTR	19	-2377.2013	4793.4120	0.0000	1.0000
F81I	15	-2385.4504	4801.5349	0.0000	1.0000
F81G	15	-2386.4196	4803.4732	0.0000	1.0000
F81IG	16	-2385.4663	4803.6522	0.0000	1.0000
K2P	12	-2394.5696	4813.5496	0.0000	1.0000
JC69I	12	-2397.1129	4818.6363	0.0000	1.0000
SYM	16	-2393.5072	4819.7341	0.0000	1.0000
JC69G	12	-2398.0849	4820.5803	0.0000	1.0000
JC69IG	13	-2397.1292	4820.7380	0.0000	1.0000
F81	14	-2406.2288	4841.0118	0.0000	1.0000
JC69	11	-2418.1491	4858.6450	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-2359.4326	4750.8652	0.2375	0.2375
GTRG	20	-2355.5878	4751.1756	0.2034	0.4409
HKYG	16	-2359.6305	4751.2611	0.1948	0.6357
GTRI	20	-2355.6695	4751.3389	0.1874	0.8231
HKYIG	17	-2359.3914	4752.7828	0.0910	0.9142
GTRIG	21	-2355.4503	4752.9006	0.0858	1.0000
K2PI	13	-2373.1923	4772.3846	0.0000	1.0000
K2PG	13	-2373.4244	4772.8489	0.0000	1.0000
K2PIG	14	-2373.1660	4774.3321	0.0000	1.0000
SYMI	17	-2372.0977	4778.1953	0.0000	1.0000
SYMG	17	-2372.3004	4778.6008	0.0000	1.0000
SYMIG	18	-2372.0606	4780.1213	0.0000	1.0000
HKY	15	-2380.5117	4791.0235	0.0000	1.0000
GTR	19	-2377.2013	4792.4027	0.0000	1.0000
F81I	15	-2385.4504	4800.9009	0.0000	1.0000
F81G	15	-2386.4196	4802.8391	0.0000	1.0000
F81IG	16	-2385.4663	4802.9326	0.0000	1.0000
K2P	12	-2394.5696	4813.1391	0.0000	1.0000
JC69I	12	-2397.1129	4818.2258	0.0000	1.0000
SYM	16	-2393.5072	4819.0145	0.0000	1.0000
JC69G	12	-2398.0849	4820.1698	0.0000	1.0000
JC69IG	13	-2397.1292	4820.2584	0.0000	1.0000
F81	14	-2406.2288	4840.4577	0.0000	1.0000
JC69	11	-2418.1491	4858.2981	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2359.4326	4825.2697	0.5264	0.5264
HKYG	16	-2359.6305	4825.6655	0.4319	0.9583
HKYIG	17	-2359.3914	4831.8375	0.0197	0.9780
K2PI	13	-2373.1923	4832.8382	0.0120	0.9900
K2PG	13	-2373.4244	4833.3025	0.0095	0.9995
K2PIG	14	-2373.1660	4839.4360	0.0004	0.9999
GTRG	20	-2355.5878	4844.1812	0.0000	1.0000
GTRI	20	-2355.6695	4844.3445	0.0000	1.0000
GTRIG	21	-2355.4503	4850.5565	0.0000	1.0000
SYMI	17	-2372.0977	4857.2500	0.0000	1.0000
SYMG	17	-2372.3004	4857.6556	0.0000	1.0000
HKY	15	-2380.5117	4860.7776	0.0000	1.0000
SYMIG	18	-2372.0606	4863.8263	0.0000	1.0000
K2P	12	-2394.5696	4868.9424	0.0000	1.0000
F81I	15	-2385.4504	4870.6550	0.0000	1.0000
F81G	15	-2386.4196	4872.5933	0.0000	1.0000
JC69I	12	-2397.1129	4874.0291	0.0000	1.0000
JC69G	12	-2398.0849	4875.9731	0.0000	1.0000
F81IG	16	-2385.4663	4877.3371	0.0000	1.0000
JC69IG	13	-2397.1292	4880.7120	0.0000	1.0000
GTR	19	-2377.2013	4880.7580	0.0000	1.0000
SYM	16	-2393.5072	4893.4189	0.0000	1.0000
F81	14	-2406.2288	4905.5616	0.0000	1.0000
JC69	11	-2418.1491	4909.4512	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/chr1_8428.phylip = 1 - 773;
 Partition Dummy = 1:591-tomodel-phylip/chr1_8428.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/chr1_8428.phylip = 1 - 773;
 Partition Dummy = 1:591-tomodel-phylip/chr1_8428.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/chr1_8428.phylip = 1 - 773;
 Partition Dummy = 1:591-tomodel-phylip/chr1_8428.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

