gambiaUTM: Gambia data

gambiaUTMR Documentation

Gambia data


This data-set was used by Diggle, Moyeed, Rowlingson, and Thomson (2002) to demonstrate how the model-based geostatistics framework of Diggle et al. (1998) could be adapted to assess the source(s) of extrabinomial variation in the data and, in particular, whether this variation was spatially structured. The malaria prevalence data set consists of measurements of the presence of malarial parasites in blood samples obtained from children in 65 villages in the Gambia. Other child- and village-level indicators include age, bed net use, whether the bed net is treated, whether or not the village belonged to the primary health care structure, and a measure of 'greenness' using a vegetation index.




A SpatVector , with column pos being the binary response for a malaria diagnosis, as well as other child-level indicators such as netuse and treated being measures of bed net use and whether the nets were treated. The column green is a village-level measure of greenness. A UTM coordinate reference system is used, where coordinates are in metres.

Source For further details on the malaria data, see Thomson et al. (1999).


Diggle, P. J., Moyeed, R. A., Rowlingson, R. and Thomson, M. (2002). Childhood Malaria in the Gambia: A case-study in model-based geostatistics. Journal of the Royal Statistical Society. Series C (Applied Statistics), 51(4): 493-506.

Diggle, P. J., Tawn, J. A. and Moyeed, R. A. (1998). Model-based geostatistics (with Discussion). Applied Statistics, 47, 299–350.

Thomson, M. C., Connor, S. J., D'Alessandro, U., Rowlingson, B., Diggle, P., Creswell, M. and Greenwood, B. (2004). Predicting malaria infection in Gambian children from satellite data and bed net use surveys: the importance of spatial correlation in the interpretation of results. American Journal of Tropical Medicine and Hygiene, 61: 2-8.


gambiaUTM = unwrap(gambiaUTM)

plot(gambiaUTM, main="gambia data")

if(require('mapmisc', quietly=TRUE)) {
  worldMap = unwrap(worldMap)
  gambia = project(worldMap[grep("Gambia", worldMap$NAME), ], crs(gambiaUTM))
  gambiaTiles = openmap(gambia, buffer=100*1000)
  plot(gambiaTiles, add=TRUE)
  plot(gambiaUTM, add=TRUE)

geostatsp documentation built on Sept. 23, 2023, 1:06 a.m.