Merging information for semiparametric density estimation

Summary.  The density ratio model specifies that the likelihood ratio of m −1 probability density functions with respect to the m th is of known parametric form without reference to any parametric model. We study the semiparametric inference problem that is related to the density ratio model by appealing to the methodology of empirical likelihood. The combined data from all the samples leads to more efficient kernel density estimators for the unknown distributions. We adopt variants of well-established techniques to choose the smoothing parameter for the density estimators proposed.     

Clustering microarray data using model-based double K-means

The microarray technology allows the measurement of expression levels of thousands of genes simultaneously. The dimension and complexity of gene expression data obtained by microarrays create challenging data analysis and management problems ranging from the analysis of images produced by microarray experiments to biological interpretation of results. Therefore, statistical and computational approaches are beginning to assume a substantial position within the molecular biology area. We consider the problem of simultaneously clustering genes and tissue samples (in general conditions) of a microarray data set. This can be useful for revealing groups of genes involved in the same molecular process as well as groups of conditions where this process takes place. The need of finding a subset of genes and tissue samples defining a homogeneous block had led to the application of double clustering techniques on gene expression data. Here, we focus on an extension of standard K-means to simultaneously cluster observations and features of a data matrix, namely double K-means introduced by Vichi (2000). We introduce this model in a probabilistic framework and discuss the advantages of using this approach. We also develop a coordinate ascent algorithm and test its performance via simulation studies and real data set. Finally, we validate the results obtained on the real data set by building resampling confidence intervals for block centroids.     

Small area estimation for longitudinal surveys

Over the last few years many studies have been carried out in Italy to identify reliable small area labour force indicators. Considering the rotated sample design of the Italian Labour Force Survey, the aim of this work is to derive a small area estimator which borrows strength from individual temporal correlation, as well as from related areas. Two small area estimators are derived as extensions of an estimation strategies proposed by Fuller (1990) for partial overlap samples. A simulation study is carried out to evaluate the gain in efficiency provided by our solutions. Results obtained for different levels of autocorrelation between repeated measurements on the same outcome and different population settings show that these estimators are always more reliable than the traditional composite one, and in some circumstances they are extremely advantageous.The present paper is financially supported by Murst-Cofin (2001) Lutilizzo di informazioni di tipo amministrativo nella stima per piccole aree e per sottoinsiemi della popolazione (National Coordinator Prof. Carlo Filippucci).     

Optimal use of two auxiliary variables in double sampling

Double sampling scheme is used when cheap auxiliary variables may be measured to improve the estimation of a finite population parameter. Several estimators for population mean, ratio of means and variance are available, when two dependent samples are drawn. However, there are few proposals for the case of independent samples. In this paper both cases of dependent and independent samples are dealt with. A general approach for estimating a finite population parameter is given, showing that all the proposed estimators are particular cases of the same general class. The minimum variance bound for any estimator in this class is provided (at the first order of approximation). Furthermore, an optimal estimator which reaches this minimum is found.     

Non-ignorable missing covariate data in survival analysis: a case-study of an International Breast Cancer Study Group trial

Summary.  Non-ignorable missing data, a serious problem in both clinical trials and observational studies, can lead to biased inferences. Quality-of-life measures have become increasingly popular in clinical trials. However, these measures are often incompletely observed, and investigators may suspect that missing quality-of-life data are likely to be non-ignorable. Although several recent references have addressed missing covariates in survival analysis, they all required the assumption that missingness is at random or that all covariates are discrete. We present a method for estimating the parameters in the Cox proportional hazards model when missing covariates may be non-ignorable and continuous or discrete. Our method is useful in reducing the bias and improving efficiency in the presence of missing data. The methodology clearly specifies assumptions about the missing data mechanism and, through sensitivity analysis, helps investigators to understand the potential effect of missing data on study results.     

Perfect samplers for mixtures of distributions

Summary. We consider the construction of perfect samplers for posterior distributions associated with mixtures of exponential families and conjugate priors, starting with a perfect slice sampler in the spirit of Mira and co-workers. The methods rely on a marginalization akin to Rao–Blackwellization and illustrate the duality principle of Diebolt and Robert. A first approximation embeds the finite support distribution on the latent variables within a continuous support distribution that is easier to simulate by slice sampling, but we later demonstrate that the approximation can be very poor. We conclude by showing that an alternative perfect sampler based on a single backward chain can be constructed. This alternative can handle much larger sample sizes than the slice sampler first proposed.     

Inference for the generalized Rayleigh distribution based on progressively censored data

In this paper, and based on a progressive type-II censored sample from the generalized Rayleigh (GR) distribution, we consider the problem of estimating the model parameters and predicting the unobserved removed data. Maximum likelihood and Bayesian approaches are used to estimate the scale and shape parameters. The Gibbs and Metropolis samplers are used to predict the life lengths of the removed units in multiple stages of the progressively censored sample. Artificial and real data analyses have been performed for illustrative purposes.     

捕获再捕获抽样估计量的模拟研究

目前,捕获再捕获抽样已广泛应用于生物科学、社会科学、生命科学和医学的调查研究中。为此,对捕获再捕获抽样的最常用的Lincoln-Petersen估计量和Chapman估计量进行模拟比较,提出新的有偏估计量,其偏差介于上述两个估计量之间。理论和模拟结果显示：新估计量有时优于其它两个估计量;在实际应用中,当上述两个估计量较难取舍时,新估计量可以作为一个更好的选择。     

An adaptive allocation sampling design for which the conventional stratified estimator is an appropriate estimator

Adaptive allocations in stratified sampling design are suitable for studying Biological and Environmental populations. Biologists tend to use the conventional stratified estimator for an adaptive allocation sampling design in spite of its adaptive nature. In this paper, we have introduced an adaptive allocation sampling design for which the conventional stratified estimator has some desired statistical properties. We conduct a simulation study on a real population. The results show that the conventional stratified estimator for the introduced adaptive allocation sampling design is more efficient than other available estimators.     

Improved attribute acceptance sampling plans based on maxima nomination sampling

This paper demonstrates the use of maxima nomination sampling (MNS) technique in design and evaluation of single AQL, LTPD, and EQL acceptance sampling plans for attributes. We exploit the effect of sample size and acceptance number on the performance of our proposed MNS plans using operating characteristic (OC) curve. Among other results, we show that MNS acceptance sampling plans with smaller sample size and bigger acceptance number perform better than commonly used acceptance sampling plans for attributes based on simple random sampling (SRS) technique. Indeed, MNS acceptance sampling plans result in OC curves which, compared to their SRS counterparts, are much closer to the ideal OC curve. A computer program is designed which can be used to specify the optimum MNS acceptance sampling plan and to show, visually, how the shape of the OC curve changes when parameters of the acceptance sampling plan vary. Theoretical results and numerical evaluations are given.