H-likelihood: problems and solutions

In recent issues of this journal it has been asserted in two papers that the use of h-likelihood is wrong, in the sense of giving unsatisfactory estimates of some parameters for binary data (Kuk and Cheng, 1999; Waddington and Thompson, 2004) or theoretically unsound (Kuk and Cheng, 1999). We wish to refute both these assertions.     

A direct derivation of the REML likelihood function

The original derivation of the widely cited form of the REML likelihood function for mixed linear models is difficult and indirect. This paper derives it directly using familiar operations with matrices and determinants.     

Nonconglomerative coherent conditional probabilities in statistical inference

Noteworthy connections among conglomerability, countable additivity and coherence are discussed in detail, reaching the conclusion that nonconglomerable conditional probabilities must not be doomed and play a significant role in statistical inference. Extended and updated version of a contributed paper presented at the International Conference on “Information Processing and Management of Uncertainty in knowledge-based systems”, IPMU 2004, Perugia, Italy.     

Data skeletons: simultaneous estimation of multiple quantiles for massive streaming datasets with applications to density estimation

We consider the problem of density estimation when the data is in the form of a continuous stream with no fixed length. In this setting, implementations of the usual methods of density estimation such as kernel density estimation are problematic. We propose a method of density estimation for massive datasets that is based upon taking the derivative of a smooth curve that has been fit through a set of quantile estimates. To achieve this, a low-storage, single-pass, sequential method is proposed for simultaneous estimation of multiple quantiles for massive datasets that form the basis of this method of density estimation. For comparison, we also consider a sequential kernel density estimator. The proposed methods are shown through simulation study to perform well and to have several distinct advantages over existing methods.     

A multivariate version of Gini's rank association coefficient

In this paper, we introduce a multivariate generalization of the population version of Gini's rank association coefficient, giving a response to this open question posed in [4]. We also study some properties of this version, present the corresponding results for the sample statistic, and provide several examples.     

Acceleration of the Multiple-Try Metropolis algorithm using antithetic and stratified sampling

The Multiple-Try Metropolis is a recent extension of the Metropolis algorithm in which the next state of the chain is selected among a set of proposals. We propose a modification of the Multiple-Try Metropolis algorithm which allows for the use of correlated proposals, particularly antithetic and stratified proposals. The method is particularly useful for random walk Metropolis in high dimensional spaces and can be used easily when the proposal distribution is Gaussian. We explore the use of quasi Monte Carlo (QMC) methods to generate highly stratified samples. A series of examples is presented to evaluate the potential of the method.     

Cox regression analysis in presence of collinearity: an application to assessment of health risks associated with occupational radiation exposure

This paper considers the analysis of time to event data in the presence of collinearity between covariates. In linear and logistic regression models, the ridge regression estimator has been applied as an alternative to the maximum likelihood estimator in the presence of collinearity. The advantage of the ridge regression estimator over the usual maximum likelihood estimator is that the former often has a smaller total mean square error and is thus more precise. In this paper, we generalized this approach for addressing collinearity to the Cox proportional hazards model. Simulation studies were conducted to evaluate the performance of the ridge regression estimator. Our approach was motivated by an occupational radiation study conducted at Oak Ridge National Laboratory to evaluate health risks associated with occupational radiation exposure in which the exposure tends to be correlated with possible confounders such as years of exposure and attained age. We applied the proposed methods to this study to evaluate the association of radiation exposure with all-cause mortality.     

Some preservation results of NBUC aging property with applications

This paper develops a new characterization of NBUC aging property, and investigates its preservation properties both under monotonic anti-star-shaped transformations and under the non-homogeneous Poisson shock models.     

Estimation of the means of the bivariate normal using moving extreme ranked set sampling with concomitant variable

The estimation of the means of the bivariate normal distribution, based on a sample obtained using a modification of the moving extreme ranked set sampling technique (MERSS) is considered. The modification involves using a concomitant random variable. Nonparametric-type methods as well as the maximum likelihood estimation are considered. The estimators obtained are compared to their counterparts based on simple random sampling (SRS). It appears that the suggested estimators are more efficient. Also, MERSS with concomitant variable is easier to use in practice than the usual ranked set sampling (RSS) with concomitant variable. The issue of robustness of the procedure is addressed. Real trees data set is used for illustration.