Table 2 Multilocus MK-test model comparison
ModelDescriptionParlog(L)2∆log(L)χ2 P valueAICcaAkaike Weightbαacαbd
All loci
M0α = 038−615.831335.90.806[0][0]
M1α ∼ (all loci)39−615.510.640.42371339.00.169[0.07][0.07]
M2α ∼ (control, immune)40−615.480.060.80651342.80.0250.040.09
Imd Pathway
M0α = 018−235.44522.10.003[0][0]
M1α ∼ (all loci)19−228.0114.860.0001511.30.695[−0.64][−0.64]
M2α ∼ (control, immune)20−226.722.580.1069513.00.302−0.20−0.91
Toll pathway
M0α = 014−181.22403.20.305[0][0]
M1α ∼ (all loci)15−178.555.340.0208402.10.520[−0.5][−0.5]
M2α ∼ (control, immune)16−177.372.360.1245404.30.175−1.09−0.24
  • a The Akaike information criterion corrected for sample size.

  • b The likelihood of the model, given the relative support for each of the models tested.

  • c Estimate of the proportion of adaptive substitutions in the control genes. Square brackets indicate where α is constrained by the model.

  • d Estimate of the proportion of adaptive substitutions in the immune genes. Square brackets indicate where α is constrained by the model.