if (getOption_IsoriX("example_maxtime") > 30) {
## We prepare the data
GNIPDataDEagg <- prepsources(data = GNIPDataDE)
## We fit the models
GermanFit <- isofit(
data = GNIPDataDEagg,
mean_model_fix = list(elev = TRUE, lat_abs = TRUE)
)
## We build the isoscapes
GermanScape <- isoscape(raster = ElevRasterDE, isofit = GermanFit)
## Saving as RDS
filename <- tempfile(fileext = ".rds") # or whatever names you want
saveRDS(GermanScape, file = filename)
## Reading RDS
GermanScape2 <- readRDS(filename)
GermanScape2
}
Run the code above in your browser using DataLab