请求版主、epoh老师指教

同求。。

对于你的问题，直接可以用locpol包中的locploSmoother函数
N <- 100
xeval <- 0:10/10
d <- data.frame(x = runif(N))
bw <- 0.125
fx <- xeval^2 - xeval + 1
## Non random
d$y <- d$x^2 - d$x + 1
lpest2 <- locPolSmootherC(d$x, d$y , xeval, bw, 2, EpaK)
lpest2$beta1  # 即为导数估计

或者用KernSmooth包中的locpoly函数，设置参数drv=1也可以估计导数

这里面的是单独非参数的一元变量的导数估计呀，半参数的貌似不可以做。