The R code included in R_code.txt allow to perform ordinal regression for Ceriodaphnia epibiont burden (example corresponding to top panel in Table 2 of the manuscript).
The Jags code included in JAGS_model.txt allows to perform a JAGS model for ordinal regression.
The file variables.pdf contains detailed information about the variables in the database database.csv.