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

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

Minimum AIC  model: HKYG
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	-679.7344	1533.2924	0.7256	0.7256
HKYI	60	-681.2492	1536.3219	0.1595	0.8852
HKYIG	61	-679.7375	1537.5046	0.0883	0.9735
HKY	59	-685.1496	1539.9781	0.0256	0.9991
GTRG	64	-678.5113	1548.0529	0.0005	0.9996
K2PG	57	-694.2016	1549.9715	0.0002	0.9998
GTRI	64	-680.0592	1551.1487	0.0001	0.9999
K2PI	57	-695.4210	1552.4104	0.0001	0.9999
GTRIG	65	-678.5117	1552.5196	0.0000	1.0000
K2PIG	58	-694.1985	1553.9911	0.0000	1.0000
GTR	63	-683.7236	1554.0788	0.0000	1.0000
K2P	56	-700.8830	1559.3661	0.0000	1.0000
SYMG	61	-692.9429	1563.9153	0.0000	1.0000
SYMI	61	-694.0932	1566.2160	0.0000	1.0000
SYMIG	62	-692.9410	1568.1806	0.0000	1.0000
SYM	60	-699.7372	1573.2978	0.0000	1.0000
F81G	59	-705.0187	1579.7163	0.0000	1.0000
F81I	59	-706.3172	1582.3132	0.0000	1.0000
F81IG	60	-705.0227	1583.8689	0.0000	1.0000
F81	58	-710.6242	1586.8427	0.0000	1.0000
JC69G	56	-716.7863	1591.1725	0.0000	1.0000
JC69I	56	-717.9720	1593.5439	0.0000	1.0000
JC69IG	57	-716.7849	1595.1381	0.0000	1.0000
JC69	55	-723.0035	1599.6949	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-679.7344	1479.4688	0.5889	0.5889
HKYIG	61	-679.7375	1481.4750	0.2160	0.8049
HKYI	60	-681.2492	1482.4983	0.1295	0.9344
GTRG	64	-678.5113	1485.0226	0.0367	0.9711
GTRIG	65	-678.5117	1487.0234	0.0135	0.9845
GTRI	64	-680.0592	1488.1184	0.0078	0.9923
HKY	59	-685.1496	1488.2993	0.0071	0.9994
GTR	63	-683.7236	1493.4472	0.0005	1.0000
K2PG	57	-694.2016	1502.4031	0.0000	1.0000
K2PIG	58	-694.1985	1504.3969	0.0000	1.0000
K2PI	57	-695.4210	1504.8421	0.0000	1.0000
SYMG	61	-692.9429	1507.8857	0.0000	1.0000
SYMIG	62	-692.9410	1509.8821	0.0000	1.0000
SYMI	61	-694.0932	1510.1863	0.0000	1.0000
K2P	56	-700.8830	1513.7661	0.0000	1.0000
SYM	60	-699.7372	1519.4743	0.0000	1.0000
F81G	59	-705.0187	1528.0374	0.0000	1.0000
F81IG	60	-705.0227	1530.0453	0.0000	1.0000
F81I	59	-706.3172	1530.6344	0.0000	1.0000
F81	58	-710.6242	1537.2485	0.0000	1.0000
JC69G	56	-716.7863	1545.5725	0.0000	1.0000
JC69IG	57	-716.7849	1547.5698	0.0000	1.0000
JC69I	56	-717.9720	1547.9439	0.0000	1.0000
JC69	55	-723.0035	1556.0070	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-679.7344	1676.4611	0.7378	0.7378
HKYI	60	-681.2492	1679.4906	0.1622	0.9000
HKYIG	61	-679.7375	1681.7504	0.0524	0.9524
HKY	59	-685.1496	1682.0083	0.0461	0.9984
K2PG	57	-694.2016	1689.5458	0.0011	0.9995
K2PI	57	-695.4210	1691.9847	0.0003	0.9998
K2PIG	58	-694.1985	1694.8227	0.0001	0.9999
GTRG	64	-678.5113	1695.1476	0.0001	1.0000
K2P	56	-700.8830	1697.6255	0.0000	1.0000
GTRI	64	-680.0592	1698.2435	0.0000	1.0000
GTR	63	-683.7236	1700.2891	0.0000	1.0000
GTRIG	65	-678.5117	1700.4317	0.0000	1.0000
SYMG	61	-692.9429	1708.1611	0.0000	1.0000
SYMI	61	-694.0932	1710.4618	0.0000	1.0000
SYMIG	62	-692.9410	1713.4407	0.0000	1.0000
SYM	60	-699.7372	1716.4665	0.0000	1.0000
F81G	59	-705.0187	1721.7465	0.0000	1.0000
F81I	59	-706.3172	1724.3434	0.0000	1.0000
F81IG	60	-705.0227	1727.0376	0.0000	1.0000
F81	58	-710.6242	1727.6743	0.0000	1.0000
JC69G	56	-716.7863	1729.4319	0.0000	1.0000
JC69I	56	-717.9720	1731.8033	0.0000	1.0000
JC69IG	57	-716.7849	1734.7124	0.0000	1.0000
JC69	55	-723.0035	1736.5832	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/chr1_32248.phylip = 1 - 197;
 Partition Dummy = 1:183-phylip/chr1_32248.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/chr1_32248.phylip = 1 - 197;
 Partition Dummy = 1:183-phylip/chr1_32248.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 183-phylip/chr1_32248.phylip = 1 - 197;
 Partition Dummy = 1:183-phylip/chr1_32248.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

