On the question of proportionality of the count of observed Scrapie cases and the size of holding
- Dankmar Böhning^{1}Email author and
- Victor J Del Rio Vilas^{2, 3}
https://doi.org/10.1186/1746-6148-5-17
© Böhning and Del Rio Vilas; licensee BioMed Central Ltd. 2009
Received: 27 November 2008
Accepted: 06 May 2009
Published: 06 May 2009
Abstract
Background
The present paper investigates the question of a suitable basic model for the number of scrapie cases in a holding and applications of this knowledge to the estimation of scrapie-affected holding population sizes and adequacy of control measures within holding. Is the number of scrapie cases proportional to the size of the holding in which case it should be incorporated into the parameter of the error distribution for the scrapie counts? Or, is there a different – potentially more complex – relationship between case count and holding size in which case the information about the size of the holding should be better incorporated as a covariate in the modeling?
Methods
We show that this question can be appropriately addressed via a simple zero-truncated Poisson model in which the hypothesis of proportionality enters as a special offset-model. Model comparisons can be achieved by means of likelihood ratio testing. The procedure is illustrated by means of surveillance data on classical scrapie in Great Britain. Furthermore, the model with the best fit is used to estimate the size of the scrapie-affected holding population in Great Britain by means of two capture-recapture estimators: the Poisson estimator and the generalized Zelterman estimator.
Results
No evidence could be found for the hypothesis of proportionality. In fact, there is some evidence that this relationship follows a curved line which increases for small holdings up to a maximum after which it declines again. Furthermore, it is pointed out how crucial the correct model choice is when applied to capture-recapture estimation on the basis of zero-truncated Poisson models as well as on the basis of the generalized Zelterman estimator. Estimators based on the proportionality model return very different and unreasonable estimates for the population sizes.
Conclusion
Our results stress the importance of an adequate modelling approach to the association between holding size and the number of cases of classical scrapie within holding. Reporting artefacts and speculative biological effects are hypothesized as the underlying causes of the observed curved relationship. The lack of adjustment for these artefacts might well render ineffective the current strategies for the control of the disease.
Keywords
Background
Surveillance efforts must adjust to the levels of occurrence of disease especially in the face of multiple threats and finite resources. There is a need to consider, among other parameters that would inform the level of priority allocated to a given disease, its prevalence and, ideally, that adjusted for any source of under-ascertainment. We define ascertainment as the definitive and complete determination of individuals with a particular trait of interest, scrapie in our case. Readily available methods for the estimation of the under-ascertainment-adjusted holding prevalence and within holding prevalence are required when planning surveillance strategies and control measures in animal health settings.
The occurrence of classical scrapie, a fatal, neurological disease of small ruminants, appears to be decreasing in Great Britain, both in the number of holdings affected and the number of sheep infected within holding [1]. Such a decreasing trend will undoubtedly lead to the reconsideration by policy makers of the overall efforts dedicated to the surveillance of the disease.
The importance of the size of the holding in the occurrence of classical scrapie is well described in the literature [2–6]. The results from these studies have been used to parameterize mathematical models describing scrapie transmission between sheep and between flocks [7]. In this particular case, a positive linear relationship was assumed between the occurrence of the disease within holding and the size of the holding. Note that some of the studies above indicated some form of non proportionality between holding size and the occurrence of scrapie [5, 6]. There is a need to describe this relationship to direct disease control measures in Great Britain. Sampling and TSE-testing of scrapie-affected flocks is one of the measures in place after the introduction of the Compulsory Scrapie Flocks Scheme (CSFS) in Great Britain in 2004 [8]. The calculation of the sample size is based on the number of adult sheep within the holding and increases proportional to the size of the holding (under the assumption of a perfect test). A non-linear relationship between the occurrence of the disease and the size of the flock might render these sample calculations ineffective for larger holdings at the right end of the size distribution. Furthermore, there is a need to inform this relationship not only when dealing with reported clinical disease, potentially affected by reporting artefacts. The count of scrapie cases arising from the TSE-testing within holding under the CSFS provides the data for the assessment of this relationship free from potential reporting biases and merits study.
Previous studies have applied capture-recapture (CRC) techniques to obtain adjusted estimates of the prevalence of scrapie-affected holdings in Great Britain [9]. Recently, two approaches [10, 11] were pursued for the incorporation of observed holding-specific variability in their CRC models: via an extension of Zelterman's (1988) estimator either as a covariate in a logistic model [10] or as a proportional term [11]. The form of this relationship influences the estimates of the scrapie-affected holding population and requires further analyses.
Hence, our objective in this paper is three-fold:
i) To model the relationship between holding size and the occurrence of disease.
ii) To study the effects of this relationship in capture-recapture estimates of the size of the scrapie-affected holding population and
iii) To discuss the consequences for within-holding sampling schemes for the detection of the disease.
Methods
Materials
Several sets of analysis were conducted and for each we used different datasets. The first analysis used year-specific disease data from the Scrapie Notifications Database (SND) (see [6] for more details). The SND collects all the suspect clinical cases of scrapie reported by farmers to the veterinary authorities. More specifically, we restricted our analyses to the number of confirmed clinical cases for the years 2002, 2003 and 2004. This provided a snapshot of the relationship between the holding characteristic of interest, its size, and the presence of clinical disease. The year-specific datasets were also used to estimate the under-ascertainment-adjusted population of scrapie-affected holdings, with clinical disease, in Great Britain per year. Note that as a by product of the application of CRC models to the SND dataset, the surveillance's sensitivity, for the detection of scrapie-affected holdings with clinical disease, was estimated at around 40–50% [10].
A second set of analysis used the list of holdings sampled and TSE-tested within the CSFS during 2005 and 2006. More specifically, animals tested and confirmed from the initial cull (IC) route (see [8] for more details). In theory, this route deals only with healthy animals, randomly selected from the eligible cull population within holdings, so if any case arises it is likely to be pre-clinical. This set of analysis would inform the relationship between the number of pre-clinical cases and the size of the holding, free from reporting artefacts potentially affecting the SND results. We extended the analyses on this dataset to estimate the overall population size of scrapie-affected holdings (in the remainder of the paper simply referred to as the population size estimate). The three regular surveillance sources for scrapie in Great Britain, the SND and the fallen stock survey targeting clinical disease and the abattoir survey targeting infection (note that the SND is an exhaustive list and the two surveys are sample-based sources), feed cases into the CSFS. As a result, any population size estimate based on the CSFS would be comprehensive and represent the overall burden of the disease, regardless of its manifestation (i.e. clinical disease or infection). The authors showed the limitations of population size estimates based on CSFS data [12]. We have conducted these analyses for illustrative purposes to show the effect of an inadequate parameterisation on the model's results.
Note that for the estimation of the relationship between counts of cases and holding size for the CSFS dataset, index cases, those detected through the regular surveillance activities (for further details see [13, 14]) that triggered the incorporation of the holdings in the CSFS, were initially included in the analyses. Due to the weight of the holdings detected through the SND in the CSFS dataset [12] we would expect similar results to those obtained from modelling the SND data alone. In a further analysis, we removed the index cases to reduce the impact of potential reporting artefacts. We extended this approach to study the relationship between the number of tested animals within a holding in the CSFS and the count of cases. Here we would expect to see a linear relationship to occur in the sense that the more animals we test the more we detect (assuming a constant within-holding prevalence across categories of holding size). A non-linear relationship between case count and number tested might reveal a deviation from the assumption of constant within-holding prevalence across holdings of different sizes.
For each year and holding, our unit of analysis, we obtained the holding size from the Census data [15].
Statistical Methods
Relationship between holding size and number of cases
which we call Model 3. Note that Model 2 is nested within Model 3.
The parameters need to be estimated and this is done by means of maximum likelihood. An observed scrapie-affected holding is defined by having a non-zero confirmed count of scrapie cases within the holding. Hence zero-counts of scrapie affected holdings cannot occur. Consequently, a zero-truncated Poisson likelihood is the basis of the inference and since all models are nested, model comparisons can be achieved using the likelihood ratio test (LRT). Log-likelihoods, the Akaike Information Criterion (AIC) and the Bayesian Information Criterion (BIC) were computed for model selection. These models were applied to the SND and CSFS datasets. The latter with the index cases included.
Modelling for CSFS data with zero counts
with α = log μ, x_{ i }= log n_{ i }and β_{1} = 1. This is the offset-model if the slope parameter is fixed to 1. More general models are possible by allowing the slope parameter to be arbitrary or adding curvature terms.
2.2.3 Adjusting estimates of the size of the population of scrapie-affected holdings
where the estimates , in the linear predictor can be found by means of a logistic regression of the binary variate B (B = 1 if the scrapie count = 2, B = 0 if the scrapie count = 1) on the covariates (here only x) in the model. Either estimator, Poisson (2.4) and the generalized Zelterman (2.5), allow for the different forms of the linear predictor λ as explained in section 2.2.2. We have run these estimators on SND year-specific data and CSFS data. Note that for the latter, the index cases, those that triggered the incorporation of the scrapie-affected holdings into the CSFS, had to be retained in the dataset.
Software
All computational analysis has been done using MINITAB version 15 and STATA version 10.0.
Results
Demographic characteristics of the SND and CSFS databases: two variables are given – the non-zero case count of scrapie cases per holding and the size of holding
Variable | mean | median | minimum | maximum | n |
---|---|---|---|---|---|
SND 2002 | |||||
Count of scrapie cases per holding | 3.278 | 1 | 1 | 41 | 144 |
Holding size | 213.2 | 99 | 2 | 4433 | 125 |
SND 2003 | |||||
Count of scrapie cases per holding | 3.418 | 2 | 1 | 30 | 134 |
Holding size | 221.7 | 102.5 | 1 | 4433 | 122 |
SND 2004 | |||||
Count of scrapie cases per holding | 2.497 | 1 | 1 | 32 | 151 |
Holding size | 192.9 | 60 | 2 | 1577 | 135 |
CSFS | |||||
Count of scrapie cases per holding | 1.614 | 1 | 1 | 18 | 251 |
Holding size | 705.3 | 502 | 4 | 4264 | 214 |
Results of zero-truncated Poisson regression modelling for the offset-model (Model 1: offset is log-size of holding), for the model treating log-size as a free covariate (Model 2) and the model with a quadratic term included (Model 3) for the years 2002, 2003, and 2004 based upon the SND data. is the scrapie-affected population estimated by means of the generalized Horvitz-Thompson estimator and n the number of holdings confirmed with scrapie by the SND in each year.
α (S.E., ^{1}Z) | β _{1} | β _{2} | log-likelihood (AIC, BIC) | ||
---|---|---|---|---|---|
2002 | (n = 144) | ||||
Model 1 | -4.3160 (0.0554, -77.92) | 1 (fixed) | - | 327 | -686.43 (1374.9, 1377.7) |
Model 2 | 0.8721 (0.1857, 4.70) | 0.0560 (0.0386, 3.43) | - | 151 | -431.81 (867.6, 873.3) |
Model 3 | -0.9608 (0.6344, -1.51) | 0.9498 (0.2768, 3.43) | -0.1002 (0.0298, -3.36) | 153 | -423.10 (852.2, 860.7) |
2003 | (n = 134) | ||||
Model 1 | -4.3845 (0.0571, -76.78) | 1 | - | 359 | -700.69 (1544.3, 1547.1) |
Model 2 | 1.2485 (0.1889, 6.61) | -0.0239 (0.0405, -0.59) | - | 139 | -399.90 (803.8, 809.4) |
Model 3 | -0.8759 (0.6186, -1.42) | 1.0482 (0.2841, 3.69) | -0.1252 (0.0324, -3.86) | 147 | -388.96 (783.9, 792.3) |
2004 | (n = 151) | ||||
Model 1 | -4.5671 (0.0653, -69.98) | 1 (fixed) | - | 363 | -439.16 (880.3, 883.2) |
Model 2 | 0.5373 (0.2810, 1.91) | 0.0566 (0.0573, 0.99) | 170 | -316.82 (637.6, 643.4) | |
Model 3 | -1.0624 (0.9407, -1.13) | 0.7776 (0.3953, 1.97) | -0.0773 (0.0412, -1.88) | 172 | -314.68 (635.3, 644.1) |
Results of generalized Zelterman regression modelling estimating population size based upon the offset-model (Model 1: offset if log-size of holding), the model treating log-size as a free covariate (Model 2) and the model with a quadratic term included (Model 3 for SND year-specific data.
Results of zero-truncated Poisson regression modelling including population size estimates (columns 4 contains the robust Zelterman estimate in brackets) for the offset-model (Model 1: offset if log-size of holding), for the model treating log-size as a free covariate (Model 2) and the model with a quadratic term included (Model 3) based upon the CSFS data (n = 214).
Model 1 | coefficient | S.E./Z = coeff./S.E. | log-likelihood (AIC, BIC) | |
---|---|---|---|---|
α | -6.5573 | 0.0785/-83.48 | 457 (1077, 0 – 12,882) | -274.38 |
β _{1} | 1 (fixed) | (550.76, 554.13) | ||
Model 2 | coefficient | S.E./Z = coeff./S.E. | log-likelihood (AIC, BIC) | |
α | -3.8869 | 0.6123/-6.35 | 197 (398, 199 – 597) | -265.92 |
β _{1} | 0.6164 | 0.0893/6.90 | (535.83, 542.57) | |
Model 3 | coefficient | S.E./Z = coeff./S.E. | log-likelihood (AIC, BIC) | |
α | 1.6045 | 1.0116/1.59 | 228 (403, 200 – 606) | -257.76 |
β _{1} | -1.2240 | 0.3430/-3.57 | (521.52, 531.61) | |
β _{2} | 0.1486 | 0.0291/5.11 |
Results of zero-inflation Poisson regression modelling for the offset-model (Model 1: offset if log-size of holding), for the model treating log-size as a free covariate (Model 2) and the model with a quadratic term included (Model 3) based upon the CSFS data without index cases in the case count (n = 214).
Model 1 | coefficient | S.E. | Z = coeff./S.E. | log-likelihood (AIC, BIC) |
---|---|---|---|---|
α | -6.0402 | 0.0981 | -61.55 | -236.27 |
β _{1} | 1 (fixed) | (476.54, 483.27) | ||
Model 2 | coefficient | S.E. | Z = coeff./S.E. | log-likelihood (AIC, BIC) |
α | -4.0355 | 0.2789 | -5.54 | -232.63 |
β _{1} | 0.7140 | 0.1043 | 6.84 | (471.27, 481.37) |
Model 3 | coefficient | S.E. | Z = coeff./S.E. | log-likelihood (AIC, BIC) |
α | 2.1013 | 1.1687 | 1.80 | -223.99 |
β _{1} | -1.3656 | 0.0333 | 5.07 | (455.99. 469.45) |
β _{2} | 0.1686 | 0.0270 | -1.24 |
Results of zero-inflation Poisson regression modelling for the offset-model (Model 1: offset if log-size of number of tested animals), for the model treating log-size as a free covariate (Model 2) and the model with a quadratic term included (Model 3) based upon the CSFS data with index cases removed from the data set (n = 174).
Model 1 | coefficient | S.E. | Z = coeff./S.E. | log-likelihood (AIC, BIC) |
---|---|---|---|---|
α | -4.3154 | 0.1091 | -39.56 | -208.00 |
β _{1} | 1 (fixed) | (420.00, 426.32) | ||
Model 2 | coefficient | S.E. | Z = coeff./S.E. | log-likelihood (AIC, BIC) |
α | -3.9731 | 1.1692 | -3.40 | -207.96 |
β _{1} | 0.9332 | 0.2275 | 4.10 | (421.92, 431.40) |
Model 3 | coefficient | S.E. | Z = coeff./S.E. | log-likelihood (AIC, BIC) |
α | -6.3097 | 5.1816 | -1.22 | -207.83 |
β _{1} | 1.9721 | 2.2259 | 0.89 | (423.67, 436.30) |
β _{2} | -0.1130 | 0.2387 | -0.47 |
Discussion
Our results show that the relationship between the holding size and the number of confirmed scrapie cases followed a curved line for the years 2002, 2003 and 2004 of SND data. The scrapie count increases with holding size with a peak around holding sizes of 100 adult sheep after which the scrapie case count decreases again. This pattern not only appears consistently for the SND in the 3 years, it also occurs in the same way for the CSFS data set. This similarity between these two sets of analysis was expected due to the large weight of the SND data (70.6%) in the CSFS dataset.
In broad terms, the observed curved relationship between the case count and the holding size can be a manifestation of the dynamics of infection within the holding or an artefact related to the reporting and/or testing of cases; or more likely, a combination of both. We can only speculate about the processes potentially responsible for the observed relationship. Under a biological approach, only the exhaustion of all susceptible animals would stop the progression of the epidemic within holding. For scrapie, where the outcome is always fatal, the exhaustion of the susceptible population is only achieved by means of depopulation or selection of resistant genotypes to the disease [18]. To relate either intervention to our results, their effects would have to be different depending on the holding size. Larger holdings might have pursued genotyping of their animals at a greater rate than smaller ones. Or, similarly, the increased number of movements into larger holdings might have facilitated the shift in their genotype profiles towards resistant types more rapidly than in smaller, more static holdings. Either option might explain the reduced number of scrapie cases observed at the right end of our data.
The effect of reporting artefacts appears plausible and allows meaningful interpretation of our SND results. We could hypothesize a proportional relationship between the holding size and the level of supervision, defined in this setting as the regular observation of the holding by the farmer. It follows that with greater supervision, all other things equal (e.g. no biological effects associated with the occurrence of disease as suggested above), we should expect, from a larger holding, a greater case count. Our results on SND data show different indicating that some of the steps in the logic above do not hold. To test the rationale above, for SND data, is not straightforward as there is no clear definition of "supervision" and even less clear approach to its quantification and measurement. Our analyses of the CSFS data could provide some proxy indication of the relationship between supervision and case count. Within the CSFS, the TSE-testing equates to the supervision in SND holdings. More specifically, our results show the proportional relationship between the number of samples tested within holding and the number of cases detected in the CSFS. This is different from the curved pattern observed with the SND data and would indicate the lack of proportionality between holding size and supervision in the field. This explanation would fit well with the large body of evidence that support the occurrence of reporting artefacts affecting the ascertainment of scrapie [3, 4, 9, 10, 12, 19]. It would appear intuitive to think that larger holdings may have more difficulty in identifying all the cases of scrapie. This reduced detection capability would account for the significant quadratic relationship in the case of the SND datasets.
It appears difficult to reconcile our results with those of previous works [2, 3, 20] where larger holdings appeared at greater risk of having scrapie. On the other hand, it was shown that the risk of scrapie followed a quadratic shape with increasing holding size for Great Britain and SND data from 1993 to 2002 [6]. Note that this previous work compares counts of holdings stratified by holding size with the occurrence of scrapie regardless of the number of cases within holding. The results presented here are not directly comparable as we focus on the count of confirmed scrapie cases within holding and hence, we are not informing any measure of risk of scrapie. However, our present results and those of [6] are consistent with the occurrence of reporting artefacts: a non-proportional supervision in larger holdings will result in fewer of them reporting suspect cases.
Our results on the zero-inflated CSFS dataset inform a pre-clinical stage and, by definition, one free from human-related artefacts. They also refer to a different disease stage, infection, from that, clinical disease, shown by the analyses on SND data. Unfortunately, two sources of biases are still pertinent to our CSFS results. The first one originates from the sampling scheme performed within the CSFS, as per EU requirements, by which sampling stops at 150 animals regardless of the holding size [8]. This sampling strategy is aimed at detecting, with an assumed within-holding prevalence of 2%, at least one infected animal. Sampling for detection of disease does not favour case counts in larger holdings [21]. The second source of biases, a less certain one, would affect the representativeness of our results on the general scrapie-affected population. In [8] a large proportion of CSFS holdings with just one confirmed case of scrapie is reported, the index case. This might reflect the active search of scrapie cases by farmers in order to benefit from the compensation schemes introduced with the CSFS in 2004. Indeed, a significant increase in the number of cases in some areas in Great Britain around the introduction of the CSFS was reported recently [22].
Scrapie-affected population
It would seem natural to expect that the case count rose linearly if not proportional with the holding size and, had this been correct, it could be used beneficially in estimating the hidden burden of the scrapie epidemic. The offset-model (proportionality hypothesis) was considered in [11] and needs to be revised on the basis of these findings here. For example, for SND data and 2002 and 2004, the proportionality model returns estimates of the scrapie-affected population twice as large as those returned by the model with the best fit, that with the quadratic term. Clearly, the appropriate choice of a model is crucial for deriving an appropriate estimate of the population size of scrapie affected holdings. It appears founded then that for the surveillance schemes of SND and CSFS the hypothesis of proportionality does not hold. It was also shown that the robust generalised Zelterman estimation based upon fitting a logistic model [10] with quadratic terms for the year-specific SND data will lead to a more realistic value of the population size. Note that the appropriate choice of the model is not only crucial for avoiding a spurious population size estimate of scrapie (see Table 2), it is also important for achieving valid variance estimates leading to trustworthy confidence intervals. As Table 3 shows, all confidence intervals computed on the basis of the – inappropriate – model 1 are meaningless since they are too large. For the robust regression model the inclusion of the quadratic term is less consequential: for the years 2002 and 2004 the population size estimates based upon model 2 and model 3 are quite close (see Table 3) whereas for the year 2003 model 2 will be preferred. We can conclude from this analysis that it is less important to consider quadratic terms in the robust generalized Zelterman approach. This can easily be explained since the robust approach ignores higher scrapie counts in the holdings and, hence, is less sensitive to model changes affecting only larger scrapie counts.
Our efforts to prevent the effect of reporting artefacts on the relationship between holding size and the case count concentrated on the analyses of CSFS data after the removal of the index case. This is also a more realistic scenario as the accumulation of clinical cases within holding is no longer possible after the introduction of the CSFS in 2004 [12]. For the latter dataset, and to increase the number of observations (holdings) in our lists, we joined the two years of data, 2005 and 2006. The joint list prevented the estimation of a biologically meaningful measure of frequency from this dataset and, hence, comparisons with previous works [12].
Practical rationale
In general terms, our results have two applications: 1) for the correct adjustment of population estimates of scrapie-affected holdings and 2) to advice on the correct relationship between holding size and the count of scrapie cases within holding. The former constitutes one of the first and most basic parameters in surveillance planning. The latter will help in the development of predictive models and, on a more practical side, in the refinement and increased understanding of targeted surveillance approaches to classical scrapie. It is important to stress at this point that our results are applicable only to the classical form of scrapie. Clinical atypical scrapie was first diagnosed in Great Britain in 2005 [23], after the introduction of the CSFS and, so far, only one holding has shown multiple cases of atypical scrapie (Del Rio Vilas, personal communication). It is worth mentioning that the testing regime applied under the CSFS scheme allows the discrimination of the two types of scrapie [8].
The implications of our results are large, if not for what they are, for what they show: either a biased, for the SND, or a truncated, for the CSFS, picture of the effects of holding size in the case count within holding. This limitation affects our understanding of the results derived from the scrapie surveillance. The current active surveillance for scrapie is a targeted one: i) at the individual level by targeting sheep older than 18 months of age where the likelihood of detecting infection is greater and ii) at the holding level by, inadvertently, targeting large holdings [24]. The latter requires careful consideration after our results. The occurrence of under-ascertainment in larger holdings by means of reduced supervision will affect the detection capability of one of the surveys, the fallen stock. The lack of opportunity to spot disease can easily be extrapolated to insufficient supervision to identify and promptly report the dead-on-farm sheep. This might contribute to explain the reduced average size of scrapie-affected holdings detected by this survey relative to those detected by the abattoir survey [24].
Reporting artefacts should not be a problem for the abattoir survey due to its random nature. Farmers can do little to influence the sampling approach of the survey. Only a biological effect, e.g. a greater proportion of resistant genotypes in larger holdings, might affect the results from this survey. In fact, if there was such a biological effect, the natural tendency of the survey to over-sample larger holdings might result in the underestimation of the prevalence of scrapie. This effect might also explain the reduced sensitivity of the abattoir survey relative to the fallen stock [25].
Conclusion
With the very low numbers of detected scrapie in recent years [1] there appears to be little rationale to support a consistent sampling scheme that focuses on efforts to identify all the scrapie cases within the CSFS holdings. The very detection of scrapie cases implies the end of their infectiousness as they are detected, either through fallen stock or regular culls, at the point of removal from the holding. It is those not detected through the truncated sampling regimes that pose a risk. Our results do not allow knowing if the reduced number of scrapie cases in CSFS holdings is the result of the truncated sampling regime or there is no such effect and large holdings, for whatever reason, present fewer cases of scrapie. Only the full TSE-testing of a representative sample of CSFS holdings would answer this question. The lessons from that exercise would allow in-depth assessment of the power of the current control measures, and in particular, of the sampling scheme to detect and eliminate all cases of scrapie from the holding.
Declarations
Authors’ Affiliations
References
- Dawson M, Del Rio Vilas VJ: The control of classical scrapie in sheep in Great Britain. In practice. 2008, 30: 330-333.View ArticleGoogle Scholar
- McLean AR, Hoek AR, Hoinville LJ, Gravenor MB: Scrapie transmission in Britain: a recipe for a mathematical model. Proc Royal Society London. 1999, 266: 2531-2538. 10.1098/rspb.1999.0956.View ArticleGoogle Scholar
- Hoinville LJ, Hoek AR, Gravenor MB, McLean AR: Descriptive epidemiology of scrapie in Great Britain: results of a postal survey. Vet Rec. 2000, 146 (16): 455-461.View ArticlePubMedGoogle Scholar
- Sivam K, Baylis M, Gravenor MB, Gubbins S, Wilesmith JW: Occurrence of scrapie in GB: results of a postal survey in 2002. Veterinary Record. 2003, 153: 782-783.PubMedGoogle Scholar
- McIntyre KM, Gubbins S, Sivam SK, Baylis M: Flock-level risk factors for scrapie in Great Britain: analysis of a 2002 anonymous postal survey. BMC Vet Res. 2006, 2: 2-25. 10.1186/1746-6148-2-2.View ArticleGoogle Scholar
- Del Rio Vilas VJ, Guitian J, Pfeiffer DU, Wilesmith JW: Analysis of data from the passive surveillance of scrapie in Great Britain between 1993 and 2002. Vet Rec. 2006, 159 (24): 799-804.PubMedGoogle Scholar
- Hagenaars TJ, Donnelly CA, Ferguson NM: Epidemiological analysis of data for scrapie in Great Britain. Epidemiology and Infection. 2006, 134: 359-367. 10.1017/S0950268805004966.PubMed CentralView ArticlePubMedGoogle Scholar
- Ortiz-Pelaez A, Del Rio Vilas VJ: Within holding prevalence of sheep classical scrapie in Great Britain. BMC Veterinary Research. 2009, 5 (1): 1-10.1186/1746-6148-5-1.PubMed CentralView ArticlePubMedGoogle Scholar
- Del Rio Vilas VJ, Sayers R, Sivam K, Pfeiffer DU, Guitian J, Wilesmith JW: A case study of capture-recapture methodology using scrapie surveillance data in Great Britain. Preventive Veterinary Medicine. 2005, 67: 303-317. 10.1016/j.prevetmed.2004.12.003.View ArticlePubMedGoogle Scholar
- Böhning D, Del Rio Vilas VJ: Estimating the hidden number of scrapie affected holdings in Great Britain using a simple, truncated count model allowing for heterogeneity. Journal of Agricultural, Biological and Environmental Statistics. 2008, 13: 1-22. 10.1198/108571108X277904.View ArticleGoogle Scholar
- Navaratna WCW, Del Rio Vilas VJ, Böhning D: Extending Zelterman's approach for robust estimation of population size from zero-truncated clustered data. Biometrical Journal. 2008, 50: 584-596. 10.1002/bimj.200710441.View ArticlePubMedGoogle Scholar
- Del Rio Vilas VJ, Böhning D: Application of one-list capture-recapture models to scrapie surveillance data in Great Britain. Preventive Veterinary Medicine. 2008, 85: 253-266. 10.1016/j.prevetmed.2008.02.003.View ArticlePubMedGoogle Scholar
- Elliott H, Gubbins S, Ryan J, Ryder S, Tongue SC, Watkins G, Wilesmith JW: Prevalence of scrapie in sheep in Great Britain estimated from abattoir surveys during 2002 and 2003. Vet Rec. 2005, 157 (14): 418-9.View ArticlePubMedGoogle Scholar
- Del Rio Vilas VJ, Ryan J, Elliott HG, Tongue SC, Wilesmith JW: Prevalence of scrapie in sheep: results from fallen stock surveys in Great Britain in 2002 and 2003. Vet Rec. 2005, 157 (23): 744-745.View ArticlePubMedGoogle Scholar
- Anon: The June Agricultural Census. 2004, [http://www.statistics.gov.uk/hub/agriculture-environment/agriculture/farm-outputs]Google Scholar
- Heijden van der PGM, Bustami R, Cruyff M, Engbersen G, van Houwelingen HC: Point and interval estimation of the population size using the truncated Poisson regression model. Statistical Modelling – An International Journal. 2003, 3: 305-322. 10.1191/1471082X03st057oa.View ArticleGoogle Scholar
- Zelterman D: Robust estimation in truncated discrete distributions with application to capture-recapture experiments. Journal of Statistical Planning and Inference. 1988, 18: 225-237. 10.1016/0378-3758(88)90007-9.View ArticleGoogle Scholar
- Dawson M, Hoinville LJ, Hosie BD, Hunter N: Guidance on the use of PrP genotyping as an aid to the control of clinical scrapie. Vet Rec. 1998, 142 (23): 623-625.PubMedGoogle Scholar
- Kuchler F, Hamm S: Animal Disease Incidence and indemnity eradication programs. Agricultural Economics. 2000, 22: 299-308. 10.1111/j.1574-0862.2000.tb00076.x.View ArticleGoogle Scholar
- Sivam K, Baylis M, Gravenor MB, Gubbins S: Descriptive analysis of the results of an anonymous postal survey of the occurrence of scrapie in Great Britain in 2002. Veterinary Record. 2006, 158: 501-506.View ArticlePubMedGoogle Scholar
- Cannon RM, Roe RT: Livestock disease surveys: a field manual for veterinarians, Australian Bureau of Animal Health, Camberra 1982.Google Scholar
- Ali M, Del Rio Vilas VJ: Studying trends in scrapie reporting in Great Britain using multivariate time series. Proceedings of the Society for Veterinary Epidemiology and Preventive Medicine. [http://www.svepm.org.uk/posters/2007/Ali.pdf]
- Konold T, Davis A, Bone G, Bracegirdle J, Everitt S, Chaplin M, Saunders GC, Cawthraw S, Simmons MM: Clinical findings in two cases of atypical scrapie in sheep: a case report. BMC Veterinary Research. 2007, 3: 2-10.1186/1746-6148-3-2.PubMed CentralView ArticlePubMedGoogle Scholar
- Birch CPD, Del Rio Vilas VJ, McDonald R, Chikukwa A: The distribution of sheep sampled for scrapie in Great Britain. proceedings of Prion Fontenay aux Roses, France: Neuroprion. 2006, 52-[http://www.neuroprion.com/pdfdocs/conferences/prion2006/abstractbook.pdf]Google Scholar
- Del Rio Vilas VJ, Hopp P, Nunes T, Ru G, Sivam K, Pelaez-Ortiz A: Explaining the heterogeneous scrapie-surveillance figures across Europe: a meta-regression approach. BMC Veterinary Research. 2007, 28 (3): 13-10.1186/1746-6148-3-13.View ArticleGoogle Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.