# NOT RUN {
data(GoettingenHealth1)
# general practitioners, psychotherapists and pharmacies
area_goe <- 1753000000
# area of Landkreis Goettingen (sqm)
area_nom <- 1267000000
# area of Landkreis Northeim (sqm)
area_gn <- area_goe+area_nom
sqrt(area_gn/pi)
# this takes some seconds
ripley(GoettingenHealth1[GoettingenHealth1$type == "phys_gen",],
"location", "lat", "lon", area = area_gn, t.max = 30000, t.sep = 300)
ripley(GoettingenHealth1[GoettingenHealth1$type == "pharm",],
"location", "lat", "lon", area = area_gn, t.max = 30000, t.sep = 300)
ripley(GoettingenHealth1[GoettingenHealth1$type == "psych",],
"location", "lat", "lon", area = area_gn, t.max = 30000, t.sep = 300)
# }
Run the code above in your browser using DataLab