Change-Point Estimation in Long Memory Nonparametric Models with Applications

We consider the estimation of a change point or discontinuity in a regression function for random design model with long memory errors. We provide several change-point estimators and investigate the consistency of the estimators. Using the fractional ARIMA process as an example of long memory process, we report a small Monte Carlo experiment to compare the performance of the estimators in finite samples. We finish by applying the method to a climatological data example.     

On the asymptotic properties of parameter estimates in a regression model with non-normally distributed errors

As pointed out in a recent paper by Amirkhalkhali and Rao (1986) (henceforth referred to as A&R), the usual assumption of normality for the error terms of a regression model isoften untenable. However, when this assumption is dropped, it may be difficult to characterize parameter estimates for the model. For example, A&R (p. 189) state that “if the regression errors are non-normal, we are not even sure of their [e.g., the generalized least squares parameter estimates1] asymptotic properties.” A partial answer, however, is given by Spall and Wall (1984), which presents an asymptotic distribution theory for Kalman filter estimates for cases where the random terms of the state space model are not necessarily Gaussian. Certain of these asymptotic distribution results are also discussed in Spall (1985) in the context of model validation (diagnostic checking)     

Optimal allocation of time points for the random effects model

In this article the problem of the optimal selection and allocation of time points in repeated measures experiments is considered. D‐ optimal designs for linear regression models with a random intercept and first order auto‐regressive serial correlations are computed numerically and compared with designs having equally spaced time points. When the order of the polynomial is known and the serial correlations are not too small, the comparison shows that for any fixed number of repeated measures, a design with equally spaced time points is almost as efficient as the D‐ optimal design. When, however, there is no prior knowledge about the order of the underlying polynomial, the best choice in terms of efficiency is a D‐ optimal design for the highest possible relevant order of the polynomial. A design with equally‐spaced time points is the second best choice     

ML and Bayes estimation in a two-phase tandem queue with a second optional service and random feedback

ABSTRACT

In this article, we consider a two-phase tandem queueing model with a second optional service and random feedback. The first phase of service is essential for all customers and after the completion of the first phase of service, any customer receives the second phase of service with probability α, feedback to the tail of the first queue with probability β if the service is not successful and leaves the system with probability 1 ? α ? β. In this model, our main purpose is to estimate the parameters of the model, traffic intensity, and mean system size, in the steady state, via maximum likelihood and Bayesian methods. Furthermore, we find asymptotic confidence intervals for mean system size. Finally, by a simulation study, we compute the confidence levels and mean length for asymptotic confidence intervals of mean system size with a nominal level 0.95.     

A Weighting Approach for GEE Analysis with Missing Data

We propose a new weighting (WT) method to handle missing categorical outcomes in longitudinal data analysis using generalized estimating equations (GEE). The proposed WT provides a valid GEE estimator when the data are missing at random (MAR), and has more stable weights and shows advantage in efficiency compared to the inverse probability weighing method in the presence of small observation probabilities. The WT estimator is similar to the stabilized weighting (SWT) estimator under mild conditions, but it is more stable and efficient than SWT when the associations of the outcome with the observation probabilities and the covariate are strong.     

Akaike Information Criterion for Selecting Variables in the Nested Error Regression Model

The Akaike Information Criterion (AIC) is developed for selecting the variables of the nested error regression model where an unobservable random effect is present. Using the idea of decomposing the likelihood into two parts of “within” and “between” analysis of variance, we derive the AIC when the number of groups is large and the ratio of the variances of the random effects and the random errors is an unknown parameter. The proposed AIC is compared, using simulation, with Mallows' C _p , Akaike's AIC, and Sugiura's exact AIC. Based on the rates of selecting the true model, it is shown that the proposed AIC performs better.     

Convergence Rate of Wavelet Expansions of Gaussian Random Processes

This article characterizes uniform convergence rate for general classes of wavelet expansions of stationary Gaussian random processes. The convergence in probability is considered.     

A General Model of Random Variation

A statistical distribution of a random variable is uniquely represented by its normal-based quantile function. For a symmetrical distribution it is S-shaped (for negative kurtosis) and inverted S-shaped (otherwise). As skewness departs from zero, the quantile function gradually transforms into a monotone convex function (positive skewness) or concave function (otherwise). Recently, a new general modeling platform has been introduced, response modeling methodology, which delivers good representation to monotone convex relationships due to its unique “continuous monotone convexity” property. In this article, this property is exploited to model the normal-based quantile function, and explored using a set of 27 distributions.     

A Multivariate Adaptive Exponentially Weighted Moving Average Control Chart

A multivariate extension of the adaptive exponentially weighted moving average (AEWMA) control chart is proposed. The new multivariate scheme can detect small and large shifts in the process mean vector effectively. The proposed scheme can be viewed as a smooth combination of a multivariate exponentially weighted moving average (MEWMA) chart and a Shewhart χ²-chart. The optimal design of the proposed chart is given according to a pre-specified in-control average run length and two shift sizes; a small and large shift each measured in terms of the non centrality parameter. The signal resistance of the newly proposed multivariate chart is also given. Comparisons among the new chart, the MEWMA chart, and the combined Shewhart-MEWMA (S-MEWMA) chart in terms of the standard and worst-case average run length profiles are presented. In addition, the three charts are compared with respect to their worst-case signal resistance values. The proposed chart gives somewhat better worst-case ARL and signal resistance values than the competing charts. It also gives better standard ARL performance especially for moderate and large shifts. The effectiveness of our proposed chart is illustrated through an example with simulated data set.