Skip to main content
Log in

Posterior Contraction Rates of Density Derivative Estimation

  • Published:
Sankhya A Aims and scope Submit manuscript

Abstract

In this paper, we study the problem of Bayesian estimation of derivatives of a density function on the unit interval. We use a finite random series prior based on B-splines and study the asymptotic properties of the posterior distribution under the setting of fixed smoothness of the true function. We obtain the posterior contraction rate under both the L 2- and \(L_{\infty }\)-distances. The rate under L 2-distance agrees with the minimax optimal rate. This result is then extended to the estimation of a multivariate density function on the unit cube and its mixed partial derivatives using tensor product B-splines.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Institutional subscriptions

Similar content being viewed by others

References

  • Bartlett MS (1963) Statistical estimation of density function. Sankhyā, Ser A 25:245–254.

  • Belitser E, Serra P (2014) Adaptive priors based on splines with random knots. Bayesian Anal 9:859–882.

  • Bhattacharya PK (1967) Estimation of a probability density function and its derivatives. Sankhyā Ser A 29:373–382.

  • Castillo I (2014) On Bayesian supremum norm contraction rates. Ann Statist 42:2058–2091.

  • Chacón JE, Duong T (2013) Data-driven density derivative estimation, with applications to nonparametric clustering and bump hunting. Electron J Stat 7:499–532.

  • Donoho DL (1988) One-sided inference about functionals of a density. Ann Statist 16:1390–1420.

  • Farrell RH (1972) On the best obtainable asymptotic rates of convergence in estimation of a density function at a point. Ann Math Statist 43:170–180.

  • Fukunaga K, Hostetler L (1975) The estimation of the gradient of a density function, with applications in pattern recognition. IEEE Trans Inform Theory 21:32–40.

  • Genovese CR, Perone-Pacifico M, Verdinelli I, Wasserman L (2016) Non-parametric inference for density modes. J Roy Statist Soc, Ser B 78:99–126.

  • Ghosal S, van der Vaart A (2001) Entropies and rates of convergence for maximum likelihood and Bayes estimation for mixtures of normal densities. Ann Statist 29:1233–1263.

  • Ghosal S, van der Vaart A (2007a) Convergence rates of posterior distributions for noniid observations. Ann Statist 35:192–223.

  • Ghosal S, van der Vaart A (2007b) Posterior convergence rates of Dirichlet mixtures at smooth densities. Ann Statist 35:697–723.

  • Ghosal S, Ghosh JK, Ramamoorthi RV (1999) Posterior consistency of Dirichlet mixtures in density estimation. Ann Statist 27:143–158.

  • Ghosal S, Ghosh JK, van der Vaart A (2000) Convergence rates of posterior distributions. Ann Statist 28:500–531.

  • Giné E, Nickl R (2011) Rates of contraction for posterior distributions in L r-metrics, \(1 \leq r \leq \infty \). Ann Statist 39:2883–2911.

  • Hall P, Yatchew A (2007) Nonparametric estimation when data on derivatives are available. Ann Statist 35:300–323.

  • Hosseinioun N, Doosti H, Niroumand HA (2011) Nonparametric estimation of a multivariate probability density for mixing sequences by the method of wavelets. Ital J Pure Appl Math 28:31–40.

  • Hosseinioun N, Doosti H, Nirumand HA (2012) Nonparametric estimation of the derivatives of a density by the method of wavelet for mixing sequences. Statist Papers 53:195–203.

  • de Jonge R, van Zanten H (2012) Adaptive estimation of multivariate functions using conditionally gaussian tensor-product spline priors. Electron J Stat 6:1984–2001.

  • Knapik BT, van der Vaart AW, van Zanten JH (2011) Bayesian inverse problems with Gaussian priors. Ann Statist 39.

  • Parzen E (1962) On the estimation of a probability density and mode. Ann Math Statist 33:1065–1076.

  • Prakasa Rao B (1996) Nonparametric estimation of the derivatives of a density by the method of wavelets. Bull Inform Cybernet 28:91–100.

  • Qiao W, Polonik W (2016) Theoretical analysis of nonparametric filament estimation. Ann Statist 44:1269–1297.

  • Ray K (2013) Bayesian inverse problems with non-conjugate priors. Electron J Stat 7:2516–2549.

  • Rivoirard V, Rousseau J (2012) Posterior concentration rates for infinite dimensional exponential families. Bayesian Anal 7:311–334.

  • Rosenblatt M (1956) Remarks on some nonparametric estimates of a density function. Ann Math Statist 27:832–837.

  • Sasaki H, Noh YK, Niu G, Sugiyama M (2016) Direct density derivative estimation. Neural Computation 28:1101–1140.

  • Schumaker L (2007) Spline Functions: Basic Theory. Cambridge University Press.

  • Schuster EF (1969) Estimation of a probability density function and its derivatives. Ann Math Statist 40:1187–1195.

  • Shen W, Ghosal S (2015) Adaptive Bayesian procedures using random series priors. Scand J Statist 42:1194–1213.

  • Shen W, Ghosal S (2016) Adaptive Bayesian density regression for high dimensional data. Bernoulli 22:396–420.

  • Shen W, Tokdar ST, Ghosal S (2013) Adaptive Bayesian multivariate density estimation with Dirichlet mixtures. Biometrika 100:623–640.

  • Silverman BW (1986) Density Estimation for Statistics and Data Analysis. CRC press.

  • Singh RS (1977a) Applications of estimators of a density and its derivatives to certain statistical problems. J Roy Statist Soc, Ser B 39:357–363.

  • Singh RS (1977b) Improvement on some known nonparametric uniformly consistent estimators of derivatives of a density. Ann Statist 5:394–399.

  • Stone CJ (1990) Large-sample inference for log-spline models. Ann Statist 18:717–741.

  • Szabo B, van der Vaart AW, van Zanten JH (2015) Frequentist coverage of adaptive nonparametric Bayesian credible sets. Ann Statist 43.

  • Tokdar ST, Ghosh JK (2007) Posterior consistency of logistic gaussian process priors in density estimation. J Statist Plann Inference 137:34–42.

  • van der Vaart A, van Zanten H (2008) Rates of contraction of posterior distributions based on Gaussian process priors. Ann Statist 36:1435–1463.

  • van der Vaart A, van Zanten JH (2009) Adaptive Bayesian estimation using a Gaussian random field with inverse gamma bandwidth. Ann Statist 37:2655–2675.

  • Yoo WW, Ghosal S (2016) Supremum norm posterior contraction and credible sets for nonparametric multivariate regression. Ann Statist 44:1069–1102.

  • Zhou S, Wolfe DA (2000) On derivative estimation in spline regression. Statist Sinica 10:93–108.

Download references

Acknowledgements

The authors would like to thank a referee for insightful comments. The research of the second author is partially supported by National Science Foundation grant number DMS-1510238.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Weining Shen.

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Shen, W., Ghosal, S. Posterior Contraction Rates of Density Derivative Estimation. Sankhya A 79, 336–354 (2017). https://doi.org/10.1007/s13171-017-0105-7

Download citation

  • Received:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s13171-017-0105-7

Keywords and phrases

AMS (2000) subject classification

Navigation