Modelling anthelmintic resistance by extending eggCounts package to allow individual efficacy

The same anthelmintic treatment can have variable efficacy on individual animals even if the parasite population is homogenously susceptible. An extension of the R package eggCounts is proposed to take individual efficacy into account using a Bayesian hierarchical model. A simulation study is conducted to compare the performance of five different methods on estimating faecal egg count reduction and its uncertainty interval. Simulation results showed the individual efficacy model offered robust inference to two different data simulation procedures with low root mean squared error on the reduction estimate and appropriate uncertainty estimates. Different methods were used to evaluate the anthelmintic resistance in a dataset from USA with sheep and cattle faecal egg counts, where a strong anthelmintic resistance was detected. Open-source statistical tools were updated to include the proposed model.