R/adj.series.grad.cos.R

In DistanceDevelopment/mrds: Mark-Recapture Distance Sampling

#' Series of the gradient of the cosine adjustment series w.r.t. the scaled 
#' distance.
#' 
#' For internal use only -- not to be called by 'mrds' or 'Distance' users 
#' directly.
#'
#' @param distance perpendicular distance vector/scalar
#' @param scaling scale parameter
#' @param adj.order the adjustment order
#' @param adj.parm vector of parameters (a_j)
#' @param adj.exp boolean, defaults to FALSE
#' 
#' @return scalar or vector containing the gradient of the cosine adjustment 
#' series for every value in \code{distance} argument
#'
#' @author Felix Petersma
adj.series.grad.cos <- function(distance, scaling = 1, adj.order, 
                                adj.parm = NULL, adj.exp = FALSE){
  
  # Check the adjustment parameters
  if(is.null(adj.parm)){
    adj.parm <- as.vector(rep(1, length(adj.order)))
  }
  
  adj.order <- as.vector(adj.order)
  
  cossum <- 0
  
  for(i in seq_along(adj.order)){
    cossum <- cossum + 
      (adj.parm[i] * -(adj.order[i] * pi * sin(adj.order[i] * pi * 
                                                 distance / scaling)))
  }
  
  # if adj.exp return exp(cossum) to keep f(x)>0  
  if(adj.exp){
    return(exp(cossum))
  }else{
    return(cossum)
  }
}