plot.RIG: Density Plot by Reciprocal Inverse Gaussian kernel

In AsyK: Kernel Density Estimation

Description

Plot density by using Reciprocal Inverse Gaussian Kernel.

Usage

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

Arguments

x	an object of class "RIG"
...	Not presently used in this implementation

Value

nothing

Author(s)

Javaria Ahmad Khan, Atif Akbar.

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

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))

AsyK documentation built on Dec. 11, 2021, 10:07 a.m.