library(nlme)
data(ralu.model)
x = c("DEPTH_F", "HLI_F", "CTI_F", "cti", "ffp")
gm_h1 <- gravity(y = "DPS", x = x, d = "DISTANCE", group = "FROM_SITE",
data = ralu.model, ln = FALSE, method="ML")
gravity.es(gm_h1)
Run the code above in your browser using DataLab