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

Input data from file "261-aggressive-phylip/5017_hg18.phylip" (29 taxa, 256 characters)

Minimum AIC  model: GTRI
Minimum AICc model: HKYG
Minimum BIC  model: HKYG


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	-769.7667	1697.0718	0.4545	0.4545
HKYI	60	-769.8870	1697.3125	0.4030	0.8574
HKYIG	61	-769.7685	1700.5267	0.0808	0.9382
GTRI	64	-765.4864	1702.5329	0.0296	0.9678
GTRG	64	-765.5682	1702.6967	0.0273	0.9951
GTRIG	65	-765.6063	1706.3705	0.0043	0.9995
HKY	59	-778.9601	1712.0426	0.0003	0.9997
K2PG	57	-783.0905	1713.5750	0.0001	0.9999
K2PI	57	-783.4183	1714.2306	0.0001	0.9999
GTR	63	-774.3150	1716.6300	0.0000	1.0000
K2PIG	58	-783.0656	1716.8723	0.0000	1.0000
SYMG	61	-780.0000	1720.9898	0.0000	1.0000
SYMI	61	-780.2256	1721.4409	0.0000	1.0000
SYMIG	62	-779.9992	1724.4751	0.0000	1.0000
F81I	59	-786.9148	1727.9520	0.0000	1.0000
F81G	59	-787.0303	1728.1830	0.0000	1.0000
K2P	56	-792.7053	1729.4910	0.0000	1.0000
F81IG	60	-787.0298	1731.5981	0.0000	1.0000
SYM	60	-788.3510	1734.2404	0.0000	1.0000
F81	58	-795.8572	1742.4555	0.0000	1.0000
JC69I	56	-799.2471	1742.5746	0.0000	1.0000
JC69G	56	-799.2917	1742.6638	0.0000	1.0000
JC69IG	57	-799.2912	1745.9763	0.0000	1.0000
JC69	55	-807.9651	1756.7302	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	64	-765.4864	1658.9727	0.2531	0.2531
GTRG	64	-765.5682	1659.1365	0.2332	0.4863
HKYG	60	-769.7667	1659.5334	0.1912	0.6775
HKYI	60	-769.8870	1659.7740	0.1695	0.8470
GTRIG	65	-765.6063	1661.2126	0.0826	0.9296
HKYIG	61	-769.7685	1661.5370	0.0702	0.9998
GTR	63	-774.3150	1674.6300	0.0001	0.9999
HKY	59	-778.9601	1675.9201	0.0001	1.0000
K2PG	57	-783.0905	1680.1810	0.0000	1.0000
K2PI	57	-783.4183	1680.8367	0.0000	1.0000
SYMG	61	-780.0000	1682.0001	0.0000	1.0000
K2PIG	58	-783.0656	1682.1312	0.0000	1.0000
SYMI	61	-780.2256	1682.4512	0.0000	1.0000
SYMIG	62	-779.9992	1683.9984	0.0000	1.0000
F81I	59	-786.9148	1691.8296	0.0000	1.0000
F81G	59	-787.0303	1692.0606	0.0000	1.0000
F81IG	60	-787.0298	1694.0596	0.0000	1.0000
SYM	60	-788.3510	1696.7019	0.0000	1.0000
K2P	56	-792.7053	1697.4106	0.0000	1.0000
F81	58	-795.8572	1707.7144	0.0000	1.0000
JC69I	56	-799.2471	1710.4942	0.0000	1.0000
JC69G	56	-799.2917	1710.5834	0.0000	1.0000
JC69IG	57	-799.2912	1712.5824	0.0000	1.0000
JC69	55	-807.9651	1725.9302	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-769.7667	1872.2440	0.5090	0.5090
HKYI	60	-769.8870	1872.4847	0.4513	0.9602
HKYIG	61	-769.7685	1877.7929	0.0318	0.9920
K2PG	57	-783.0905	1882.2561	0.0034	0.9954
K2PI	57	-783.4183	1882.9118	0.0025	0.9978
HKY	59	-778.9601	1885.0856	0.0008	0.9987
GTRI	64	-765.4864	1885.8641	0.0006	0.9992
GTRG	64	-765.5682	1886.0278	0.0005	0.9997
K2PIG	58	-783.0656	1887.7515	0.0002	1.0000
GTRIG	65	-765.6063	1891.6492	0.0000	1.0000
K2P	56	-792.7053	1895.9406	0.0000	1.0000
GTR	63	-774.3150	1897.9762	0.0000	1.0000
SYMG	61	-780.0000	1898.2559	0.0000	1.0000
SYMI	61	-780.2256	1898.7070	0.0000	1.0000
F81I	59	-786.9148	1900.9950	0.0000	1.0000
F81G	59	-787.0303	1901.2260	0.0000	1.0000
SYMIG	62	-779.9992	1903.7994	0.0000	1.0000
F81IG	60	-787.0298	1906.7702	0.0000	1.0000
JC69I	56	-799.2471	1909.0242	0.0000	1.0000
JC69G	56	-799.2917	1909.1134	0.0000	1.0000
SYM	60	-788.3510	1909.4125	0.0000	1.0000
F81	58	-795.8572	1913.3347	0.0000	1.0000
JC69IG	57	-799.2912	1914.6575	0.0000	1.0000
JC69	55	-807.9651	1920.9149	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/5017_hg18.phylip = 1 - 256;
 Partition Dummy = 1:261-aggressive-phylip/5017_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/5017_hg18.phylip = 1 - 256;
 Partition Dummy = 1:261-aggressive-phylip/5017_hg18.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 (HKYG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/5017_hg18.phylip = 1 - 256;
 Partition Dummy = 1:261-aggressive-phylip/5017_hg18.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;


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

