Search

Small-angle X-ray scattering (SAXS) is a technique for obtaining low-resolution structural information about biological macromolecules, by exposing a dilute solution to a high-intensity X-ray beam and capturing the resulting scattering pattern on a two-dimensional detector. The two-dimensional pattern is reduced to a one-dimensional curve through radial averaging; that is, by averaging across annuli on the detector plane. Subsequent analysis of structure relies on these one-dimensional data. This paper reviews the technique of SAXS and investigates autocorrelation structure in the detector plane and in the radial averages. Across a range of experimental conditions and molecular types, spatial autocorrelation in the detector plane is present and is well-described by a stationary kernel convolution model. The corresponding autocorrelation structure for the radial averages is non-stationary. Implications of the autocorrelation structure for inference about macromolecular structure are discussed.

Many environmental data sets have a continuous domain, in time and/or space, and complex features that may be poorly modelled with a stationary (in space and time) Gaussian process (GP). We adapt stochastic volatility modelling to this context, resulting in a stochastic heteroscedastic process (SHP), which is unconditionally stationary and non-Gaussian. Conditional on a latent GP, the SHP is a heteroscedastic GP with non-stationary (in space and time) covariance structure. The realizations from SHP are versatile and can represent spatial inhomogeneities. The unconditional correlation functions of SHP form a rich isotropic class that can allow for a smoothed nugget effect. We apply an importance sampling strategy to implement pseudo maximum likelihood parameter estimation for the SHP. To predict the process at unobserved locations, we develop a plug-in best predictor. We extend the single-realization SHP model to handle replicates across time of SHP realizations in space. Empirical results with simulated data show that SHP is nearly as efficient as a stationary GP in out-of-sample prediction when the true process is a stationary GP, and outperforms a stationary GP substantially when the true process is SHP. The SHP methodology is applied to enhanced vegetation index data and US NO3 deposition data for illustration.

We consider bootstrap construction and calibration of prediction intervals for nonGaussian autoregressions. In particular, we address the question of prediction conditioned on the last p observations of the process, for which we offer an exact simulation technique and an approximate bootstrap approach. In simulations for a variety of first-order autoregressions, we compare various nonparametric prediction intervals and find that calibration gives reasonably narrow prediction intervals with the lowest coverage probability mean squared error among the methods used.