Stepwise local influence in generalized autoregressive conditional heteroskedasticity models

Detection of outliers or influential observations is an important work in statistical modeling, especially for the correlated time series data. In this paper we propose a new procedure to detect patch of influential observations in the generalized autoregressive conditional heteroskedasticity (GARCH) model. Firstly we compare the performance of innovative perturbation scheme, additive perturbation scheme and data perturbation scheme in local influence analysis. We find that the innovative perturbation scheme give better result than other two schemes although this perturbation scheme may suffer from masking effects. Then we use the stepwise local influence method under innovative perturbation scheme to detect patch of influential observations and uncover the masking effects. The simulated studies show that the new technique can successfully detect a patch of influential observations or outliers under innovative perturbation scheme. The analysis based on simulation studies and two real data sets show that the stepwise local influence method under innovative perturbation scheme is efficient for detecting multiple influential observations and dealing with masking effects in the GARCH model.     

Estimation of regression parameters in the presence of outliers in the response

In this paper, we consider a regression model with non-spherical covariance structure and outliers in the response. The generalized least squares estimator obtained from the full data set is generally not used in the presence of outliers and an estimator based on only the non-outlying observations is preferred. Here we propose as an estimator a convex combination of the full set and the deleted set estimators and compare its performance with the other two.     

Effect of individual observations on the Box–Cox transformation

In this paper, we consider the influence of individual observations on inferences about the Box–Cox power transformation parameter from a Bayesian point of view. We compare Bayesian diagnostic measures with the ‘forward’ method of analysis due to Riani and Atkinson. In particular, we look at the effect of omitting observations on the inference by comparing particular choices of transformation using the conditional predictive ordinate and the k _d measure of Pettit and Young. We illustrate the methods using a designed experiment. We show that a group of masked outliers can be detected using these single deletion diagnostics. Also, we show that Bayesian diagnostic measures are simpler to use to investigate the effect of observations on transformations than the forward search method.     

Asymptotic normality of numbers of observations in random regions determined by order statistics

In this paper, we study the joint limiting behaviour of numbers of observations that fall into regions determined by order statistics and Borel sets. We show that suitably centred and normed versions of these numbers are asymptotically multivariate normal under some conditions. We consider two cases: one where the population distribution function is discontinuous and the other where it is continuous and the order statistics are extreme. Finally, we compare results obtained for the two cases with their analogues for absolutely continuous distribution function and central-order statistics.     

Residuals for log-Burr XII regression models in survival analysis

In this paper, we compare three residuals to assess departures from the error assumptions as well as to detect outlying observations in log-Burr XII regression models with censored observations. These residuals can also be used for the log-logistic regression model, which is a special case of the log-Burr XII regression model. For different parameter settings, sample sizes and censoring percentages, various simulation studies are performed and the empirical distribution of each residual is displayed and compared with the standard normal distribution. These studies suggest that the residual analysis usually performed in normal linear regression models can be straightforwardly extended to the modified martingale-type residual in log-Burr XII regression models with censored data.     

Information in an observation in robust designs

We consider the Information contained 1n each observation in a given design robust with respect to the estlmability of parameters and against the unavailability of observations. We compare the observations in various 1-, 2- and 3- dimensional designs on the basis of their informations.     

βExpectation and β-content tolerance intervals for dependent observations

Methods of constructing exact tolerance intervals (β-expectation and β-content) for independent observations are well known. For the case of dependent observations, obtaining exact results is not possible. In this article we provide an approximate method of constructing β-expectation tolerance intervals via a Taylor series expansion. Examples of independent observations are considered to compare the intervals constructed with those obtained by the exact method. For the case of non-stationary type processes we have proposed a method of constructing approximate β-content tolerance intervals. Once again an example is given to illustrate the results.     

Variance estimation in nonparametric regression with jump discontinuities

Variance estimation is an important topic in nonparametric regression. In this paper, we propose a pairwise regression method for estimating the residual variance. Specifically, we regress the squared difference between observations on the squared distance between design points, and then estimate the residual variance as the intercept. Unlike most existing difference-based estimators that require a smooth regression function, our method applies to regression models with jump discontinuities. Our method also applies to the situations where the design points are unequally spaced. Finally, we conduct extensive simulation studies to evaluate the finite-sample performance of the proposed method and compare it with some existing competitors.     

Predictive Accuracy Gain From Disaggregate Sampling in ARIMA Models

We compare the forecast accuracy of autoregressive integrated moving average (ARIMA) models based on data observed with high and low frequency, respectively. We discuss how, for instance, a quarterly model can be used to predict one quarter ahead even if only annual data are available, and we compare the variance of the prediction error in this case with the variance if quarterly observations were indeed available. Results on the expected information gain are presented for a number of ARIMA models including models that describe the seasonally adjusted gross national product (GNP) series in the Netherlands. Disaggregation from annual to quarterly GNP data has reduced the variance of short-run forecast errors considerably, but further disaggregation from quarterly to monthly data is found to hardly improve the accuracy of monthly forecasts.     

A Comparison of Three-Level Orthogonal Arrays in the Presence of Different Correlation Structures in Observations

When the experimenter suspects that there might be a quadratic relation between the response variable and the explanatory parameters, a design with at least three points must be employed to establish and explore this relation (second-order design). Orthogonal arrays (OAs) with three levels are often used as second-order response surface designs. Generally, we assume that the data are independent observations; however, there are many situations where this assumption may not be sustainable. In this paper, we want to compare three-level OAs with 18, 27, and 36 runs under the presence of three specific forms of correlation in observations. The aim is to derive the best designs that can be efficiently used for response surface modeling.     

Comparison of record data and random observations based on statistical evidence

“Science looks to Statistics for an objective measure of the strength of evidence in a given body of observations. The Law of Likelihood explains that the strength of statistical evidence for one hypothesis over another is measured by their likelihood ratio” (Blume, 2002). In this paper, we compare probabilities of weak and strong misleading evidence based on record data with those based on the same number of iid observations from the original distribution. We shall also use a criterion defined as a combination of probabilities of weak and strong misleading evidence to do the above comparison. We also give numerical results of a simulated comparison.     

Local Influence in Generalized Estimating Equations

Abstract.  We investigate the influence of subjects or observations on regression coefficients of generalized estimating equations (GEEs) using local influence. The GEE approach does not require the full multivariate distribution of the response vector. We extend the likelihood displacement to a quasi-likelihood displacement, and propose local influence diagnostics under several perturbation schemes. An illustrative example in GEEs is given and we compare the results using the local influence and deletion methods.     

Model selection of generalized partially linear models with missing covariates

In this paper, a generalized partially linear model (GPLM) with missing covariates is studied and a Monte Carlo EM (MCEM) algorithm with penalized-spline (P-spline) technique is developed to estimate the regression coefficients and nonparametric function, respectively. As classical model selection procedures such as Akaike's information criterion become invalid for our considered models with incomplete data, some new model selection criterions for GPLMs with missing covariates are proposed under two different missingness mechanism, say, missing at random (MAR) and missing not at random (MNAR). The most attractive point of our method is that it is rather general and can be extended to various situations with missing observations based on EM algorithm, especially when no missing data involved, our new model selection criterions are reduced to classical AIC. Therefore, we can not only compare models with missing observations under MAR/MNAR settings, but also can compare missing data models with complete-data models simultaneously. Theoretical properties of the proposed estimator, including consistency of the model selection criterions are investigated. A simulation study and a real example are used to illustrate the proposed methodology.     

Influence functions of the Spearman and Kendall correlation measures

Nonparametric correlation estimators as the Kendall and Spearman correlation are widely used in the applied sciences. They are often said to be robust, in the sense of being resistant to outlying observations. In this paper we formally study their robustness by means of their influence functions and gross-error sensitivities. Since robustness of an estimator often comes at the price of an increased variance, we also compute statistical efficiencies at the normal model. We conclude that both the Spearman and Kendall correlation estimators combine a bounded and smooth influence function with a high efficiency. In a simulation experiment we compare these nonparametric estimators with correlations based on a robust covariance matrix estimator.     

Confidence Intervals Following Group Sequential Tests in Clinical Trails with Multivariate Observations

Several authors ( e.g. Kim and DeMets, 1987a, Biometrics) have developed methods for estimation following group sequential tests in clinical trials when each patient has only one response. In many long-term clinical trials, the subjects enter the study sequentially and yield repeated measurements or other types of multivariate observations at successive follow-up visits. Typically, investigators want to compare a parameter of interest such as the slope over time in a repeated measures trial etc. In this article, we propose an exact confidence interval for these parameters in a repeated measures trial, and compare it with a naive confidence interval using Monte Carlo simulation. This method is illustrated with a real example for bone density measurements.     

Group sequential x2 statistic for testing the difference of two mean vectors in clinical trials

In this study we discuss the group sequential procedures for comparing two treatments based on multivariate observations in clinical trials. Also we suppose that a response vector on each of two treatments has a multivariate normal distribution with unknown covariance matrix. Then we propose a group sequential x² statistic in order to carry out repeated significance test for hypothesis of no difference between two population mean vectors. In order to realize the group sequential test where average sample number is reduced, we propose another modified group sequential x² statistic by extension of Jennison and Turnbull ( 1991 ). After construction of repeated confidence boundaries for making the repeated significance test, we compare two group sequential procedures based on two statistics regarding the average sample number and the power of the test in the simulations.     

Sequential Hypothesis Tests for Identifying the Minimum Dose with a Threshold Effect

We propose a procedure to identify a lowest dose having greater effect than a threshold dose under the assumption of monotonicity of dose mean response in dose response test. So, we use statistics based on contrasts among sample means and apply a group sequential procedure to our procedure to identify effectively the dose. If we can identify the dose at an early step in the sequential test, since we can terminate the procedure with a few observations, the procedure is useful from an economical point of view. In a simulation studies, we compare the superiority among these procedures based on three contrasts.     

Approximate maximum likelihood predictors of future failure times of shifted exponential distributions under multiple type II censoring

Suppose we consider a general multiple type II censored sample (some middle observations being censored) from a shifted exponential distribution. The maximum likelihood prediction method does not admit explicit solutions. We introduce a simple approximation to one of prediction likelihood equations and derive approximate predictors of missing failure times. We compute their mean square prediction errors by simulation and compare them with the best linear predictors. Further, we present two real examples to illustrate this method of prediction.AMS Subject Classification (2000): 62G30, 62M20, 62F99     

Accounting for plasma levels below detection limits in a one-compartment zero-order absorption pharmacokinetics model

We provide a simple method for fitting a one-compartment, zero-order absorption pharmacokinetics model in the presence of observations below the detection limit. This method may be extended to more complex pharmacokinetics models. We demonstrate, using a small simulation study, that the method provides accurate parameter estimates over a range of detection limits and we compare it to an ad hoc midpoint method. An applied example is provided from a pharmacokinetic investigation of a nicotine nasal spray.     

A discrete probabilistic model for analyzing pairwise comparison matrices

The pairwise comparison matrix is often used for the estimation of the priorities in the analytic hierarchy process. In this paper, we propose an estimation method based on the discrete probabilistic expression of each choice. Moreover, we show numerical examples to compare our method with commonly used ones. As a result, it is shown that, using a robust divergence measure for the estimation, the proposed method can extract the priorities more stably even if some outlying observations are included.