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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: K2P


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	-761.8624	1675.9716	0.3497	0.3497
HKY	59	-763.7045	1676.4616	0.2737	0.6234
HKYI	60	-762.1565	1676.5597	0.2606	0.8840
HKYIG	61	-761.9099	1679.2888	0.0666	0.9506
K2P	56	-771.2678	1682.1719	0.0158	0.9663
K2PG	57	-769.8071	1682.3619	0.0143	0.9806
K2PI	57	-770.0791	1682.9061	0.0109	0.9915
K2PIG	58	-769.8118	1685.5100	0.0030	0.9945
GTR	63	-762.3150	1686.6300	0.0017	0.9962
GTRG	64	-760.7294	1686.7683	0.0016	0.9978
GTRI	64	-761.0750	1687.4594	0.0011	0.9989
SYM	60	-768.8093	1689.8653	0.0003	0.9992
GTRIG	65	-760.8683	1690.3852	0.0003	0.9995
SYMG	61	-767.4812	1690.4314	0.0003	0.9998
SYMI	61	-767.7632	1690.9954	0.0002	1.0000
SYMIG	62	-767.5558	1693.8315	0.0000	1.0000
F81	58	-786.4259	1718.7383	0.0000	1.0000
F81G	59	-785.0971	1719.2468	0.0000	1.0000
F81I	59	-785.3597	1719.7720	0.0000	1.0000
JC69	55	-792.7664	1722.0846	0.0000	1.0000
F81IG	60	-785.1068	1722.4603	0.0000	1.0000
JC69G	56	-791.5376	1722.7115	0.0000	1.0000
JC69I	56	-791.8004	1723.2372	0.0000	1.0000
JC69IG	57	-791.5446	1725.8369	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-761.8624	1643.7249	0.3707	0.3707
HKYI	60	-762.1565	1644.3130	0.2763	0.6470
HKY	59	-763.7045	1645.4090	0.1597	0.8068
HKYIG	61	-761.9099	1645.8197	0.1301	0.9368
GTRG	64	-760.7294	1649.4589	0.0211	0.9579
GTRI	64	-761.0750	1650.1500	0.0149	0.9728
GTR	63	-762.3150	1650.6300	0.0117	0.9846
GTRIG	65	-760.8683	1651.7365	0.0068	0.9913
K2PG	57	-769.8071	1653.6141	0.0026	0.9940
K2PI	57	-770.0791	1654.1582	0.0020	0.9960
K2P	56	-771.2678	1654.5356	0.0017	0.9976
K2PIG	58	-769.8118	1655.6235	0.0010	0.9986
SYMG	61	-767.4812	1656.9623	0.0005	0.9991
SYMI	61	-767.7632	1657.5263	0.0004	0.9995
SYM	60	-768.8093	1657.6186	0.0004	0.9998
SYMIG	62	-767.5558	1659.1115	0.0002	1.0000
F81G	59	-785.0971	1688.1941	0.0000	1.0000
F81I	59	-785.3597	1688.7193	0.0000	1.0000
F81	58	-786.4259	1688.8518	0.0000	1.0000
F81IG	60	-785.1068	1690.2136	0.0000	1.0000
JC69G	56	-791.5376	1695.0751	0.0000	1.0000
JC69	55	-792.7664	1695.5329	0.0000	1.0000
JC69I	56	-791.8004	1695.6009	0.0000	1.0000
JC69IG	57	-791.5446	1697.0891	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2P	56	-771.2678	1859.6614	0.4716	0.4716
HKY	59	-763.7045	1861.5236	0.1859	0.6574
K2PG	57	-769.8071	1862.4028	0.1197	0.7772
K2PI	57	-770.0791	1862.9470	0.0912	0.8684
HKYG	60	-761.8624	1863.5025	0.0691	0.9375
HKYI	60	-762.1565	1864.0907	0.0515	0.9890
K2PIG	58	-769.8118	1868.0752	0.0070	0.9960
HKYIG	61	-761.9099	1869.2603	0.0039	0.9999
SYM	60	-768.8093	1877.3963	0.0001	1.0000
SYMG	61	-767.4812	1880.4029	0.0000	1.0000
SYMI	61	-767.7632	1880.9669	0.0000	1.0000
GTR	63	-762.3150	1881.3965	0.0000	1.0000
GTRG	64	-760.7294	1883.8884	0.0000	1.0000
GTRI	64	-761.0750	1884.5794	0.0000	1.0000
SYMIG	62	-767.5558	1886.2151	0.0000	1.0000
GTRIG	65	-760.8683	1889.8290	0.0000	1.0000
JC69	55	-792.7664	1896.9957	0.0000	1.0000
JC69G	56	-791.5376	1900.2009	0.0000	1.0000
JC69I	56	-791.8004	1900.7267	0.0000	1.0000
F81	58	-786.4259	1901.3035	0.0000	1.0000
F81G	59	-785.0971	1904.3088	0.0000	1.0000
F81I	59	-785.3597	1904.8340	0.0000	1.0000
JC69IG	57	-791.5446	1905.8779	0.0000	1.0000
F81IG	60	-785.1068	1909.9912	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 183-phylip/chr7_10443.phylip = 1 - 288;
 Partition Dummy = 1:183-phylip/chr7_10443.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 183-phylip/chr7_10443.phylip = 1 - 288;
 Partition Dummy = 1:183-phylip/chr7_10443.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 (K2P)]

BEGIN MRBAYES;
 Charset 183-phylip/chr7_10443.phylip = 1 - 288;
 Partition Dummy = 1:183-phylip/chr7_10443.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0);
END;


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

