Variations in body condition and PC1 P were assessed with general linear models (LMs) initially including terms for cluster membership, body size (SVL) and host stage (classified as in Jackson et al. [18] ), with deletion of non-significant terms. For analyses of individual macroparasite species, LMs with normal errors or equivalent generalized linear models (GLMs) with Poisson errors or negative binomial errors were used as appropriate (again with deletion of non-significant terms). A Poisson distribution was used as a surrogate distribution in the case of ordinal data for P. serrata. In the case of negative binomial errors, the dispersion parameter k was estimated from that of the entire host population and any significant results reported are also consistent with an analysis of log-transformed data in an LM.