Measuring agreement between two statistics with applications to age standardization

Summary.  We detail a general method for measuring agreement between two statistics. An application is two ratios of directly standardized rates which differ only by the choice of the standard. If the statistics have a high value for the coefficient of agreement then the expected squared difference between the statistics is small relative to the variance of the average of the two statistics, and inferences vary little by changing statistics. The estimation of a coefficient of agreement between two statistics is not straightforward because there is only one pair of observed values, each statistic calculated from the data. We introduce estimators of the coefficient of agreement for two statistics and discuss their use, especially as applied to functions of standardized rates.     

A Comparison of Partially Adaptive and Reweighted Least Squares Estimation

The small sample performance of least median of squares, reweighted least squares, least squares, least absolute deviations, and three partially adaptive estimators are compared using Monte Carlo simulations. Two data problems are addressed in the paper: (1) data generated from non-normal error distributions and (2) contaminated data. Breakdown plots are used to investigate the sensitivity of partially adaptive estimators to data contamination relative to RLS. One partially adaptive estimator performs especially well when the errors are skewed, while another partially adaptive estimator and RLS perform particularly well when the errors are extremely leptokur-totic. In comparison with RLS, partially adaptive estimators are only moderately effective in resisting data contamination; however, they outperform least squares and least absolute deviation estimators.     

La médiane simpliciale d'Oja: Existence,unicité et stabilité

Oja (1983) examined various ways of measuring location, scatter, skewness, and kurtosis for multivariate distributions. Among other measures of location, he introduced a generalised median known in this paper under the name of the Oja median. In our study of the existence of that median, we show that Oja's definition can only be applied to distributions having a mean. In dimension d θ 2, we establish that the usual method of extension breaks down, which raises the question of the validity of the concept as a notion of median. Two fundamental theoretical properties of that median are also considered: uniqueness and consistency.     

A Weighted Least Squares Approach to Levene's Test of Homogeneity of Variance

In 1960 Levene suggested a potentially robust test of homogeneity of variance based on an ordinary least squares analysis of variance of the absolute values of mean-based residuals. Levene's test has since been shown to have inflated levels of significance when based on the F-distribution, and tests a hypothesis other than homogeneity of variance when treatments are unequally replicated, but the incorrect formulation is now standard output in several statistical packages. This paper develops a weighted least squares analysis of variance of the absolute values of both mean-based and median-based residuals. It shows how to adjust the residuals so that tests using the F -statistic focus on homogeneity of variance for both balanced and unbalanced designs. It shows how to modify the F -statistics currently produced by statistical packages so that the distribution of the resultant test statistic is closer to an F-distribution than is currently the case. The weighted least squares approach also produces component mean squares that are unbiased irrespective of which variable is used in Levene's test. To complete this aspect of the investigation the paper derives exact second-order moments of the component sums of squares used in the calculation of the mean-based test statistic. It shows that, for large samples, both ordinary and weighted least squares test statistics are equivalent; however they are over-dispersed compared to an F variable.     

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     

Online Control Charts for Process Averages Based on Repeated Median Filters

Two types of estimates of process level, namely repeated median estimates (Siegel, 1982 Siegel , A. F. ( 1982 ). Robust regression using repeated medians . Biometrika 69 : 242 – 244 .[Crossref], [Web of Science ®] , [Google Scholar]) and full online estimates (Gather et al., 2006 Gather , U. , Schettlinger , K. , Fried , R. ( 2006 ). Online signal extraction by robust linear regression . Computational Statistics 21 : 33 – 51 .[Crossref], [Web of Science ®] , [Google Scholar]) based on repeated median filters, are used to develop control charts. The distributional properties of the estimates are studied using simulation and these are found to closely follow normal distribution. The repeated median being robust against outliers with asymptotically 50% breakdown value and having small standard deviation is found to be useful as a basis for monitoring process averages. The control charts using repeated median estimates have been recommended for general use.     

Generalized Laplacian Distributions and Autoregressive Processes

Generalized Laplacian distribution is considered. A new distribution called geometric generalized Laplacian distribution is introduced and its properties are studied. First- and higher-order autoregressive processes with these stationary marginal distributions are developed and studied. Simulation studies are conducted and trajectories of the process are obtained for selected values of the parameters. Various areas of application of these models are discussed.     

The Joint Distribution of the Sum and the Maximum of IID Exponential Random Variables

The authors establish the joint distribution of the sum X and the maximum Y of IID exponential random variables. They derive exact formuli describing the random vector (X, Y), including its joint PDF, CDF, and other characteristics; marginal and conditional distributions; moments and related parameters; and stochastic representations leading to further properties of infinite divisibility and self-decomposability. The authors also discuss parameter estimation and include an example from climatology that illustrates the modeling potential of this new bivariate model.     

Normal Deviation and Poisson Approximation of a Security Market by the Geometric Markov Renewal Processes

We consider the geometric Markov renewal processes (GMRP) as a model for a security market. Normal deviations of the geometric Markov renewal processes for ergodic averaging and double averaging schemes are derived. We introduce Poisson averaging scheme for the geometric Markov renewal processes. European call option pricing formulas for GMRP are presented.     

Inference Under Right Censoring in a Discrete Setup

In this article, we consider the right random censoring scheme in a discrete setup when the lifetime and censoring variables are independent and have geometric distributions with means 1/θ₁ and 1/θ₂, respectively. We first obtain the Maximum Likelihood and Method of Moment estimators of the unknown parameters. We also find the Bayes and Posterior Regret Gamma Minimax estimators of the parameters for the two cases when the prior distributions are dependent and independent, assuming a squared error loss function. We then discuss the Proportional Hazard model, and obtain Maximum Likelihood estimators of the unknown parameters and derive the Bayes estimators assuming squared error loss using Markov Chain Monte Carlo methods.