In this paper, we examine the validity of non-parametric spatial bootstrap as a procedure to quantify errors in estimates of N-point correlation functions. We do this by means of a small simulation study with simple point process models and estimating the two-point correlation functions and their errors. The coverage of confidence intervals obtained using bootstrap is compared with those obtained from assuming Poisson errors. The bootstrap procedure considered here is adapted for use with spatial (i.e. dependent) data. In particular, we describe a marked point bootstrap where, instead of resampling points or blocks of points, we resample marks assigned to the data points. These marks are numerical values that are based on the statistic of interest. We describe how the marks are defined for the two- and three-point correlation functions. By resampling marks, the bootstrap samples retain more of the dependence structure present in the data. Furthermore, this method of bootstrap can be performed much quicker than some other bootstrap methods for spatial data, making it a more practical method with large datasets. We find that with clustered point datasets, confidence intervals obtained using the marked point bootstrap has empirical coverage closer to the nominal level than the confidence intervals obtained using Poisson errors. The bootstrap errors were also found to be closer to the true errors for the clustered point datasets.

Thanks to the recent availability of large surveys, there has been renewed interest in third-order correlation statistics. Measures of third-order clustering are sensitive to the structure of filaments and voids in the universe and are useful for studying large-scale structure. Thus, statistics of these third-order measures can be used to test and constrain parameters in cosmological models. Third-order measures such as the three-point correlation function are now commonly estimated for galaxy surveys. Studies of third-order clustering of absorption systems will complement these analyses. We define a statistic, which we denote K, that measures third-order clustering of a data set of point observations and focus on estimating this statistic for an absorber catalog. The statistic K can be considered a third-order version of the second-order Ripley K-function and allows one to study the abundance of various configurations of point triplets. In particular, configurations consisting of point triplets that lie close to a straight line can be examined. Studying third-order clustering of absorbers requires consideration of the absorbers as a three-dimensional process, observed on QSO lines of sight that extend radially in three-dimensional space from Earth. Since most of this three-dimensional space is not probed by the lines of sight, edge corrections become important. We use an analytical form of edge correction weights and construct an estimator of the statistic K for use with an absorber catalog. We show that with these weights, ratio-unbiased estimates of K can be obtained. Results from a simulation study also verify unbiasedness and provide information on the decrease of standard errors with increasing number of lines of sight.

We examine the three-dimensional clustering of C IV absorption-line systems, using an extensive catalog of QSO heavy-element absorbers drawn from the literature. We measure clustering by a volume-weighted integral of the correlation function called the reduced second-moment measure, and include information from both along and across QSO lines of sight, thus enabling a full determination of the three-dimensional clustering of absorbers, as well as a comparison of line- and cross-line-of-sight clustering properties. Here we present the three-dimensional reduced second-moment estimator for a three-dimensional point process probed by one-dimensional lines of sight, and apply our algorithm to a sample of 345 C IV absorbers with median redshift of 2.2, drawn from the spectra of 276 QSOs. We confirm the existence of significant clustering on comoving scales up to 100 Mpc, and find that the additional cross-line-of-sight information strengthens the evidence for clustering on scales from 100 Mpc to 150 Mpc. There is no evidence of absorber clustering along or across lines of sight for scales from 150 Mpc to 300 Mpc. We show that with a 300-times larger catalog, such as that to be compiled by the Sloan Digital Sky Survey (100,000 QSOs), use of the full three-dimensional estimator and cross-line-of-sight information will substantially increase clustering sensitivity. We find that standard errors are reduced by an additional factor of 2 to 20 on scales of 30 to 200 Mpc, effectively increasing the sample size by an extra factor of 4 to 400 at large distances.

Motivated by the study of an important data set for understanding the large-scale structure of the universe, this work considers the estimation of the reduced second moment function, or K-function, of a stationary point process observed over a large number of segments of possibly varying lengths. Theory and simulation are used to compare the behavior of isotropic and rigid motion correction estimators and some modifications of these estimators. These results generally support the use of modified versions of the rigid motion correction. When applied to a catalog of astronomical objects known as absorbers, the proposed methods confirm results from earlier analyses of the absorber catalog showing clear evidence of clustering up to 50 Mpc and marginal evidence for clustering of matter on spatial scales beyond 100 Mpc, which is beyond the distance at which clustering of matter is now generally accepted to exist.