Risk mapping of Rinderpest sero-prevalence in Central and Southern Somalia based on spatial and network risk factors
© Ortiz-Pelaez et al; licensee BioMed Central Ltd. 2010
Received: 20 July 2009
Accepted: 28 April 2010
Published: 28 April 2010
In contrast to most pastoral systems, the Somali livestock production system is oriented towards domestic trade and export with seasonal movement patterns of herds/flocks in search of water and pasture and towards export points. Data from a rinderpest survey and other data sources have been integrated to explore the topology of a contact network of cattle herds based on a spatial proximity criterion and other attributes related to cattle herd dynamics. The objective of the study is to integrate spatial mobility and other attributes with GIS and network approaches in order to develop a predictive spatial model of presence of rinderpest.
A spatial logistic regression model was fitted using data for 562 point locations. It includes three statistically significant continuous-scale variables that increase the risk of rinderpest: home range radius, herd density and clustering coefficient of the node of the network whose link was established if the sum of the home ranges of every pair of nodes was equal or greater than the shortest distance between the points. The sensitivity of the model is 85.1% and the specificity 84.6%, correctly classifying 84.7% of the observations. The spatial autocorrelation not accounted for by the model is negligible and visual assessment of a semivariogram of the residuals indicated that there was no undue amount of spatial autocorrelation. The predictive model was applied to a set of 6176 point locations covering the study area. Areas at high risk of having serological evidence of rinderpest are located mainly in the coastal districts of Lower and Middle Juba, the coastal area of Lower Shabele and in the regions of Middle Shabele and Bay. There are also isolated spots of high risk along the border with Kenya and the southern area of the border with Ethiopia.
The identification of point locations and areas with high risk of presence of rinderpest and their spatial visualization as a risk map will be useful for informing the prioritization of disease surveillance and control activities for rinderpest in Somalia. The methodology applied here, involving spatial and network parameters, could also be applied to other diseases and/or species as part of a standardized approach for the design of risk-based surveillance activities in nomadic pastoral settings.
Somalia's livestock production sector accounts for at least 40% of the gross domestic product (GDP) with 55% of the human population being directly involved in the rearing of livestock . Pastoral movements across Somalia's borders for the purposes of grazing livestock and livestock-related trade have occurred for centuries , resulting in seasonal movement patterns of herds/flocks in search of water and pasture. Moreover, in contrast to most pastoral systems, which are normally aimed at household subsistence, the Somali livestock production system is oriented towards domestic trade and export . Livestock are shipped to various countries in the Arabian Peninsula, and trekked or transported to markets in Kenya, Djibouti, and Ethiopia . In 2007 some 1,639,625 heads of cattle, sheep, goats and camels were exported through Bossasso port and a total of 1,633,793 through Berbera port , the two major export markets in the country.
Despite the collapse of the Somali government in 1991 and the lack of continued delivery of public services, an informal system of stateless order, social trust and an informal economy have allowed rural and urban populations to survive unfavourable economic and political circumstances . But the Somali livestock industry is therefore now even more vulnerable to the introduction of bans by importing countries for two reasons. One being the stringent measures for livestock trade specified under the Sanitary and Phytosanitary (SPS) Agreement of the World Trade Organization http://www.wto.org, of which Somalia is not a member, and the other the poor standards of veterinary services and the absence of control measures for fighting trans-boundary animal diseases. The high mobility of the livestock population poses an additional challenge for the control and establishment of credible certification systems for the major trans-boundary diseases occurring in Somalia. An understanding of the aggregation/dispersion mechanisms and contact structure of the Somali livestock potentially could assist in setting up appropriate spatial risk-based surveillance activities and control measures that may lead to the establishment of an internationally accepted certification system for this nomadic pastoral livestock production system.
Data on the geographical patterns of human and animal population distributions are now commonly considered in epidemiological studies. High density of susceptible populations has been shown to be a key factor in the transmission of infectious diseases in humans . In the case of animal populations, spatial proximity is closely linked to the transmission of many infectious diseases . The geo-referencing of animals kept in herds through the point location of the farm is relatively simple and cost-effective . However in settings where mobility is a major feature of the husbandry systems the geographical location of animals is no longer a discrete entity hence alternative approaches are required for the collection and analysis of their spatial data.
Network data are not usually linked to geographical data, and social network analysis rarely considers the spatial configuration of the links  apart from the pure visualization of the networks. Definitions of links within networks based on spatial criteria such as the distance between pairs of nodes, usually estimated using the Euclidean distance ("as the crow flies") are rarely available in the literature. For example, Webb  made the assumption in her analysis that there was a link between two farms if their postcodes were less than 25 km apart, as a proxy for the catchment area of farm animal veterinary practices. Dent et al.  considered poultry premises to be linked if they were within 3 km distance from each other.
The purpose of epidemiological modelling of spatial data is to explain or predict the occurrence of disease [12–18]. The production of risk maps can guide decision makers provided that the underlying assumptions of uncertainty and variability are exposed . Using the data from a previous rinderpest survey conducted in the Central and Southern regions of Somalia, it was possible to explore the topology of a contact network of cattle herds based on a spatial proximity criterion and other attributes related to cattle herd dynamics. The objective of the current study is to integrate spatial mobility and other attributes with GIS and network approaches in order to investigate their effect on disease presence and develop a predictive spatial model of presence of rinderpest. The risk map generated by the predictive model can be used to inform the design of risk-based surveillance activities of rinderpest by identifying high-risk point locations where to prioritize disease surveillance and control activities. Surveillance efforts could be targeted at high-risk areas which would lead to more effective use of the scarce veterinary resources in Somalia, and ultimately, contribute to the establishment of an internationally credible certification system for livestock export by nomadic pastoral systems. The standardisation of the proposed methodology for cattle will serve as a model for other livestock species (e.g. small ruminants and camelids) and other diseases.
At selected sampling sites, a minimum number of 15 eligible animals were sampled per herd and questionnaires were administered to the livestock owners of the sampled herds in order to collect data on species, herd/flock size and location during each climatic season over the two years prior to the date the survey was conducted in that point and livestock markets (primary and secondary) where the pastoralists sold their animals during the same period. A total of 9216 serum samples were collected from cattle aged 1 to 3 years at 562 sampling sites. In addition, 1071 sera were collected from cattle older than 3 years at 58 sampling sites. All serum samples collected were tested for the presence of RP antibodies using a RP Competitive Enzyme Linked Immunosorbent Assay (C-ELISA) directed against the H protein of the virus . All samples having a percentage of inhibition (PI) above 50% were considered positive in accordance with the OIE recommendation for testing large numbers of sera.
Attributes for the 562 point locations were generated through spatial queries and analysis in ArcGIS 9.2. (© ESRI), using as unique identifier a point location reference number linked to the projected coordinates. Seven attributes were extracted for each point location. (1) The type of point (categorical): villages, towns , water points , primary/secondary markets , and grazing points . (2) Proximity to roads (binary): less than 500 meters from a main road. (3) Proximity to trade routes (binary): the locations of three major trade routes were recorded in southern and central Somalia by Terra Nuova  by use of a radio collar GPS fitted to an animal within the traded herds being trekked from secondary markets to export points. The proximity was based on whether a point's estimated home range as described above intersected the trade routes. (4) Speed at trade route (numeric): the speed of the movement at each trade route point was calculated for segments of time and distance between points along the trade routes derived from the spatial join of collar (cattle) and drover (handler) GPS data and then aggregated based on drover GPS locations. If the point was linked to a trade route, speed in km/h was included. (5) Home range radius (numeric): radius of the home range assuming circular shape (method described above). (6) Herd density (numeric): herd density values assigned to the survey points through intersection of the points with a distribution map of herd density generated through moving average neighbourhood analysis  from livestock census data . (7) Herd size (numeric): herd size values assigned to the survey points through intersection of the points with a distribution map of herd sizes .
Spatial logistic model
where β0 is the intercept, x i the covariates, β i the regression coefficients, ξ i normally distributed non-correlated residuals and S i the spatial component (spatial autoregressive coefficient), which is assumed to have zero mean, variance σ2 and an isotropic covariance function ρ (d ij θ), where θ is the range coefficient and measures the decay of the spatial autocorrelation and d ij is the Euclidean distance between each pair of points.
A generalized linear mixed model with a logit link function was fitted. The outcome variable was the count of rinderpest positive samples expressed as a binomial distribution to account for the variation in the number of animals tested. The following covariates were introduced in the multivariable model as fixed effects: type of point, next to road, linked to trade route, speed at trade route, herd density, herd size, home range radius, degree, betweenness, closeness, inclusion in the giant weak component and clustering coefficient. An exponential spatial covariance structure (exp (-d ij /θ)) was included as a random effect. The analysis was implemented using restricted pseudo likelihood with the GLIMMIX procedure of SAS 9.1 (© 2002-2003 by SAS Institute Inc., Cary, NC, USA). Linearity of the effects was tested by categorising continuous variables into three-category variables. The statistical significance of the fixed effect variables was determined using Type III tests. The ability of the binomial model to predict the outcome was estimated using as cut-off point the mean of the fitted probabilities predicted by the model in order to transform them into a binary outcome . Predicted probabilities below the cut-off value were considered as predicting absence (0) and probabilities above the cut-off value were considered predicting presence (1).
Once the final model was fitted, a semivariogram of the residuals was generated using the ArcGIS 9.2 Geostatistical Analyst extension (© ESRI) for assessing the potential presence of spatial autocorrelation in the model residuals. The logistic regression model equation was applied to the set of 6176 point locations using their specific covariate patterns. Attributes for the significant variables of the final model were extracted for the set of point locations using the same methods as for the survey points. The predicted prevalence using the final model equation for each point location was exported to ArcGIS 9.2 (© ESRI) for visualization of the resulting risk map. A kernel smoothed intensity map was produced with 50 km bandwidth and 5 km grid cells using the Spatial Analyst extension of ArcGIS 9.2 (© ESRI).
Coefficients of the final binomial generalized mixed model (GLM) with standard errors, odds ratios and 95% confidence intervals
Odds Ratio (95% CI)
P < 0.001
Home range radius (in Km)
P < 0.001
P < 0.001
Clustering coefficient (cc1)
P < 0.001
Predicted outcomes from the final binomial generalized mixed model (GLM) compared to the observed data
The network contains 6738 nodes and 319,887 links with a density of 1.4% and an average degree of 224.68 (range: 1-605). The average distance between reachable pairs is 7.7 links with the maximum distance between two points being 24. There is a giant weak component containing 2950 nodes (43.78%). Within the network the 562 points of the survey had an average degree of 169 (STD: 118 Range: 2-605), average betweenness of 0.001 (STD: 0.004 Range: 0-0.05), average closeness of 0.12 (STD: 0.023 Range: 0.08-0.17) and average clustering coefficient 0.67 (STD: 0.01 Range: 0.3-1).
This study has applied spatial and network analysis techniques to a combination of observed and estimated data with the objective to produce a predictive model for the presence of serological evidence of rinderpest in the southern and central regions of Somalia. Incorporating the network framework into GIS has been attempted before mainly in geography and social sciences. For example, Faust et al.  showed how social relations can be determined by geographical location and other spatially-related attributes. In the context of disease transmission, Rothenberg et al.  demonstrated that social and geographical distance can support transmission of disease. Following the same rationale the main hypothesis tested in this study was whether the presence of rinderpest is associated with the network centrality related to its geographical location and/or to other attributes of the geographical locations. The results have shown the significance of both types of factors. One of these is the degree of interactions between herds due to overlapping home ranges in neighbouring areas, expressed by the clustering coefficient of the nodes derived from the spatial contact network. The other is the home range radius or area where nomadic herds graze throughout the year and the herd density are both significant risk factors. The combination of the three significant factors reveals the importance for the presence of rinderpest of local neighbourhood and geographic compactness interacting in densely populated areas, an effect that has been reported before for human diseases such as STDs and HIV-AIDS [35, 6].
Network parameters that take into account paths between all possible pairs of nodes such as betweenness and closeness did not appear to be significantly associated with the presence of disease. Large intermediaries or globally relevant nodes in the network represented by high betweenness and closeness are not roles that increase the risk of disease presence in the population of the survey, nor is the membership in within-network structures such as the giant weak component. Although two nodes/points far apart can be linked via a path of overlapping home ranges, being reachable by many other nodes in the network, in between many other shorter paths, or included in a connected subpopulation of nodes, the results suggest that the amount of overlap between neighbouring herd home ranges is the only network feature that substantially increases the risk of rinderpest. The spread of rinderpest to distant areas may follow other paths unaccounted for in this study. For example, the effect of markets has not been addressed in this analysis. It would be advised to investigate additional potential risk factors such as presence of or distance to local or regional markets. Nor has the possible risk attributed to long-distance movements via trade routes been assessed given the presence in the dataset of only three trade routes linked to each other and located in the southern districts of the study area. Incorporating more comprehensive data on putative risk factors such as the movement of the herds to markets and borders with neighbouring countries would allow testing these hypotheses.
The inclusion of variables with a strong spatial component in the multivariable analysis required the application of an analytical method that could account for the dependence of the network data due to both the inherent nature of the network data and the spatial definition of the link. Thus, the general assumption of dependence within the network data analyzed is confirmed by both the network parameterization and the spatial nature of the link. In fact, the spatial component of the binomial generalized linear mixed model, the spatial autoregressive coefficient, is similar to the local and interaction effects in network models  in what Anselin  called the family of autocorrelation models. In both cases the specifications of the chosen weight matrix are important for testing the statistical significance of a particular parameter. In this case, the spatial correlation function accounts at the same time for the similarity of home range radius and herd density between proximal points and the similarity of clustering coefficient of points linked in the network.
The internal validation of the final model showed a very high predictive ability. Its application to a different dataset has produced risk maps of point locations and areas based on the model estimates. In general, the location of the high-risk points agrees with the observed points positive to rinderpest in the original dataset, which suggests that the areas affected by rinderpest are characterised by large home ranges, large herd density and strong local interaction between neighbouring nomadic herds. However there is an area in Gedo region close to the border with Ethiopia and Somalia where the risk map predicts points of high prevalence whereas the kernel smoothed intensity map does not. This is due to the presence in that area of a small number of points with very large home ranges and high herd density that increase the risk at point level but dilute the risk at area level. This is a methodological constraint rather than an incoherent result given the difficulty of obtaining stable risk estimates if intensities vary substantially across an area .
Both, the point and area risk maps, identified a high risk area that the survey sampled at low intensity and did not detect rinderpest. It is located in the northern part of the Bay region close to the border with Bakool. These areas of predicted high risk and for which disease data is sparse or absent and were not well covered by the survey should therefore be prioritised for surveillance purposes.
The risk maps have identified areas in the Central and Southern regions of Somalia where veterinary authorities should concentrate their efforts for implementation of surveillance activities. It is recommended to validate these results with local experts and through further surveys to ascertain the rinderpest status in areas of predicted high risk where the disease is not currently known to be present. In this respect similar serological studies subsequently conducted in the study area detected rinderpest antibodies in the high-risk areas identified under this study. However, the observed sero-prevalence in these areas was observed to decline over time [38–40]. Similar sero-prevalence studies conducted in neighbouring areas of Kenya and Ethiopia in subsequent years also showed a very low sero-prevalence or the absence of rinderpest antibodies . The progressive decline in rinderpest antibody prevalence is consistent with virus extinction. The verification of freedom from infection from the Somali Eco-System is currently in progress.
This methodology can also be used for other diseases and/or species as part of a standardized approach for the design of risk-based surveillance activities in nomadic pastoral settings. The detection of high risk areas and hot spots will assist veterinary authorities in strengthening the veterinary surveillance in such areas in order to apply control measures at the local level and to minimize the potential spread of the disease to other regions. The application of disease control measures based on geographical factors and applied at local/compartmental level could be a starting point for the establishment of an internationally recognised certification system in settings where large-scale surveillance and control activities cannot be implemented.
This study demonstrates the potential for the integration of different data sources and analytical methods for the development of a spatial and network predictive model of the presence of rinderpest in central and southern Somalia. The identification of point locations and areas with high risk of presence of rinderpest and their spatial visualization as a risk map can assist in prioritizing disease surveillance and control activities for rinderpest in Somalia.
Although considerable progress on RP eradication has been achieved in the study area since this study was conducted, the value and application of the applied methodology remains of interest for the investigation of livestock diseases in nomadic and semi-nomadic pastoral systems. This methodology can be used as part of a standardized approach for developing new decision-support tools for the design of risk-based surveillance activities in nomadic pastoral settings where the particular husbandry system, the lack of infrastructure and the inadequacy of the available veterinary services prevent the implementation of comprehensive control policies.
The paper originates from activities carried out within the frame of the "Improvement and Diversification of Somali Livestock Trade and Marketing", project funded by the EC and Royal Danish Embassy and implemented by Terra Nuova and ILRI.
The content of this publication is the sole responsibility of the authors and in no way reflects the views of the European Union or/and the Royal Danish Embassy.
- Food and Agriculture Organization (FAO): Somalia, towards a livestock sector strategy. Final Report. 2004, FAO-World Bank Cooperative Programme. European Union, Report No.: 04/001 IC-SOMGoogle Scholar
- Lindley A: Somalia Country Study. A part of the report on Informal Remittance Systems in Africa, Caribbean and Pacific (ACP) countries (Ref: RO2CS008). ESRC Centre on Migration, Policy and Society (COMPAS). Oxford, UK 2005.Google Scholar
- Abdullahi AM: Pastoral Production Systems in Africa. A Study of Nomadic Household Economy and Livestock Marketing in Central Somalia. Farming Systems and Resource Economics in the Tropics (Germany). Wissenschaftsverlag Vauk. Kiel (Germany), 8: 1990, 266.Google Scholar
- Food Security Analysis Unit (FSAU): Livestock export data for Somalia. Food Security Analysis Unit-Somalia, Food and Agriculture Organization. 2008, [http://www.fsausomali.org/fileadmin/uploads/1273_01.pdf]Google Scholar
- Little PD: Somalia: Economy Without State. Bloomington and Indianapolis: Indiana University Press, 2003, 206.Google Scholar
- Rothenberg R, Muth SQ, Malone S, Potterat JJ, Woodhouse DE: Social and geographic distance in HIV risk. Sex Transm Dis. 2005, 32 (8): 506-12. 10.1097/01.olq.0000161191.12026.ca.PubMedView ArticleGoogle Scholar
- Pfeiffer DU, Robinson TP, Stevenson M, Stevens K, Rogers DJ, Clements ACA: Spatial Analysis in Epidemiology. Oxford University Press. Oxford, UK, 2008, 142.View ArticleGoogle Scholar
- Durr PA, Froggatt AEA: How best to geo-reference farms? A case study from Cornwall, England. Prev Vet Med. 2002, 56 (1): 51-62. 10.1016/S0167-5877(02)00123-X.PubMedView ArticleGoogle Scholar
- Faust K, Entwistel B, Rindfuss RR, Walsh SJ, Sawangdee Y: Spatial arrangements of social and economic networks among villages in Nang Rong District, Thailand. Social Networks. 1999, 21: 311-337. 10.1016/S0378-8733(99)00014-3.View ArticleGoogle Scholar
- Webb CR: Farm animal networks: unravelling the contact structure of the British sheep population. Prev Vet Med. 2005, 68 (1): 3-17. 10.1016/j.prevetmed.2005.01.003.PubMedView ArticleGoogle Scholar
- Dent JE, Kao RR, Kiss IZ, Hyder K, Arnold M: Contact structures in the poultry industry in Great Britain: Exploring transmission routes for a potential avian influenza virus epidemic. BMC Veterinary Research. 2008, 4: 27-10.1186/1746-6148-4-27.PubMed CentralPubMedView ArticleGoogle Scholar
- Pfeiffer DU, Duchateau L, Kruska RL, Ushewokunze-Obatolu U, Perry BD: A spatially predictive logistic regression model for occurrence of theileriosis outbreaks in Zimbabwe. Epidemiologie et santé animale. 1997, 31-32: 1212.1-3.Google Scholar
- Baylis M, Mellor PS, Wittman EJ, Rogers DJ: Prediction of areas around the Mediterranean at risk of bluetongue by modelling the distribution of its vector using satellite imaging. The Veterinary Record. 2001, 149: 639-643. 10.1136/vr.149.21.639.PubMedView ArticleGoogle Scholar
- McKenzie JS, Morris RS, Pfeiffer DU, Dymond JR: Application of remote sensing to enhance the control of wildlife-associated Mycobacterium bovis infection. PE&RS - Photogrammetric Engineering & Remote Sensing. 2002, 68 (2): 153-159.Google Scholar
- Tum S, Puotinen ML, Copeman DB: A geographic information systems model for mapping risk of fasciolosis in cattle and buffaloes in Cambodia. Veterinary Parasitology. 2004, 122: 141-149. 10.1016/j.vetpar.2004.03.016.PubMedView ArticleGoogle Scholar
- Pfeiffer DU, Minh PQ, Martin V, Epprecht M, Otte MJ: An analysis of the spatial and temporal patterns of highly pathogenic avian influenza occurrence in Vietnam using national surveillance data. The Veterinary Journal. 2007, 174: 302-309. 10.1016/j.tvjl.2007.05.010.PubMedView ArticleGoogle Scholar
- Gilbert M, Xiao X, Pfeiffer DU, Epprecht M, Boles S, Czarnecki C, Chaitaweesub P, Kalpravidh W, Minh PQ, Otte MJ, Martin V, Slingenbergh J: Mapping H5N1 highly pathogenic avian influenza risk in Southeast Asia. PNAS March 25. 2008, 105 (12): 4769-4774. 10.1073/pnas.0710581105.View ArticleGoogle Scholar
- Boender GJ, Nodelijk G, Hagenaars TJ, Elbers ARW, de Jong MCM: Local spread of classical swine fever upon virus introduction into The Netherlands: Mapping of areas at high risk. BMC Veterinary Research. 2008, 4: 9-10.1186/1746-6148-4-9.PubMed CentralPubMedView ArticleGoogle Scholar
- Pfeiffer DU: Geographical information science and spatial analysis in animal health. GIS and spatial analysis in veterinary science. Edited by: Durr P, Gatrell A. CABI publishing, Wallingford, Oxfordshire, UK, 2004, 303.Google Scholar
- Anon: Somali Livestock Statistics 1988/1989. Ministry of Livestock, Forestry and Range - Department of Planning and Statistics, Mogadishu, Somalia
- Tempia S: The dynamic of Rinderpest in nomadic pastoral systems: the Somali surveillance example. PhD thesis. Colorado State University. Fort Collins-Colorado
- Anderson J, McKey JA, Butcher RN: The use of monoclonal antibodies in competitive ELISA for the detection of antibodies to rinderpest and peste des petits ruminants. The Sero-Monitoring of Rinderpest Through Africa - Phase I. Proceedings of the Final Research Coordination Meeting of IAEA Rinderpest Control Projects. International Atomic Energy Agency, Vienna 1991.Google Scholar
- Africa Data Sampler: A Geo-Referenced Database for All African Countries. World Resources Institute in Collaboration with World Conservation Monitoring Centre and Padco, Inc. CD-ROM 2000.Google Scholar
- UNDOS: Somalia Topomaps. United Nations Development Office for Somalia (UNDOS), Nairobi, CD-ROM 2003.Google Scholar
- White GC, Garrott RA: Analysis of wildlife radiotracking data. Academic Press, New York. New York, USA 1990.Google Scholar
- Carrat F, Valleron AJ: Epidemiologic mapping using the "Kriging" method: application to an influenza-like illness epidemic in France. American Journal of Epidemiology. 1992, 135: 1293-1300.PubMedGoogle Scholar
- Olivier MA, Webster R: Kriging: a method of interpolation for geographical information systems. International Journal of Geographical Information Systems. 1990, 4: 313-332. 10.1080/02693799008941549.View ArticleGoogle Scholar
- Stein ML: Interpolation of Spatial Data - Some Theory for Kriging. Springer-Verlag, Inc., New York 1999.View ArticleGoogle Scholar
- Webster R, Olivier MA, Muir KR, Mann JR: Kriging the local risk of a rare diseases from a register of diagnosis. Geographical Analysis. 1994, 26: 168-185. 10.1111/j.1538-4632.1994.tb00318.x.View ArticleGoogle Scholar
- Lessard P, L'Eplattenier R, Norval RAI, Kundert K, Dolan TT, Croze H: Geographical information systems for studying the epidemiology of cattle disease caused by Theileria parva. The Veterinary Record. 1990, 126: 255-262.PubMedGoogle Scholar
- Cressie NAC: Statistics for spatial data. John Wiley & Sons, New York, USA 1993.Google Scholar
- Terra Nuova: Improvement and Diversification of Somali Livestock Trade and Marketing. Fourth Progress Report. August 2006 to January 2007. Terra Nuova, Nairobi, Kenya 2007.Google Scholar
- Terra Nuova: PanAfrican Control for the Epizootics - Somali Component. Rinderpest Cross-Sectional Sero-Survey in Central and Southern Somalia-Technical Report. Terra Nuova, Nairobi. 2003Google Scholar
- Hardin J, Hilbe J: Generalized Linear Models and Extensions. Stata Press, Texas. USA 2001.Google Scholar
- Zenilman JM, Ellish N, Fresia A, Glass G: The geography of sexual partnerships in Baltimore: Applications of core theory dynamics using a geographic information system. Sex Transm Dis. 1999, 26: 75-81. 10.1097/00007435-199902000-00002.PubMedView ArticleGoogle Scholar
- Leenders R Th. : Modeling social influence through network autocorrelation: constructing the weight matrix. Social Networks. 2002, 24: 21-47. 10.1016/S0378-8733(01)00049-1.View ArticleGoogle Scholar
- Anselin L: Specification tests on the structure of interaction in spatial econometric models. Papers of the Regional Science Association. 1984, 54: 165-182. 10.1007/BF01940131.View ArticleGoogle Scholar
- Terra Nuova: Follow-up investigation of the rinderpest situation in central and southern Somalia. Somali Animal Health Services Project - Technical Report, Terra Nuova, Nairobi 2005.Google Scholar
- Terra Nuova: Participatory disease search and purposive sero-surveillance for rinderpest in central and southern Somalia. Somali Animal Health Services Project - Technical Report, Terra Nuova, Nairobi 2006.Google Scholar
- Terra Nuova: Rinderpest cross-sectional sero-survey in the Somali Eco-System: Somalia component. Somali Animal Health Services Project - Technical Report, Terra Nuova, Nairobi 2006.Google Scholar
- African Union-Interafrican Bureau for Animal Resources (AU-IBAR): Rinderpest situation in the Somali Eco-System. AU-IBAR, Nairobi 2006.Google Scholar