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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKY


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
HKYG	16	-809.5909	1652.6363	0.2942	0.2942
HKYI	16	-809.5918	1652.6380	0.2940	0.5882
HKY	15	-811.2115	1653.7031	0.1726	0.7608
HKYIG	17	-809.5948	1654.8304	0.0982	0.8590
GTRG	20	-807.0369	1656.3441	0.0461	0.9051
GTRI	20	-807.0511	1656.3725	0.0454	0.9505
GTR	19	-808.4369	1656.9223	0.0345	0.9850
GTRIG	21	-807.0425	1658.5891	0.0150	1.0000
K2PI	13	-824.5697	1676.1049	0.0000	1.0000
K2PG	13	-824.6177	1676.2009	0.0000	1.0000
K2P	12	-826.3694	1677.5642	0.0000	1.0000
K2PIG	14	-824.6282	1678.3735	0.0000	1.0000
F81I	15	-824.0462	1679.3724	0.0000	1.0000
F81G	15	-824.0567	1679.3934	0.0000	1.0000
F81	14	-825.6607	1680.4385	0.0000	1.0000
SYMI	17	-822.5439	1680.7286	0.0000	1.0000
SYMG	17	-822.5974	1680.8356	0.0000	1.0000
F81IG	16	-824.0801	1681.6147	0.0000	1.0000
SYM	16	-824.1119	1681.6783	0.0000	1.0000
SYMIG	18	-822.6055	1683.0496	0.0000	1.0000
JC69I	12	-838.0784	1700.9821	0.0000	1.0000
JC69G	12	-838.1252	1701.0758	0.0000	1.0000
JC69	11	-839.7219	1702.1404	0.0000	1.0000
JC69IG	13	-838.1324	1703.2304	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-809.5909	1651.1817	0.2766	0.2766
HKYI	16	-809.5918	1651.1835	0.2764	0.5530
HKY	15	-811.2115	1652.4231	0.1487	0.7017
HKYIG	17	-809.5948	1653.1897	0.1014	0.8031
GTRG	20	-807.0369	1654.0739	0.0651	0.8683
GTRI	20	-807.0511	1654.1023	0.0642	0.9325
GTR	19	-808.4369	1654.8738	0.0437	0.9762
GTRIG	21	-807.0425	1656.0850	0.0238	1.0000
K2PI	13	-824.5697	1675.1394	0.0000	1.0000
K2PG	13	-824.6177	1675.2354	0.0000	1.0000
K2P	12	-826.3694	1676.7388	0.0000	1.0000
K2PIG	14	-824.6282	1677.2565	0.0000	1.0000
F81I	15	-824.0462	1678.0924	0.0000	1.0000
F81G	15	-824.0567	1678.1134	0.0000	1.0000
SYMI	17	-822.5439	1679.0878	0.0000	1.0000
SYMG	17	-822.5974	1679.1948	0.0000	1.0000
F81	14	-825.6607	1679.3214	0.0000	1.0000
F81IG	16	-824.0801	1680.1602	0.0000	1.0000
SYM	16	-824.1119	1680.2238	0.0000	1.0000
SYMIG	18	-822.6055	1681.2109	0.0000	1.0000
JC69I	12	-838.0784	1700.1567	0.0000	1.0000
JC69G	12	-838.1252	1700.2504	0.0000	1.0000
JC69	11	-839.7219	1701.4439	0.0000	1.0000
JC69IG	13	-838.1324	1702.2649	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKY	15	-811.2115	1711.9537	0.6547	0.6547
HKYG	16	-809.5909	1714.6811	0.1674	0.8221
HKYI	16	-809.5918	1714.6828	0.1673	0.9893
HKYIG	17	-809.5948	1720.6577	0.0084	0.9978
K2P	12	-826.3694	1724.3633	0.0013	0.9991
K2PI	13	-824.5697	1726.7326	0.0004	0.9995
K2PG	13	-824.6177	1726.8286	0.0004	0.9999
GTR	19	-808.4369	1730.2793	0.0001	0.9999
K2PIG	14	-824.6282	1732.8184	0.0000	1.0000
GTRG	20	-807.0369	1733.4480	0.0000	1.0000
GTRI	20	-807.0511	1733.4764	0.0000	1.0000
F81	14	-825.6607	1734.8833	0.0000	1.0000
F81I	15	-824.0462	1737.6230	0.0000	1.0000
F81G	15	-824.0567	1737.6440	0.0000	1.0000
GTRIG	21	-807.0425	1739.4279	0.0000	1.0000
F81IG	16	-824.0801	1743.6595	0.0000	1.0000
SYM	16	-824.1119	1743.7231	0.0000	1.0000
JC69	11	-839.7219	1745.0996	0.0000	1.0000
SYMI	17	-822.5439	1746.5558	0.0000	1.0000
SYMG	17	-822.5974	1746.6629	0.0000	1.0000
JC69I	12	-838.0784	1747.7812	0.0000	1.0000
JC69G	12	-838.1252	1747.8749	0.0000	1.0000
SYMIG	18	-822.6055	1752.6477	0.0000	1.0000
JC69IG	13	-838.1324	1753.8581	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr9_3566.phylip = 1 - 391;
 Partition Dummy = 1:591-tomodel-phylip/chr9_3566.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr9_3566.phylip = 1 - 391;
 Partition Dummy = 1:591-tomodel-phylip/chr9_3566.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr9_3566.phylip = 1 - 391;
 Partition Dummy = 1:591-tomodel-phylip/chr9_3566.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

