## 1- Random field generation
n.site <- 30
locations <- matrix(runif(2*n.site, 0, 10), ncol = 2)
colnames(locations) <- c("lon", "lat")
data <- rmaxstab(60, locations, cov.mod = "whitmat", sill = 1, range =
3, smooth = 1)
## 2- Fit a max-stable processes
schlather <- fitmaxstab(data, locations, "whitmat",
fit.marge = FALSE)
## 3- Plot the extremal coefficient
extcoeff(schlather)Run the code above in your browser using DataLab