# Pharmacies in a 500 m buffer:
data(health1)
# Health service locations (physicians and pharmacies)
data(health2)
# 420 regions
health2_city <- health2[health2$city == 1,]
# only regions in the city of Goettingen (city=1)
pharmacy <- health1[health1$type == "pharm",]
# Only pharmacies
pharm_city <- dist.buf(health2_city, "region", "lat", "lon", pharmacy,
"location", "lat", "lon", bufdist = 500)
# Pharmacies in a 500 m buffer from the statistical districts in Goettingen
Run the code above in your browser using DataLab