R/ksmooth2.R

In KarimOualkacha/SmoothMSC: The smooth methylation status calling (SMSC) procedure

ksmooth2 <-
function (x, y, kernel = c("box", "normal"), bandwidth = 0.5, 
                    range.x = range(x), n.points = max(100L, length(x)), x.points,band) 
{
  if (missing(y) || is.null(y)) 
    stop("numeric y must be supplied.\nFor density estimation use density()")
  kernel <- match.arg(kernel)
  krn <- switch(kernel, box = 1L, normal = 2L)
  x.points <- if (missing(x.points)) 
    seq.int(range.x[1L], range.x[2L], length.out = n.points)
  else {
    n.points <- length(x.points)
    sort(x.points)
  }
  ord <- order(x)
  .Call("ksmooth2", x[ord], y[ord], x.points, krn, band)
}