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

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

Minimum AIC  model: SYMG
Minimum AICc model: SYMG
Minimum BIC  model: SYMG


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
SYMG	61	-4627.0339	9385.9554	0.6782	0.6782
GTRIG	65	-4623.5935	9388.4615	0.1937	0.8719
SYMIG	62	-4627.8616	9389.9483	0.0921	0.9641
GTRG	64	-4626.4558	9391.8301	0.0359	1.0000
SYMI	61	-4669.8398	9471.5672	0.0000	1.0000
GTRI	64	-4667.0588	9473.0362	0.0000	1.0000
HKYG	60	-4676.1633	9481.8828	0.0000	1.0000
HKYIG	61	-4675.8238	9483.5352	0.0000	1.0000
K2PG	57	-4688.2166	9499.0315	0.0000	1.0000
K2PIG	58	-4688.5510	9502.0134	0.0000	1.0000
HKYI	60	-4716.6423	9562.8407	0.0000	1.0000
K2PI	57	-4727.9321	9578.4624	0.0000	1.0000
F81IG	60	-4784.0811	9697.7184	0.0000	1.0000
F81G	59	-4787.3446	9701.9199	0.0000	1.0000
JC69IG	57	-4800.6923	9723.9828	0.0000	1.0000
JC69G	56	-4805.6728	9731.6365	0.0000	1.0000
F81I	59	-4825.5005	9778.2317	0.0000	1.0000
JC69I	56	-4843.2905	9806.8720	0.0000	1.0000
GTR	63	-4835.7091	9807.9869	0.0000	1.0000
SYM	60	-4841.6857	9812.9276	0.0000	1.0000
HKY	59	-4873.2188	9873.6683	0.0000	1.0000
K2P	56	-4888.0469	9896.3848	0.0000	1.0000
F81	58	-4980.4519	10085.8154	0.0000	1.0000
JC69	55	-4999.2076	10116.4048	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
SYMG	61	-4627.0339	9376.0679	0.5491	0.5491
GTRIG	65	-4623.5935	9377.1869	0.3138	0.8630
SYMIG	62	-4627.8616	9379.7231	0.0883	0.9513
GTRG	64	-4626.4558	9380.9115	0.0487	1.0000
SYMI	61	-4669.8398	9461.6796	0.0000	1.0000
GTRI	64	-4667.0588	9462.1176	0.0000	1.0000
HKYG	60	-4676.1633	9472.3266	0.0000	1.0000
HKYIG	61	-4675.8238	9473.6476	0.0000	1.0000
K2PG	57	-4688.2166	9490.4333	0.0000	1.0000
K2PIG	58	-4688.5510	9493.1020	0.0000	1.0000
HKYI	60	-4716.6423	9553.2845	0.0000	1.0000
K2PI	57	-4727.9321	9569.8642	0.0000	1.0000
F81IG	60	-4784.0811	9688.1623	0.0000	1.0000
F81G	59	-4787.3446	9692.6891	0.0000	1.0000
JC69IG	57	-4800.6923	9715.3846	0.0000	1.0000
JC69G	56	-4805.6728	9723.3456	0.0000	1.0000
F81I	59	-4825.5005	9769.0010	0.0000	1.0000
GTR	63	-4835.7091	9797.4181	0.0000	1.0000
JC69I	56	-4843.2905	9798.5811	0.0000	1.0000
SYM	60	-4841.6857	9803.3715	0.0000	1.0000
HKY	59	-4873.2188	9864.4376	0.0000	1.0000
K2P	56	-4888.0469	9888.0939	0.0000	1.0000
F81	58	-4980.4519	10076.9039	0.0000	1.0000
JC69	55	-4999.2076	10108.4152	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
SYMG	61	-4627.0339	9663.8539	0.9849	0.9849
SYMIG	62	-4627.8616	9672.2270	0.0150	0.9999
GTRG	64	-4626.4558	9682.8510	0.0001	1.0000
GTRIG	65	-4623.5935	9683.8442	0.0000	1.0000
SYMI	61	-4669.8398	9749.4657	0.0000	1.0000
HKYG	60	-4676.1633	9755.3949	0.0000	1.0000
K2PG	57	-4688.2166	9759.3481	0.0000	1.0000
HKYIG	61	-4675.8238	9761.4337	0.0000	1.0000
GTRI	64	-4667.0588	9764.0571	0.0000	1.0000
K2PIG	58	-4688.5510	9766.7347	0.0000	1.0000
HKYI	60	-4716.6423	9836.3528	0.0000	1.0000
K2PI	57	-4727.9321	9838.7791	0.0000	1.0000
F81G	59	-4787.3446	9971.0396	0.0000	1.0000
F81IG	60	-4784.0811	9971.2306	0.0000	1.0000
JC69IG	57	-4800.6923	9984.2995	0.0000	1.0000
JC69G	56	-4805.6728	9987.5427	0.0000	1.0000
F81I	59	-4825.5005	10047.3514	0.0000	1.0000
JC69I	56	-4843.2905	10062.7781	0.0000	1.0000
SYM	60	-4841.6857	10086.4398	0.0000	1.0000
GTR	63	-4835.7091	10094.6398	0.0000	1.0000
HKY	59	-4873.2188	10142.7881	0.0000	1.0000
K2P	56	-4888.0469	10152.2909	0.0000	1.0000
F81	58	-4980.4519	10350.5366	0.0000	1.0000
JC69	55	-4999.2076	10367.8944	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/1062_hg18.phylip = 1 - 827;
 Partition Dummy = 1:261-aggressive-phylip/1062_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=Fixed(Equal);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/1062_hg18.phylip = 1 - 827;
 Partition Dummy = 1:261-aggressive-phylip/1062_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=Fixed(Equal);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/1062_hg18.phylip = 1 - 827;
 Partition Dummy = 1:261-aggressive-phylip/1062_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=Fixed(Equal);
END;


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

