Robust location estimation with missing data

In a missing data setting, we have a sample in which a vector of explanatory variables ${\bf x}_i$ is observed for every subject i, while scalar responses $y_i$ are missing by happenstance on some individuals. In this work we propose robust estimators of the distribution of the responses assuming missing at random (MAR) data, under a semiparametric regression model. Our approach allows the consistent estimation of any weakly continuous functional of the response's distribution. In particular, strongly consistent estimators of any continuous location functional, such as the median, L‐functionals and M‐functionals, are proposed. A robust fit for the regression model combined with the robust properties of the location functional gives rise to a robust recipe for estimating the location parameter. Robustness is quantified through the breakdown point of the proposed procedure. The asymptotic distribution of the location estimators is also derived. The proofs of the theorems are presented in Supplementary Material available online. The Canadian Journal of Statistics 41: 111–132; 2013 © 2012 Statistical Society of Canada     

A Mechanical Interpretation of the Coefficient of Rank Correlation and other Analogies

Periodically, the pyramid or “chain letter” scheme is offered to Americans under the guise of a business dealership. Recently, the FTC ordered Glen Turner's “Dare to be Great” firm to repay 44 million dollars to participants. In order to demonstrate that the potential gains are misrepresented by promoters, a probability model of the pyramid scheme is developed. The major implications are that the vast majority of participants have less than a ten percent chance of recouping their initial investment when a small profit is achieved as soon as they recruit three people and that, on the average, half of the participants will recruit no one else and lose all their money.     

长三角科技资源配置能力与城市化进程的协调耦合关系研究

科技资源配置能力已成为长江地区城市化进程的重要推力。因此,在构建科技资源与城市化耦合系统评价指标体系的基础上,利用耦合协调度模型综合评价长江区域科技资源配置与城市化的交互耦合关系。结果显示:从时间维度看,研究期内长三角科技资源配置能力与城市化进程之间的耦合关系基本上处于磨合阶段,耦合协调度以低度协调耦合状态为主并逐步过渡到中度协调耦合状态;从空间维度看,长三角各城市的科技资源配置与城市化耦合协调度在极度协调耦合、高度协调耦合、中度协调耦合、低度协调耦合四个阶段均有分布,层级分明且大部分城市尚处于低度协调耦合状态。     

基于灰色关联的西部煤炭产业集中度影响因素研究

在概述西部原煤产量及其矿井类型与数量变化趋势的基础上,利用2003-2010年相关时序数据测度了西部煤炭产业在这一时期的集中度。随后,运用灰色关联理论,分析了期初集中度、生产扩张能力、进入壁垒、市场容量等相关因素与西部煤炭产业集中度的灰色关联程度。研究表明:期初集中度与市场容量是影响西部煤炭产业集中度诸多因素中最为重要的两个因素。在进一步分析产生这一结果深层原因的基础上,提出了提高西部煤炭产业集中度的政策建议。     

Mean estimation with data missing at random for functional covariables

In a missing-data setting, we want to estimate the mean of a scalar outcome, based on a sample in which an explanatory variable is observed for every subject while responses are missing by happenstance for some of them. We consider two kinds of estimates of the mean response when the explanatory variable is functional. One is based on the average of the predicted values and the second one is a functional adaptation of the Horvitz–Thompson estimator. We show that the infinite dimensionality of the problem does not affect the rates of convergence by stating that the estimates are root-n consistent, under missing at random (MAR) assumption. These asymptotic features are completed by simulated experiments illustrating the easiness of implementation and the good behaviour on finite sample sizes of the method. This is the first paper emphasizing that the insensitiveness of averaged estimates, well known in multivariate non-parametric statistics, remains true for an infinite-dimensional covariable. In this sense, this work opens the way for various other results of this kind in functional data analysis.     

On the use of curve fitting to model the error of the Cornish-Fisher expansion of the pearson type-VI distribution

The Cornish-Fisher expansion of the Pearson type VI distribution is known to be reasonably accurate when both degrees of freedom are relatively large (say greater than or equal to 5). However, when either or both degrees of freedom are less than 5, the accuracy of the computed percentage point begins to suffer; in some cases severely. To correct for this, the error surface in the degrees of freedom plane is modeled by least squares curve fitting for selected levels of tail probability (.025, .05, and .10) which can be used to adjust the percentage point obtained from the usual Cornish-Fisher expansion. This adjustment procedure produces a computing algorithm that computes percentage points of the Pearson type VI distribution at the above probability levels, accurate to at least + 1 in 3 digits in approximately 11 milliseconds per subroutine call on an IBM 370/145. This adjusted routine is valid for both degrees of freedom greater than or equal to 1.     

Variable Selection in Joint Location and Scale Models of the Skew-t-Normal Distribution

Variable selection is an important issue in all regression analysis, and in this article, we investigate the simultaneous variable selection in joint location and scale models of the skew-t-normal distribution when the dataset under consideration involves heavy tail and asymmetric outcomes. We propose a unified penalized likelihood method which can simultaneously select significant variables in the location and scale models. Furthermore, the proposed variable selection method can simultaneously perform parameter estimation and variable selection in the location and scale models. With appropriate selection of the tuning parameters, we establish the consistency and the oracle property of the regularized estimators. These estimators are compared by simulation studies.     

Tests of hypotheses based on ranks in the general linear model

A unified approach is developed for testing hypotheses in the general linear model based on the ranks of the residuals. It complements the nonparametric estimation procedures recently reported in the literature. The testing and estimation procedures together provide a robust alternative to least squares. The methods are similar in spirit to least squares so that results are simple to interpret. Hypotheses concerning a subset of specified parameters can be tested, while the remaining parameters are treated as nuisance parameters. Asymptotically, the test statistic is shown to have a chi-square distribution under the null hypothesis. This result is then extended to cover a sequence of contiguous alternatives from which the Pitman efficacy is derived. The general application of the test requires the consistent estimation of a functional of the underlying distribution and one such estimate is furnished.     

Simulation of the p-generalized Gaussian distribution

In this paper, we introduce the p-generalized polar methods for the simulation of the p-generalized Gaussian distribution. On the basis of geometric measure representations, the well-known Box–Muller method and the Marsaglia–Bray rejecting polar method for the simulation of the Gaussian distribution are generalized to simulate the p-generalized Gaussian distribution, which fits much more flexibly to data than the Gaussian distribution and has already been applied in various fields of modern sciences. To prove the correctness of the p-generalized polar methods, we give stochastic representations, and to demonstrate their adequacy, we perform a comparison of six simulation techniques w.r.t. the goodness of fit and the complexity. The competing methods include adapted general methods and another special method. Furthermore, we prove stochastic representations for all the adapted methods.