AsyK (version 1.5.5)

plot.RIG: Density Plot by Reciprocal Inverse Gaussian kernel

Description

Plot density by using Reciprocal Inverse Gaussian Kernel.

Usage

# S3 method for RIG
plot(x, ...)

Arguments

x

an object of class "RIG"

…

Not presently used in this implementation

Value

nothing

References

Scaillet, O. 2004. Density estimation using inverse and reciprocal inverse Gaussian kernels. Nonparametric Statistics, 16, 217-226.

See Also

To examine RIG estimated values for density see RIG and for Mean Squared Error mse. Similarly, for plot of Laplace kernel plot.Laplace.

Examples

# NOT RUN {
y <- rexp(200, 1)
h <- 0.79 * IQR(y) * length(y) ^ (-1/5)
xx <- seq(min(y) + 0.05, max(y), length = 200)
den <- RIG(x = xx, y = y, k = 200, h = h)
plot(den, type = "l")


##other details can also be added
y <- rexp(200, 1)
h <- 0.79 * IQR(y) * length(y) ^ (-1/5)
den <- RIG(x = xx, y = y, k = 200, h = h)
plot(den, type = "s", ylab = "Density Function", lty = 1, xlab = "Time")

## To add true density along with estimated
d1 <- density(y, bw = h)
lines(d1, type = "p", col = "red")
legend("topright", c("Real Density", "Density by RIG Kernel"),
col = c("red", "black"), lty = c(1, 2))
# }