makecumsumSurface: Create cumulative sum probability surface

In cjcampbell/IsoModAT: Isotope Origin Clustering and Assignment Tools

Create cumulative sum probability surface

Description

Converts normalized probability surface (e.g. one layer output of isotopeAssignmentModel function) to cumulative sum surfaces, i.e., one where the new value of a given cell is equal to the sum of all old values less than or equal to the old value of the cell.

Usage

makecumsumSurface(indivraster, rescale = FALSE, rename = FALSE)

Arguments

indivraster	Normalized probability surface RasterLayer
rescale	Rescale between 0 and 1? Defaults to FALSE.
rename	Character value to append to raster name (e.g. "_odds"). Defaults to FALSE.

Value

Returns RasterLayer rescaled to Cumulative Sum values.

See Also

cumsumAtSamplingLocation

Examples

# Generate example probability surfaces.
myiso <- rasterFromXYZ(isoscape)
myiso_sd <- rasterFromXYZ(isoscape_sd)
df <- data.frame(
        ID = c(-100, -80, -50),
        isotopeValue = c(-100, -80, -50),
        SD_indv = rep(5, 3)
        )
assignmentModels <- isotopeAssignmentModel(
        ID = df$ID,
        isotopeValue = df$isotopeValue,
        SD_indv = df$SD_indv,
        precip_raster = myiso,
        precip_SD_raster = myiso_sd,
        nClusters = FALSE
        )

# Convert to cumulative sum surface.
cumulative_sum_surface <- stack(
     lapply( unstack( assignmentModels ), makecumsumSurface )
     )
plot(cumulative_sum_surface)

cjcampbell/IsoModAT documentation built on Feb. 21, 2024, 10:27 a.m.