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

Input data from file "261-aggressive-phylip/2321_hg18.phylip" (29 taxa, 568 characters)

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


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	-2489.0010	5122.5428	0.6710	0.6710
GTRIG	65	-2488.6263	5124.3442	0.2726	0.9436
SYMIG	62	-2494.9173	5129.3038	0.0228	0.9665
SYMG	61	-2496.4431	5129.8348	0.0175	0.9840
HKYG	60	-2498.3673	5131.1724	0.0090	0.9929
HKYIG	61	-2497.9605	5132.8696	0.0038	0.9968
K2PG	57	-2503.5298	5134.0244	0.0022	0.9989
K2PIG	58	-2503.1392	5135.7243	0.0009	0.9999
GTRI	64	-2497.6792	5139.8992	0.0001	1.0000
SYMI	61	-2502.9361	5142.8209	0.0000	1.0000
K2PI	57	-2510.0302	5147.0251	0.0000	1.0000
HKYI	60	-2506.4882	5147.4144	0.0000	1.0000
GTR	63	-2536.6936	5215.3872	0.0000	1.0000
F81G	59	-2542.5894	5217.1157	0.0000	1.0000
SYM	60	-2541.7644	5217.9667	0.0000	1.0000
JC69G	56	-2546.9089	5218.3109	0.0000	1.0000
F81IG	60	-2542.2416	5218.9210	0.0000	1.0000
JC69IG	57	-2546.5772	5220.1191	0.0000	1.0000
HKY	59	-2545.7401	5223.4172	0.0000	1.0000
K2P	56	-2549.7177	5223.9286	0.0000	1.0000
F81I	59	-2548.6436	5229.2242	0.0000	1.0000
JC69I	56	-2553.2764	5231.0459	0.0000	1.0000
F81	58	-2587.9577	5305.3615	0.0000	1.0000
JC69	55	-2591.9644	5305.9601	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2489.0010	5106.0021	0.6349	0.6349
GTRIG	65	-2488.6263	5107.2525	0.3397	0.9746
SYMIG	62	-2494.9173	5113.8345	0.0126	0.9873
SYMG	61	-2496.4431	5114.8862	0.0075	0.9947
HKYG	60	-2498.3673	5116.7346	0.0030	0.9977
HKYIG	61	-2497.9605	5117.9210	0.0016	0.9994
K2PG	57	-2503.5298	5121.0597	0.0003	0.9997
K2PIG	58	-2503.1392	5122.2783	0.0002	0.9999
GTRI	64	-2497.6792	5123.3584	0.0001	1.0000
SYMI	61	-2502.9361	5127.8723	0.0000	1.0000
HKYI	60	-2506.4882	5132.9765	0.0000	1.0000
K2PI	57	-2510.0302	5134.0604	0.0000	1.0000
GTR	63	-2536.6936	5199.3872	0.0000	1.0000
F81G	59	-2542.5894	5203.1787	0.0000	1.0000
SYM	60	-2541.7644	5203.5288	0.0000	1.0000
F81IG	60	-2542.2416	5204.4831	0.0000	1.0000
JC69G	56	-2546.9089	5205.8177	0.0000	1.0000
JC69IG	57	-2546.5772	5207.1544	0.0000	1.0000
HKY	59	-2545.7401	5209.4802	0.0000	1.0000
K2P	56	-2549.7177	5211.4355	0.0000	1.0000
F81I	59	-2548.6436	5215.2872	0.0000	1.0000
JC69I	56	-2553.2764	5218.5528	0.0000	1.0000
F81	58	-2587.9577	5291.9155	0.0000	1.0000
JC69	55	-2591.9644	5293.9289	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-2503.5298	5368.5606	0.9241	0.9241
K2PIG	58	-2503.1392	5374.1214	0.0573	0.9814
HKYG	60	-2498.3673	5377.2619	0.0119	0.9933
SYMG	61	-2496.4431	5379.7556	0.0034	0.9967
K2PI	57	-2510.0302	5381.5613	0.0014	0.9981
HKYIG	61	-2497.9605	5382.7904	0.0008	0.9989
SYMIG	62	-2494.9173	5383.0460	0.0007	0.9995
GTRG	64	-2489.0010	5383.8978	0.0004	1.0000
GTRIG	65	-2488.6263	5389.4904	0.0000	1.0000
SYMI	61	-2502.9361	5392.7417	0.0000	1.0000
HKYI	60	-2506.4882	5393.5038	0.0000	1.0000
GTRI	64	-2497.6792	5401.2542	0.0000	1.0000
JC69G	56	-2546.9089	5448.9765	0.0000	1.0000
K2P	56	-2549.7177	5454.5943	0.0000	1.0000
JC69IG	57	-2546.5772	5454.6553	0.0000	1.0000
F81G	59	-2542.5894	5459.3639	0.0000	1.0000
JC69I	56	-2553.2764	5461.7116	0.0000	1.0000
SYM	60	-2541.7644	5464.0561	0.0000	1.0000
F81IG	60	-2542.2416	5465.0104	0.0000	1.0000
HKY	59	-2545.7401	5465.6654	0.0000	1.0000
F81I	59	-2548.6436	5471.4723	0.0000	1.0000
GTR	63	-2536.6936	5472.9408	0.0000	1.0000
JC69	55	-2591.9644	5532.7455	0.0000	1.0000
F81	58	-2587.9577	5543.7585	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/2321_hg18.phylip = 1 - 568;
 Partition Dummy = 1:261-aggressive-phylip/2321_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/2321_hg18.phylip = 1 - 568;
 Partition Dummy = 1:261-aggressive-phylip/2321_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 (K2PG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/2321_hg18.phylip = 1 - 568;
 Partition Dummy = 1:261-aggressive-phylip/2321_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

