R/unionP.R

In SPATIAL-Lab/isorig: Rescale isoscape using known-origin tissue isotope values, and assign origin of unknown samples with isotopic data

unionP <- function(pdR){
  if(class(pdR) != "RasterStack"){
    stop("input probability density map (pdR) should be RasterLayer")
  }
  result <- (1-pdR[[1]])
  for(i in 2:nlayers(pdR)){
    result <- overlay(result, pdR[[i]], fun = function(x,y){return(x*(1-y))})
  }
  plot(1-result)
  title("Union Probability")
  return(1-result)
}