
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
21:50 04/04/2011

Input data from file "591-tomodel-phylip/chr6_3274.phylip" (7 taxa, 743 characters)

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	20	-3091.6000	6224.3634	0.6644	0.6644
GTRIG	21	-3091.6083	6226.4981	0.2285	0.8929
HKYG	16	-3097.9390	6228.6273	0.0788	0.9717
HKYIG	17	-3097.9386	6230.7214	0.0277	0.9993
GTRI	20	-3098.7184	6238.6002	0.0005	0.9999
HKYI	16	-3104.7929	6242.3352	0.0001	1.0000
K2PG	13	-3109.6100	6245.7194	0.0000	1.0000
SYMG	17	-3105.6691	6246.1824	0.0000	1.0000
K2PIG	14	-3109.6097	6247.7963	0.0000	1.0000
SYMIG	18	-3105.6691	6248.2829	0.0000	1.0000
K2PI	13	-3116.4183	6259.3360	0.0000	1.0000
SYMI	17	-3112.2731	6259.3903	0.0000	1.0000
GTR	19	-3143.3153	6325.6818	0.0000	1.0000
HKY	15	-3150.0489	6330.7580	0.0000	1.0000
K2P	12	-3160.8139	6346.0551	0.0000	1.0000
SYM	16	-3156.9108	6346.5708	0.0000	1.0000
F81G	15	-3168.6643	6367.9888	0.0000	1.0000
F81IG	16	-3168.6704	6370.0902	0.0000	1.0000
F81I	15	-3172.9497	6376.5597	0.0000	1.0000
JC69G	12	-3182.7235	6389.8745	0.0000	1.0000
JC69IG	13	-3182.7253	6391.9499	0.0000	1.0000
JC69I	12	-3186.7878	6398.0031	0.0000	1.0000
F81	14	-3211.9293	6452.4355	0.0000	1.0000
JC69	11	-3225.9246	6474.2104	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-3091.6000	6223.1999	0.6677	0.6677
GTRIG	21	-3091.6083	6225.2165	0.2436	0.9113
HKYG	16	-3097.9390	6227.8780	0.0644	0.9757
HKYIG	17	-3097.9386	6229.8772	0.0237	0.9994
GTRI	20	-3098.7184	6237.4367	0.0005	0.9999
HKYI	16	-3104.7929	6241.5858	0.0001	1.0000
K2PG	13	-3109.6100	6245.2200	0.0000	1.0000
SYMG	17	-3105.6691	6245.3383	0.0000	1.0000
K2PIG	14	-3109.6097	6247.2194	0.0000	1.0000
SYMIG	18	-3105.6691	6247.3381	0.0000	1.0000
SYMI	17	-3112.2731	6258.5461	0.0000	1.0000
K2PI	13	-3116.4183	6258.8367	0.0000	1.0000
GTR	19	-3143.3153	6324.6306	0.0000	1.0000
HKY	15	-3150.0489	6330.0977	0.0000	1.0000
K2P	12	-3160.8139	6345.6277	0.0000	1.0000
SYM	16	-3156.9108	6345.8215	0.0000	1.0000
F81G	15	-3168.6643	6367.3285	0.0000	1.0000
F81IG	16	-3168.6704	6369.3409	0.0000	1.0000
F81I	15	-3172.9497	6375.8994	0.0000	1.0000
JC69G	12	-3182.7235	6389.4471	0.0000	1.0000
JC69IG	13	-3182.7253	6391.4506	0.0000	1.0000
JC69I	12	-3186.7878	6397.5757	0.0000	1.0000
F81	14	-3211.9293	6451.8586	0.0000	1.0000
JC69	11	-3225.9246	6473.8492	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3097.9390	6301.6491	0.8208	0.8208
K2PG	13	-3109.6100	6305.1591	0.1419	0.9628
HKYIG	17	-3097.9386	6308.2591	0.0301	0.9929
K2PIG	14	-3109.6097	6311.7692	0.0052	0.9981
HKYI	16	-3104.7929	6315.3570	0.0009	0.9990
GTRG	20	-3091.6000	6315.4139	0.0008	0.9998
K2PI	13	-3116.4183	6318.7757	0.0002	1.0000
GTRIG	21	-3091.6083	6322.0411	0.0000	1.0000
SYMG	17	-3105.6691	6323.7201	0.0000	1.0000
GTRI	20	-3098.7184	6329.6507	0.0000	1.0000
SYMIG	18	-3105.6691	6330.3307	0.0000	1.0000
SYMI	17	-3112.2731	6336.9280	0.0000	1.0000
HKY	15	-3150.0489	6399.2582	0.0000	1.0000
K2P	12	-3160.8139	6400.9561	0.0000	1.0000
GTR	19	-3143.3153	6412.2338	0.0000	1.0000
SYM	16	-3156.9108	6419.5927	0.0000	1.0000
F81G	15	-3168.6643	6436.4890	0.0000	1.0000
F81IG	16	-3168.6704	6443.1120	0.0000	1.0000
JC69G	12	-3182.7235	6444.7754	0.0000	1.0000
F81I	15	-3172.9497	6445.0599	0.0000	1.0000
JC69IG	13	-3182.7253	6451.3896	0.0000	1.0000
JC69I	12	-3186.7878	6452.9040	0.0000	1.0000
F81	14	-3211.9293	6516.4084	0.0000	1.0000
JC69	11	-3225.9246	6524.5669	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 591-tomodel-phylip/chr6_3274.phylip = 1 - 743;
 Partition Dummy = 1:591-tomodel-phylip/chr6_3274.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 (GTRG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr6_3274.phylip = 1 - 743;
 Partition Dummy = 1:591-tomodel-phylip/chr6_3274.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 BIC model (HKYG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr6_3274.phylip = 1 - 743;
 Partition Dummy = 1:591-tomodel-phylip/chr6_3274.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

