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

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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
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
GTRIG	65	-4662.5135	9464.4453	0.9072	0.9072
GTRG	64	-4666.3838	9469.8904	0.0596	0.9668
HKYIG	61	-4670.4796	9471.2258	0.0306	0.9974
HKYG	60	-4674.0678	9476.1268	0.0026	1.0000
K2PIG	58	-4696.0472	9515.5498	0.0000	1.0000
GTRI	64	-4690.6364	9518.3957	0.0000	1.0000
K2PG	57	-4699.4607	9520.1162	0.0000	1.0000
HKYI	60	-4696.2562	9520.5037	0.0000	1.0000
SYMIG	62	-4694.2718	9521.0906	0.0000	1.0000
SYMG	61	-4696.5878	9523.4423	0.0000	1.0000
K2PI	57	-4718.3490	9557.8928	0.0000	1.0000
SYMI	61	-4717.0082	9564.2831	0.0000	1.0000
F81IG	60	-4825.4204	9778.8320	0.0000	1.0000
F81G	59	-4826.6372	9778.9953	0.0000	1.0000
JC69G	56	-4843.6210	9806.1811	0.0000	1.0000
JC69IG	57	-4843.0491	9807.2930	0.0000	1.0000
F81I	59	-4846.5713	9818.8634	0.0000	1.0000
JC69I	56	-4863.4501	9845.8393	0.0000	1.0000
GTR	63	-4860.9060	9856.6444	0.0000	1.0000
HKY	59	-4870.0649	9865.8506	0.0000	1.0000
K2P	56	-4899.1605	9917.2602	0.0000	1.0000
SYM	60	-4894.6981	9917.3875	0.0000	1.0000
F81	58	-5018.6598	10160.7750	0.0000	1.0000
JC69	55	-5039.3635	10195.4154	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-4662.5135	9455.0271	0.9285	0.9285
GTRG	64	-4666.3838	9460.7676	0.0526	0.9811
HKYIG	61	-4670.4796	9462.9592	0.0176	0.9987
HKYG	60	-4674.0678	9468.1355	0.0013	1.0000
K2PIG	58	-4696.0472	9508.0944	0.0000	1.0000
GTRI	64	-4690.6364	9509.2729	0.0000	1.0000
HKYI	60	-4696.2562	9512.5124	0.0000	1.0000
SYMIG	62	-4694.2718	9512.5436	0.0000	1.0000
K2PG	57	-4699.4607	9512.9214	0.0000	1.0000
SYMG	61	-4696.5878	9515.1756	0.0000	1.0000
K2PI	57	-4718.3490	9550.6980	0.0000	1.0000
SYMI	61	-4717.0082	9556.0164	0.0000	1.0000
F81IG	60	-4825.4204	9770.8407	0.0000	1.0000
F81G	59	-4826.6372	9771.2744	0.0000	1.0000
JC69G	56	-4843.6210	9799.2420	0.0000	1.0000
JC69IG	57	-4843.0491	9800.0982	0.0000	1.0000
F81I	59	-4846.5713	9811.1426	0.0000	1.0000
JC69I	56	-4863.4501	9838.9002	0.0000	1.0000
GTR	63	-4860.9060	9847.8120	0.0000	1.0000
HKY	59	-4870.0649	9858.1298	0.0000	1.0000
SYM	60	-4894.6981	9909.3963	0.0000	1.0000
K2P	56	-4899.1605	9910.3210	0.0000	1.0000
F81	58	-5018.6598	10153.3196	0.0000	1.0000
JC69	55	-5039.3635	10188.7270	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYIG	61	-4670.4796	9760.9129	0.5350	0.5350
HKYG	60	-4674.0678	9761.2047	0.4623	0.9973
GTRIG	65	-4662.5135	9772.5187	0.0016	0.9989
GTRG	64	-4666.3838	9773.3748	0.0011	1.0000
K2PG	57	-4699.4607	9791.3372	0.0000	1.0000
K2PIG	58	-4696.0472	9791.3947	0.0000	1.0000
HKYI	60	-4696.2562	9805.5816	0.0000	1.0000
SYMG	61	-4696.5878	9813.1293	0.0000	1.0000
SYMIG	62	-4694.2718	9815.3818	0.0000	1.0000
GTRI	64	-4690.6364	9821.8800	0.0000	1.0000
K2PI	57	-4718.3490	9829.1138	0.0000	1.0000
SYMI	61	-4717.0082	9853.9701	0.0000	1.0000
F81G	59	-4826.6372	10059.4592	0.0000	1.0000
F81IG	60	-4825.4204	10063.9099	0.0000	1.0000
JC69G	56	-4843.6210	10072.7732	0.0000	1.0000
JC69IG	57	-4843.0491	10078.5139	0.0000	1.0000
F81I	59	-4846.5713	10099.3273	0.0000	1.0000
JC69I	56	-4863.4501	10112.4314	0.0000	1.0000
HKY	59	-4870.0649	10146.3145	0.0000	1.0000
GTR	63	-4860.9060	10155.5346	0.0000	1.0000
K2P	56	-4899.1605	10183.8523	0.0000	1.0000
SYM	60	-4894.6981	10202.4655	0.0000	1.0000
F81	58	-5018.6598	10436.6198	0.0000	1.0000
JC69	55	-5039.3635	10457.3738	0.0000	1.0000



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

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

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


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

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


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/712_hg18.phylip = 1 - 977;
 Partition Dummy = 1:261-aggressive-phylip/712_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

