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

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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
Minimum BIC  model: GTRIG


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	-2652.8715	5447.8788	0.9761	0.9761
GTRG	64	-2657.7746	5455.3007	0.0239	1.0000
GTRI	64	-2669.0060	5477.7635	0.0000	1.0000
HKYIG	61	-2672.8444	5478.3273	0.0000	1.0000
HKYG	60	-2677.2991	5484.8792	0.0000	1.0000
HKYI	60	-2684.3802	5499.0414	0.0000	1.0000
SYMIG	62	-2688.0657	5511.1342	0.0000	1.0000
K2PIG	58	-2695.0172	5515.6198	0.0000	1.0000
SYMG	61	-2692.7112	5518.0610	0.0000	1.0000
K2PG	57	-2699.5738	5522.3952	0.0000	1.0000
SYMI	61	-2700.6998	5534.0382	0.0000	1.0000
K2PI	57	-2706.8483	5536.9442	0.0000	1.0000
F81IG	60	-2731.3761	5593.0332	0.0000	1.0000
F81G	59	-2735.5943	5599.1185	0.0000	1.0000
GTR	63	-2731.1626	5599.6990	0.0000	1.0000
F81I	59	-2742.3225	5612.5749	0.0000	1.0000
JC69IG	57	-2751.4342	5626.1159	0.0000	1.0000
HKY	59	-2750.0061	5627.9420	0.0000	1.0000
JC69G	56	-2754.8973	5630.7108	0.0000	1.0000
JC69I	56	-2762.2442	5645.4045	0.0000	1.0000
SYM	60	-2766.0935	5662.4679	0.0000	1.0000
K2P	56	-2774.3873	5669.6909	0.0000	1.0000
F81	58	-2807.3535	5740.2925	0.0000	1.0000
JC69	55	-2828.0068	5774.6049	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-2652.8715	5435.7430	0.9802	0.9802
GTRG	64	-2657.7746	5443.5493	0.0198	1.0000
GTRI	64	-2669.0060	5466.0121	0.0000	1.0000
HKYIG	61	-2672.8444	5467.6887	0.0000	1.0000
HKYG	60	-2677.2991	5474.5983	0.0000	1.0000
HKYI	60	-2684.3802	5488.7605	0.0000	1.0000
SYMIG	62	-2688.0657	5500.1314	0.0000	1.0000
K2PIG	58	-2695.0172	5506.0344	0.0000	1.0000
SYMG	61	-2692.7112	5507.4225	0.0000	1.0000
K2PG	57	-2699.5738	5513.1476	0.0000	1.0000
SYMI	61	-2700.6998	5523.3997	0.0000	1.0000
K2PI	57	-2706.8483	5527.6967	0.0000	1.0000
F81IG	60	-2731.3761	5582.7523	0.0000	1.0000
GTR	63	-2731.1626	5588.3253	0.0000	1.0000
F81G	59	-2735.5943	5589.1887	0.0000	1.0000
F81I	59	-2742.3225	5602.6450	0.0000	1.0000
JC69IG	57	-2751.4342	5616.8683	0.0000	1.0000
HKY	59	-2750.0061	5618.0122	0.0000	1.0000
JC69G	56	-2754.8973	5621.7946	0.0000	1.0000
JC69I	56	-2762.2442	5636.4883	0.0000	1.0000
SYM	60	-2766.0935	5652.1870	0.0000	1.0000
K2P	56	-2774.3873	5660.7747	0.0000	1.0000
F81	58	-2807.3535	5730.7070	0.0000	1.0000
JC69	55	-2828.0068	5766.0136	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRIG	65	-2652.8715	5738.0112	0.8278	0.8278
GTRG	64	-2657.7746	5741.1672	0.1708	0.9986
HKYIG	61	-2672.8444	5751.3557	0.0010	0.9997
HKYG	60	-2677.2991	5753.6150	0.0003	1.0000
GTRI	64	-2669.0060	5763.6299	0.0000	1.0000
HKYI	60	-2684.3802	5767.7772	0.0000	1.0000
K2PIG	58	-2695.0172	5775.7506	0.0000	1.0000
K2PG	57	-2699.5738	5778.2135	0.0000	1.0000
SYMIG	62	-2688.0657	5788.4487	0.0000	1.0000
SYMG	61	-2692.7112	5791.0895	0.0000	1.0000
K2PI	57	-2706.8483	5792.7626	0.0000	1.0000
SYMI	61	-2700.6998	5807.0667	0.0000	1.0000
F81IG	60	-2731.3761	5861.7690	0.0000	1.0000
F81G	59	-2735.5943	5863.5551	0.0000	1.0000
F81I	59	-2742.3225	5877.0115	0.0000	1.0000
GTR	63	-2731.1626	5881.2928	0.0000	1.0000
JC69IG	57	-2751.4342	5881.9342	0.0000	1.0000
JC69G	56	-2754.8973	5882.2102	0.0000	1.0000
HKY	59	-2750.0061	5892.3786	0.0000	1.0000
JC69I	56	-2762.2442	5896.9039	0.0000	1.0000
K2P	56	-2774.3873	5921.1903	0.0000	1.0000
SYM	60	-2766.0935	5931.2038	0.0000	1.0000
F81	58	-2807.3535	6000.4232	0.0000	1.0000
JC69	55	-2828.0068	6021.7789	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/6894_hg18.phylip = 1 - 773;
 Partition Dummy = 1:261-aggressive-phylip/6894_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/6894_hg18.phylip = 1 - 773;
 Partition Dummy = 1:261-aggressive-phylip/6894_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 (GTRIG)]

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


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

