
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
00:44 02/17/2011

Input data from file "183-phylip/chr3_18148.phylip" (29 taxa, 197 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	57	-771.3874	1704.3432	0.5204	0.5204
K2PG	57	-771.6284	1704.8252	0.4089	0.9293
K2PIG	58	-771.6279	1708.8500	0.0547	0.9840
HKYI	60	-769.4556	1712.7348	0.0078	0.9918
HKYG	60	-769.7882	1713.4000	0.0056	0.9974
HKYIG	61	-769.4552	1716.9401	0.0010	0.9984
SYMI	61	-769.5564	1717.1424	0.0009	0.9992
SYMG	61	-769.8258	1717.6811	0.0007	0.9999
SYMIG	62	-769.8853	1722.0691	0.0001	1.0000
GTRI	64	-767.3608	1725.7520	0.0000	1.0000
GTRG	64	-767.7407	1726.5118	0.0000	1.0000
K2P	56	-785.6944	1728.9887	0.0000	1.0000
GTRIG	65	-767.3378	1730.1718	0.0000	1.0000
HKY	59	-783.2347	1736.1482	0.0000	1.0000
SYM	60	-783.3886	1740.6008	0.0000	1.0000
GTR	63	-780.7866	1748.2048	0.0000	1.0000
JC69I	56	-807.1931	1771.9861	0.0000	1.0000
JC69G	56	-807.4880	1772.5760	0.0000	1.0000
JC69IG	57	-807.4874	1776.5431	0.0000	1.0000
F81I	59	-805.4879	1780.6547	0.0000	1.0000
F81G	59	-805.8386	1781.3561	0.0000	1.0000
F81IG	60	-805.9466	1785.7167	0.0000	1.0000
JC69	55	-820.7857	1795.2593	0.0000	1.0000
F81	58	-819.0143	1803.6228	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2PI	57	-771.3874	1656.7748	0.3196	0.3196
K2PG	57	-771.6284	1657.2569	0.2512	0.5708
HKYI	60	-769.4556	1658.9112	0.1098	0.6806
K2PIG	58	-771.6279	1659.2558	0.0924	0.7731
HKYG	60	-769.7882	1659.5765	0.0788	0.8518
HKYIG	61	-769.4552	1660.9105	0.0404	0.8922
SYMI	61	-769.5564	1661.1127	0.0365	0.9288
SYMG	61	-769.8258	1661.6515	0.0279	0.9567
GTRI	64	-767.3608	1662.7217	0.0163	0.9730
GTRG	64	-767.7407	1663.4815	0.0112	0.9842
SYMIG	62	-769.8853	1663.7706	0.0097	0.9938
GTRIG	65	-767.3378	1664.6756	0.0062	1.0000
K2P	56	-785.6944	1683.3887	0.0000	1.0000
HKY	59	-783.2347	1684.4693	0.0000	1.0000
SYM	60	-783.3886	1686.7773	0.0000	1.0000
GTR	63	-780.7866	1687.5732	0.0000	1.0000
JC69I	56	-807.1931	1726.3861	0.0000	1.0000
JC69G	56	-807.4880	1726.9760	0.0000	1.0000
JC69IG	57	-807.4874	1728.9748	0.0000	1.0000
F81I	59	-805.4879	1728.9758	0.0000	1.0000
F81G	59	-805.8386	1729.6772	0.0000	1.0000
F81IG	60	-805.9466	1731.8932	0.0000	1.0000
JC69	55	-820.7857	1751.5714	0.0000	1.0000
F81	58	-819.0143	1754.0286	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	57	-771.3874	1843.9174	0.5415	0.5415
K2PG	57	-771.6284	1844.3995	0.4255	0.9671
K2PIG	58	-771.6279	1849.6817	0.0303	0.9974
HKYI	60	-769.4556	1855.9035	0.0014	0.9988
HKYG	60	-769.7882	1856.5687	0.0010	0.9997
HKYIG	61	-769.4552	1861.1859	0.0001	0.9998
SYMI	61	-769.5564	1861.3882	0.0001	0.9999
SYMG	61	-769.8258	1861.9269	0.0001	1.0000
K2P	56	-785.6944	1867.2481	0.0000	1.0000
SYMIG	62	-769.8853	1867.3293	0.0000	1.0000
GTRI	64	-767.3608	1872.8467	0.0000	1.0000
GTRG	64	-767.7407	1873.6065	0.0000	1.0000
GTRIG	65	-767.3378	1878.0838	0.0000	1.0000
HKY	59	-783.2347	1878.1783	0.0000	1.0000
SYM	60	-783.3886	1883.7695	0.0000	1.0000
GTR	63	-780.7866	1894.4151	0.0000	1.0000
JC69I	56	-807.1931	1910.2455	0.0000	1.0000
JC69G	56	-807.4880	1910.8354	0.0000	1.0000
JC69IG	57	-807.4874	1916.1174	0.0000	1.0000
F81I	59	-805.4879	1922.6849	0.0000	1.0000
F81G	59	-805.8386	1923.3863	0.0000	1.0000
F81IG	60	-805.9466	1928.8854	0.0000	1.0000
JC69	55	-820.7857	1932.1476	0.0000	1.0000
F81	58	-819.0143	1944.4544	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 183-phylip/chr3_18148.phylip = 1 - 197;
 Partition Dummy = 1:183-phylip/chr3_18148.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 183-phylip/chr3_18148.phylip = 1 - 197;
 Partition Dummy = 1:183-phylip/chr3_18148.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 183-phylip/chr3_18148.phylip = 1 - 197;
 Partition Dummy = 1:183-phylip/chr3_18148.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

