assignR (version 2.2.3)

unionP: Union probability of origin

Description

Probabilities that at least one individual came from each location in the assignment area (union of probabilities)

Usage

unionP(pdR)

Value

RasterLayer.

Arguments

pdR

RasterStack or RasterBrick of probability density maps, e.g., as produced by pdRaster.

Examples

# load North America boundary and isoscape
data("naMap")
data("d2h_lrNA")

# load hydrogen isotope data for human hair in North America
d = subOrigData(group = "Modern human", mask = naMap, niter = 100)

# rescale from environmental isoscape to tissue isoscape
r = calRaster(known = d, isoscape = d2h_lrNA, mask = naMap)

# four unknown-origin examples
id = c("A", "B", "C", "D")
d2H = c(-110, -90, -105, -102)
un = data.frame(id,d2H)

# posterior probabilities
pp = pdRaster(r, unknown = un, mask = naMap)

# probability that one or more individuals are from a given location
unionP(pp)