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"))
MAD(grid1, grid2, strata, eval = "original")
MAD(grid1, grid2, strata, eval = "multiple")
if (FALSE) {
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)
MAD(veg_obs1, veg_pre1, veg_mask1, eval = "multiple")
}
par(old.par)
Run the code above in your browser using DataLab