if (FALSE) {
## 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)"))
}
Run the code above in your browser using DataLab