Vitale | R Documentation |
"Vitale"
This class deals with Vitale (1975) Bernstein Polynomial approximation as described in Leblanc (2009). The polynomial estimator is computed using the provided data samples. Using this polynomial estimator, the methods implemented in the class can be used to compute densities, values of the distribution function, quantiles, sample the distribution and obtain graphical representations.
Objects can be created by using the generator function vitale
.
dataPointsCache
:a numeric vector containing points within the [lower.limit,upper.limit]
interval
densityCache
:a numeric vector containing the density for each point in dataPointsCache
distributionCache
:a numeric vector used to cache the values of the distribution function. This slot is included to improve the performance of the methods when multiple calculations of the distribution function are used
dataPoints
:a numeric vector containing data samples within the [lower.limit,upper.limit]
interval. These data samples are used to obtain the kernel estimator
m
:the order of the polynomial approximation
lower.limit
:a numeric value for the lower limit of the bounded interval for the data
upper.limit
:a numeric value for the upper limit of the bounded interval for the data
See "density"
for details
See "distribution"
for details
See "quantile"
for details
See "rsample"
for details
See "plot"
for details
See "getdataPointsCache"
for details
See "getdensityCache"
for details
See "getdistributionCache"
for details
See "getdataPoints"
for details
See "getm"
for details
Guzman Santafe, Borja Calvo and Aritz Perez
Vitale, R. A. (1975). A Bernstein polynomial approach to density function estimation. Statistical Inference and Related Topics, 2, 87-99.
Leblanc, A. (2010). A bias-reduced approach to density estimation using Bernstein polynomials. Journal of Nonparametric Statistics, 22(4), 459-475.
# create the model model <- vitale(dataPoints = tuna.r, m = 25) # examples of usual functions density(model,0.5) distribution(model,0.5,discreteApproximation=FALSE) # graphical representation hist(tuna.r,freq=FALSE,main="Tuna Data") lines(model, col="red",lwd=2) # graphical representation using ggplot2 graph <- gplot(model, show=TRUE, includePoints=TRUE)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.