Maximum likelihood estimation in a Weibull regression model with type-1 censoring: a Monte Carlo study

Results of the Monte Carlo study of the performance of a maximum likelihood estimation in a Weibull parametric regression model with two explanatory variables are presented. One simulation run contained 1000 samples censored on the average by the amount of 0-30%. Each simulatedsample was generated in a form of two-factor two-level balanced experiment. The confidence intervals were computed using the large-sample normal approximation via the matrix of observed information. For small sample sizes the estimates of the scale parameter b of the loglifetime were significantly negatively biased, which resulted in a poor quality of confidence intervals for b and the low-level quantiles. All estimators improved their quality when the nominal value of b decreased. A moderate amount of censoring improved the quality of point and confidence estimation. The reparametrization b 7 produced rather accurate confidence intervals. Exact confidence intervals for b in case of non-censoring were obtained using the pivotal quantity b/b.     

POOLED SLICING METHODS VERSUS SLICING METHODS

We consider the semiparametric regression model introduced by [1] Duan, N. and Li, K. C. 1991. Slicing regression: a link-free regression method. The Annals of Statistics, 19: 505–530. [Crossref], [Web of Science ®] , [Google Scholar]. The dependent variable y is linked to the index x′ β through an unknown link function. [1] Duan, N. and Li, K. C. 1991. Slicing regression: a link-free regression method. The Annals of Statistics, 19: 505–530. [Crossref], [Web of Science ®] , [Google Scholar] and [2] Li, K. C. 1991. Sliced inverse regression for dimension reduction, with discussions. Journal of the American Statistical Association, 86: 316–342. [Taylor & Francis Online], [Web of Science ®] , [Google Scholar] present Slicing methods (the Sliced Inverse Regression methods SIR-I, SIR-II and SIR_α) in order to estimate the direction of the unknown slope parameter β. These methods are computationally simple and fast but depend on the choice of an arbitrary slicing fixed by the user. When the sample size is small, the number and the position of slices have an influence on the estimated direction. In this paper, we suggest to use the corresponding Pooled Slicing methods: PSIR-I (proposed by [3] Aragon, Y. and Saracco, J. 1997. Sliced Inverse Regression (SIR): an appraisal of small sample alternatives to slicing. Computational Statistics, 12: 109–130. [Web of Science ®] , [Google Scholar]), PSIR-II and PSIR_α. These methods combine the results from a number of slicings. We compare the sample behaviour of Slicing and Pooled Slicing methods on simulations. We also propose a practical choice of α in SIR_α and PSIR_α methods.     

Use of the Correlation Coefficient with Normal Probability Plots

The use of the correlation coefficient is suggested as a technique for summarizing and objectively evaluating the information contained in probability plots. Goodness-of-fit tests are constructed using this technique for several commonly used plotting positions for the normal distribution. Empirical sampling methods are used to construct the null distribution for these tests, which are then compared on the basis of power against certain nonnormal alternatives. Commonly used regression tests of fit are also included in the comparisons. The results indicate that use of the plotting position p_i = (i - .375)/(n + .25) yields a competitive regression test of fit for normality.     

Book reviews

Donald J. Koosis: Statistics. A Self-Teaching Guide, Fourth Edition. Wiley 1997, ISBN 0-471-14688-9

Christopher C. Heyde: Quasi-Likelihood and Its Applications, Springer Series in Statistics, 1997, pp. 235, [ISBN 0-387-98225-61

Jeffrey S. Simonoff: Smoothing Methods in Statistics, Springer Series in Statistics, 1996, pp. 338

Andr´e I. Khuri, Thomas Mathew and Birmal K. Sinha: Statistical Tests for Mixed Linear Models, Wiley Series in Probability and Statistics, 1998, pp. 352, PSBN 0-471-1 5653-11

Ronald Christensen: Log-Linear Models and Logistic Regression, Springer Texts in Statistics, 1997, pp. 483, [ISBN 0-387-98247-71

E. L. Lehmann: Testing Statistical Hypotheses, Springer Texts in Statistics, 1997, pp. 600, [ISBN 0-387-94919-41

S. R. Searle, G. Casella and Ch. E. McCulloch: Variance Components, Wiley and Sons, 1992, pp. 496, [ISBN 0-471-62162-51

B. J. T. Morgan: Analysis of Quantal Response Data, Chapman & Hall, 1992, pp. 51 1

Carl D. Huberty: Applied Discriminant Analysis, Wiley and Sons, 1994, pp. 490

C. R. Rao and H. Toutenburg: Linear Models. Least Squares and Alternatives, Springer Series in Statistics, 1995, 188 pp., [ISBN 0-387-94562-81     

Nonparametric comparison of curves with dependent errors

Suppose that data {(x _l,i,n , y _l,i,n ): l?=?1, …, k; i?=?1, …, n} are observed from the regression models: Y _l,i,n ?=?m _l (x _l,i,n )?+?? _l,i,n , l?=?1, …, k, where the regression functions {m _l } _l=1 ^k are unknown and the random errors {? _l,i,n } are dependent, following an MA(∞) structure. A new test is proposed for testing the hypothesis H ₀: m ₁?=?·?·?·?=?m _k , without assuming that {m _l } _l=1 ^k are in a parametric family. The criterion of the test derives from a Crámer-von-Mises-type functional based on different distances between {[mcirc]} _l and {[mcirc]} _s , l?≠?s, l, s?=?1, …, k, where {[mcirc] _l } _l=1 ^k are nonparametric Gasser–Müller estimators of {m _l } _l=1 ^k . A generalization of the test to the case of unequal design points, with different sample sizes {n _l } _l=1 ^k and different design densities {f _l } _l=1 ^k , is also considered. The asymptotic normality of the test statistic is obtained under general conditions. Finally, a simulation study and an analysis with real data show a good behavior of the proposed test.     

A Comparison of Multiple Outlier Detection Methods for Regression Data

The problem of outliers in statistical data has attracted many researchers for a long time. Consequently, numerous outlier detection methods have been proposed in the statistical literature. However, no consensus has emerged as to which method is uniformly better than the others or which one is recommended for use in practical situations. In this article, we perform an extensive comparative Monte Carlo simulation study to assess the performance of the multiple outlier detection methods that are either recently proposed or frequently cited in the outlier detection literature. Our simulation experiments include a wide variety of realistic and challenging regression scenarios. We give recommendations on which method is superior to others under what conditions.     

Monotone B-Spline Smoothing for a Generalized Linear Model Response

Various methods have been proposed for smoothing under the monotonicity constraint. We review the literature and implement an approach of monotone smoothing with B-splines for a generalized linear model response. The approach is expressed as a quadratic programming problem and is easily solved using the statistical software R. In a simulation study, we find that the approach performs better than other approaches with much faster computation time. The approach can also be used for smoothing under other shape constraints or mixed constraints. Supplementary materials of the appendices and R code to implement the developed approach is available online.     

A note on super exogeneity in linear regression models

This note considers how hypotheses of invariance and super exogeneity may be formulated and tested in elliptical linear regression models. It is demonstrated that for jointly elliptical random variables super exogeneity will only hold under normality.