### small example with earthquake data
data(quakes)
set.seed(2)
quakes <- quakes[ sample(1:nrow(quakes), 50), ]
quakes.geogam <- geoGAM(response = "mag",
covariates = c("depth", "stations"),
data = quakes,
seed = 2,
max.stop = 5,
cores = 1)
summary(quakes.geogam)
summary(quakes.geogam, what = "path")
plot(quakes.geogam)
plot(quakes.geogam, what = "path")
Run the code above in your browser using DataLab