On Sampling Designs for Integral Estimation of a Random Process

Abstract

In this paper the problem of finding exactly optimal sampling designs for estimating the weighted integral of a stochastic process with a product covariance structure (R(s,t)=u(s)v(t), s<t) is discussed. The sampling designs for certain standard processes belonging to the product class are calculated. An asymptotic solution to the design problem also follows as a consequence.     

Estimated estimating equations: semiparametric inference for clustered and longitudinal data

Summary.  We introduce a flexible marginal modelling approach for statistical inference for clustered and longitudinal data under minimal assumptions. This estimated estimating equations approach is semiparametric and the proposed models are fitted by quasi-likelihood regression, where the unknown marginal means are a function of the fixed effects linear predictor with unknown smooth link, and variance–covariance is an unknown smooth function of the marginal means. We propose to estimate the nonparametric link and variance–covariance functions via smoothing methods, whereas the regression parameters are obtained via the estimated estimating equations. These are score equations that contain nonparametric function estimates. The proposed estimated estimating equations approach is motivated by its flexibility and easy implementation. Moreover, if data follow a generalized linear mixed model, with either a specified or an unspecified distribution of random effects and link function, the model proposed emerges as the corresponding marginal (population-average) version and can be used to obtain inference for the fixed effects in the underlying generalized linear mixed model, without the need to specify any other components of this generalized linear mixed model. Among marginal models, the estimated estimating equations approach provides a flexible alternative to modelling with generalized estimating equations. Applications of estimated estimating equations include diagnostics and link selection. The asymptotic distribution of the proposed estimators for the model parameters is derived, enabling statistical inference. Practical illustrations include Poisson modelling of repeated epileptic seizure counts and simulations for clustered binomial responses.     

Controlling the size of multivariate outlier tests with the MCD estimator of scatter

Multivariate outlier detection requires computation of robust distances to be compared with appropriate cut-off points. In this paper we propose a new calibration method for obtaining reliable cut-off points of distances derived from the MCD estimator of scatter. These cut-off points are based on a more accurate estimate of the extreme tail of the distribution of robust distances. We show that our procedure gives reliable tests of outlyingness in almost all situations of practical interest, provided that the sample size is not much smaller than 50. Therefore, it is a considerable improvement over all the available MCD procedures, which are unable to provide good control over the size of multiple outlier tests for the data structures considered in this paper.     

Estimation of a two-equation panel model with mixed continuous and ordered categorical outcomes and missing data

Summary.  The paper proposes an estimation approach for panel models with mixed continuous and ordered categorical outcomes based on generalized estimating equations for the mean and pseudoscore equations for the covariance parameters. A numerical study suggests that efficiency can be gained in the mean parameter estimators by using individual covariance matrices in the estimating equations for the mean parameters. The approach is applied to estimate the returns to occupational qualification in terms of income and perceived job security in a 9-year period based on the German Socio-Economic Panel. To compensate for missing data, a combined multiple imputation–weighting approach is adopted.     

A NOTE ON SAMPLING DESIGNS FOR RANDOM PROCESSES WITH NO QUADRATIC MEAN DERIVATIVE

Several authors have previously discussed the problem of obtaining asymptotically optimal design sequences for estimating the path of a stochastic process using intricate analytical techniques. In this note, an alternative treatment is provided for obtaining asymptotically optimal sampling designs for estimating the path of a second order stochastic process with known covariance function. A simple estimator is proposed which is asymptotically equivalent to the full‐fledged best linear unbiased estimator and the entire asymptotics are carried out through studying this estimator. The current approach lends an intuitive statistical perspective to the entire estimation problem.     

Modelling association between two irregularly observed spatiotemporal processes by using maximum covariance analysis

Summary.  Climatic phenomena such as the El-Niño–southern oscillation and the north Atlantic oscillation are results of complex interactions between atmospheric and oceanic processes. Understanding the interactions has enabled scientists to give early warning of the forthcoming phenomena, thereby reducing damage caused by them. Statistical methods have played an important role in revealing effects of these phenomena on different regions of the world. One such method is maximum covariance analysis (MCA). Two apparent weaknesses are associated with MCA. Firstly, it tends to produce estimates with a low signal-to-noise ratio, especially when the sample size is small. Secondly, there has been no objective way of incorporating incomplete records, which are frequently encountered in climatology and oceanographic data-bases. We introduce an MCA which incorporates a smoothing procedure on the estimates. The introduction of the smoothing procedure is shown to improve the signal-to-noise ratio on the estimates. The estimation of smoothing parameters is carried out by using a penalized likelihood approach, which makes the inclusion of incomplete records quite straightforward. The methodology is applied to investigate the association between Irish winter precipitation and sea surface temperature anomalies around the world. The results show relationships between Irish precipitation anomalies and the El-Niño–southern oscillation and the north Atlantic oscillation phenomena.     

Parameterizations and Fitting of Bi-directed Graph Models to Categorical Data

Abstract.  We discuss two parameterizations of models for marginal independencies for discrete distributions which are representable by bi-directed graph models, under the global Markov property. Such models are useful data analytic tools especially if used in combination with other graphical models. The first parameterization, in the saturated case, is also known as thenation multivariate logistic transformation, the second is a variant that allows, in some (but not all) cases, variation-independent parameters. An algorithm for maximum likelihood fitting is proposed, based on an extension of the Aitchison and Silvey method.     

A note on sequential ML estimates and their asymptotic covariances

A marginal and sequential maximum likelihood estimation method is described which can be used instead of full information maximum likelihood estimation if the latter method is unfeasible. It is shown that the sequential procedure yields strongly consistent and asymptotically normal estimates under relatively general regularity conditions. It is shown that the covariance matrix of the sequential ML estimator does not coincide with the inverse of the Fisher information matrix. Hence, the corrected covariance matrix is derived. The application of the sequential procedure to the multivariate probit model with dichotomous, ordered categorical, single-sided censored and double-sided censored endogenous variables is included. This research was partially supported by a dissertation grant of theStudienstiftung des Deutschen Volkes. Comments and suggestions on earlier drafts by Gerhard Arminger, Giorgio Calzolari, Bernd Kortzen and an anonymous referee are gratefully acknowledged.     

A simple procedure for testing linear hypotheses about the parameters of a nonlinear model using weighted least squares

Suppose the same nonlinear function involving k parameters is fit to each of t populations. Suppose further it is of interest to compare a specific parameter of the models across the populations. Such comparisons can be expressed as linear hypotheses about the parameters of the nonlinear models. A weighted linear least squares (WLLS) procedure is proposed to test these linear hypotheses. The advantages and disadvantages of the WLLS procedure are discussed. This procedure is also compared to a nonlinear least squares procedure for testing these hypotheses in nonlinear models.