If you want to do the density estimation with adaptive kernel method, you can use the "akj" function in the "quantreg" package.

Hope this helps