R/RemoveTrend.R

In detrendeR: A Graphical User Interface (GUI) to Visualize and Analyze Dendrochronological Data

RemoveTrend <- function(rw,
                        method = c("Spline", "Spline%", "ModNegExp", "Mean"),
                        BandwidthPerc = 2 / 3,
                        Bandwidth = 32,
                        P = 1 / 2) {
  switch(method,
    "Spline" = {
      curve <- .SPLINE(rw, bandwidth = Bandwidth, p = P)
      return(curve)
    },
    "Spline%" = {
      curve <- .SPLINE(rw,
        bandwidth = length(rw) * BandwidthPerc,
        p = P
      )
      return(curve)
    },
    "ModNegExp" = {
      curve <- .NegExp(rw)
      return(curve)
    },
    "Mean" = {
      series <- rw[!is.na(rw)]
      curve <- rep(mean(series), length(series))
      rw[!is.na(rw)] <- curve
      return(rw)
    }
  )

  # BandwidthPerc <- as.numeric(BandwidthPerc)
  #
  # if (method == "Spline") {
  #   curve <- SPLINE(rw, bandwidth = Bandwidth, p = P)
  #   return(curve)
  # }
  # if (method == "Spline%") {
  #   curve <- SPLINE(rw, bandwidth = length(rw) * BandwidthPerc, p = P)
  #   return(curve)
  # }
  # if (method == "ModNegExp") {
  #   curve <- NegExp(rw)
  #   return(curve)
  # }
  # if (method == "Mean") {
  #   series <- rw[!is.na(rw)]
  #   curve <- rep(mean(series), length(series))
  #   rw[!is.na(rw)] = curve
  #   return(rw)
  # }
}