On the pade tablé and its relationship to the r and sarrays and arm a modeling

The Box-Jenkins method is a popular and important technique for modeling and forecasting of time series. Unfortunately the problem of determining the appropriate ARMA forecasting model (or indeed if an ARMA model holds) is a major drawback to the use of the Box-Jenkins methodology. Gray et al. (1978) and Woodward and Gray (1979) have proposed methods of estimating p and qin ARMA modeling based on the R and Sarrays that circumvent some of these modeling difficulties.

In this paper we generalize the R and S arrays by showing a relationship to Padé approximunts and then show that these arrays have a much wider application than in just determining model order. Particular non-ARMA models can be identified as well. This includes certain processes that consist of deterministic functions plus ARMA noise, indeed we believe that the combined R and S arrays are the best overall tool so fur developed for the identification of general 2nd order (not just stationary) time scries models.     

Critical values of Mandel’s h and k, the Grubbs and the Cochran test statistic

According to ISO 5725-2 (1994), measurement results obtained in an interlaboratory experiment are inspected for consistency by plotting Mandel’s h and k statistics and for outliers by application of the Grubbs test and the Cochran test. Critical values of these statistics for significance levels α=5% and α=1% and for some numbers p of laboratories and n of repeated measurements in the laboratories are supplied in ISO 5725-2 without reference to methods for their calculation. In this paper, exact formulae for the critical values of Mandel’s h and k and approximate formulae for the critical values of the Single Grubbs test, the Double Grubbs test and the Cochran test are derived.     

A note on the equality of the OLSE and the BLUE of the parametric function in the general Gauss–Markov model

In this note we consider the equality of the ordinary least squares estimator (OLSE) and the best linear unbiased estimator (BLUE) of the estimable parametric function in the general Gauss–Markov model. Especially we consider the structures of the covariance matrix V for which the OLSE equals the BLUE. Our results are based on the properties of a particular reparametrized version of the original Gauss–Markov model.      

Generalized Concept of Relative Risk and Wider Applications of the Proportional Hazards Model and the Kaplan–Meier Estimator

The concepts of relative risk and hazard ratio are generalized for ordinary ordinal and continuous response variables, respectively. Under the generalized concepts, the Cox proportional hazards model with the Breslow's and Efron's methods can be regarded as generalizations of the Mantel–Haenszel estimator for dealing with broader types of covariates and responses. When ordinal responses can be regarded as discretized observations of a hypothetical continuous variable, the estimated relative risks from the Cox model reflect the associations between the responses and covariates. Examples are given to illustrate the generalized concepts and wider applications of the Cox model and the Kaplan–Meier estimator.     

Is the Selected Population the Best?—Location and Scale Parameter Cases

We observe X ₁,…,X _k , where X _i has density f(x,θ _i ) possessing monotone likelihood ratio. The best population corresponds to the largest θ _i . We select the population corresponding to the largest X _i . The goal is to attach the best possible p-value to the inference: the selected population has the uniquely largest θ _i . Gutmann and Maymin (1987 Gutmann , S. , Maymin , Z. ( 1987 ). Is the selected population the best? Ann. Statist . 15 : 456 – 461 .[Crossref], [Web of Science ®] , [Google Scholar]) considered the location parameter case and derived the supremum of the error probability by conditioning on S, the index of the largest X _i . Using this conditioning approach, Kannan and Panchapakesan (2009 Kannan , N. , Panchapakesan , S. ( 2009 ). Does the selected normal population have the smallest variance? Amer. J. Math. Management Sci . 29 : To appear . [Google Scholar]) considered the problem for the gamma family. We consider here a unified approach to both the location and scale parameter cases, and obtain the supremum of the error probability without using conditioning.     

Stockouts and Restocking: Monitoring the Retailer from the Supplier’s Perspective

Suppliers and retailers typically do not have identical incentives to avoid stockouts (lost sales due to the lack of product availability on the shelf). Thus, the supplier needs to monitor the retailer’s restocking efforts with the available data. We empirically assess stockout levels using only shipment and sales data that is readily available to the supplier. The model distinguishes between store stockouts (zero inventory in the store) and shelf stockouts (an empty shelf but some inventory in other parts of the store), thereby identifying the cause of the stockout to be either a supply chain or a restocking issue. We find that, as suspected by the supplier, the average stockout rate is much higher than published averages. In addition, stockout rates vary widely between stores. Moreover, almost all stockouts are shelf stockouts. The model identifies stores that may have restocking issues.     

Mixture Models and the Krätzel Integral Transform

Continuous mixture Weibull models arise in many areas of sciences such as reliability studies, communications theory, etc. Due to its wide applicability, we introduce a class of continuous mixture Weibull models which is a combination of Weibull and generalized gamma distributions. Some characteristics of the distribution are obtained. It is seen that Krätzel integral enters into the model naturally, and then the model can be called as a Krätzel density. Applications of the density function related to fading channels and ultrasonic backscatter signals modeling are discussed. A real data analysis is given to illustrate the use of this distribution.     

The Frisch–Waugh–Lovell theorem for the lasso and the ridge regression

The Frisch–Waugh–Lovell (FWL) (partitioned regression) theorem is essential in regression analysis. This is partly because it is quite useful to derive theoretical results. The lasso regression and the ridge regression, both of which are penalized least-squares regressions, have become popular statistical techniques. This article describes that the FWL theorem remains valid for these penalized least-squares regressions. More precisely, we demonstrate that the covariates corresponding to unpenalized regression parameters in these penalized least-squares regression can be projected out. Some other results related to the FWL theorem in such penalized least-squares regressions are also presented.     

Understanding and Addressing the Unbounded “Likelihood” Problem

The joint probability density function, evaluated at the observed data, is commonly used as the likelihood function to compute maximum likelihood estimates. For some models, however, there exist paths in the parameter space along which this density-approximation likelihood goes to infinity and maximum likelihood estimation breaks down. In all applications, however, observed data are really discrete due to the round-off or grouping error of measurements. The “correct likelihood” based on interval censoring can eliminate the problem of an unbounded likelihood. This article categorizes the models leading to unbounded likelihoods into three groups and illustrates the density-approximation breakdown with specific examples. Although it is usually possible to infer how given data were rounded, when this is not possible, one must choose the width for interval censoring, so we study the effect of the round-off on estimation. We also give sufficient conditions for the joint density to provide the same maximum likelihood estimate as the correct likelihood, as the round-off error goes to zero.     

On the existence and uniqueness of the maximum likelihood estimates of the parameters of Birnbaum–Saunders distribution based on Type-I,Type-II and hybrid censored samples

The Birnbaum–Saunders (BS) distribution is a positively skewed distribution and is a common model for analysing lifetime data. In this paper, we discuss the existence and uniqueness of the maximum likelihood estimates (MLEs) of the parameters of BS distribution based on Type-I, Type-II and hybrid censored samples. The line of proof is based on the monotonicity property of the likelihood function. We then describe the numerical iterative procedure for determining the MLEs of the parameters, and point out briefly some recently developed simple methods of estimation in the case of Type-II censoring. Some graphical illustrations of the approach are given for three real data from the reliability literature. Finally, for illustrative purpose, we also present an example in which the MLEs do not exist.     

Linear sufficiency and completeness in the context of estimating the parametric function in the general Gauss–Markov model

In this paper we consider linear sufficiency and linear completeness in the context of estimating the estimable parametric function K^′β${\mathbf{K}}^{\prime }\mathit{\beta }$ under the general Gauss–Markov model {y,Xβ,σ²V}$\{\mathbf{y},\mathbf{X}\mathit{\beta },{\sigma }^2\mathbf{V}\}$. We give new characterizations for linear sufficiency, and define and characterize linear completeness in a case of estimation of K^′β${\mathbf{K}}^{\prime }\mathit{\beta }$. Also, we consider a predictive approach for obtaining the best linear unbiased estimator of K^′β${\mathbf{K}}^{\prime }\mathit{\beta }$, and subsequently, we give the linear analogues of the Rao–Blackwell and Lehmann–Scheffé Theorems in the context of estimating K^′β${\mathbf{K}}^{\prime }\mathit{\beta }$.     

Discussion on the paper “On simulation and properties of the stable law” by L. Devroye and L. James

We congratulate the authors for the interesting paper. The reading has been really pleasant and instructive. We discuss briefly only some of the interesting results given in Devroye and James (Stat Methods Appl 2014) with particular attention to evolution problems. The contribution of the results collected in the paper is useful in a more wide class of applications in many areas of applied mathematics.     

Harrison–Stevens Forecasting and the Multiprocess Dynamic Linear Model

This is an expository article. The Harrison–Stevens forecasting algorithm using the multiprocess dynamic linear model is a robust method for forecasting in a nonstationary time series. The purpose of this article is to help statisticians become familiar with the method.     

A Cautionary Note on the Use of the Grønnesby and Borgan Goodness-of-Fit Test for the Cox Proportional Hazards Model

Grønnesby and Borgan (1996, Lifetime Data Analysis 2, 315–328) propose an omnibus goodness-of-fit test for the Cox proportional hazards model. The test is based on grouping the subjects by their estimated risk score and comparing the number of observed and a model based estimated number of expected events within each group. We show, using extensive simulations, that even for moderate sample sizes the choice of number of groups is critical for the test to attain the specified size. In light of these results we suggest a grouping strategy under which the test attains the correct size even for small samples. The power of the test statistic seems to be acceptable when compared to other goodness-of-fit tests.     

Open Access and Promotion and Tenure Evaluation Plans at the University of Wisconsin–Eau Claire

Department and program evaluation plans at the University of Wisconsin–Eau Claire were examined to see if these documents provide evidence that could be used to justify supporting the publication of peer-reviewed open access articles toward tenure and promotion. In an earlier study, the authors reveal that faculty members at the University of Wisconsin–Eau Claire are more unaware of open access publishing than their counterparts at larger universities. These findings dovetail with other studies that show that faculty members are reluctant to publish in open access journals because of concerns about the quality of those journals. The existing body of scholarship suggests that tenure-line faculty fear publishing in open access journals because it could adversely impact their chances of promotion and tenure. The authors of this current study sought to determine if department and program evaluation plans could influence negative perceptions faculty have of open access journals. The implications of this study for librarians, scholarly communication professionals, tenure-line faculty, departments, and programs are addressed.     

A powerful and interpretable alternative to the Jarque–Bera test of normality based on 2nd-power skewness and kurtosis,using the Rao's score test on the APD family

We introduce the 2nd-power skewness and kurtosis, which are interesting alternatives to the classical Pearson's skewness and kurtosis, called 3rd-power skewness and 4th-power kurtosis in our terminology. We use the sample 2nd-power skewness and kurtosis to build a powerful test of normality. This test can also be derived as Rao's score test on the asymmetric power distribution, which combines the large range of exponential tail behavior provided by the exponential power distribution family with various levels of asymmetry. We find that our test statistic is asymptotically chi-squared distributed. We also propose a modified test statistic, for which we show numerically that the distribution can be approximated for finite sample sizes with very high precision by a chi-square. Similarly, we propose a directional test based on sample 2nd-power kurtosis only, for the situations where the true distribution is known to be symmetric. Our tests are very similar in spirit to the famous Jarque–Bera test, and as such are also locally optimal. They offer the same nice interpretation, with in addition the gold standard power of the regression and correlation tests. An extensive empirical power analysis is performed, which shows that our tests are among the most powerful normality tests. Our test is implemented in an R package called PoweR.