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

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

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


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	-2665.3271	5480.5153	0.6631	0.6631
HKYG	60	-2673.0030	5482.7565	0.2162	0.8793
HKYIG	61	-2672.6748	5484.6982	0.0819	0.9612
GTRG	64	-2669.4888	5486.1924	0.0388	1.0000
SYMIG	62	-2683.7036	5509.3658	0.0000	1.0000
K2PG	57	-2690.2082	5509.4436	0.0000	1.0000
K2PIG	58	-2689.3349	5510.2597	0.0000	1.0000
SYMG	61	-2687.7409	5514.8304	0.0000	1.0000
HKYI	60	-2696.8800	5530.5105	0.0000	1.0000
GTRI	64	-2694.5347	5536.2841	0.0000	1.0000
SYMI	61	-2706.7003	5552.7491	0.0000	1.0000
K2PI	57	-2716.4282	5561.8836	0.0000	1.0000
F81IG	60	-2778.0923	5692.9353	0.0000	1.0000
F81G	59	-2782.1516	5698.4675	0.0000	1.0000
JC69IG	57	-2788.5964	5706.2200	0.0000	1.0000
JC69G	56	-2792.5912	5711.6586	0.0000	1.0000
F81I	59	-2800.2374	5734.6391	0.0000	1.0000
JC69I	56	-2810.0151	5746.5065	0.0000	1.0000
GTR	63	-2802.1634	5748.9075	0.0000	1.0000
HKY	59	-2813.2645	5760.6935	0.0000	1.0000
SYM	60	-2824.4940	5785.7386	0.0000	1.0000
K2P	56	-2831.6103	5789.6967	0.0000	1.0000
F81	58	-2914.7243	5961.0386	0.0000	1.0000
JC69	55	-2926.5688	5977.0743	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-2665.3271	5460.6542	0.8723	0.8723
HKYG	60	-2673.0030	5466.0059	0.0601	0.9324
GTRG	64	-2669.4888	5466.9776	0.0369	0.9693
HKYIG	61	-2672.6748	5467.3496	0.0307	1.0000
SYMIG	62	-2683.7036	5491.4072	0.0000	1.0000
K2PG	57	-2690.2082	5494.4163	0.0000	1.0000
K2PIG	58	-2689.3349	5494.6697	0.0000	1.0000
SYMG	61	-2687.7409	5497.4817	0.0000	1.0000
HKYI	60	-2696.8800	5513.7599	0.0000	1.0000
GTRI	64	-2694.5347	5517.0693	0.0000	1.0000
SYMI	61	-2706.7003	5535.4005	0.0000	1.0000
K2PI	57	-2716.4282	5546.8563	0.0000	1.0000
F81IG	60	-2778.0923	5676.1847	0.0000	1.0000
F81G	59	-2782.1516	5682.3032	0.0000	1.0000
JC69IG	57	-2788.5964	5691.1927	0.0000	1.0000
JC69G	56	-2792.5912	5697.1824	0.0000	1.0000
F81I	59	-2800.2374	5718.4747	0.0000	1.0000
GTR	63	-2802.1634	5730.3268	0.0000	1.0000
JC69I	56	-2810.0151	5732.0303	0.0000	1.0000
HKY	59	-2813.2645	5744.5291	0.0000	1.0000
SYM	60	-2824.4940	5768.9881	0.0000	1.0000
K2P	56	-2831.6103	5775.2205	0.0000	1.0000
F81	58	-2914.7243	5945.4486	0.0000	1.0000
JC69	55	-2926.5688	5963.1376	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2673.0030	5718.6419	0.9406	0.9406
HKYIG	61	-2672.6748	5724.1962	0.0585	0.9991
GTRIG	65	-2665.3271	5734.3432	0.0004	0.9995
K2PG	57	-2690.2082	5734.4205	0.0004	0.9998
GTRG	64	-2669.4888	5736.4560	0.0001	1.0000
K2PIG	58	-2689.3349	5738.8845	0.0000	1.0000
SYMIG	62	-2683.7036	5752.4644	0.0000	1.0000
SYMG	61	-2687.7409	5754.3283	0.0000	1.0000
HKYI	60	-2696.8800	5766.3959	0.0000	1.0000
GTRI	64	-2694.5347	5786.5477	0.0000	1.0000
K2PI	57	-2716.4282	5786.8605	0.0000	1.0000
SYMI	61	-2706.7003	5792.2471	0.0000	1.0000
F81IG	60	-2778.0923	5928.8207	0.0000	1.0000
F81G	59	-2782.1516	5930.7286	0.0000	1.0000
JC69IG	57	-2788.5964	5931.1969	0.0000	1.0000
JC69G	56	-2792.5912	5932.9760	0.0000	1.0000
F81I	59	-2800.2374	5966.9001	0.0000	1.0000
JC69I	56	-2810.0151	5967.8239	0.0000	1.0000
HKY	59	-2813.2645	5992.9545	0.0000	1.0000
GTR	63	-2802.1634	5995.5946	0.0000	1.0000
K2P	56	-2831.6103	6011.0141	0.0000	1.0000
SYM	60	-2824.4940	6021.6241	0.0000	1.0000
F81	58	-2914.7243	6189.6634	0.0000	1.0000
JC69	55	-2926.5688	6194.7206	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/10368_hg18.phylip = 1 - 498;
 Partition Dummy = 1:261-aggressive-phylip/10368_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/10368_hg18.phylip = 1 - 498;
 Partition Dummy = 1:261-aggressive-phylip/10368_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 (HKYG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/10368_hg18.phylip = 1 - 498;
 Partition Dummy = 1:261-aggressive-phylip/10368_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 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);
END;


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

