Biased Bootstrap Methods for Reducing the Effects of Contamination

Contamination of a sampled distribution, for example by a heavy-tailed distribution, can degrade the performance of a statistical estimator. We suggest a general approach to alleviating this problem, using a version of the weighted bootstrap. The idea is to 'tilt' away from the contaminated distribution by a given (but arbitrary) amount, in a direction that minimizes a measure of the new distribution's dispersion. This theoretical proposal has a simple empirical version, which results in each data value being assigned a weight according to an assessment of its influence on dispersion. Importantly, distance can be measured directly in terms of the likely level of contamination, without reference to an empirical measure of scale. This makes the procedure particularly attractive for use in multivariate problems. It has several forms, depending on the definitions taken for dispersion and for distance between distributions. Examples of dispersion measures include variance and generalizations based on high order moments. Practicable measures of the distance between distributions may be based on power divergence, which includes Hellinger and Kullback–Leibler distances. The resulting location estimator has a smooth, redescending influence curve and appears to avoid computational difficulties that are typically associated with redescending estimators. Its breakdown point can be located at any desired value ε∈ (0, ½) simply by 'trimming' to a known distance (depending only on ε and the choice of distance measure) from the empirical distribution. The estimator has an affine equivariant multivariate form. Further, the general method is applicable to a range of statistical problems, including regression.     

蛛网图评价法的计算机实现

论述了蛛网图评价法的思想，给出了利用Ｆｏｘｐｒｏ和Ｖｉｓｉｕａｌ Ｂａｓｉｃ混合编程的方法进行的程序设计。     

Some sampling effects of pairwise correlated observations on likelihood ratio tests for the difference between two means

We study the effects of the inclusion of pairs of correlated observations in a sample on likelihood ratio tests for the difference in two means. In particular, we assess how the inclusion of correlated data pairs (e.g., such as data inadvertently collected from sib-pairs) affects the sample size requirements necessary for the implementation of a Likelihood Ratio (LR) test for the difference between two means. Our results suggest that correlated data pairs beneficially or adversely effect sample size requirements for an LR test to a degree functionally related to the mixture parameters dictating their relative frequencies in the larger sample on which the test will be performed, the strength of the correlation between the observations, and the size of imbalances in the sample with respect to the number of observations in each group. The relevance and implications of the results for genetic and epidemiologic research are discussed.     

Estimation of the smaller and larger of two exponential location parameters

Assume independent random samples are drawn from two populations which are exponentially distributed with unknown location parameters and a common known scale parameter. We want to estimate the maximum and the minimum of the unknowo location paremeters. In this paper several estimators are proposed which are better than the natural estimations in terms of absolute bias and /or meaqn squared error.     

Influence of human immunodeficiency virus infection on neurological impairment: an analysis of longitudinal binary data with informative drop-out

A study to investigate the human immunodeficiency virus (HIV) status on the course of neurological impairment, conducted by the HIV Center at Columbia University, followed a cohort of HIV positive and negative gay men for 5 years and assessed the presence or absence of neurological impairment every 6 months. Almost half of the subjects dropped out before the end of the study for reasons that might have been related to the missing neurological data. We propose likelihood-based methods for analysing such binary longitudinal data under informative and non-informative drop-out. A transition model is assumed for the binary response, and several models for the drop-out processes are considered which are functions of the response variable (neurological impairment). The likelihood ratio test is used to compare models with informative and non-informative drop-out mechanisms. Using simulations, we investigate the percentage bias and mean-squared error (MSE) of the parameter estimates in the transition model under various assumptions for the drop-out. We find evidence for informative drop-out in the study, and we illustrate that the bias and MSE for the parameters of the transition model are not directly related to the observed drop-out or missing data rates. The effect of HIV status on the neurological impairment is found to be statistically significant under each of the models considered for the drop-out, although the regression coefficient may be biased in certain cases. The presence and relative magnitude of the bias depend on factors such as the probability of drop-out conditional on the presence of neurological impairment and the prevalence of neurological impairment in the population under study.     

Estimation of parameters in incomplete data models defined by dynamical systems

Parametric incomplete data models defined by ordinary differential equations (ODEs) are widely used in biostatistics to describe biological processes accurately. Their parameters are estimated on approximate models, whose regression functions are evaluated by a numerical integration method. Accurate and efficient estimations of these parameters are critical issues. This paper proposes parameter estimation methods involving either a stochastic approximation EM algorithm (SAEM) in the maximum likelihood estimation, or a Gibbs sampler in the Bayesian approach. Both algorithms involve the simulation of non-observed data with conditional distributions using Hastings–Metropolis (H–M) algorithms. A modified H–M algorithm, including an original local linearization scheme to solve the ODEs, is proposed to reduce the computational time significantly. The convergence on the approximate model of all these algorithms is proved. The errors induced by the numerical solving method on the conditional distribution, the likelihood and the posterior distribution are bounded. The Bayesian and maximum likelihood estimation methods are illustrated on a simulated pharmacokinetic nonlinear mixed-effects model defined by an ODE. Simulation results illustrate the ability of these algorithms to provide accurate estimates.     

A Continuous Representation of the Family of Stable Law Distributions

Conventional parametric representations of stable law distributions do not allow all members of the family to be obtained as continuous limits of the parameters. Model building (or simulation) using such representations will be numerically unstable near such limits in consequence. Existing tables are not satisfactory near such limits as interpolation cannot be carried out. We show that these difficulties are overcome by using a new shifted Cartesian representation which characterizes the entire stable law family in a completely continuous way. Standardization is still possible with this representation so that tabulation, using just two bounded parameters, can be carried out. Its use is illustrated in a non-regular threshold estimation problem involving stable distributions which are discontinuous limits in conventional representations.     

一种基于冗余小波变换的DT网格运动估计和运动补偿方法

网格划分和网格顶点的运动估计是基于不规则网格的视频压缩技术的关键。为了进一步提高网格运动估计和运动补偿的效果,在综合比较现有冗余小波变换域运动估计方法和适用规则网格的EMRMC算法的基础上,提出了一种新的基于不规则网格的运动估计和运动补偿算法,即在冗余小波变换域提取特征点和运动潜在区,网格顶点的运动估计采用结合运动潜在区的在时域进行块匹配的运动估计和运动补偿方法,而运动补偿则通过三角形仿射变换完成。同时还给出了冗余小波变换域提取运动潜在区的计算模板。理论分析和实验结果表明,该算法在补偿效果方面较前两种方法得到了改进。     

完全检查下预防维修模型建立及案例分析

为了解决维修间隔期确定的难题,建立了维修间隔期和总停机时间之间关系的目标函数,以及时间延迟维修模型.根据故障记录数据和预防维修活动的检查数据,采用最大似然法估计有关缺陷发生率和时间延迟分布等参数.在所建立目标函数和估计参数的基础上,计算出最优的维修间隔期.     

Adaptive interval estimation in one-way random effects models

Two-stage procedures are introduced to control the width and coverage (validity) of confidence intervals for the estimation of the mean, the between groups variance component and certain ratios of the variance components in one-way random effects models. The procedures use the pilot sample data to estimate an “optimal” group size and then proceed to determine the number of groups by a stopping rule. Such sampling plans give rise to unbalanced data, which are consequently analyzed by the harmonic mean method. Several asymptotic results concerning the proposed procedures are given along with simulation results to assess their performance in moderate sample size situations. The proposed procedures were found to effectively control the width and probability of coverage of the resulting confidence intervals in all cases and were also found to be robust in the presence of missing observations. From a practical point of view, the procedures are illustrated using a real data set and it is shown that the resulting unbalanced designs tend to require smaller sample sizes than is needed in a corresponding balanced design where the group size is arbitrarily pre-specified.