Improving the Cavalieri estimator under non-equidistant sampling and dropouts

Mads Stehr, Markus Kiderlen

Abstract

Motivated by the stereological problem of volume estimation from parallel section profiles, the so-called Newton-Cotes integral estimators based on random sampling nodes are analyzed. These estimators generalize the classical Cavalieri estimator and its variant for non-equidistant sampling nodes, the generalized Cavalieri estimator, and have typically a substantially smaller variance than the latter. The present paper focuses on the following points in relation to Newton-Cotes estimators: the treatment of dropouts, the construction of variance estimators, and, finally, their application in volume estimation of convex bodies.

Dropouts are eliminated points in the initial stationary point process of sampling nodes, modeled by independent thinning. Among other things, exact representations of the variance are given in terms of the thinning probability and increments of the initial points under two practically relevant sampling models. The paper presents a general estimation procedure for the variance of Newton-Cotes estimators based on the sampling nodes in a bounded interval. Finally, the findings are illustrated in an application of volume estimation for three-dimensional convex bodies with sufficiently smooth boundaries.


Keywords
Cavalieri estimator; dropouts; Newton-Cotes quadrature; numerical integration with random nodes; stationary point process; variance estimation; weakly (m,p)-piecewise smooth function

Full Text:

PDF


DOI: 10.5566/ias.2422

Creative Commons License
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.

Image Analysis & Stereology
EISSN 1854-5165 (Electronic version)
ISSN 1580-3139 (Printed version)