# NOT RUN {
y <- rexp(100, 1)
h <- 0.79 * IQR(y) * length(y) ^ (-1/5)
den<-RIG(y,200,h)
plot(den, type = "s", ylab = "Density Function", lty = 1, xlab = "Time")
d1 <- density(y, bw=h) #To add true density along with estimated
lines(d1,type="p",col="red")
legend("topright", c("Real Density", "Density by RIG Kernel"), col=c("red", "black"), lty=c(1,2))
# }
Run the code above in your browser using DataLab