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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	64	-1599.7379	3346.7351	0.6706	0.6706
GTRIG	65	-1599.7798	3349.4667	0.1711	0.8418
HKYG	60	-1606.8300	3350.4490	0.1047	0.9465
HKYIG	61	-1606.8391	3353.0666	0.0283	0.9748
K2PG	57	-1612.7458	3354.5530	0.0135	0.9882
SYMG	61	-1608.3374	3356.0633	0.0063	0.9946
K2PIG	58	-1612.7548	3357.1352	0.0037	0.9983
SYMIG	62	-1608.3465	3358.6929	0.0017	1.0000
GTR	63	-1611.4026	3367.4288	0.0000	1.0000
GTRI	64	-1611.4084	3370.0761	0.0000	1.0000
HKY	59	-1618.5542	3371.3099	0.0000	1.0000
HKYI	60	-1618.5601	3373.9091	0.0000	1.0000
K2P	56	-1624.9268	3376.3627	0.0000	1.0000
SYM	60	-1620.2969	3377.3828	0.0000	1.0000
K2PI	57	-1624.9325	3378.9265	0.0000	1.0000
SYMI	61	-1620.3026	3379.9937	0.0000	1.0000
F81G	59	-1635.0900	3404.3814	0.0000	1.0000
JC69G	56	-1639.4936	3405.4963	0.0000	1.0000
F81IG	60	-1635.0991	3406.9873	0.0000	1.0000
JC69IG	57	-1639.5026	3408.0667	0.0000	1.0000
F81	58	-1646.5202	3424.6659	0.0000	1.0000
JC69	55	-1651.2725	3426.5132	0.0000	1.0000
F81I	59	-1646.5260	3427.2534	0.0000	1.0000
JC69I	56	-1651.2781	3429.0653	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1599.7379	3327.4759	0.7026	0.7026
GTRIG	65	-1599.7798	3329.5595	0.2479	0.9505
HKYG	60	-1606.8300	3333.6600	0.0319	0.9824
HKYIG	61	-1606.8391	3335.6781	0.0116	0.9941
SYMG	61	-1608.3374	3338.6748	0.0026	0.9967
K2PG	57	-1612.7458	3339.4915	0.0017	0.9984
SYMIG	62	-1608.3465	3340.6929	0.0009	0.9993
K2PIG	58	-1612.7548	3341.5096	0.0006	1.0000
GTR	63	-1611.4026	3348.8053	0.0000	1.0000
GTRI	64	-1611.4084	3350.8168	0.0000	1.0000
HKY	59	-1618.5542	3355.1085	0.0000	1.0000
HKYI	60	-1618.5601	3357.1201	0.0000	1.0000
SYM	60	-1620.2969	3360.5938	0.0000	1.0000
K2P	56	-1624.9268	3361.8536	0.0000	1.0000
SYMI	61	-1620.3026	3362.6052	0.0000	1.0000
K2PI	57	-1624.9325	3363.8650	0.0000	1.0000
F81G	59	-1635.0900	3388.1800	0.0000	1.0000
F81IG	60	-1635.0991	3390.1983	0.0000	1.0000
JC69G	56	-1639.4936	3390.9872	0.0000	1.0000
JC69IG	57	-1639.5026	3393.0052	0.0000	1.0000
F81	58	-1646.5202	3409.0404	0.0000	1.0000
F81I	59	-1646.5260	3411.0520	0.0000	1.0000
JC69	55	-1651.2725	3412.5450	0.0000	1.0000
JC69I	56	-1651.2781	3414.5562	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1612.7458	3579.3811	0.9259	0.9259
K2PIG	58	-1612.7548	3585.6078	0.0412	0.9670
HKYG	60	-1606.8300	3586.1754	0.0310	0.9980
HKYIG	61	-1606.8391	3592.4021	0.0014	0.9994
SYMG	61	-1608.3374	3595.3988	0.0003	0.9997
GTRG	64	-1599.7379	3596.8256	0.0002	0.9999
K2P	56	-1624.9268	3597.5347	0.0001	1.0000
SYMIG	62	-1608.3465	3601.6255	0.0000	1.0000
GTRIG	65	-1599.7798	3603.1179	0.0000	1.0000
HKY	59	-1618.5542	3603.4153	0.0000	1.0000
K2PI	57	-1624.9325	3603.7546	0.0000	1.0000
HKYI	60	-1618.5601	3609.6355	0.0000	1.0000
SYM	60	-1620.2969	3613.1092	0.0000	1.0000
GTR	63	-1611.4026	3613.9464	0.0000	1.0000
SYMI	61	-1620.3026	3619.3292	0.0000	1.0000
GTRI	64	-1611.4084	3620.1666	0.0000	1.0000
JC69G	56	-1639.4936	3626.6683	0.0000	1.0000
JC69IG	57	-1639.5026	3632.8949	0.0000	1.0000
F81G	59	-1635.0900	3636.4868	0.0000	1.0000
F81IG	60	-1635.0991	3642.7137	0.0000	1.0000
JC69	55	-1651.2725	3644.0174	0.0000	1.0000
JC69I	56	-1651.2781	3650.2372	0.0000	1.0000
F81	58	-1646.5202	3653.1386	0.0000	1.0000
F81I	59	-1646.5260	3659.3588	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/7317_hg18.phylip = 1 - 497;
 Partition Dummy = 1:261-aggressive-phylip/7317_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/7317_hg18.phylip = 1 - 497;
 Partition Dummy = 1:261-aggressive-phylip/7317_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 (K2PG)]

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

