R/subsetByLOCID.R

In apear9/SSNdesign: An R Package for Optimal and Adaptive Sampling Designs on Stream Networks

subsetByLOCID <- function(ssn, path, obs, preds){
  
  ssn <- subsetSSN(ssn, path, subset = ssn@obspoints@SSNPoints[[1]]@point.data$locID %in% obs)
  if(!is.null(preds)){
    ssn <- subsetPreds(ssn, locID %in% preds)
    createDistMat(ssn, "preds", T, T)
  } else {
    createDistMat(ssn, o.write = TRUE)
  }
  return(ssn)
  
}