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

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

Minimum AIC  model: SYMG
Minimum AICc model: SYMG
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
SYMG	61	-2670.6437	5484.2985	0.8766	0.8766
K2PG	57	-2678.3895	5488.9438	0.0859	0.9625
GTRG	64	-2670.0170	5491.3393	0.0259	0.9885
HKYG	60	-2676.5636	5493.4042	0.0092	0.9977
K2PIG	58	-2681.1382	5497.1304	0.0014	0.9992
SYMIG	62	-2676.3718	5498.5041	0.0007	0.9999
HKYIG	61	-2679.4758	5501.9628	0.0001	1.0000
GTRIG	65	-2682.3938	5518.8888	0.0000	1.0000
SYMI	61	-2688.8237	5520.6586	0.0000	1.0000
GTRI	64	-2688.5444	5528.3941	0.0000	1.0000
K2PI	57	-2699.5740	5531.3128	0.0000	1.0000
HKYI	60	-2697.9132	5536.1033	0.0000	1.0000
JC69IG	57	-2774.5948	5681.3544	0.0000	1.0000
F81IG	60	-2773.8056	5687.8883	0.0000	1.0000
JC69G	56	-2788.1195	5705.7293	0.0000	1.0000
F81G	59	-2787.5595	5712.6770	0.0000	1.0000
JC69I	56	-2792.0040	5713.4984	0.0000	1.0000
GTR	63	-2790.7549	5730.0349	0.0000	1.0000
F81I	59	-2797.7262	5733.0103	0.0000	1.0000
SYM	60	-2796.6072	5733.4914	0.0000	1.0000
K2P	56	-2807.1926	5743.8755	0.0000	1.0000
HKY	59	-2805.4587	5748.4755	0.0000	1.0000
JC69	55	-2913.0746	5952.9797	0.0000	1.0000
F81	58	-2912.6472	5960.1484	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
SYMG	61	-2670.6437	5463.2874	0.8881	0.8881
GTRG	64	-2670.0170	5468.0340	0.0827	0.9709
K2PG	57	-2678.3895	5470.7790	0.0210	0.9918
HKYG	60	-2676.5636	5473.1272	0.0065	0.9983
SYMIG	62	-2676.3718	5476.7437	0.0011	0.9994
K2PIG	58	-2681.1382	5478.2764	0.0005	0.9999
HKYIG	61	-2679.4758	5480.9517	0.0001	1.0000
GTRIG	65	-2682.3938	5494.7877	0.0000	1.0000
SYMI	61	-2688.8237	5499.6475	0.0000	1.0000
GTRI	64	-2688.5444	5505.0888	0.0000	1.0000
K2PI	57	-2699.5740	5513.1479	0.0000	1.0000
HKYI	60	-2697.9132	5515.8263	0.0000	1.0000
JC69IG	57	-2774.5948	5663.1896	0.0000	1.0000
F81IG	60	-2773.8056	5667.6113	0.0000	1.0000
JC69G	56	-2788.1195	5688.2389	0.0000	1.0000
F81G	59	-2787.5595	5693.1190	0.0000	1.0000
JC69I	56	-2792.0040	5696.0080	0.0000	1.0000
GTR	63	-2790.7549	5707.5097	0.0000	1.0000
SYM	60	-2796.6072	5713.2144	0.0000	1.0000
F81I	59	-2797.7262	5713.4523	0.0000	1.0000
K2P	56	-2807.1926	5726.3851	0.0000	1.0000
HKY	59	-2805.4587	5728.9174	0.0000	1.0000
JC69	55	-2913.0746	5936.1491	0.0000	1.0000
F81	58	-2912.6472	5941.2944	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-2678.3895	5701.3443	0.9835	0.9835
SYMG	61	-2670.6437	5710.0327	0.0128	0.9962
K2PIG	58	-2681.1382	5712.8867	0.0031	0.9993
HKYG	60	-2676.5636	5715.8275	0.0007	1.0000
GTRG	64	-2670.0170	5726.9143	0.0000	1.0000
SYMIG	62	-2676.3718	5727.5340	0.0000	1.0000
HKYIG	61	-2679.4758	5727.6970	0.0000	1.0000
K2PI	57	-2699.5740	5743.7132	0.0000	1.0000
SYMI	61	-2688.8237	5746.3928	0.0000	1.0000
GTRIG	65	-2682.3938	5757.7130	0.0000	1.0000
HKYI	60	-2697.9132	5758.5266	0.0000	1.0000
GTRI	64	-2688.5444	5763.9691	0.0000	1.0000
JC69IG	57	-2774.5948	5893.7549	0.0000	1.0000
F81IG	60	-2773.8056	5910.3116	0.0000	1.0000
JC69G	56	-2788.1195	5914.7592	0.0000	1.0000
JC69I	56	-2792.0040	5922.5283	0.0000	1.0000
F81G	59	-2787.5595	5931.7743	0.0000	1.0000
F81I	59	-2797.7262	5952.1076	0.0000	1.0000
K2P	56	-2807.1926	5952.9054	0.0000	1.0000
SYM	60	-2796.6072	5955.9147	0.0000	1.0000
GTR	63	-2790.7549	5962.3451	0.0000	1.0000
HKY	59	-2805.4587	5967.5728	0.0000	1.0000
JC69	55	-2913.0746	6158.6244	0.0000	1.0000
F81	58	-2912.6472	6175.9047	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/6344_hg18.phylip = 1 - 422;
 Partition Dummy = 1:261-aggressive-phylip/6344_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/6344_hg18.phylip = 1 - 422;
 Partition Dummy = 1:261-aggressive-phylip/6344_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 (K2PG)]

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

