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

Input data from file "183-phylip/chr9_393.phylip" (29 taxa, 226 characters)

Minimum AIC  model: GTRG
Minimum AICc model: HKYG
Minimum BIC  model: K2PG


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	60	-758.2807	1680.9250	0.5242	0.5242
HKYI	60	-759.1040	1682.5717	0.2301	0.7543
HKYIG	61	-757.9881	1684.0982	0.1073	0.8616
K2PG	57	-765.9212	1685.1996	0.0618	0.9234
K2PI	57	-766.6574	1686.6719	0.0296	0.9531
GTRG	64	-753.9122	1687.5015	0.0196	0.9726
K2PIG	58	-765.7004	1688.3828	0.0126	0.9852
GTRI	64	-754.7769	1689.2307	0.0082	0.9934
GTRIG	65	-753.7172	1691.0593	0.0033	0.9967
SYMG	61	-761.9750	1692.0719	0.0020	0.9987
SYMI	61	-762.7440	1693.6099	0.0009	0.9997
SYMIG	62	-761.8393	1695.6050	0.0003	1.0000
HKY	59	-776.4325	1713.5156	0.0000	1.0000
F81G	59	-776.6605	1713.9716	0.0000	1.0000
JC69G	56	-782.6689	1715.1130	0.0000	1.0000
F81I	59	-777.4059	1715.4624	0.0000	1.0000
K2P	56	-783.3299	1716.4350	0.0000	1.0000
JC69I	56	-783.3789	1716.5328	0.0000	1.0000
F81IG	60	-776.4685	1717.3006	0.0000	1.0000
JC69IG	57	-782.4973	1718.3518	0.0000	1.0000
GTR	63	-771.6691	1719.1160	0.0000	1.0000
SYM	60	-778.9027	1722.1690	0.0000	1.0000
F81	58	-793.6965	1744.3750	0.0000	1.0000
JC69	55	-799.5583	1745.3519	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-753.9122	1635.8245	0.3096	0.3096
HKYG	60	-758.2807	1636.5613	0.2142	0.5238
GTRIG	65	-753.7172	1637.4343	0.1384	0.6623
GTRI	64	-754.7769	1637.5537	0.1304	0.7927
HKYIG	61	-757.9881	1637.9763	0.1056	0.8983
HKYI	60	-759.1040	1638.2080	0.0940	0.9923
K2PG	57	-765.9212	1645.8425	0.0021	0.9944
SYMG	61	-761.9750	1645.9499	0.0020	0.9963
K2PI	57	-766.6574	1647.3147	0.0010	0.9973
K2PIG	58	-765.7004	1647.4008	0.0009	0.9983
SYMI	61	-762.7440	1647.4880	0.0009	0.9992
SYMIG	62	-761.8393	1647.6786	0.0008	1.0000
GTR	63	-771.6691	1669.3383	0.0000	1.0000
HKY	59	-776.4325	1670.8650	0.0000	1.0000
F81G	59	-776.6605	1671.3210	0.0000	1.0000
F81I	59	-777.4059	1672.8118	0.0000	1.0000
F81IG	60	-776.4685	1672.9370	0.0000	1.0000
JC69G	56	-782.6689	1677.3379	0.0000	1.0000
SYM	60	-778.9027	1677.8054	0.0000	1.0000
K2P	56	-783.3299	1678.6599	0.0000	1.0000
JC69I	56	-783.3789	1678.7577	0.0000	1.0000
JC69IG	57	-782.4973	1678.9947	0.0000	1.0000
F81	58	-793.6965	1703.3930	0.0000	1.0000
JC69	55	-799.5583	1709.1166	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-765.9212	1840.8130	0.3999	0.3999
HKYG	60	-758.2807	1841.7934	0.2449	0.6448
K2PI	57	-766.6574	1842.2852	0.1915	0.8363
HKYI	60	-759.1040	1843.4401	0.1075	0.9438
K2PIG	58	-765.7004	1845.7918	0.0332	0.9770
HKYIG	61	-757.9881	1846.6289	0.0218	0.9988
SYMG	61	-761.9750	1854.6026	0.0004	0.9992
GTRG	64	-753.9122	1854.7387	0.0004	0.9996
SYMI	61	-762.7440	1856.1406	0.0002	0.9998
GTRI	64	-754.7769	1856.4680	0.0002	0.9999
SYMIG	62	-761.8393	1859.7518	0.0000	1.0000
GTRIG	65	-753.7172	1859.7691	0.0000	1.0000
JC69G	56	-782.6689	1868.8879	0.0000	1.0000
K2P	56	-783.3299	1870.2098	0.0000	1.0000
JC69I	56	-783.3789	1870.3077	0.0000	1.0000
HKY	59	-776.4325	1872.6765	0.0000	1.0000
F81G	59	-776.6605	1873.1325	0.0000	1.0000
JC69IG	57	-782.4973	1873.9652	0.0000	1.0000
F81I	59	-777.4059	1874.6234	0.0000	1.0000
F81IG	60	-776.4685	1878.1691	0.0000	1.0000
SYM	60	-778.9027	1883.0375	0.0000	1.0000
GTR	63	-771.6691	1884.8320	0.0000	1.0000
JC69	55	-799.5583	1897.2461	0.0000	1.0000
F81	58	-793.6965	1901.7840	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr9_393.phylip = 1 - 226;
 Partition Dummy = 1:183-phylip/chr9_393.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr9_393.phylip = 1 - 226;
 Partition Dummy = 1:183-phylip/chr9_393.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 (K2PG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr9_393.phylip = 1 - 226;
 Partition Dummy = 1:183-phylip/chr9_393.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

