An Example Of The Use Of Andrews’ Plots To Detect Time Variations In Model Parameters And To Detect Outlying Observations

Andrews (1972) introduced a method of plotting highdimensional data in two dimensions. This method is exploited as a graphical tool for the examination of changes over time in the parameters of a time series model. An example using a Fourier series model is given to illustrate the method. It is also shown how outlying observations in the data can be found.

AMS (MOS) Subject Classifivations: 62M10, 62H30     

A tale of two countries: The Craig–Sakamoto–Matusita theorem

This paper has a long history. Junjiro Ogawa and I have been friends and collaborators for many years. We started this joint paper in the early 1990s and prepared a draft in September 1995. However, we were trying to improve and simplify some of the proofs, and we revised the paper in January 1998. But circumstances were such that we did not submit the paper, but planned to find further simplifications. Junjiro's untimely death in 2003 brought it forth again. Junjiro always liked this topic, and I am pleased to submit this joint paper in his memory. This paper contains both an historical summary as well as several new results. (Ingram Olkin)     

The Likelihood Ratio Test with the Box–Cox Transformation for the Normal Mixture Problem: Power and Sample Size Study

Abstract

Through simulation and regression, we study the alternative distribution of the likelihood ratio test in which the null hypothesis postulates that the data are from a normal distribution after a restricted Box–Cox transformation and the alternative hypothesis postulates that they are from a mixture of two normals after a restricted (possibly different) Box–Cox transformation. The number of observations in the sample is called N. The standardized distance between components (after transformation) is D = (μ₂ ? μ₁)/σ, where μ₁ and μ₂ are the component means and σ² is their common variance. One component contains the fraction π of observed, and the other 1 ? π. The simulation results demonstrate a dependence of power on the mixing proportion, with power decreasing as the mixing proportion differs from 0.5. The alternative distribution appears to be a non-central chi-squared with approximately 2.48 + 10N ^?0.75 degrees of freedom and non-centrality parameter 0.174N(D ? 1.4)² × [π(1 ? π)]. At least 900 observations are needed to have power 95% for a 5% test when D = 2. For fixed values of D, power, and significance level, substantially more observations are necessary when π ≥ 0.90 or π ≤ 0.10. We give the estimated powers for the alternatives studied and a table of sample sizes needed for 50%, 80%, 90%, and 95% power.     

Notes on Kolassa’s Method for Estimating the Power of Wilcoxon–Mann–Whitney Test

The Kolassa method implemented in the nQuery Advisor software has been widely used for approximating the power of the Wilcoxon–Mann–Whitney (WMW) test for ordered categorical data, in which Edgeworth approximation is used to estimate the power of an unconditional test based on the WMW U statistic. When the sample size is small or when the sizes in the two groups are unequal, Kolassa’s method may yield quite poor approximation to the power of the conditional WMW test that is commonly implemented in statistical packages. Two modifications of Kolassa’s formula are proposed and assessed by simulation studies.     

The Topp–Leone odd log-logistic family of distributions

We introduce a new class of continuous distributions named the Topp–Leone odd log-logistic family, which extends the one-parameter distribution pioneered by Topp and Leone [A family of J-shaped frequency functions. J Amer Statist Assoc. 1955;50:209–219]. We study some of its mathematical properties and describe two special cases. Further, we propose a regression model based on the new Topp–Leone odd log-logistic Weibull distribution. The usefulness and flexibility of the proposed family are illustrated by means of three real data sets.     

The mixture of left–right truncated normal distributions

This paper introduces the mixture of left–right truncated normal distributions, from the spreads between bid and ask prices, as a statistical model for handle non-normality of asset price returns. It has been proved that there is only one maximum for the likelihood function of the new model.     

The Risk of James–Stein and Lasso Shrinkage

This article compares the mean-squared error (or ?₂ risk) of ordinary least squares (OLS), James–Stein, and least absolute shrinkage and selection operator (Lasso) shrinkage estimators in simple linear regression where the number of regressors is smaller than the sample size. We compare and contrast the known risk bounds for these estimators, which shows that neither James–Stein nor Lasso uniformly dominates the other. We investigate the finite sample risk using a simple simulation experiment. We find that the risk of Lasso estimation is particularly sensitive to coefficient parameterization, and for a significant portion of the parameter space Lasso has higher mean-squared error than OLS. This investigation suggests that there are potential pitfalls arising with Lasso estimation, and simulation studies need to be more attentive to careful exploration of the parameter space.     

The Weibull–Dagum distribution: Properties and applications

ABSTRACT

In this article, we define a new lifetime model called the Weibull–Dagum distribution. The proposed model is based on the Weibull–G class. It can also be defined by a simple transformation of the Weibull random variable. Its density function is very flexible and can be symmetrical, left-skewed, right-skewed, and reversed-J shaped. It has constant, increasing, decreasing, upside-down bathtub, bathtub, and reversed-J shaped hazard rate. Various structural properties are derived including explicit expressions for the quantile function, ordinary and incomplete moments, and probability weighted moments. We also provide explicit expressions for the Rényi and q-entropies. We derive the density function of the order statistics as a mixture of Dagum densities. We use maximum likelihood to estimate the model parameters and illustrate the potentiality of the new model by means of a simulation study and two applications to real data. In fact, the proposed model outperforms the beta-Dagum, McDonald–Dagum, and Dagum models in these applications.     

The Rayleigh–Lindley model: properties and applications

In this paper, the Rayleigh–Lindley (RL) distribution is introduced, obtained by compounding the Rayleigh and Lindley discrete distributions, where the compounding procedure follows an approach similar to the one previously studied by Adamidis and Loukas in some other contexts. The resulting distribution is a two-parameter model, which is competitive with other parsimonious models such as the gamma and Weibull distributions. We study some properties of this new model such as the moments and the mean residual life. The estimation was approached via EM algorithm. The behavior of these estimators was studied in finite samples through a simulation study. Finally, we report two real data illustrations in order to show the performance of the proposed model versus other common two-parameter models in the literature. The main conclusion is that the model proposed can be a valid alternative to other competing models well established in the literature.     

The Poisson–exponential distribution: a Bayesian approach

In this paper, we proposed a new two-parameter lifetime distribution with increasing failure rate. The new distribution arises on a latent complementary risk scenario. The properties of the proposed distribution are discussed, including a formal proof of its density function and an explicit algebraic formulae for its quantiles and survival and hazard functions. Also, we have discussed inference aspects of the model proposed via Bayesian inference by using Markov chain Monte Carlo simulation. A simulation study investigates the frequentist properties of the proposed estimators obtained under the assumptions of non-informative priors. Further, some discussions on models selection criteria are given. The developed methodology is illustrated on a real data set.     

The Jackknife estimate of covariance of two Kaplan–Meier integrals with covariables

In this paper the Jackknife estimate of covariance of two Kaplan–Meier integrals with covariates is introduced. Its strong consistency is established under mild conditions. Several applications of the estimator are discussed.     

Comparison of Statistical Methods for Pretest–Posttest Designs in Terms of Type I Error Probability and Statistical Power

The pretest–posttest design is widely used to investigate the effect of an experimental treatment in biomedical research. The treatment effect may be assessed using analysis of variance (ANOVA) or analysis of covariance (ANCOVA). The normality assumption for parametric ANOVA and ANCOVA may be violated due to outliers and skewness of data. Nonparametric methods, robust statistics, and data transformation may be used to address the nonnormality issue. However, there is no simultaneous comparison for the four statistical approaches in terms of empirical type I error probability and statistical power. We studied 13 ANOVA and ANCOVA models based on parametric approach, rank and normal score-based nonparametric approach, Huber M-estimation, and Box–Cox transformation using normal data with and without outliers and lognormal data. We found that ANCOVA models preserve the nominal significance level better and are more powerful than their ANOVA counterparts when the dependent variable and covariate are correlated. Huber M-estimation is the most liberal method. Nonparametric ANCOVA, especially ANCOVA based on normal score transformation, preserves the nominal significance level, has good statistical power, and is robust for data distribution.     

The Little Magazine Collection at the University of Wisconsin–Madison

Abstract

Serials Review visits the University of Wisconsin–Madison Memorial Library to hear about the Little Magazine Collection. This collection consists primarily of literary journals, but it also contains magazines, which reflect trends in art, society, and politics. Readers of Serials Review will recognize this collection as the source of a regular column, “Little Magazine Interview Index,” compiled by Yvonne Schofer and Barbara Richards. The author discusses the history of the collection and its current development. Serials Review 2002; 28:153–155.     

Power comparison of the Kolmogorov–Smirnov test under ranked set sampling and simple random sampling

Many studies have been used to compare the power of several goodness-of-fit (GOF) tests under simple random sampling (SRS) and ranked set sampling (RSS). In our study, a different design procedure and ranking process in RSS are thoroughly investigated. A simulation study is conducted to compare the power of the Kolmogorov–Smirnov test under SRS and RSS with different sets and cycle sizes for several distributions. Level-2 sampling design and partially rank-ordered sets are used. Also, we benefited from auxiliary variables in the ranking process. Finally, results are presented with tables and figures. Under these conditions we show that the RSS has better performance against the SRS in finite population.     

The precision of regression-type estimator for incremental cost–effectiveness ratio

The estimation of incremental cost–effectiveness ratio (ICER) has received increasing attention recently. It is expressed in terms of the ratio of the change in costs of a therapeutic intervention to the change in the effects of the intervention. Despite the intuitive interpretation of ICER as an additional cost per additional benefit unit, it is a challenge to estimate the distribution of a ratio of two stochastically dependent distributions. A vast literature regarding the statistical methods of ICER has developed in the past two decades, but none of these methods provide an unbiased estimator. Here, to obtain the unbiased estimator of the cost–effectiveness ratio (CER), the zero intercept of the bivariate normal regression is assumed. In equal sample sizes, the Iman–Conover algorithm is applied to construct the desired variance–covariance matrix of two random bivariate samples, and the estimation then follows the same approach as CER to obtain the unbiased estimator of ICER. The bootstrapping method with the Iman–Conover algorithm is employed for unequal sample sizes. Simulation experiments are conducted to evaluate the proposed method. The regression-type estimator performs overwhelmingly better than the sample mean estimator in terms of mean squared error in all cases.     

Gambler’s Ruin: The Duration of Play

We study the gambler’s ruin problem with a general distribution of the payoffs in each game. Assuming the expected value of the payoff distribution is negative, so that eventual ruin occurs with probability 1, we are interested in the distribution of the duration to ruin, also known as the first-passage time distribution. A generating function for this distribution is obtained. Exact expressions for the expected value and variance of this distribution, as well as asymptotic expressions for the case of large initial wealth, are derived.     

The alpha–beta skew normal distribution: properties and applications

This paper introduces a new class of skew distributions by extending the alpha skew normal distribution proposed by Elal-Olivero [Elal-Olivero, D. Alpha-skew-normal distribution. Proyecciones. 2010;29:224–240]. Statistical properties of the new family are studied in details. In particular, explicit expressions for the moments and the shape parameters including the skewness and the kurtosis coefficients and the moment generating function are derived. The problem of estimating parameters on the basis of a random sample coming from the new class of distribution is considered. To examine the performance of the obtained estimators, a Monte Carlo simulation study is conducted. Flexibility and usefulness of the proposed family of distributions are illustrated by analysing three real data sets.     

The inverted Nadarajah–Haghighi distribution: estimation methods and applications

The inverted (or inverse) distributions are sometimes very useful to explore additional properties of the phenomenons which non-inverted distributions cannot. We introduce a new inverted model called the inverted Nadarajah–Haghighi distribution which exhibits decreasing and unimodal (right-skewed) density while the hazard rate shapes are decreasing and upside-down bathtub. Our main focus is the estimation (from both frequentist and Bayesian points of view) of the unknown parameters along with some mathematical properties of the new model. The Bayes estimators and the associated credible intervals are obtained using Markov Chain Monte Carlo techniques under squared error loss function. The gamma priors are adopted for both scale and shape parameters. The potentiality of the distribution is analysed by means of two real data sets. In fact, it is found to be superior in its ability to sufficiently model the data as compared to the inverted Weibull, inverted Rayleigh, inverted exponential, inverted gamma, inverted Lindley and inverted power Lindley models.