## disease mapping example
data(leukDat)
data(GDRmap) # this is the map of the former GDR
mix.leuk <- mixalg(obs = "oleuk", pop.at.risk = "eleuk",
data = leukDat, family = "poisson")
mix.leuk
plot(GDRmap, col = mix.leuk@classification)
Run the code above in your browser using DataLab