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

Input data from file "591-tomodel-phylip/chr20_1391.phylip" (7 taxa, 368 characters)

Minimum AIC  model: K2PI
Minimum AICc model: K2PI
Minimum BIC  model: K2PI


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
K2PI	13	-985.3871	1997.8025	0.3878	0.3878
K2PG	13	-985.7617	1998.5517	0.2666	0.6544
K2PIG	14	-985.7618	2000.7133	0.0905	0.7449
SYMI	17	-982.6528	2001.0542	0.0763	0.8212
SYMG	17	-982.7849	2001.3184	0.0669	0.8880
HKYI	16	-984.4827	2002.5153	0.0367	0.9248
HKYG	16	-984.8524	2003.2547	0.0254	0.9502
SYMIG	18	-982.7849	2003.5297	0.0221	0.9723
GTRI	20	-981.4917	2005.4042	0.0087	0.9810
GTRG	20	-981.4998	2005.4203	0.0086	0.9896
HKYIG	17	-984.8525	2005.4536	0.0085	0.9980
GTRIG	21	-981.8526	2008.3757	0.0020	1.0000
K2P	12	-1000.0354	2024.9498	0.0000	1.0000
SYM	16	-996.6068	2026.7634	0.0000	1.0000
JC69I	12	-1001.0212	2026.9213	0.0000	1.0000
JC69G	12	-1001.4337	2027.7464	0.0000	1.0000
HKY	15	-998.9838	2029.3313	0.0000	1.0000
JC69IG	13	-1001.4337	2029.8957	0.0000	1.0000
GTR	19	-994.9030	2029.9898	0.0000	1.0000
F81I	15	-1000.1543	2031.6723	0.0000	1.0000
F81G	15	-1000.5565	2032.4767	0.0000	1.0000
F81IG	16	-1000.5565	2034.6629	0.0000	1.0000
JC69	11	-1013.1710	2049.0836	0.0000	1.0000
F81	14	-1012.1312	2053.4522	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2PI	13	-985.3871	1996.7743	0.3447	0.3447
K2PG	13	-985.7617	1997.5234	0.2370	0.5818
SYMI	17	-982.6528	1999.3056	0.0972	0.6790
K2PIG	14	-985.7618	1999.5235	0.0872	0.7662
SYMG	17	-982.7849	1999.5698	0.0852	0.8514
HKYI	16	-984.4827	2000.9655	0.0424	0.8938
SYMIG	18	-982.7849	2001.5698	0.0313	0.9252
HKYG	16	-984.8524	2001.7049	0.0293	0.9545
GTRI	20	-981.4917	2002.9835	0.0155	0.9699
GTRG	20	-981.4998	2002.9996	0.0153	0.9853
HKYIG	17	-984.8525	2003.7050	0.0108	0.9960
GTRIG	21	-981.8526	2005.7052	0.0040	1.0000
K2P	12	-1000.0354	2024.0709	0.0000	1.0000
SYM	16	-996.6068	2025.2136	0.0000	1.0000
JC69I	12	-1001.0212	2026.0425	0.0000	1.0000
JC69G	12	-1001.4337	2026.8675	0.0000	1.0000
GTR	19	-994.9030	2027.8059	0.0000	1.0000
HKY	15	-998.9838	2027.9677	0.0000	1.0000
JC69IG	13	-1001.4337	2028.8675	0.0000	1.0000
F81I	15	-1000.1543	2030.3087	0.0000	1.0000
F81G	15	-1000.5565	2031.1130	0.0000	1.0000
F81IG	16	-1000.5565	2033.1130	0.0000	1.0000
JC69	11	-1013.1710	2048.3420	0.0000	1.0000
F81	14	-1012.1312	2052.2624	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	13	-985.3871	2047.5794	0.5800	0.5800
K2PG	13	-985.7617	2048.3285	0.3988	0.9787
K2PIG	14	-985.7618	2054.2367	0.0208	0.9995
HKYI	16	-984.4827	2063.4948	0.0002	0.9997
HKYG	16	-984.8524	2064.2342	0.0001	0.9999
SYMI	17	-982.6528	2065.7430	0.0001	0.9999
SYMG	17	-982.7849	2066.0072	0.0001	1.0000
HKYIG	17	-984.8525	2070.1424	0.0000	1.0000
K2P	12	-1000.0354	2070.9679	0.0000	1.0000
SYMIG	18	-982.7849	2071.9153	0.0000	1.0000
JC69I	12	-1001.0212	2072.9395	0.0000	1.0000
JC69G	12	-1001.4337	2073.7645	0.0000	1.0000
JC69IG	13	-1001.4337	2079.6725	0.0000	1.0000
GTRI	20	-981.4917	2081.1451	0.0000	1.0000
GTRG	20	-981.4998	2081.1612	0.0000	1.0000
HKY	15	-998.9838	2086.5889	0.0000	1.0000
SYM	16	-996.6068	2087.7429	0.0000	1.0000
GTRIG	21	-981.8526	2087.7749	0.0000	1.0000
F81I	15	-1000.1543	2088.9299	0.0000	1.0000
F81G	15	-1000.5565	2089.7343	0.0000	1.0000
JC69	11	-1013.1710	2091.3309	0.0000	1.0000
F81IG	16	-1000.5565	2095.6423	0.0000	1.0000
GTR	19	-994.9030	2102.0595	0.0000	1.0000
F81	14	-1012.1312	2106.9756	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr20_1391.phylip = 1 - 368;
 Partition Dummy = 1:591-tomodel-phylip/chr20_1391.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr20_1391.phylip = 1 - 368;
 Partition Dummy = 1:591-tomodel-phylip/chr20_1391.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr20_1391.phylip = 1 - 368;
 Partition Dummy = 1:591-tomodel-phylip/chr20_1391.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

