R/obs.R

In conservation-decisions/smsPOMDP: POMDP problem solver: When to stop managing or surveying cryptic threatened species ?

#' @export
obs <- function(p0, pm, d0, dm, ds, V, Cm, Cs){
  #tests the inputs
  stopifnot(p0>=0,p0<=1) #checks if p0 is a probability
  stopifnot(pm>=0,pm<=1) #checks if pm is a probability
  stopifnot(d0>=0,d0<=1) #checks if d0 is a probability
  stopifnot(dm>=0,dm<=1) #checks if dm is a probability
  stopifnot(ds>=0,ds<=1) #checks if ds is a probability
  stopifnot(V>=0, Cm >= 0, Cs >= 0) #checks if values and costs are positif
  
  manage <- matrix(c(dm,0,1-dm,1), nrow = 2)
  survey <- matrix(c(ds,0,1-ds,1), nrow = 2)
  nothing <- matrix(c(d0,0,1-d0,1), nrow = 2)
  output <- array(c(manage, survey, nothing), dim = c(2,2,3))
  return(output)
}