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

Input data from file "261-aggressive-phylip/10566_hg18.phylip" (29 taxa, 782 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	64	-3170.7409	6481.0857	0.8586	0.8586
GTRIG	65	-3171.3651	6484.7134	0.1400	0.9986
HKYG	60	-3181.9126	6493.9777	0.0014	1.0000
GTRI	64	-3181.3126	6502.2291	0.0000	1.0000
HKYIG	61	-3184.8871	6502.2798	0.0000	1.0000
HKYI	60	-3190.4171	6510.9867	0.0000	1.0000
SYMG	61	-3197.5166	6527.5387	0.0000	1.0000
SYMIG	62	-3197.7002	6530.2655	0.0000	1.0000
K2PG	57	-3206.8584	6536.8494	0.0000	1.0000
K2PIG	58	-3205.7851	6537.0363	0.0000	1.0000
SYMI	61	-3205.7537	6544.0129	0.0000	1.0000
K2PI	57	-3213.8685	6550.8695	0.0000	1.0000
GTR	63	-3222.5188	6582.2687	0.0000	1.0000
HKY	59	-3230.3545	6588.5151	0.0000	1.0000
F81G	59	-3242.1291	6612.0643	0.0000	1.0000
F81IG	60	-3242.6813	6615.5152	0.0000	1.0000
SYM	60	-3244.9794	6620.1113	0.0000	1.0000
F81I	59	-3250.3517	6628.5096	0.0000	1.0000
K2P	56	-3254.8217	6630.4489	0.0000	1.0000
JC69G	56	-3256.3880	6633.5815	0.0000	1.0000
JC69IG	57	-3256.6325	6636.3976	0.0000	1.0000
JC69I	56	-3264.1053	6649.0162	0.0000	1.0000
F81	58	-3289.4628	6704.3917	0.0000	1.0000
JC69	55	-3303.3810	6725.2469	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-3170.7409	6469.4818	0.8348	0.8348
GTRIG	65	-3171.3651	6472.7302	0.1645	0.9993
HKYG	60	-3181.9126	6483.8251	0.0006	1.0000
GTRI	64	-3181.3126	6490.6252	0.0000	1.0000
HKYIG	61	-3184.8871	6491.7743	0.0000	1.0000
HKYI	60	-3190.4171	6500.8342	0.0000	1.0000
SYMG	61	-3197.5166	6517.0331	0.0000	1.0000
SYMIG	62	-3197.7002	6519.4004	0.0000	1.0000
K2PIG	58	-3205.7851	6527.5702	0.0000	1.0000
K2PG	57	-3206.8584	6527.7168	0.0000	1.0000
SYMI	61	-3205.7537	6533.5073	0.0000	1.0000
K2PI	57	-3213.8685	6541.7369	0.0000	1.0000
GTR	63	-3222.5188	6571.0375	0.0000	1.0000
HKY	59	-3230.3545	6578.7090	0.0000	1.0000
F81G	59	-3242.1291	6602.2582	0.0000	1.0000
F81IG	60	-3242.6813	6605.3626	0.0000	1.0000
SYM	60	-3244.9794	6609.9587	0.0000	1.0000
F81I	59	-3250.3517	6618.7035	0.0000	1.0000
K2P	56	-3254.8217	6621.6433	0.0000	1.0000
JC69G	56	-3256.3880	6624.7760	0.0000	1.0000
JC69IG	57	-3256.6325	6627.2650	0.0000	1.0000
JC69I	56	-3264.1053	6640.2107	0.0000	1.0000
F81	58	-3289.4628	6694.9256	0.0000	1.0000
JC69	55	-3303.3810	6716.7621	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-3181.9126	6763.5364	0.8925	0.8925
GTRG	64	-3170.7409	6767.8405	0.1037	0.9962
GTRIG	65	-3171.3651	6775.7507	0.0020	0.9982
HKYIG	61	-3184.8871	6776.1474	0.0016	0.9998
HKYI	60	-3190.4171	6780.5455	0.0002	1.0000
GTRI	64	-3181.3126	6788.9839	0.0000	1.0000
K2PG	57	-3206.8584	6793.4425	0.0000	1.0000
K2PIG	58	-3205.7851	6797.9578	0.0000	1.0000
SYMG	61	-3197.5166	6801.4063	0.0000	1.0000
K2PI	57	-3213.8685	6807.4627	0.0000	1.0000
SYMIG	62	-3197.7002	6808.4354	0.0000	1.0000
SYMI	61	-3205.7537	6817.8805	0.0000	1.0000
HKY	59	-3230.3545	6853.7584	0.0000	1.0000
GTR	63	-3222.5188	6864.7343	0.0000	1.0000
F81G	59	-3242.1291	6877.3076	0.0000	1.0000
K2P	56	-3254.8217	6882.7072	0.0000	1.0000
F81IG	60	-3242.6813	6885.0739	0.0000	1.0000
JC69G	56	-3256.3880	6885.8398	0.0000	1.0000
SYM	60	-3244.9794	6889.6700	0.0000	1.0000
JC69IG	57	-3256.6325	6892.9907	0.0000	1.0000
F81I	59	-3250.3517	6893.7529	0.0000	1.0000
JC69I	56	-3264.1053	6901.2746	0.0000	1.0000
F81	58	-3289.4628	6965.3132	0.0000	1.0000
JC69	55	-3303.3810	6973.1641	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 261-aggressive-phylip/10566_hg18.phylip = 1 - 782;
 Partition Dummy = 1:261-aggressive-phylip/10566_hg18.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 261-aggressive-phylip/10566_hg18.phylip = 1 - 782;
 Partition Dummy = 1:261-aggressive-phylip/10566_hg18.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 261-aggressive-phylip/10566_hg18.phylip = 1 - 782;
 Partition Dummy = 1:261-aggressive-phylip/10566_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

