Skip to main content
  • Research article
  • Open access
  • Published:

Development of a high number, high coverage dog rabies vaccination programme in Sri Lanka



Rabies is estimated to cause 59,000 deaths and economic losses of US$8.6 billion every year. Despite several years of rabies surveillance and awareness programmes, increased availability of post-exposure prophylaxis vaccinations and dog population control, the disease still remains prevalent in Sri Lanka. This study reports the roll-out of a high number, high coverage canine rabies vaccination campaign in Sri Lanka, providing estimates for the vaccination coverage achieved, analysing the local dog demographics, and identifying barriers of attendance to static vaccination clinics.


A mass dog vaccination campaign was undertaken in Negombo, Sri Lanka. The campaign was composed of static point and door-to-door vaccination stages, with a final survey of vaccination coverage. A large volume of data on the distribution, health, and signalment of vaccinated dogs was collected through a mobile phone application. A logistic regression model was developed to investigate which socio-spatial and dog-related factors influenced attendance of owners to static vaccination points.


The campaign vaccinated over 7800 dogs achieving a vaccination coverage of 75.8%. A dog:human ratio of 1:17 was estimated. Most dogs were owned, and the dog population was mostly male, adult, and non-sterilized. Unawareness, unavailability and handling problems were the most common reasons given by owners to explain failure to attend a static vaccination point. The regression analysis showed that increasing distance to a static point, in addition to young age and poor health of the dog, were associated with a decrease in the likelihood of attendance to a static vaccination points.


This study demonstrates the feasibility of high number, high coverage vaccination campaigns in Sri Lanka. The information on dog ecology and barriers of attendance to static point vaccination clinics will facilitate development of future vaccination campaigns.

Peer Review reports


Rabies is still a prevalent and underreported disease in many developing countries, causing 59,000 deaths each year, and economic losses amounting to 8.6 billion USD annually [1]. Ninety-nine percent of human rabies cases can be traced back to bites from rabies infected dogs [2]. Infected patients develop fatal encephalitis unless they are treated shortly after the bite with post-exposure prophylaxis (PEP), a treatment not widely available in many rabies-endemic regions [3, 4]. Since dog bites are the main source of human infections [2], mass dog vaccination campaigns represent the most effective course of action to reduce rabies incidence [1, 5], and have been shown to reduce human infection rates in multiple settings [6,7,8]. The World Health Organization recommends a minimum annual vaccination coverage of 70% of the dog population [9, 10] in order to achieve herd immunity, reduce rabies’ incidence, and minimize the burden of the disease in both dog and human populations [5].

In order to achieve this goal, different dog vaccination strategies are used depending on the local dog ecology, dog ownership structure, and the resources available. Static point (SP) vaccination approaches are commonly used, as they are easy to establish and are an efficient way to vaccinate a large number of dogs with limited personnel. However, reliance on SP clinics only often leads to a failure to vaccinate a sufficient proportion of the population, since they depend on high ownership levels and other socio-economic and cultural factors which may influence dog owners’ attendance [11]. Although more logistically challenging and costly, “door-to-door” (D2D) approaches based on visiting all households in a community and vaccinating free-roaming and owned dogs are able to achieve high coverages very effectively [12, 13]. For this reason, many vaccination campaigns combine SP and D2D stages to maximise vaccination coverage in a feasible and cost-effective manner [14, 15].

Rabies-related human deaths in Sri Lanka have steadily diminished since the establishment of the rabies control program in 1975. Despite more than 30 years of efforts, the disease still represents a serious concern for the island, especially in regions with economies largely sustained by tourism. The latest surveys in 2014 estimated a dog vaccination coverage of 48% [16], which is considered to be too low to rapidly eliminate the disease. Rabies remains prevalent among the Sri Lankan canine population, resulting in around thirty human deaths each year [16]. Consequently, PEP expenditure represents an important economic burden, amounting to more than 300,000 prescriptions annually [16, 17]. Nationwide efforts to reduce the prevalence of the disease continue to be based on treatment availability, dog vaccination, and stray dog population control [16,17,18]. However, dog vaccination and neutering methods in Sri Lanka have been poorly reported, with scarce data available on dog demographics and vaccination coverage. The lack of a working template for high numbers, high coverage campaigns is one of the factors hindering the development of efficient vaccination operations that could be applied in many of regions of Sri Lanka to enable the nation to reach sufficient dog vaccination coverage. This report describes a mass dog vaccination campaign carried out in Negombo by Mission Rabies [19] and the Dogstar Foundation [20], which was able to vaccinate a large number of dogs, obtaining a high vaccination coverage, and whose design could be implemented in other regions of the country to reduce the incidence of rabies. The objectives of this study included: the estimation of the vaccination coverage achieved, the analysis of the demographics of the local dog population, and the identification of barriers of attendance to SP.


Study area

Negombo is the biggest city in the Gampaha District, located in the western coast of Sri Lanka. Negombo also constitutes one of the biggest hubs of the country, with a population of 142,136 inhabitants [21].

Mission rabies 2016 vaccination campaign: period and course of the campaign

A pilot D2D campaign in 2015 was used to demonstrate the feasibility and effectiveness of working protocols. The 2016 campaign was performed between June and September, and covered most of the city, whose area was divided into 33 wards (Fig. 1), denominated with numbers 1 to 34 and missing number 21. Based on previous information regarding ownership levels from the 2015 pilot campaign [19] and reports from another municipality in Gampaha District [22], a SP stage was included in the campaign. A total of 146 static vaccination points were set up daily in different areas of the city (Fig. 1), from June 15th to September 1st, giving coloured collars to the animals administered vaccination (Nobivac® Rabies, MSD Animal Health). The SP efforts were complemented by a D2D stage, from September 12th to 28th, in which the staff members covered the entire area of each working zone, vaccinating and marking (with paint and/or a collar) any dogs found on the street and knocking on household doors to offer vaccination to owned dogs. Free-roaming dogs which did not approach the teams willingly were caught using lightweight Balinese-style nets, and restrained securely. All staff members involved in the vaccination process were trained to do so causing the minimum distress to the animals. A survey to assess the vaccination coverage was performed from September 13th to 29th, recording the presence or absence of paint marks and collars on any dogs seen while travelling on every traversable road of each ward using a tuk tuk. Surveys were usually carried the day after the D2D vaccination, otherwise surveys were carried out maximum 2 days after. In 2 wards where the coverage results were below 70%, the vaccination teams were sent for a second round of D2D vaccination and subsequent survey. Survey entries gathered in wards 13, 24, 25, 26, 27, 28, 29, 30 and 34 were collected by a different surveyor, since the original surveying supervisor was not available. Due to the lack of training, entries from these 9 wards were not included in the coverage assessment. Field data during all stages was collected using the Mission Rabies App [23], a web-based platform created for simple management and streamlined entry of field data. The app collected global positioning system (GPS) coordinates and timestamps for each dog vaccinated automatically and offered a path-tracking tool allowing the staff to check their spatial coverage in real time. The app also requested additional relevant data of the dog such as sex, age, health status or neuter status in addition to other information such as the team responsible for the vaccination, the provider of any previous vaccine, and the opinion of owners on matters such as neutering procedures and the reason why they did not bring their dog to a SP.

Fig. 1
figure 1

Topographical division of Negombo in 33 wards. The division into wards was performed according to the working zone shapefile provided by Mission Rabies. The location of the 146 static vaccination points is marked as dots. The location of Negombo within Sri Lanka is shown on the right. Background map tiles by Stamen Design, under CC BY 3.0. Data by OpenStreetMap, under ODbL

Data sources

Collection tools and data collected

Field data collected using the Mission Rabies App [23] was compiled into datasets and used for this study. A shapefile containing the specifics of the polygons used by Mission Rabies to divide Negombo in the 33 wards was also provided.

Other sources

Data regarding additional geospatial variables for the regression analysis was obtained from publicly available sources. Weather data for Negombo during the campaign was obtained from the World Weather Online database [24], averaging the tri-hourly measurements for temperature and precipitation during the working hours (6:00–15:00) for each day of the SP campaign. A raster file depicting population density in 2015 (100 m resolution) was obtained from WorldPop [25] and used to extract population density at each D2D entry coordinates. Since the raster did not cover the entirety of the working area, entries located outside of the raster’s data grid were applied a buffer averaging the values found in a 550 m radius. In the absence of suitable poverty geodata for Negombo, a raster file depicting the number of underweight children under five years old (in the year 2000) with a 2.5-min resolution was obtained from the Socioeconomic Data and Applications Center (SEDAC) [26]. Hunger data for each D2D entry coordinate was extracted from the raster and used as proxy data for poverty. A shapefile containing land types (in 2012) based on the national 7-class classification scheme was obtained from the World Bank Group [27] and used to determine land type for every D2D entry coordinate. The R package ggmap [28] was used to plot all maps, using background tiling sourced from Stamen Design (using data by OpenStreetMap [29]), available under CC BY 3.0 license.

Data analysis

Data manipulation and analysis was performed using the R statistical software environment version 3.4.3 [30].. The use of specific R packages for particular purposes is detailed below:

Identification and management of GPS outliers

The GPS data was checked to detect any erroneous entries caused by coordinate inaccuracy by automated recording systems [31]. These inaccuracies were mostly represented by coordinates ending up outside the visited ward, or in the ocean. These were considered to occur at random, and a system for outlier coordinate detection was created. This system was irrelevant for the SP stage since coordinates for all the vaccination points were known. For the D2D locations, due to the lack of reference coordinates, this discrimination was performed by spatial clustering. Using the dbscan function from the fpc package [32], clusters of vaccination coordinates were determined. The function required 2 parameters to be specified in advance: the size of the epsilon neighbourhood was set at 0.005, and the minimum number of neighbours was set to 3. Clusters containing < 20% of the total points for a given day and team responsible were considered outliers. The detection system was complemented by a visual screening of the coordinates. At the end of the process, 836 entries were marked as outliers and not included in the regression analysis due to the need for spatial accuracy, however they were included in the vaccination coverage analysis after adjusting their ward according to Mission Rabies working schedules. The assumption that the outliers occurred at random was checked by comparing the distribution of the variables included in the regression model between the outlier non-outlier entries.

Estimation of coverage by ward

Vaccination coverages for each ward were calculated based on the number of dogs marked with a collar or paint sighted during the survey out of the total number of dogs sighted. The 95% binomial confidence interval (CI) was calculated using the binom.test function from base R [30], which carries out an exact binomial test. The over function from the rgeos package [33] was used to match the coordinates from the data with the ward shapefile, which was imported using the rgdal package [34].

Analysis of dog demographics

The dog population in Negombo was estimated using the Chapman estimator [35] for mark and recapture. The dog population density was calculated using the area for the Mission Rabies working zone (30 km2). The 95% CI for the population size was calculated using the ciChapman function from the recapr package [36], with the default bootstrap method. Data obtained from the D2D stage was used to study the dog demographics on sex, age, ownership and neuter status of the dogs, as it contained information on both the owned and stray populations. In order to determine any relationship between the dog-related variables (sex, age, ownership status, vaccination status and neuter status), the Chi-Squared Test for Independence was used, with a significance level (α) of 0.05. In cases where any of the expected frequencies were lower than 5, the Fisher’s Exact Test was used instead.

Analysis of owner opinions on neutering and failure to attend SP

The answers given by the owners when asked “are you opposed to dog sterilization, and if so why?” and “why didn’t you attend to the SP campaign?” were compiled into a frequency table. Infrequent answers were included under the “Other” category.

Development of the logistic regression model

A logistic regression model analysing the effect of several dog-related and geospatial factors on attendance of owners to SP was built using the glm (generalised linear model) function. Attendance of a dog to SP was defined as the outcome variable, supplied as TRUE/FALSE. For the model, only entries regarding owned dogs found during the D2D stage were considered (4310 entries), filtering out any coordinate outliers (675 entries). The RANN package [37] was used to determine the closest static vaccination point from each owner’s household. Packages rgdal [34] and raster [38] were used to manipulate shapefiles and raster files, respectively. The linearity of continuous predictive variables was tested using the Box-Tidwell Transformation test [39], and variables “distance to SP” and “precipitation” were transformed into categorical variables as linearity could not be assumed. “Distance to SP” was divided into four quantiles: [2.1–126], (126–194], (194–279] and (279–1100] (distance in metres), whereas the “precipitation” variable was categorized into slight (< 0.5 mm/h), moderate (0.5–4 mm/h), and heavy (> 4 mm/h) precipitation, according to the UK meteorological office [40]. In order to increase the counts of young animals, “juvenile” and “puppy” dogs were merged into “young”, as opposed to “adult” dogs (over one year of age).

Best model selection

A series of different models were built by a mixture of forward and backwards selection using combinations of the relevant explanatory variables and their interactions. After filtering unowned dogs and outlier entries, the dataset was composed of 3635 entries. To determine the performance of the proposed regression models, this data was randomly partitioned into “training” and “testing” subsets, in a 7:3 ratio, using the vtreat package [41]. This partition allowed to generate the model from the “training” subset, and test its predictive power on the “testing” subset using the caret package [42]. To aid the variable selection process, 5-fold cross validation was used [43], partitioning the “training” data into 5 sets. This process produced goodness of fit estimations (area under the receiving operator curve (AUC)) averaged from the 5 sequential analyses, that allowed for the comparison of the different models in order to select the best. This decision was taken based on two parameters: the AUC, and the Akaike information criterion (AIC) [44]. Once the best model was chosen, its predictive power was determined using the initial “testing” subset by estimating the AUC using the package ROCR [45]. Using the final model, ward and vaccination team were included as random effects through mixed-effects regression. Since their influence on the odds ratios was negligible and for the sake of simplicity, the fixed effects model was chosen.


Number of dogs vaccinated

During the campaign, a total of 7804 dogs were vaccinated. 4382 dogs were vaccinated at the 146 different SP across the city. During the D2D vaccination stage, a further 3422 dogs were vaccinated. 115 dogs (3.2% of all encountered) were not vaccinated during the D2D stage, the most common reasons being lack of owner consent (61%) and handling difficulties (22.9%). Out of the 5177 dogs seen during the D2D campaign, 1019 (19.7%) had been vaccinated during the SP campaign.

Estimation of the vaccination coverage

The average vaccination coverage across the 24 accurately surveyed wards was estimated as 75.8% (95% CI: 73.6–77.9%), with 1194 dogs out of the 1576 sighted being marked as vaccinated. Vaccination coverages above 70% were estimated in 23 out of the 24 wards (Fig. 2). Table 1 includes a summary of the vaccination efforts during the SP and D2D stages, including vaccination coverages per ward. Dogs vaccinated at SP represented 62.1% of the total dogs found to be previously vaccinated. The remaining had been vaccinated by a private veterinarian (25.4%), in a Dogstar clinic (5.4%) or by another source (7%).

Fig. 2
figure 2

Plot of vaccination coverage achieved in each ward. The 95% Confidence Interval (CI) is represented as vertical bars. Wards in which the coverage surpassed 70% are coloured in green, while wards under 70% are coloured in red. Wards whose coverage and 95% CI lower bound surpass 70% are coloured in blue

Table 1 Summary of vaccination numbers and estimated coverage per ward

General dog demographics

Using the Chapman estimation, a dog population of 8370 was estimated (95% CI: 8138 – 8608). The dog population density was estimated to be 279 dogs/km2 in Negombo. Based on a human population density of 4737 inhabitants/km2 [21], these estimates produced a dog:human ratio of 1:17. There was a modestly lower proportion of female dogs (42.2%) compared to males (56.3%). Most of the dogs found during the D2D campaign were considered adults (over one year of age) (87.2%), with only a small proportion classified as juvenile (3 to 12 months old) (5.9%) or puppies (less than 3 months old) (6.9%). A positive relationship was found between female sex and young age (p < 0.001). The great majority of the dogs (83.2%) were recorded as owned. More female dogs were considered stray, compared to male dogs (p < 0.001). A higher proportion of young (under 1 year of age) stray animals was also found, (p < 0.001) compared to a larger, adult-owned population. Almost a third (31.6%) of the dogs seen had been previously vaccinated (at the SP stage, in Dogstar clinics, or by other veterinarians), while the rest (68.3%) were unvaccinated. Young animals were seldom vaccinated (p < 0.001), and male dogs were more likely to be already vaccinated (p < 0.001). In regards to the confinement status, 37% of the dogs were found free-roaming, 39.5% were found chained or leashed to a wall or post, while 18.7% were confined in a kennel or inside a household. More details on dog population based on sex, age, ownership, vaccination, neutering, and confinement can be found in Fig. 3. Proportion tables for all the relationships drawn in this section are provided in Additional file 1: Figure S1.

Fig. 3
figure 3

Descriptive tables on dog demographics. The tables show the counts and proportions of dogs in each category, including the number of dogs whose category was not recorded

From all the dogs brought to SP, 1961 (44.8%) had never been vaccinated before or their vaccine was overdue. 2421 (55.2%) were found to be already vaccinated, having been immunized at a Dogstar clinic (251, 10.4%), by a municipal veterinarian (2059, 85.6%), or by a private veterinarian (89, 3.7%).

Analysis of dog neuter status

Only a minority (1377, 26.6%) of the dogs encountered during the D2D campaign were neutered. The neuter status of 186 dogs was not recorded. A negative relationship between ownership and neutering (p 0.001) was detected, as dogs owned were typically entire. Female animals were more likely neutered (p < 0.001), since male dogs were much more frequently owned and entire. Young animals were seldom found neutered (p < 0.001). Entire dogs represented a 71.5% of the owned population. When owners were asked whether they were interested in getting their pet sterilised, 23.7% stated that they would like to have their dog neutered, while 69.4% reported that they did not want their dog to be neutered. Table 2 shows their responses when asked why they opposed neutering. Considering neutering as an unnecessary procedure was the most common answer (42.6%), followed by a desire to breed the dog (13.4%).

Table 2 Reasons given by dog owners against the sterilization of their dogs

Identification of barriers of attendance to SP

The barriers of attendance to SP were investigated by a regression analysis. The results from univariable logistic regression analysis performed to assess each of the variables considered can be found in Additional file 2: Figure S2. The final model represented a compromise between predictive power (greater AUC), and simplicity (less predictors). This model showed that the odds of attendance to SP decrease with increasing distance from household to SP, in young dogs, and in dogs with poor body condition (under/overweight). Odds of attendance were higher in male and/or neutered dogs. Infant hunger showed a very slight negative association with attendance to SP, although not statistically significant (p 0.065). The odd ratios for predictors related to unknown categories (unknown sex and unknown age) were not interpreted due to their ambiguous and non-defined nature. This analysis is represented in Fig. 4, showing the odds ratio and CI for each of the predictive variables. A complete table including this information in a numerical format and the P values for each predictor can be found in Table 3. The model presented a moderate predictive power (AUC = 0.642).

Fig. 4
figure 4

Graphical representation of the chosen regression model. The 95% Confidence Interval (CI) is represented as horizontal bars. The value for the odds ratio is indicated above the 95% CI. A positive relationship between the variable and the response (attendance) is coloured in red (odds ratio > 1), while a negative relationship is coloured in blue (odds ratio < 1). Baseline categories included

Table 3 Numerical results of the regression model

Table 4 shows the most common responses when owners of unvaccinated dogs were asked why they did not take their dogs to a SP. Most common reasons reported were unawareness (37.9%), unavailability (29.3%), and problems handling the dog (10.2%). The dog being too young (5.3%) and excessive distance to the SP (4.5%) were also responses commonly given. Responses regarding rabies vaccination as unnecessary or harmful were uncommon.

Table 4 Reasons given by dog owners explaining failure to attend to a Static Point


The vaccination effort described in this report demonstrates the feasibility of mass dog vaccination campaigns with high coverage in Sri Lanka. A total of 7804 dogs were vaccinated during the course of the campaign, 3422 at SP locations and 4382 during the D2D campaign, over the 33 wards of Negombo.

The coverage obtained for the Negombo working area achieved the 70% minimum coverage recommended by the WHO [9], and surpassed the 48% vaccination coverage estimated for all of Sri Lanka in previous campaigns [16]. However, reliable survey data is lacking for 9 of the wards. In future campaigns, this can be improved by increasing the number of people trained to assess vaccination coverage in case of unavailability of the primary supervisor. Vaccination coverage over 70% was achieved in one out of the two wards where the initial survey was not satisfactory, after returning to both for a second round of vaccination and survey. The information gathered during the campaign has provided extensive insights on the dog population of Negombo.

The estimated dog:human ratio from Negombo (1:17) is much lower than the ratio reported in Mirigama in 2000 (1:4.6) [22], the ratio estimated by Knobel et al. in 2005 for urban settings in the Asian continent (1:7.5) [46], and the nation-wide estimates of 2014 (1:6.7) [16]. This disparity could be due to the large timespan between the studies and this report, the different human population densities, and the recent increase of population control efforts, provided by organizations such as the Dogstar Foundation [20], which regularly arranges mobile sterilization units. The findings on dog demographics show many similarities with previous studies on the dog population of Mirigama [22] and Colombo [47]. Male dogs were overrepresented in the owned dog populations of the three cities, representing also a bigger proportion of the vaccinated population in the Mirigama study. As in Negombo, this study also reported high proportions of owned dogs among the studied population, and lower proportions of young dogs (less than 1 year of age). The Colombo study also reported higher proportions of neutered females, compared to males. Conversely, while owned dogs in Negombo were seldom neutered, high percentages of neutered owned dogs were described in Colombo. In conclusion, very similar demographics could be found from the studies in the three cities, and although it may not be representative across the nation, the characterization of the dog population in each setting is a vital step for the development of successful vaccination campaigns.

As a further step, the regression model was able to identify several factors associated with a failure to attend a static point, notably the distance between household and the SP, and the young age of the dog. It is hoped that this information can be used to inform the roll out of future vaccination strategies in Sri Lanka. In a similar study using the same technique for a campaign in Malawi, Mazeri et al. [11] found that there was a decreasing likelihood of attending a SP with increasing distance from home to said SP, in dogs of young age, and in animals with poor health. These results parallel the findings from the Negombo study. Distance was also the sixth most common reason dog owners gave for not attending a SP, and the deterring effect of distance between domicile and healthcare facilities has been reported in previous studies [48,49,50,51]. The young age of the dog appears as the fifth most common quoted reason for not attending the SP. According to the regression analysis, odds of dogs under one year of age to be brought to a SP were 3.7 times lower compared to adult dogs. This finding has been observed in other studies where a misleading “too young to vaccinate” attitude prevails, leading to dogs remaining unimmunised in the population for several months [52]. The reduction in the odds of attendance with dogs with poor Body Condition Score (BCS) might be linked to owner care as a factor. In addition, the most common reasons the owners gave for not attending a SP were unawareness, unavailability and dog handling problems. These three reasons were also the most common answers to that same question in studies performed in different settings [11, 53,54,55,56]. However, answers conveying distrust in the campaign or the procedure, and/or lack of rabies disease knowledge, which were common in some of these studies, represented only a small percentage of the reasons provided by Negombian owners. This might be related to the high level of rabies awareness of the Sri Lankan population [18]. Increased focus into awareness prior to the vaccination and the arrangement of educational campaigns on responsible dog ownership and handling techniques could potentially allow more owners, willing to vaccinate their dogs, to attend SP campaigns. Increasing the number of SP deployed across the region might increase attendance rates by reducing distance to households. In addition, efforts to increase community awareness about the fact that puppy vaccination is a safe and necessary procedure are also required in order to effectively decrease the number of susceptible dogs in the population.

The 2016 vaccination campaign provided large amounts of data on the local dog population, enabling to investigate the local dog demographics. Moreover, the GPS coordinate recordings allowed the inclusion of additional variables of geospatial nature, with the goal of enriching the regression analysis and increasing its explanatory power. The unexplained variation reducing the predictability of the model could be due to the poor resolution of the geographic information data available for Sri Lanka, or due to the lack thereof, preventing the inclusion of other potentially explanatory variables, such as poverty data. Vaccination coverage results from 9 wards were unusable due to data collection unreliability, losing the ability to completely assess the breadth of the city-wide campaign. D2D data regarding the vaccination status of dogs shows similar percentages of successfully vaccinated dogs in these 9 wards (65%, CI 63.4–66.3) compared to the rest of the wards (72%, CI 69.2–75.3), and dog demographics were nearly indistinguishable between them. For these reasons, we expect the vaccination efforts to be successful and homogeneous throughout the working area, even though vaccination coverages could not be produced for these 9 wards due to inaccurate data. We acknowledge this fact as a limitation of our study, and suggest increasing overall staff competency appointing and training additional back-up supervisors for the different campaign stages. The vaccination coverages reported in this study might be underestimates, as collars and paint, used to mark vaccinated dogs, could have been lost or worn off respectively due to the timespan between stages. Lastly, any survey method chosen to assess vaccination coverage includes the possibility of missing some dogs in the population. In this study, this was minimized by the use of small vehicles and the ability to observe the urban layout through the use of the Mission Rabies App, identifying inspection routes dynamically to ensure thorough coverage and preventing path retracing.


This report describes a high numbers, high coverage mass dog vaccination campaign in Negombo, Sri Lanka. Over 7800 dogs were vaccinated, achieving a coverage of 75.8%. The data collected during the campaign allowed to study the local dog population and identify factors influencing attendance to static vaccination points. This campaign serves as a working template that can be used to develop dog vaccination strategies in other parts of the country. In addition, the insights obtained on dog ecology and barriers of attendance can be used to increase the cost-efficiency of said campaigns, facilitate their development, and improve their coverage.

Availability of data and materials

The vaccination and survey datasets generated and analysed during the current study are not publicly available. This data is held by Mission Rabies. GIS datasets are available from the respective sources as described in the “Materials and Methods: Data sources: Other sources” section.



Akaike information criterion


Area under the receiving operator curve


Body condition score


Confidence Interval




Global positioning system


Post-exposure prophylaxis


Static point


World Health Organization


  1. Hampson K, Coudeville L, Lembo T, Sambo M, Kieffer A, Attlan M, Barrat J, Blanton J, Briggs D, Cleaveland S, Costa P, Freuling C, Hiby E, Knopf L, Leanes F, Meslin F, Metlin A, Miranda M, Müller T, Nel L, Recuenco S, Rupprecht C, Schumacher C, Taylor L, Vigilato M, Zinsstag J, Dushoff J. Estimating the global burden of endemic canine rabies. PLOS Negl Trop Dis. 2015;9(4):e0003709.

    Article  Google Scholar 

  2. World Health Organization, Rabies fact sheet. (2017). .

  3. Hemachudha T, Ugolini G, Wacharapluesadee S, Sungkarat W, Shuangshoti S, Laothamatas J. Human rabies: neuropathogenesis, diagnosis, and management. Lancet Neurol. 2013;12:498–513.

    Article  Google Scholar 

  4. Fooks A, Banyard A, Horton D, Johnson N, McElhinney L, Jackson A. Current status of rabies and prospects for elimination. Lancet. 2014;384:1389–99.

    Article  Google Scholar 

  5. Hampson K, Dushoff J, Cleaveland S, Haydon D, Kaare M, Packer C, Dobson A. Transmission dynamics and prospects for the elimination of canine rabies. PLoS Biol. 2009;7:e1000053.

    Article  Google Scholar 

  6. Cleaveland S, Kaare M, Knobel D, Laurenson MK. Canine vaccination – providing broader benefits for disease control. Vet Microbiol. 2006;117:43–50.

    Article  CAS  Google Scholar 

  7. Lembo T, Hampson K, Kaare M, Ernest E, Knobel D, Kazwala R, Haydon D, Cleaveland S. The feasibility of canine rabies elimination in Africa: dispelling doubts with Data. PLoS Negl Trop Dis. 2010;4:e626.

    Article  Google Scholar 

  8. Morters M, McKinley T, Horton D, Cleaveland S, Schoeman J, Restif O, Whay H, Goddard A, Fooks A, Damriyasa I, Wood J. Achieving population-level immunity to rabies in free-roaming dogs in Africa and Asia. PLoS Negl Trop Dis. 2014;8:e3160.

    Article  Google Scholar 

  9. WHO Expert Consultation on Rabies: second report (2013). Accessed 16 Nov 2018.

  10. Coleman P, Dye C. Immunization coverage required to prevent outbreaks of dog rabies. Vaccine. 1996;14:185–6.

    Article  CAS  Google Scholar 

  11. Mazeri S, Gibson A, Meunier N, Bronsvoort B, Handel I, Mellanby R, Gamble L. Barriers of attendance to dog rabies static point vaccination clinics in Blantyre, Malawi. PLoS Negl Trop Dis. 2018;12:e0006159.

    Article  Google Scholar 

  12. Tenzin T, Ahmed R, Debnath N, Ahmed G, Yamage M. Free-roaming dog population estimation and status of the dog population management and rabies control program in Dhaka City, Bangladesh. PLoS Negl Trop Dis. 2015;9:e0003784.

    Article  Google Scholar 

  13. Gibson A, Ohal P, Shervell K, Handel I, Bronsvoort B, Mellanby R, Gamble L. Vaccinate-assess-move method of mass canine rabies vaccination utilising mobile technology data collection in Ranchi, India. BMC Infect Dis. 2015;15(1).

  14. Léchenne M, Oussiguere A, Naissengar K, Mindekem R, Mosimann L, Rives G, Hattendorf J, Moto D, Alfaroukh I, Zinsstag J. Operational performance and analysis of two rabies vaccination campaigns in N’Djamena, Chad. Vaccine. 2016;34:571–7.

    Article  Google Scholar 

  15. Gibson A, Handel I, Shervell K, Roux T, Mayer D, Muyila S, Maruwo G, Nkhulungo E, Foster R, Chikungwa P, Chimera B, Bronsvoort B, Mellanby R, Gamble L. The vaccination of 35,000 dogs in 20 working days using combined static point and door-to-door methods in Blantyre, Malawi. PLoS Negl Trop Dis. 2016;10:0004824.

    Google Scholar 

  16. Health Statistical Bulletin. (2015). Accessed 16 Nov 2018.

  17. Abela-Ridder B, Lionel Harischandra P, Gunesekera A, Janakan N, Gongal G. Sri Lanka takes action towards a target of zero rabies death by 2020. WHO South-East Asia. J Public Health. 2016;5(2):113.

    Google Scholar 

  18. Kumarapeli V, Awerbuch-Friedlander T. Human rabies focusing on dog ecology — a challenge to public health in Sri Lanka. Acta Trop. 2009;112:33–7.

    Article  Google Scholar 

  19. Mission Rabies Sri Lanka vaccination campaigns. (2017). Accessed 16 Nov 2018.

  20. Dogstar Foundation. (2018). Accessed 16 Nov 2018.

  21. Census of Population and Housing of Sri Lanka (2012). .

  22. Matter H, Wandeler A, Neuenschwander B, Harischandra L, Meslin F. Study of the dog population and the rabies control activities in the Mirigama area of Sri Lanka. Acta Trop. 2000;75:95–108.

    Article  CAS  Google Scholar 

  23. Gibson A, Mazeri S, Lohr F, Mayer D, Bailey J, Wallace R, Handel I, Shervell K, Bronsvoort B, Mellanby R, Gamble L. One million dog vaccinations recorded on mHealth innovation used to direct teams in numerous rabies control campaigns. PLOS ONE. 2018;13:e0200942.

    Article  Google Scholar 

  24. World Weather Online. (2016). Accessed 16 Nov 2018.

  25. WorldPop. (2015). Accessed 16 Nov 2018.

  26. Socioeconomic Data and Applications Center. (2018). Accessed 16 Nov 2018.

  27. World Bank Group. (2012). Accessed 16 Nov 2018.

  28. Kahle D, Wickham H. Ggmap: spatial visualization with ggplot2. R Journal. 2013;5:144–61.

    Article  Google Scholar 

  29. Open Street Map. (2018). Accessed 16 Nov 2018.

  30. R Core Team (2017). R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria.

  31. Ordoñez C, Martínez J, Rodríguez-Pérez J, Reyes A. Detection of outliers in GPS measurements by using functional-Data analysis. J Surv Eng. 2011;137:150–5.

    Article  Google Scholar 

  32. Hennig, C. (2018). Fpc: flexible procedures for clustering. R package version 2.1–11.

  33. Bivand, R. and Rundel, C. (2017). Rgeos: Interface to geometry engine - open source ('GEOS'). R package version 0.3–26.

  34. Bivand, R., Keitt, T., and Rowlingson, B. (2018). Rgdal: bindings for the 'Geospatial' Data abstraction library. R package version 1.2–18.

  35. Chapman D. The estimation of biological populations. Ann Math Stat. 1954;25:1–15.

    Article  Google Scholar 

  36. Tyers, M. (2016). Recapr: estimating, testing, and simulating abundance in a mark-recapture. R package version 0.3.9002.

  37. Sunil Arya, David Mount, Samuel E. Kemp and Gregory Jefferis (2018). RANN: fast nearest neighbour search (wraps ANN library) using L2 metric. R package version 2.6.

  38. Robert J. Hijmans (2017). Raster: geographic Data analysis and modeling. R package version 2.6–7.

  39. Box G, Tidwell P. Transformation of the independent variables. Technometrics. 1962;4:531.

    Article  Google Scholar 

  40. National Meteorological Library and Archive Fact sheet 3. – Water in the atmosphere (2012). Accessed 16 Nov 2018.

  41. Mount J, Zumel N. (2018). Vtreat: a statistically sound 'data.frame' processor/conditioner. R package version 1.3.0.

  42. Max Kuhn. Contributions from Jed Wing, Steve Weston, Andre Williams, Chris Keefer, Allan Engelhardt, Tony Cooper, Zachary Mayer, Brenton Kenkel, the R Core Team, Michael Benesty, Reynald Lescarbeau, Andrew Ziem, Luca Scrucca, Yuan Tang, Can Candan and Tyler Hunt (2018). Caret: classification and regression training. R package version 6.0–80.

  43. Stone M. Cross-Validatory choice and assessment of statistical predictions. J R Stat Soc Ser B Methodol. 1974;36:111–47.

    Google Scholar 

  44. Akaike H. A new look at the statistical model identification. IEEE Trans Autom Control. 1974;19:716–23.

    Article  Google Scholar 

  45. Sing T, Sander O, Beerenwinkel N, Lengauer T. ROCR: visualizing classifier performance in R. Bioinformatics. 2005;21:7881.

    Article  Google Scholar 

  46. Knobel D, Cleaveland S, Coleman P, Fevre E, Meltzer M, Miranda M, Shaw A, Zinsstag J, Meslin F. Re-evaluating the burden of rabies in Africa and Asia. Bull World Health Organ. 2005;83:360–8.

    PubMed  PubMed Central  Google Scholar 

  47. Seneviratne M, Subasinghe D, Watson P. A survey of pet feeding practices of dog owners visiting a veterinary practice in Colombo, Sri Lanka. Vet Med Sci. 2016;2:106–16.

    Article  Google Scholar 

  48. Stock R. Distance and the utilization of health facilities in rural Nigeria. Soc Sci Med. 1983;17:563–70.

    Article  CAS  Google Scholar 

  49. Müller I, Smith T, Mellor S, Rare L, Genton B. The effect of distance from home on attendance at a small rural health Centre in Papua New Guinea. Int J Epidemiol. 1998;27:878–84.

    Article  Google Scholar 

  50. Buor D. Analysing the primacy of distance in the utilization of health services in the Ahafo-Ano south district, Ghana. Int J Health Plann Manag. 2003;18:293–311.

    Article  Google Scholar 

  51. Kumar S, Dansereau E, Murray C. Does distance matter for institutional delivery in rural India? Appl Econ. 2014;46:4091–103.

    Article  Google Scholar 

  52. Morters M, McNabb S, Horton D, Fooks A, Schoeman J, Whay H, Wood J, Cleaveland S. Effective vaccination against rabies in puppies in rabies endemic regions. Vet Rec. 2015;177:150.

    Article  CAS  Google Scholar 

  53. Muthiani Y, Traoré A, Mauti S, Zinsstag J, Hattendoorf J. Low coverage of central point vaccination against dog rabies in Bamako, Mali. Prev Vet Med. 2014;120:203–9.

    Article  Google Scholar 

  54. Bardosh K, Sambo M, Sikana L, Hampson K, Welburn S. Eliminating rabies in Tanzania? Local understandings and responses to mass dog vaccination in Kilombero and Ulanga districts. PLoS Negl Trop Dis. 2014;8:e0002935.

    Article  Google Scholar 

  55. Castillo-Neyra R, Brown J, Borrini K, Arevalo C, Levy M, Buttenheim A, Hunter G, Becerra V, Behrman J, Paz-Soldan V. Barriers to dog rabies vaccination during an urban rabies outbreak: qualitative findings from Arequipa, Peru. PLoS Negl Trop Dis. 2017;11:e0005460.

    Article  Google Scholar 

  56. Mulipukwa C, Mudenda B, Mbewe A. Insights and efforts to control rabies in Zambia: evaluation of determinants and barriers to dog vaccination in Nyimba district. PLoS Negl Trop Dis. 2017;11:e0005946.

    Article  Google Scholar 

Download references


We are grateful to Dogs Trust Worldwide for funding the work of Mission Rabies and Dogstar in Sri Lanka. We are grateful to MSD Animal Health for donating all Nobivac® Rabies vaccines used for the project. We thank all the volunteering personnel and the local veterinarians who participated in the execution of the campaign.


Barend M. deC. Bronsvoort was supported through BBSRC through the Institute Strategic Programme funding (BB/J004235/1 and BB/P013740/1). The funders had no role in study design, data collection and analysis, decision to publish or preparation of the manuscript.

Author information

Authors and Affiliations



ADG, LG, JLBB, SG and MG developed, carried out and coordinated the vaccination campaign, supervised the data collection and compiled the datasets used in the study. LG acquired the funding for the campaign and study. CSS, ADG and SM conceptualized the study. CSS and SM performed the statistical analyses described in the methodology. BMB, IGH, RJM and SM provided supervision during the analyses. CSS created the original draft of the article. CSS, IGH, RJM and SM reviewed and edited the manuscript. All authors mentioned read and approved the final manuscript.

Corresponding author

Correspondence to Carlos Sánchez-Soriano.

Ethics declarations

Ethics approval and consent to participate

The analysis of the data gathered in the Mission Rabies 2016 Negombo vaccination campaign was approved by the Veterinary Ethical Review Committee at the University of Edinburgh, with reference 64/15. The vaccination program and data collection constituted a non-research public health campaign, and were approved by the Municipal Council of Negombo and the Commissioner. Informed verbal consent was obtained from the person presenting the dog for vaccination. When an owner could not be identified, dogs were vaccinated in accordance with the Government Public Health protocol.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Figure S1.

Proportions of the dog population according to different sets of variables. Each 2 by 2 table describes the proportions of the dog population according to each pair of variables. The P value obtained for the test of independence between both variables is shown.

Additional file 2: Figure S2.

Results from the univariable regression analysis. Variables whose P value was under 0.2 were considered for the model and included in the models created mostly by forward selection. The age category was divided into adult (over one year of age) and young (less than one year) for the final model due to low counts of puppy dogs.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Sánchez-Soriano, C., Gibson, A.D., Gamble, L. et al. Development of a high number, high coverage dog rabies vaccination programme in Sri Lanka. BMC Infect Dis 19, 977 (2019).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: