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

Input data from file "261-aggressive-phylip/10381_hg18.phylip" (29 taxa, 1049 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	-4441.1009	9018.6571	0.9539	0.9539
GTRIG	65	-4443.7476	9026.2236	0.0217	0.9756
SYMG	61	-4448.4375	9026.5386	0.0185	0.9941
SYMIG	62	-4448.4627	9028.8484	0.0058	1.0000
K2PG	57	-4459.5296	9039.7312	0.0000	1.0000
HKYIG	61	-4455.4253	9040.5142	0.0000	1.0000
K2PIG	58	-4460.5676	9044.0484	0.0000	1.0000
HKYG	60	-4461.4430	9050.2948	0.0000	1.0000
GTRI	64	-4545.2607	9226.9768	0.0000	1.0000
SYMI	61	-4557.3719	9244.4074	0.0000	1.0000
GTR	63	-4557.9829	9250.1527	0.0000	1.0000
SYM	60	-4561.6422	9250.6932	0.0000	1.0000
HKYI	60	-4562.0553	9251.5196	0.0000	1.0000
K2PI	57	-4568.3996	9257.4712	0.0000	1.0000
K2P	56	-4572.7552	9263.9459	0.0000	1.0000
HKY	59	-4574.7262	9274.6112	0.0000	1.0000
JC69G	56	-4611.1291	9340.6938	0.0000	1.0000
JC69IG	57	-4616.2566	9353.1853	0.0000	1.0000
F81IG	60	-4615.3146	9358.0380	0.0000	1.0000
F81G	59	-4616.8941	9358.9470	0.0000	1.0000
JC69I	56	-4714.4659	9547.3672	0.0000	1.0000
JC69	55	-4717.7481	9551.6997	0.0000	1.0000
F81I	59	-4721.7825	9568.7238	0.0000	1.0000
F81	58	-4725.5382	9573.9896	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-4441.1009	9010.2018	0.9580	0.9580
GTRIG	65	-4443.7476	9017.4952	0.0250	0.9829
SYMG	61	-4448.4375	9018.8749	0.0125	0.9955
SYMIG	62	-4448.4627	9020.9254	0.0045	1.0000
HKYIG	61	-4455.4253	9032.8506	0.0000	1.0000
K2PG	57	-4459.5296	9033.0591	0.0000	1.0000
K2PIG	58	-4460.5676	9037.1352	0.0000	1.0000
HKYG	60	-4461.4430	9042.8859	0.0000	1.0000
GTRI	64	-4545.2607	9218.5215	0.0000	1.0000
SYMI	61	-4557.3719	9236.7438	0.0000	1.0000
GTR	63	-4557.9829	9241.9659	0.0000	1.0000
SYM	60	-4561.6422	9243.2843	0.0000	1.0000
HKYI	60	-4562.0553	9244.1107	0.0000	1.0000
K2PI	57	-4568.3996	9250.7991	0.0000	1.0000
K2P	56	-4572.7552	9257.5104	0.0000	1.0000
HKY	59	-4574.7262	9267.4525	0.0000	1.0000
JC69G	56	-4611.1291	9334.2583	0.0000	1.0000
JC69IG	57	-4616.2566	9346.5132	0.0000	1.0000
F81IG	60	-4615.3146	9350.6291	0.0000	1.0000
F81G	59	-4616.8941	9351.7882	0.0000	1.0000
JC69I	56	-4714.4659	9540.9318	0.0000	1.0000
JC69	55	-4717.7481	9545.4963	0.0000	1.0000
F81I	59	-4721.7825	9561.5651	0.0000	1.0000
F81	58	-4725.5382	9567.0765	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-4459.5296	9315.5279	0.9301	0.9301
SYMG	61	-4448.4375	9321.1661	0.0555	0.9856
K2PIG	58	-4460.5676	9324.5596	0.0102	0.9958
GTRG	64	-4441.1009	9327.3597	0.0025	0.9983
SYMIG	62	-4448.4627	9328.1722	0.0017	0.9999
HKYIG	61	-4455.4253	9335.1418	0.0001	1.0000
GTRIG	65	-4443.7476	9339.6087	0.0000	1.0000
HKYG	60	-4461.4430	9340.2215	0.0000	1.0000
K2PI	57	-4568.3996	9533.2679	0.0000	1.0000
K2P	56	-4572.7552	9535.0236	0.0000	1.0000
GTRI	64	-4545.2607	9535.6794	0.0000	1.0000
SYMI	61	-4557.3719	9539.0349	0.0000	1.0000
SYM	60	-4561.6422	9540.6199	0.0000	1.0000
HKYI	60	-4562.0553	9541.4462	0.0000	1.0000
GTR	63	-4557.9829	9554.1682	0.0000	1.0000
HKY	59	-4574.7262	9559.8324	0.0000	1.0000
JC69G	56	-4611.1291	9611.7715	0.0000	1.0000
JC69IG	57	-4616.2566	9628.9820	0.0000	1.0000
F81G	59	-4616.8941	9644.1682	0.0000	1.0000
F81IG	60	-4615.3146	9647.9647	0.0000	1.0000
JC69	55	-4717.7481	9818.0539	0.0000	1.0000
JC69I	56	-4714.4659	9818.4449	0.0000	1.0000
F81I	59	-4721.7825	9853.9450	0.0000	1.0000
F81	58	-4725.5382	9854.5008	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/10381_hg18.phylip = 1 - 1049;
 Partition Dummy = 1:261-aggressive-phylip/10381_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/10381_hg18.phylip = 1 - 1049;
 Partition Dummy = 1:261-aggressive-phylip/10381_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/10381_hg18.phylip = 1 - 1049;
 Partition Dummy = 1:261-aggressive-phylip/10381_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

