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

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

Minimum AIC  model: HKYIG
Minimum AICc model: HKYIG
Minimum BIC  model: HKYIG


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
HKYIG	61	-4658.9039	9449.1345	0.8381	0.8381
GTRIG	65	-4655.8956	9452.4233	0.1619	0.9999
HKYG	60	-4669.7326	9468.4800	0.0001	1.0000
GTRG	64	-4666.6043	9471.5057	0.0000	1.0000
K2PG	57	-4691.1816	9504.4760	0.0000	1.0000
HKYI	60	-4688.5999	9506.2147	0.0000	1.0000
K2PIG	58	-4691.6526	9507.7131	0.0000	1.0000
GTRI	64	-4686.3755	9511.0480	0.0000	1.0000
SYMG	61	-4694.3031	9519.9330	0.0000	1.0000
SYMIG	62	-4694.3071	9522.2587	0.0000	1.0000
SYMI	61	-4714.6579	9560.6425	0.0000	1.0000
K2PI	57	-4722.7307	9567.5744	0.0000	1.0000
F81IG	60	-4732.2884	9593.5916	0.0000	1.0000
F81G	59	-4744.0394	9614.7872	0.0000	1.0000
JC69IG	57	-4752.2690	9626.6510	0.0000	1.0000
HKY	59	-4754.2781	9635.2646	0.0000	1.0000
GTR	63	-4753.0243	9642.0164	0.0000	1.0000
JC69G	56	-4764.5789	9648.9814	0.0000	1.0000
F81I	59	-4761.6167	9649.9419	0.0000	1.0000
K2P	56	-4769.4965	9658.8165	0.0000	1.0000
SYM	60	-4777.1591	9683.3330	0.0000	1.0000
JC69I	56	-4783.0385	9685.9005	0.0000	1.0000
F81	58	-4824.9818	9774.3715	0.0000	1.0000
JC69	55	-4845.0281	9807.5961	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYIG	61	-4658.9039	9439.8077	0.7294	0.7294
GTRIG	65	-4655.8956	9441.7913	0.2705	0.9999
HKYG	60	-4669.7326	9459.4652	0.0000	1.0000
GTRG	64	-4666.6043	9461.2087	0.0000	1.0000
K2PG	57	-4691.1816	9496.3631	0.0000	1.0000
HKYI	60	-4688.5999	9497.1999	0.0000	1.0000
K2PIG	58	-4691.6526	9499.3053	0.0000	1.0000
GTRI	64	-4686.3755	9500.7510	0.0000	1.0000
SYMG	61	-4694.3031	9510.6063	0.0000	1.0000
SYMIG	62	-4694.3071	9512.6143	0.0000	1.0000
SYMI	61	-4714.6579	9551.3157	0.0000	1.0000
K2PI	57	-4722.7307	9559.4615	0.0000	1.0000
F81IG	60	-4732.2884	9584.5768	0.0000	1.0000
F81G	59	-4744.0394	9606.0788	0.0000	1.0000
JC69IG	57	-4752.2690	9618.5381	0.0000	1.0000
HKY	59	-4754.2781	9626.5561	0.0000	1.0000
GTR	63	-4753.0243	9632.0485	0.0000	1.0000
JC69G	56	-4764.5789	9641.1579	0.0000	1.0000
F81I	59	-4761.6167	9641.2334	0.0000	1.0000
K2P	56	-4769.4965	9650.9929	0.0000	1.0000
SYM	60	-4777.1591	9674.3182	0.0000	1.0000
JC69I	56	-4783.0385	9678.0769	0.0000	1.0000
F81	58	-4824.9818	9765.9637	0.0000	1.0000
JC69	55	-4845.0281	9800.0563	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYIG	61	-4658.9039	9730.8958	0.9994	0.9994
HKYG	60	-4669.7326	9745.7814	0.0006	1.0000
GTRIG	65	-4655.8956	9751.9671	0.0000	1.0000
GTRG	64	-4666.6043	9766.6126	0.0000	1.0000
K2PG	57	-4691.1816	9768.3634	0.0000	1.0000
K2PIG	58	-4691.6526	9776.0775	0.0000	1.0000
HKYI	60	-4688.5999	9783.5160	0.0000	1.0000
SYMG	61	-4694.3031	9801.6943	0.0000	1.0000
GTRI	64	-4686.3755	9806.1548	0.0000	1.0000
SYMIG	62	-4694.3071	9808.4743	0.0000	1.0000
K2PI	57	-4722.7307	9831.4618	0.0000	1.0000
SYMI	61	-4714.6579	9842.4038	0.0000	1.0000
F81IG	60	-4732.2884	9870.8929	0.0000	1.0000
F81G	59	-4744.0394	9887.6230	0.0000	1.0000
JC69IG	57	-4752.2690	9890.5384	0.0000	1.0000
HKY	59	-4754.2781	9908.1003	0.0000	1.0000
JC69G	56	-4764.5789	9908.3863	0.0000	1.0000
K2P	56	-4769.4965	9918.2213	0.0000	1.0000
F81I	59	-4761.6167	9922.7776	0.0000	1.0000
GTR	63	-4753.0243	9932.6805	0.0000	1.0000
JC69I	56	-4783.0385	9945.3053	0.0000	1.0000
SYM	60	-4777.1591	9960.6343	0.0000	1.0000
F81	58	-4824.9818	10042.7359	0.0000	1.0000
JC69	55	-4845.0281	10062.5127	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (HKYIG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/7854_hg18.phylip = 1 - 873;
 Partition Dummy = 1:261-aggressive-phylip/7854_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (HKYIG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/7854_hg18.phylip = 1 - 873;
 Partition Dummy = 1:261-aggressive-phylip/7854_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best BIC model (HKYIG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/7854_hg18.phylip = 1 - 873;
 Partition Dummy = 1:261-aggressive-phylip/7854_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

