Highly variable relationships between shedding (MAP) and MAP-specific cellular (LPT) and humoral (ELISA) immune responses. We calculate correlation coefficients (using Spearman rank test) between measurements of humoral and cellular immunity (panel A), between shedding and cellular immunity (panel B), and shedding and humoral immunity (panel C). Because not all data were available at all time points, missing values were generated using linear interpolations (see Materials and methods for detail). Numerical values for the estimated correlation coefficients are given in Additional file 1. The histograms of the correlation coefficients were similar with fewer correlations being statistically significant if we restrict our analysis to data where all 3 variables (LPT, ELISA, MAP) were measured at the same time (results not shown).