R/simPrey.R

In HFpostdoc/sarModel: Sarracenia micro-ecosystem model

#' @title simPrey
#' 
#' Description simPrey
#' 
#' @param days days
#' @param prey prey
#' @param prey prey
#' @param perturb perturb
#' @param perturb perturb
#' @export simPrey
#' @return simPrey
simPrey <- function(days = 30, prey.mass = 10,prey.rate = 3,perturb.mass = 10, perturb.rate = 3){

    fw <- rep(c(rep(0,prey.rate - 1),prey.mass),ceiling(days / prey.rate))[1:days]
    fw.perturb <- rep(c(rep(0,perturb.rate - 1),perturb.mass),ceiling(days / perturb.rate))[1:days]
    fw + fw.perturb

}