## Not run:
# ## load zambia data and map
# data("ForestHealth")
# data("BeechBnd")
#
# fm <- bayesx(defoliation ~ stand + fertilized +
# humus + moisture + alkali + ph + soil +
# sx(age) + sx(inclination) + sx(canopy) +
# sx(year) + sx(elevation),
# family = "cumlogit", method = "REML", data = ForestHealth)
#
# summary(fm)
# plot(fm, term = c("sx(age)", "sx(inclination)",
# "sx(canopy)", "sx(year)", "sx(elevation)"))
# ## End(Not run)
Run the code above in your browser using DataLab