old.par <- par(no.readonly = TRUE)
grid1 <- rast(system.file("external/GRID1_INT.rst", package = "diffeR"))
grid2 <- rast(system.file("external/GRID2_INT.rst", package = "diffeR"))
strata <- rast(system.file("external/strata_int.rst", package = "diffeR"))
MADscatterplot(grid1, grid2, strata)
veg_obs1 <- rast(system.file("external/veg_obs1.rst", package = "diffeR"))
veg_pre1 <- rast(system.file("external/veg_pre1.rst", package = "diffeR"))
veg_mask1 <- rast(system.file("external/veg_mask1.rst", package = "diffeR"))
MADscatterplot(veg_obs1, veg_pre1, veg_mask1)
par(old.par)
Run the code above in your browser using DataLab