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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: GTRG


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
GTRG	64	-1563.6652	3269.2203	0.7751	0.7751
GTRIG	65	-1563.6747	3271.6972	0.2247	0.9998
HKYG	60	-1577.1278	3286.3949	0.0001	0.9999
HKYIG	61	-1577.1372	3288.8392	0.0000	1.0000
SYMG	61	-1577.5369	3289.6386	0.0000	1.0000
SYMIG	62	-1577.5464	3292.0910	0.0000	1.0000
K2PG	57	-1587.9073	3300.7255	0.0000	1.0000
K2PIG	58	-1587.9168	3303.1460	0.0000	1.0000
GTR	63	-1589.2803	3318.0007	0.0000	1.0000
GTRI	64	-1589.2866	3320.4630	0.0000	1.0000
F81G	59	-1598.7098	3327.1414	0.0000	1.0000
F81IG	60	-1598.7192	3329.5777	0.0000	1.0000
HKY	59	-1601.3818	3332.4855	0.0000	1.0000
HKYI	60	-1601.3883	3334.9159	0.0000	1.0000
SYM	60	-1602.0222	3336.1838	0.0000	1.0000
JC69G	56	-1607.8431	3338.2036	0.0000	1.0000
SYMI	61	-1602.0285	3338.6218	0.0000	1.0000
JC69IG	57	-1607.8525	3340.6160	0.0000	1.0000
K2P	56	-1611.9452	3346.4077	0.0000	1.0000
K2PI	57	-1611.9515	3348.8140	0.0000	1.0000
F81	58	-1622.0130	3371.3385	0.0000	1.0000
F81I	59	-1622.0194	3373.7606	0.0000	1.0000
JC69	55	-1631.1591	3382.4498	0.0000	1.0000
JC69I	56	-1631.1653	3384.8480	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1563.6652	3255.3305	0.7328	0.7328
GTRIG	65	-1563.6747	3257.3494	0.2671	0.9999
HKYG	60	-1577.1278	3274.2556	0.0001	1.0000
HKYIG	61	-1577.1372	3276.2744	0.0000	1.0000
SYMG	61	-1577.5369	3277.0738	0.0000	1.0000
SYMIG	62	-1577.5464	3279.0927	0.0000	1.0000
K2PG	57	-1587.9073	3289.8146	0.0000	1.0000
K2PIG	58	-1587.9168	3291.8336	0.0000	1.0000
GTR	63	-1589.2803	3304.5607	0.0000	1.0000
GTRI	64	-1589.2866	3306.5732	0.0000	1.0000
F81G	59	-1598.7098	3315.4195	0.0000	1.0000
F81IG	60	-1598.7192	3317.4384	0.0000	1.0000
HKY	59	-1601.3818	3320.7636	0.0000	1.0000
HKYI	60	-1601.3883	3322.7766	0.0000	1.0000
SYM	60	-1602.0222	3324.0445	0.0000	1.0000
SYMI	61	-1602.0285	3326.0570	0.0000	1.0000
JC69G	56	-1607.8431	3327.6863	0.0000	1.0000
JC69IG	57	-1607.8525	3329.7051	0.0000	1.0000
K2P	56	-1611.9452	3335.8904	0.0000	1.0000
K2PI	57	-1611.9515	3337.9031	0.0000	1.0000
F81	58	-1622.0130	3360.0261	0.0000	1.0000
F81I	59	-1622.0194	3362.0387	0.0000	1.0000
JC69	55	-1631.1591	3372.3182	0.0000	1.0000
JC69I	56	-1631.1653	3374.3307	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-1563.6652	3543.2205	0.5157	0.5157
HKYG	60	-1577.1278	3544.1525	0.3236	0.8393
K2PG	57	-1587.9073	3546.2167	0.1153	0.9546
GTRIG	65	-1563.6747	3549.7377	0.0198	0.9745
HKYIG	61	-1577.1372	3550.6696	0.0124	0.9869
SYMG	61	-1577.5369	3551.4690	0.0083	0.9952
K2PIG	58	-1587.9168	3552.7340	0.0044	0.9997
SYMIG	62	-1577.5464	3557.9862	0.0003	1.0000
JC69G	56	-1607.8431	3579.5901	0.0000	1.0000
F81G	59	-1598.7098	3580.8182	0.0000	1.0000
JC69IG	57	-1607.8525	3586.1072	0.0000	1.0000
HKY	59	-1601.3818	3586.1623	0.0000	1.0000
F81IG	60	-1598.7192	3587.3354	0.0000	1.0000
K2P	56	-1611.9452	3587.7942	0.0000	1.0000
GTR	63	-1589.2803	3587.9525	0.0000	1.0000
HKYI	60	-1601.3883	3592.6735	0.0000	1.0000
SYM	60	-1602.0222	3593.9414	0.0000	1.0000
K2PI	57	-1611.9515	3594.3051	0.0000	1.0000
GTRI	64	-1589.2866	3594.4632	0.0000	1.0000
SYMI	61	-1602.0285	3600.4522	0.0000	1.0000
JC69	55	-1631.1591	3619.7238	0.0000	1.0000
F81	58	-1622.0130	3620.9264	0.0000	1.0000
JC69I	56	-1631.1653	3626.2345	0.0000	1.0000
F81I	59	-1622.0194	3627.4374	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/6928_hg18.phylip = 1 - 664;
 Partition Dummy = 1:261-aggressive-phylip/6928_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=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/6928_hg18.phylip = 1 - 664;
 Partition Dummy = 1:261-aggressive-phylip/6928_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=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/6928_hg18.phylip = 1 - 664;
 Partition Dummy = 1:261-aggressive-phylip/6928_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=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

