“…The approach can in principle be extended to include estimation of the threshold, however in many cases the threshold model may be a poor description of the under-reporting mechanism which could, for example, be related to more than one covariate. Oliveira et al (2017) presents an alternative to this approach, which treats the binary under-reporting indicator I i,t,s as unobserved and therefore random. The classification of the data is characterised by I i,t,s ∼ Bernoulli(π i,t,s ), such that π i,t,s is the probability of any data point suffering from under-reporting, which is potentially informed by covariates.…”