MODELING THE LIFETIME OF LONGITUDINAL ELEMENTS

In the study of the stochastic behaviour of the lifetime of an element as a function of its length, it is often observed that the failure time (or lifetime) decreases as the length increases. In probabilistic terms, such an idea can be expressed as follows. Let T be the lifetime of a specimen of length x, so the survival function, which denotes the probability that an element of length x survives till time t, will be given by S_T (t, x) = P(T > t/α(x), where α(x) is a monotonically decreasing function. In particular, it is often assumed that T has a Weibull distribution. In this paper, we propose a generalization of this Weibull model by assuming that the distribution of T is Generalized gamma (GG). Since the GG model contains the Weibull, Gamma and Lognormal models as special and limiting cases, a GG regression model is an appropriate tool for describing the size effect on the lifetime and for selecting among the embedded models. Maximum likelihood estimates are obtained for the GG regression model with α(x) = cx^b . As a special case this provide an alternative to the usual approach to estimation for the GG distribution which involves reparametrization. Related parametric inference issues are addressed and illustrated using two experimental data sets. Some discussion of censored data is also provided.     

Performance of Confidence Intervals in Regression Models with Unbalanced One-Fold Nested Error Structures

Abstract

In this article we consider the problem of constructing confidence intervals for a linear regression model with unbalanced nested error structure. A popular approach is the likelihood-based method employed by PROC MIXED of SAS. In this article, we examine the ability of MIXED to produce confidence intervals that maintain the stated confidence coefficient. Our results suggest that intervals for the regression coefficients work well, but intervals for the variance component associated with the primary level cannot be recommended. Accordingly, we propose alternative methods for constructing confidence intervals on the primary level variance component. Computer simulation is used to compare the proposed methods. A numerical example and SAS code are provided to demonstrate the methods.     

Efficiency and Validity Analyses of Two-Stage Estimation Procedures and Derived Testing Procedures in Quantitative Linear Models with AR(1) Errors

Abstract

In a quantitative linear model with errors following a stationary Gaussian, first-order autoregressive or AR(1) process, Generalized Least Squares (GLS) on raw data and Ordinary Least Squares (OLS) on prewhitened data are efficient methods of estimation of the slope parameters when the autocorrelation parameter of the error AR(1) process, ρ, is known. In practice, ρ is generally unknown. In the so-called two-stage estimation procedures, ρ is then estimated first before using the estimate of ρ to transform the data and estimate the slope parameters by OLS on the transformed data. Different estimators of ρ have been considered in previous studies. In this article, we study nine two-stage estimation procedures for their efficiency in estimating the slope parameters. Six of them (i.e., three noniterative, three iterative) are based on three estimators of ρ that have been considered previously. Two more (i.e., one noniterative, one iterative) are based on a new estimator of ρ that we propose: it is provided by the sample autocorrelation coefficient of the OLS residuals at lag 1, denoted r(1). Lastly, REstricted Maximum Likelihood (REML) represents a different type of two-stage estimation procedure whose efficiency has not been compared to the others yet. We also study the validity of the testing procedures derived from GLS and the nine two-stage estimation procedures. Efficiency and validity are analyzed in a Monte Carlo study. Three types of explanatory variable x in a simple quantitative linear model with AR(1) errors are considered in the time domain: Case 1, x is fixed; Case 2, x is purely random; and Case 3, x follows an AR(1) process with the same autocorrelation parameter value as the error AR(1) process. In a preliminary step, the number of inadmissible estimates and the efficiency of the different estimators of ρ are compared empirically, whereas their approximate expected value in finite samples and their asymptotic variance are derived theoretically. Thereafter, the efficiency of the estimation procedures and the validity of the derived testing procedures are discussed in terms of the sample size and the magnitude and sign of ρ. The noniterative two-stage estimation procedure based on the new estimator of ρ is shown to be more efficient for moderate values of ρ at small sample sizes. With the exception of small sample sizes, REML and its derived F-test perform the best overall. The asymptotic equivalence of two-stage estimation procedures, besides REML, is observed empirically. Differences related to the nature, fixed or random (uncorrelated or autocorrelated), of the explanatory variable are also discussed.     

LOCAL INFLUENCE IN NONLINEAR REPRODUCTIVE DISPERSION MODELS

Nonlinear reproductive dispersion models (NRDM, Jorgensen 1997) include a wider range of distributions and nonlinear models such as the possibility of correlated errors and nonlinear hypotheses dropping the exponential family assumption. Based on the generalized Cook distance and the conformal normal curvature of Poon & Poon (1999), local influence of minor perturbations on the data set is investigated for NRDM. Two examples are used to illustrate our results.     

GENERALIZED SIGNED-RANK ESTIMATORS FOR AUTOREGRESSION PARAMETERS

ABSTRACT

A new class of weighted signed-rank-based estimates for estimating the parameter vector of an autoregressive time series is considered. The Wilcoxon signed-rank estimate and the GR-estimates of Terpstra et al. (GR-Estimates for an Autoregressive Time Series. Statistics and Probability Letters 2001, 51, 165–172; Generalized Rank Estimates for an Autoregressive Time Series: A U-Statistic Approach. Statistical Inference for Stochastic Processes 2001, 4, 155–179) can be viewed as special cases of the so-called GSR-estimates. Asymptotic linearity properties are derived for the GSR-estimates. Based on these properties, and a symmetry assumption, the GSR-estimates are shown to be asymptotically normal at rate n ^1/2. The theory of U-Statistics along with a characterization of weak dependence that is inherent in stationary AR(p) models are the primary tools used to obtain the results. Tests of hypotheses as well as standard errors for confidence interval procedures can be based on such results. An efficiency study indicates that, for an appropriately chosen set of weights, the GSR-estimate is more efficient than the GR-estimate. Furthermore, the GSR-estimate has an added advantage in that an intercept term can be estimated simultaneously; unlike the GR-estimate. Two examples and a small simulation study are used to illustrate the computational and robust aspects of the GSR-estimates.     

RELIABILITY PROPERTIES OF SERIES AND PARALLEL SYSTEMS FROM BIVARIATE EXPONENTIAL MODELS

ABSTRACT

In this paper we study the classification of the generalized mixtures of two or three exponential distributions, in the ILR and DLR classes, and consequently in the IFR and DFR classes. We apply these results to classify the aging of the series and parallel systems, in accord with some common bivariate exponential models of their components.     

Out-of-Sample Performance of Discrete-Time Spot Interest Rate Models

We provide a comprehensive analysis of the out-of-sample performance of a wide variety of spot rate models in forecasting the probability density of future interest rates. Although the most parsimonious models perform best in forecasting the conditional mean of many financial time series, we find that the spot rate models that incorporate conditional heteroscedasticity and excess kurtosis or heavy tails have better density forecasts. Generalized autoregressive conditional heteroscedasticity significantly improves the modeling of the conditional variance and kurtosis, whereas regime switching and jumps improve the modeling of the marginal density of interest rates. Our analysis shows that the sophisticated spot rate models in the existing literature are important for applications involving density forecasts of interest rates.     

Quadratic estimators of covariance components in a multivariate mixed linear model

It is known that the Henderson Method III (Biometrics 9:226–252, 1953) is of special interest for the mixed linear models because the estimators of the variance components are unaffected by the parameters of the fixed factor (or factors). This article deals with generalizations and minor extensions of the results obtained for the univariate linear models. A MANOVA mixed model is presented in a convenient form and the covariance components estimators are given on finite dimensional linear spaces. The results use both the usual parametric representations and the coordinate-free approach of Kruskal (Ann Math Statist 39:70–75, 1968) and Eaton (Ann Math Statist 41:528–538, 1970). The normal equations are generalized and it is given a necessary and sufficient condition for the existence of quadratic unbiased estimators for covariance components in the considered model.     

Generalized inferences on the common mean vector of several multivariate normal populations

The hypothesis testing and confidence region are considered for the common mean vector of several multivariate normal populations when the covariance matrices are unknown and possibly unequal. A generalized confidence region is derived using the concepts of generalized method based on the generalized p$p$-value. The generalized confidence region is illustrated with two numerical examples. The merits of the proposed method are numerically compared with those of existing methods with respect to their expected area or expected d-dimensional volumes and coverage probabilities under different scenarios.