# NOT RUN {
raster_2005 <- raster::raster(system.file("external/RASTER_2005.RST", package="intensity.analysis"))
raster_2010 <- raster::raster(system.file("external/RASTER_2010.RST", package="intensity.analysis"))
raster_2012 <- raster::raster(system.file("external/RASTER_2012.RST", package="intensity.analysis"))
raster.layers <- list(raster_2005, raster_2010, raster_2012)
time.points <- c("2005","2010","2012")
categories <- c("Water","Trees","Impervious")
crosstabulation <- multicrosstab(raster.layers, time.points, categories)
IIA.output <- IIA(crosstabulation, time.points)
filename <- file.path(normalizePath(tempdir(), winslash = "/"), "IIA.csv")
IIA2csv(IIA.output,time.points, filename)
# }
Run the code above in your browser using DataLab