swh:1:snp:16c54c84bc54885e783d4424d714e5cc82f479a1
Tip revision: d4080a3a276c7de207e13dbcd1535f15d2b7f7d1 authored by Roger Koenker on 27 October 2020, 18:10:02 UTC
version 5.75
version 5.75
Tip revision: d4080a3
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)
}