Analysis of incomplete durations with application to contraceptive use

Models for analysing incomplete durations obtained from cross-sectional surveys are presented. The aim of the paper is to develop a framework for analysing the incomplete duration of episodes in progress at the time of the survey by formulating generalized linear models and fitting and assessing them by using standard statistical packages. The maximum quasi-likelihood method is used for model fitting. The choice of the distribution and the diagnostic procedures are discussed. Simulated data from two distributions (the Weibull and log-logistic distributions) are used to evaluate the methodology developed and to assess model misspecifications. A data set on the current use of the contraceptive pill from a cross-sectional survey in Egypt is analysed.     

A modified likelihood ratio test for homogeneity in finite mixture models

Testing for homogeneity in finite mixture models has been investigated by many researchers. The asymptotic null distribution of the likelihood ratio test (LRT) is very complex and difficult to use in practice. We propose a modified LRT for homogeneity in finite mixture models with a general parametric kernel distribution family. The modified LRT has a χ-type of null limiting distribution and is asymptotically most powerful under local alternatives. Simulations show that it performs better than competing tests. They also reveal that the limiting distribution with some adjustment can satisfactorily approximate the quantiles of the test statistic, even for moderate sample sizes.     

Higher order moments of order statistics from INID symmetric random variables

In this paper we present analogues of Balakrishnan's (1989) relations that relate the triple and quadruple moments of order statistics from independent and nonidentically distributed (I.NI.D.) random variables from a symmetric distribution to those of the folded distribution. We then apply these results, along with the corresponding recurrence relations for the exponential distribution derived recently by Childs (2003), to study the robustness of the Winsorized variance.     

A Construction of Lancaster Probabilities with Margins in the Multidimensional Meixner Class

The well-known Meixner class (Meixner, 1934) of probabilities on R has been extended recently to R ^d (Pommeret, 1996). This generalized Meixner class corresponds to the simple quadratic natural exponential families characterized by Casalis (1996). Following Lancaster (1975), the present paper offers a characterization of the joint probability of a randomvector ( X, Y ) such that the two variables X and Y on R ^d belong to the multidimensional Meixner class and fulfil a bi-orthogonality condition involving orthogonal polynomials. The joint probabilities, called Lancaster probabilities, are characterized by two sequences of orthogonal polynomials with respect to the margins and a sequence of expectations of products. Some multivariate probabilities are studied, namely the Poisson-Gaussian and the gamma-Gaussian.     

Bootstrap tests: how many bootstraps?

In practice, bootstrap tests must use a finite number of bootstrap samples. This means that the outcome of the test will depend on the sequence of random numbers used to generate the bootstrap samples, and it necessarily results in some loss of power. We examine the extent of this power loss and propose a simple pretest procedure for choosing the number of bootstrap samples so as to minimize experimental randomness. Simulation experiments suggest that this procedure will work very well in practice.     

Linearity of regression for non-adjacent record values

Let X₁,X₂,… be a sequence of iid random variables having a continuous distribution; by R₁,R₂,… denote the corresponding record values. All the distributions allowing linearity of regressions either E(R_m+k|R_m) or E(R_m|R_m+k) are identified.     

Recent and classical goodness-of-fit tests for the Poisson distribution

We give a critical synopsis of classical and recent tests for Poissonity, our emphasis being on procedures which are consistent against general alternatives. Two classes of weighted Cramér–von Mises type test statistics, based on the empirical probability generating function process, are studied in more detail. Both of them generalize already known test statistics by introducing a weighting parameter, thus providing more flexibility with regard to power against specific alternatives. In both cases, we prove convergence in distribution of the statistics under the null hypothesis in the setting of a triangular array of rowwise independent and identically distributed random variables as well as consistency of the corresponding test against general alternatives. Therefore, a sound theoretical basis is provided for the parametric bootstrap procedure, which is applied to obtain critical values in a large-scale simulation study. Each of the tests considered in this study, when implemented via the parametric bootstrap method, maintains a nominal level of significance very closely, even for small sample sizes. The procedures are applied to four well-known data sets.     

Designing and integrating composite networks for monitoring multivariate gaussian pollution fields

Networks of ambient monitoring stations are used to monitor environmental pollution fields such as those for acid rain and air pollution. Such stations provide regular measurements of pollutant concentrations. The networks are established for a variety of purposes at various times so often several stations measuring different subsets of pollutant concentrations can be found in compact geographical regions. The problem of statistically combining these disparate information sources into a single 'network' then arises. Capitalizing on the efficiencies so achieved can then lead to the secondary problem of extending this network. The subject of this paper is a set of 31 air pollution monitoring stations in southern Ontario. Each of these regularly measures a particular subset of ionic sulphate, sulphite, nitrite and ozone. However, this subset varies from station to station. For example only two stations measure all four. Some measure just one. We describe a Bayesian framework for integrating the measurements of these stations to yield a spatial predictive distribution for unmonitored sites and unmeasured concentrations at existing stations. Furthermore we show how this network can be extended by using an entropy maximization criterion. The methods assume that the multivariate response field being measured has a joint Gaussian distribution conditional on its mean and covariance function. A conjugate prior is used for these parameters, some of its hyperparameters being fitted empirically.     

Asymptotic Bias in Parameter Estimation of AR-Processes Using Recursive Least Squares with Exponential Forgetting

The recursive least squares technique is often extended with exponential forgetting as a tool for parameter estimation in time-varying systems. The distribution of the resulting parameter estimates is, however, unknown when the forgetting factor is less than one. In this paper an approximative expression for bias of the recursively obtained parameter estimates in a time-invariant AR( n_a ) process with arbitrary noise is given, showing that the bias is non-zero and giving bounds on the approximation errors. Simulations confirm the approximation expressions.     

Testing for No Effect by Cosine Series Methods

The properties of three lack-of-fit tests that are related to non-parametric cosine regression analysis are examined in the context of testing for a constant mean function. Analytic power comparisons of these tests vs a most powerful test are made using intermediate asymptotic relative efficiency. In particular, a data-driven test is produced which is asymptotically as efficient as the most powerful test over a class of alternatives. A small scale simulation experiment is conducted to ascertain the extent that the large sample comparisons are applicable to finite samples.