Modelling of connected processes

The problem of comparing, contrasting and combining information from different sets of data is an enduring one in many practical applications of statistics. A specific problem of combining information from different sources arose in integrating information from three different sets of data generated by three different sampling campaigns at the input stage as well as at the output stage of a grey-water treatment process. For each stage, a common process trend function needs to be estimated to describe the input and output material process behaviours. Once the common input and output process models are established, it is required to estimate the efficiency of the grey-water treatment method. A synthesized tool for modelling different sets of process data is created by assembling and organizing a number of existing techniques: (i) a mixed model of fixed and random effects, extended to allow for a nonlinear fixed effect, (ii) variogram modelling, a geostatistical technique, (iii) a weighted least squares regression embedded in an iterative maximum-likelihood technique to handle linear/nonlinear fixed and random effects and (iv) a formulation of a transfer-function model for the input and output processes together with a corresponding nonlinear maximum-likelihood method for estimation of a transfer function. The synthesized tool is demonstrated, in a new case study, to contrast and combine information from connected process models and to determine the change in one quality characteristic, namely pH, of the input and output materials of a grey-water filtering process.     

Transect-spacing design of ice cores on the Antarctic continent

Chemical analyses of ice cores, drilled deep into an ice sheet, provide a historical record of the earth's atmosphere that dates back as far as 400,000–500,000 years. Although the atmosphere mixes quite well, it is recognized that spatial variability associated with ice-core locations should be allowed for. In this article, spatial statistical methodology is applied to the design question of finding the best spacing of ice-core locations on a partial transect of Antarctica.     

AN ALGORITHM FOR THE DESIGN OF 2? FACTORIAL EXPERIMENTS ON CONTINUOUS PROCESSES

Saunders & Eccleston (1992) presented an approach to the design of 2-level factorial experiments for continuous processes. It determined sets of contrasts between the observations that could be well estimated, and then selected a design so that those contrasts estimated the parameters of interest. This paper shows that a well-estimated contrast must have a large number of changes of sign or level, and also be ‘paired’ in a particular sense. It develops an algorithm which constructs designs that must have a large number of changes of sign, evenly spread among the contrasts and optimal or near optimal. When such designs exist they are often preferable to those produced by the reverse foldover algorithm of Cheng & Steinberg (1991).     

Childhood mortality spatial distribution in Ethiopia

The risk of a child dying before completing five years of age is highest in Sub-Saharan African countries. But Child mortality rates have shown substantial decline in Ethiopia. For this study, the 2000, 2005 and 2011 Ethiopian Demographic Survey (EDHS) was used. Generalized linear mixed model with spatial covariance structure was adapted. The model allowed for spatial correlation, and leads to the more realistic estimate for under-five mortality risk factors. The analysis showed that the risk of under-five mortality shows decline in years. But, some regions showed increase in years. The study highlights the need to implement better education for family planning and child care to improve the under-five mortality situation in some administrative areas.     

Comparisons of changepoint estimators

Statistical comparisons of several nonparametric changepoint estimators which are currently offered in the literature are undertaken. The criteria for comparison are Pitman nearness, probability concentration and mean squared error. Recommendations for small, moderate and large sample cases are given.     

A spatiotemporal model for Mexico City ozone levels

Summary.  We consider hourly readings of concentrations of ozone over Mexico City and propose a model for spatial as well as temporal interpolation and prediction. The model is based on a time-varying regression of the observed readings on air temperature. Such a regression requires interpolated values of temperature at locations and times where readings are not available. These are obtained from a time-varying spatiotemporal model that is coupled to the model for the ozone readings. Two location-dependent harmonic components are added to account for the main periodicities that ozone presents during a given day and that are not explained through the covariate. The model incorporates spatial covariance structure for the observations and the parameters that define the harmonic components. Using the dynamic linear model framework, we show how to compute smoothed means and predictive values for ozone. We illustrate the methodology on data from September 1997.     

Quotient spaces and statistical models

The purpose of this paper is to draw attention to the widespread occurrence of quotient spaces in statistical work. Quotient spaces are intrinsic to probability distributions, residuals and interaction in linear models, covariance functions and variograms of stochastic processes, etc. The theme is that explicit recognition of the quotient space can offer surprising conceptual simplification. The advantages of working directly with the quotient space are hard to describe in general. As the examples demonstrate, the answer lies partly in directness of approach.     

Estimation of causal continuous‐time autoregressive moving average random fields

We estimate model parameters of Lévy‐driven causal continuous‐time autoregressive moving average random fields by fitting the empirical variogram to the theoretical counterpart using a weighted least squares (WLS) approach. Subsequent to deriving asymptotic results for the variogram estimator, we show strong consistency and asymptotic normality of the parameter estimator. Furthermore, we conduct a simulation study to assess the quality of the WLS estimator for finite samples. For the simulation, we utilize numerical approximation schemes based on truncation and discretization of stochastic integrals and we analyze the associated simulation errors in detail. Finally, we apply our results to real data of the cosmic microwave background.     

A Kernel Variogram Estimator for Clustered Data

Abstract.  The variogram provides an important method for measuring the dependence of attribute values between spatial locations. Suppose that the nature of the sampling process leads to the presence of clustered data; it would be advisable to use a variogram estimator that aims to adjust for clustering of samples. In this setting, the use of a non-parametric weighted estimator, obtained by considering an inverse weight to a given neighbourhood density combined with the kernel method, seems to have a satisfactory behaviour in practice. This paper pursues a theoretical study of the cluster robust estimator, by proving that it is asymptotically unbiased as well as consistent and by providing criteria for selection of the bandwidth parameter and the neighbourhood radius. Numerical studies are also included to illustrate the performance of the considered estimator and the suggested approaches.     

On asymptotic properties of the mark variogram estimator of a marked point process

The mark variogram [Cressie, 1993. Statistics for Spatial Data. Wiley, New York] is a useful tool to analyze data from marked point processes. In this paper, we investigate the asymptotic properties of its estimator. Our main findings are that the sample mark variogram is a consistent estimator for the true mark variogram and is asymptotically normal under some mild conditions. These results hold for both the geostatistical marking case (i.e., the case where the marks and points are independent) and the non-geostatistical marking case (i.e., the case where the marks and points are dependent). As an application we develop a general test for spatial isotropy and study our methodology through a simulation study and an application to a data set on long leaf pine trees.