R/locf2func.R

In Bioconductor/bioDist: Different distance measures

locf2func <- function(x, ..., maxk = 100) {
    if(!require(locfit))
		stop("can only be used if locfit library is available")
    for(i in 1:10) {
		l <- try(locfit(~x, ..., maxk = i * maxk))
		if(!is(l, "try-error"))
			break;
	}
    den1 <- preplot(l,...)
    xvals <- den1$xev[[1]]
    yvals <- den1$trans(den1$fit)
    f <- function(w) approx(xvals, yvals, w, yleft=0, yright=0)$y
    class(f) <- "dfun"
    f
}