https://github.com/cran/quantreg
Revision f9636959970240346858a0b1e698208977cc912c authored by Roger Koenker on 14 August 2013, 00:00:00 UTC, committed by Gabor Csardi on 14 August 2013, 00:00:00 UTC
1 parent c4c8cfa
Raw File
Tip revision: f9636959970240346858a0b1e698208977cc912c authored by Roger Koenker on 14 August 2013, 00:00:00 UTC
version 5.02
Tip revision: f963695
lprq.R
lprq <- function(x, y, h, tau = .5, m = 50)
{
	## A toy routine to do locally polynomial quantile regression
        xx <- seq(min(x),max(x),length=m)
        fv <- xx
        dv <- xx
        for(i in 1:length(xx)) {
                z <- x - xx[i]
                wx <- dnorm(z/h)
                r <- rq(y~z, weights=wx, tau=tau, ci=FALSE)
                fv[i] <- r$coef[1.]
                dv[i] <- r$coef[2.]
        }
        list(xx = xx, fv = fv, dv = dv)
}

back to top