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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: GTRG


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	64	-3012.9331	6167.8025	0.6748	0.6748
GTRIG	65	-3012.4334	6169.2627	0.3252	1.0000
HKYIG	61	-3028.6683	6191.9432	0.0000	1.0000
HKYG	60	-3030.5903	6193.3602	0.0000	1.0000
GTRI	64	-3032.3160	6206.5684	0.0000	1.0000
SYMG	61	-3037.5120	6209.6306	0.0000	1.0000
SYMIG	62	-3037.1344	6211.3105	0.0000	1.0000
K2PIG	58	-3045.5041	6218.3580	0.0000	1.0000
K2PG	57	-3048.9677	6222.8825	0.0000	1.0000
HKYI	60	-3048.0846	6228.3489	0.0000	1.0000
SYMI	61	-3054.2751	6243.1568	0.0000	1.0000
K2PI	57	-3062.3309	6249.6089	0.0000	1.0000
F81G	59	-3099.8479	6329.4565	0.0000	1.0000
F81IG	60	-3099.0718	6330.3232	0.0000	1.0000
JC69G	56	-3111.3727	6345.2974	0.0000	1.0000
JC69IG	57	-3110.5730	6346.0929	0.0000	1.0000
GTR	63	-3106.8238	6353.1325	0.0000	1.0000
F81I	59	-3119.7685	6369.2979	0.0000	1.0000
HKY	59	-3120.6872	6371.1353	0.0000	1.0000
JC69I	56	-3131.4618	6385.4756	0.0000	1.0000
SYM	60	-3127.6358	6387.4512	0.0000	1.0000
K2P	56	-3135.8529	6394.2579	0.0000	1.0000
F81	58	-3190.9374	6509.2247	0.0000	1.0000
JC69	55	-3202.5610	6525.2870	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-3012.9331	6153.8662	0.6225	0.6225
GTRIG	65	-3012.4334	6154.8667	0.3775	1.0000
HKYIG	61	-3028.6683	6179.3366	0.0000	1.0000
HKYG	60	-3030.5903	6181.1805	0.0000	1.0000
GTRI	64	-3032.3160	6192.6321	0.0000	1.0000
SYMG	61	-3037.5120	6197.0240	0.0000	1.0000
SYMIG	62	-3037.1344	6198.2687	0.0000	1.0000
K2PIG	58	-3045.5041	6207.0081	0.0000	1.0000
K2PG	57	-3048.9677	6211.9355	0.0000	1.0000
HKYI	60	-3048.0846	6216.1692	0.0000	1.0000
SYMI	61	-3054.2751	6230.5501	0.0000	1.0000
K2PI	57	-3062.3309	6238.6619	0.0000	1.0000
F81G	59	-3099.8479	6317.6957	0.0000	1.0000
F81IG	60	-3099.0718	6318.1435	0.0000	1.0000
JC69G	56	-3111.3727	6334.7453	0.0000	1.0000
JC69IG	57	-3110.5730	6335.1459	0.0000	1.0000
GTR	63	-3106.8238	6339.6475	0.0000	1.0000
F81I	59	-3119.7685	6357.5371	0.0000	1.0000
HKY	59	-3120.6872	6359.3745	0.0000	1.0000
JC69I	56	-3131.4618	6374.9235	0.0000	1.0000
SYM	60	-3127.6358	6375.2715	0.0000	1.0000
K2P	56	-3135.8529	6383.7058	0.0000	1.0000
F81	58	-3190.9374	6497.8748	0.0000	1.0000
JC69	55	-3202.5610	6515.1220	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-3012.9331	6441.5632	0.9294	0.9294
GTRIG	65	-3012.4334	6447.0590	0.0595	0.9889
HKYG	60	-3030.5903	6450.8965	0.0087	0.9977
HKYIG	61	-3028.6683	6453.5478	0.0023	1.0000
K2PIG	58	-3045.5041	6467.7335	0.0000	1.0000
K2PG	57	-3048.9677	6468.1656	0.0000	1.0000
SYMG	61	-3037.5120	6471.2352	0.0000	1.0000
SYMIG	62	-3037.1344	6476.9752	0.0000	1.0000
GTRI	64	-3032.3160	6480.3291	0.0000	1.0000
HKYI	60	-3048.0846	6485.8852	0.0000	1.0000
K2PI	57	-3062.3309	6494.8920	0.0000	1.0000
SYMI	61	-3054.2751	6504.7613	0.0000	1.0000
F81G	59	-3099.8479	6582.9164	0.0000	1.0000
JC69G	56	-3111.3727	6586.4802	0.0000	1.0000
F81IG	60	-3099.0718	6587.8595	0.0000	1.0000
JC69IG	57	-3110.5730	6591.3760	0.0000	1.0000
F81I	59	-3119.7685	6622.7577	0.0000	1.0000
GTR	63	-3106.8238	6622.8493	0.0000	1.0000
HKY	59	-3120.6872	6624.5951	0.0000	1.0000
JC69I	56	-3131.4618	6626.6584	0.0000	1.0000
K2P	56	-3135.8529	6635.4407	0.0000	1.0000
SYM	60	-3127.6358	6644.9875	0.0000	1.0000
F81	58	-3190.9374	6758.6002	0.0000	1.0000
JC69	55	-3202.5610	6762.3616	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_6343.phylip = 1 - 662;
 Partition Dummy = 1:183-phylip/chr9_6343.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 183-phylip/chr9_6343.phylip = 1 - 662;
 Partition Dummy = 1:183-phylip/chr9_6343.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 (GTRG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr9_6343.phylip = 1 - 662;
 Partition Dummy = 1:183-phylip/chr9_6343.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;


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

