On Rates of Convergence for Posterior Distributions Under Misspecification

We extend the approach of Walker (2003 Walker , S. ( 2003 ). On sufficient conditions for Bayesian consistency . Biometrika 90 ( 2 ): 482 – 488 .[Crossref], [Web of Science ®] , [Google Scholar]); (2004 Walker , S. ( 2004 ). New approaches to Bayesian consistency . Ann. Statist. 32 ( 5 ): 2028 – 2043 .[Crossref], [Web of Science ®] , [Google Scholar]) to the case of misspecified models. A sufficient condition for establishing rates of convergence is given based on a key identity involving martingales, which does not require construction of tests. We also show roughly that the result obtained by using tests can also be obtained by our approach, which demonstrates the potential wider applicability of this method.     

The Modified Pulak and Al-Sultan's Model for Determining the Optimum Process Parameters

Economic selection of process parameters has been an important topic in modern statistical process control. The optimum process parameters setting have a major effect on the expected profit/cost per item. There are some concerns on the problem of setting process parameters. Boucher and Jafari (1991 Boucher , T. O. , Jafari , M. A. ( 1991 ). The optimum target value for single filling operations with quality sampling plans . J. Qual. Technol. 23 : 44 – 47 . [CSA] [Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) first considered the attribute single sampling plan applied in the selection of process target. Pulak and Al-Sultan (1996 Pulak , M. F. S. , Al-Sultan , K. S. ( 1996 ). The optimum targeting for a single filling operation with rectifying inspection . Omega 24 : 727 – 733 . [CSA] [CROSSREF] [Crossref], [Web of Science ®] , [Google Scholar]) extended Boucher and Jafari's model and presented the rectifying inspection plan for determining the optimum process mean. In this article, we further propose a modified Pulak and Al-Sultan model for determining the optimum process mean and standard deviation under the rectifying inspection plan with the average outgoing quality limit (AOQL) protection. Taguchi's (1986 Taguchi , G. ( 1986 ). Introduction to Quality Engineering . Asian Productivity Organization . [Google Scholar]) symmetric quadratic quality loss function is adopted for evaluating the product quality. By solving the modified model, we can obtain the optimum process parameters with the maximum expected profit per item and the specified quality level can be reached.     

From Ordinary to Shrinkage Square-Root Estimators

Consider a skewed population. Suppose an intelligent guess could be made about an interval that contains the population mean. There may exist biased estimators with smaller mean squared error than the arithmetic mean within such an interval. This article indicates when it is advisable to shrink the arithmetic mean towards a guessed interval using root estimators. The goal is to obtain an estimator that is better near the average of natural origins. An estimator proposed. This estimator contains the Thompson (1968 Thompson , J. R. ( 1968 ). Accuracy borrowing in the estimation of the mean by shrinkage towards an interval . J. Amer. Statist. Assoc. 63 : 953 – 963 . [CSA] [CROSSREF] [Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) ordinary shrinkage estimator, the Jenkins et al. (1973 Jenkins , O. C. , Ringer , L. J. , Hartley , H. O. ( 1973 ). Root estimators . J Amer. Statist. Assoc. 68 : 414 – 419 . [CSA] [CROSSREF] [Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) square-root estimator, and the arithmetic sample mean as special cases. The bias and the mean squared error of the proposed more general estimator is compared with the three special cases. Shrinkage coefficients that yield minimum mean squared error estimators are obtained. The proposed estimator is considerably more efficient than the three special cases. This remains true for highly skewed populations. The merits of the proposed shrinkage square-root estimator are supported by the results of numerical and simulation studies.     

Lag Length Selection in DF-GLS Unit Root Tests

We reinvestigate the empirical problem of lag length selection in unit root tests when using the augmented Dickey–Fuller test based on GLS-detrending. We extend the Ng and Perron (1995 Ng , S. , Perron , P. ( 1995 ). Unit root tests in ARMA models with data-dependent methods for the selection of the truncation lag . Journal of American Statistical Association 90 : 268 – 281 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) work on this issue by applying the finite sample critical values calculated using the formulae proposed by Cheung and Lai (1995 Cheung , Y. W. , Lai , K. S. ( 1995 ). Lag order and critical values of a modified Dickey–Fuller test . Oxford Bulletin of Business and Economics 57 : 411 – 418 .[Crossref] , [Google Scholar]). Unlike Ng and Perron (2001 Ng , S. , Perron , P. (2001). Lag length selection and the construction of unit root tests with good size and power. Econometrica 69:1519–1554.[Crossref], [Web of Science ®] , [Google Scholar]) we find through simulation studies that the method of selecting lag length using the sequential t-test in the ADF regression of GLS-detrended series performs the best in most cases.     

Online Control Charts for Process Averages Based on Repeated Median Filters

Two types of estimates of process level, namely repeated median estimates (Siegel, 1982 Siegel , A. F. ( 1982 ). Robust regression using repeated medians . Biometrika 69 : 242 – 244 .[Crossref], [Web of Science ®] , [Google Scholar]) and full online estimates (Gather et al., 2006 Gather , U. , Schettlinger , K. , Fried , R. ( 2006 ). Online signal extraction by robust linear regression . Computational Statistics 21 : 33 – 51 .[Crossref], [Web of Science ®] , [Google Scholar]) based on repeated median filters, are used to develop control charts. The distributional properties of the estimates are studied using simulation and these are found to closely follow normal distribution. The repeated median being robust against outliers with asymptotically 50% breakdown value and having small standard deviation is found to be useful as a basis for monitoring process averages. The control charts using repeated median estimates have been recommended for general use.     

Regression Analysis with Covariates Missing at Random: A Piece-wise Nonparametric Model for Missing Covariates

Statistical analysis for the regression model f _β(y | x, z) with missing values in the covariate vector X requires modeling of the covariate distribution g(x | z). Likelihood methods, including Ibrahim (1990 Ibrahim , J. G. ( 1990 ). Incomplete data in generalized linear models . J. Amer. Statist. Assoc. 85 : 765 – 769 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), Chen (2004 Chen , H. Y. (2004). Nonparametric and semiparametric models for missing covariates in parametric regression. J. Amer. Statist. Assoc. 99:1176–1189.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), and Zhao (2005 Zhao , Y. ( 2005 ). Design and Efficient Estimation in Regression Analysis with Missing Data in Two-Phase Studies. Ph.D. thesis , University of Waterloo . [Google Scholar]), need either X or Z to be discrete. This article considers extending the likelihood methods to deal with cases where both X and Z may be continuous. We propose modeling the covariate distribution g(x | z) using a piece-wise nonparametric model, then a maximum likelihood estimate (MLE) of β can be computed following the maximum likelihood estimating procedure of Chen (2004 Chen , H. Y. (2004). Nonparametric and semiparametric models for missing covariates in parametric regression. J. Amer. Statist. Assoc. 99:1176–1189.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) or Zhao (2005 Zhao , Y. ( 2005 ). Design and Efficient Estimation in Regression Analysis with Missing Data in Two-Phase Studies. Ph.D. thesis , University of Waterloo . [Google Scholar]). The resulting estimation method is easy to implement and the asymptotic properties of the MLE follow under certain conditions. Extensive simulation studies for different models indicate that the proposed method is acceptable for practical implementation. A real data example is used to illustrate the method.     

Two Statistics for Unconditional Exact Tests in Matched-Pairs Design

Abstract

In a recent article Hsueh et al. (Hsueh, H.-M., Liu, J.-P., Chen, J. J. (2001 Hsueh, H.-M., Liu, J.-P. and Chen, J. J. 2001. Unconditional exact tests for equivalence or noninferiority for paired binary endpoints. Biometrics, 57: 478–483. [Crossref], [PubMed], [Web of Science ®] , [Google Scholar]). Unconditional exact tests for equivalence or noninferiority for paired binary endpoints. Biometrics 57:478–483.) considered unconditional exact tests for paired binary endpoints. They suggested two statistics one of which is based on the restricted maximum-likelihood estimator. Properties of these statistics and the related tests are treated in this article.     

Parametric Nonlinear Regression with Endogenous Switching

Based on the insightful work of Olsen (1980 Olsen , R. J. ( 1980 ). A least squares correction for selectivity bias . Econometrica 48 : 1815 – 1820 .[Crossref], [Web of Science ®] , [Google Scholar]) for the linear context, a generic and unifying framework is developed that affords a simple extension of the classical method of Heckman (1974 Heckman , J. ( 1974 ). Shadow prices, market wages, and labor supply . Econometrica 42 : 679 – 694 .[Crossref], [Web of Science ®] , [Google Scholar], 1976 Heckman , J. ( 1976 ). The common structure of statistical models of truncation sample selection and limited dependent variables and a simple estimator for such models . Annals of Economic and Social Measurement 5 : 475 – 492 . [Google Scholar], 1978 Heckman , J. ( 1978 ). Dummy endogenous variables in a simultaneous equation system . Econometrica 46 : 931 – 959 .[Crossref], [Web of Science ®] , [Google Scholar], 1979 Heckman , J. ( 1979 ). Sample selection bias as a specification error . Econometrica 47 : 153 – 161 .[Crossref], [Web of Science ®] , [Google Scholar]) to a broad class of nonlinear regression models involving endogenous switching and its two most common incarnations, endogenous sample selection and endogenous treatment effects. The approach should be appealing to applied researchers for three reasons. First, econometric applications involving endogenous switching abound. Secondly, the approach requires neither linearity of the regression function nor full parametric specification of the model. It can, in fact, be applied under the minimal parametric assumptions—i.e., specification of only the conditional means of the outcome and switching variables. Finally, it is amenable to relatively straightforward estimation methods. Examples of applications of the method are discussed.     

GMM in linear regression for longitudinal data with multiple covariates measured with error

Griliches and Hausman 5 Griliches, Z. and Hausman, J. A. 1986. Errors in variables in panel data. J. Econometrics, 32: 93–118. [Crossref], [Web of Science ®] [Google Scholar] and Wansbeek 11 Wansbeek, T. J. 2001. GMM estimation in panel data models with measurement error. J. Econometrics, 104: 259–268. [Crossref], [Web of Science ®] [Google Scholar] proposed using the generalized method of moments (GMM) to obtain consistent estimators in linear regression models for longitudinal data with measurement error in one covariate, without requiring additional validation or replicate data. For usefulness of this methodology, we must extend it to the more realistic situation where more than one covariate are measured with error. Such an extension is not straightforward, since measurement errors across different covariates may be correlated. By a careful construction of the measurement error correlation structure, we are able to extend Wansbeek's GMM and show that the extended Griliches and Hausman's GMM is equivalent to the extended Wansbeek's GMM. For illustration, we apply the extended GMM to data from two medical studies, and compare it with the naive method and the method assuming only one covariate having measurement error.     

Refined Inference on Long Memory in Realized Volatility

There is an emerging consensus in empirical finance that realized volatility series typically display long range dependence with a memory parameter (d) around 0.4 (Andersen et al., 2001 Andersen , T. G. , Bollerslev , T. , Diebold , F. X. , Labys , P. ( 2001 ). The distribution of realized exchange rate volatility . Journal of the American Statistical Association 96 ( 453 ): 42 – 55 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]; Martens et al., 2004 Martnes , M. , Van Dijk , D. , De Pooter , M. ( 2004 ). Modeling and forecasting S&P 500 volatility: Long memory, structural breaks and nonlinearity. Tinbergen Institute Discussion Paper 2004-067/4 . [Google Scholar]). The present article provides some illustrative analysis of how long memory may arise from the accumulative process underlying realized volatility. The article also uses results in Lieberman and Phillips (2004 Lieberman , O. , Phillips , P. C. B. ( 2004 ). Expansions for the distribution of the maximum likelihood estimator of the fractional difference parameter . Econometric Theory 20 ( 3 ): 464 – 484 . [Google Scholar], 2005 Lieberman , O. , Phillips , P. C. B. ( 2005 ). Expansions for approximate maximum likelihood estimators of the fractional difference parameter . The Econometrics Journal 8 : 367 – 379 . [Google Scholar]) to refine statistical inference about d by higher order theory. Standard asymptotic theory has an O(n ^?1/2) error rate for error rejection probabilities, and the theory used here refines the approximation to an error rate of o(n ^?1/2). The new formula is independent of unknown parameters, is simple to calculate and user-friendly. The method is applied to test whether the reported long memory parameter estimates of Andersen et al. (2001 Andersen , T. G. , Bollerslev , T. , Diebold , F. X. , Labys , P. ( 2001 ). The distribution of realized exchange rate volatility . Journal of the American Statistical Association 96 ( 453 ): 42 – 55 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) and Martens et al. (2004 Martnes , M. , Van Dijk , D. , De Pooter , M. ( 2004 ). Modeling and forecasting S&P 500 volatility: Long memory, structural breaks and nonlinearity. Tinbergen Institute Discussion Paper 2004-067/4 . [Google Scholar]) differ significantly from the lower boundary (d = 0.5) of nonstationary long memory, and generally confirms earlier findings.     

An Inverse-Probability-Weighted Approach to the Estimation of Successive Duration Times under Left-Truncation and Right-Censoring

In incident cohort studies, survival data often include subjects who have had an initiate event at recruitment and may potentially experience two successive events (first and second) during the follow-up period. Since the second duration process becomes observable only if the first event has occurred, left-truncation and dependent censoring arise if the two duration times are correlated. To confront the two potential sampling biases, Chang and Tzeng (2006 Chang , S.-H. , Tzeng , S.-J. (2006). Noparametric estimation of sojourn time distributions for truncated serial event data- a weight-adjusted approach. Lifetime Data Anal. 53–67.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) provided an inverse-probability-weighted (IPW) approach for estimating the joint probability function of successive duration times. In this note, an alternative IPW approach is proposed. A simulation study is conducted to compare the two IPW approaches.     

A Similarity Measure Based on Species Proportions

ABSTRACT

There are several indices for measuring the similarity of two populations, including the ratio of the number of shared species to the number of distinct species (Jaccard's index) and the conditional probability of observing a shared species (Smith et al., 1996 Smith , W. , Solow , A. R. , Preston , P. E. ( 1996 ). An estimator of species overlap using a modified beta-binomial model. Biometrics 52 : 1472 – 1477 . [CSA] [CROSSREF] [Crossref], [Web of Science ®] , [Google Scholar]). However, these indices only take into account the number of species and species proportions of shared species. In this article, we propose a new similarity index which includes the species proportions of both the shared and non shared species in each population, and also propose a Nonparametric Maximum Likelihood Estimator (NPMLE) for this index. Bootstrap and delta methods are used to evaluate the standard errors of the NPMLE. Based on a loss function, we also compare a class of nonparametric estimators for the proposed index in various situations.     

Adaptive Test for Periodicity in Self-Exciting Threshold Autoregressive Models

This article is devoted to the study of the periodicity testing problem in a self-exciting threshold autoregressive (SETAR) model. The local asymptotic normality (LAN) property is shown via the adapted sufficient conditions due to Swensen (1985 Swensen , A. R. ( 1985 ). The asymptotic distribution of the likelihood ratio for autoregressive time series with a regression trend . Journal of Multivariate Analysis 16 : 54 – 70 .[Crossref], [Web of Science ®] , [Google Scholar]). Moreover, the LAN of the central sequence is established. First, we consider the case where the innovation density is specified and we obtain a parametric local asymptotic test. Second, we construct an adaptive test in the case where this density is unspecified but symmetric. The performances of these established tests are shown via simulation studies.     

Statistical Properties of Threshold Models

This article focuses the attention on the Self Exciting Threshold Autoregressive Moving Average model (SETARMA) proposed in Tong (1983 Tong , H. ( 1983 ). Threshold Models in Nonlinear Time Series Analysis . London : Springer-Verlag .[Crossref] , [Google Scholar]). The stochastic structure of the model is discussed and different specifications are presented. Starting from one of them, we give sufficient conditions for the weak stationarity of the model that are discussed and critically compared to other results given in literature. In particular, after showing that the SETARMA model belongs to the class of the Random Coefficients Autoregressive models, widely discussed in Nicholls and Quinn (1982 Nicholls , D. F. , Quinn , B. G. (1982). Random Coefficients Autoregressive Models. An Introduction . New York : Springer-Verlag.[Crossref] , [Google Scholar]), we give some issues on the weak stationarity of its stochastic structure that are more general than those given in the existing literature and appear not affected by the moving average component.     

A Note on the Poisson Likelihood Ratio Test Statistic for Kulldorff's Scan Methods

Methods based on scan statistics are widely used in health-related applications to detect clusters of disease. The most common methods are based on the Bernoulli and Poisson models. Kulldorff (1997 Kulldorff , M. ( 1997 ). A spatial scan statistic . Communications in Statistics—Theory and Methods 26 : 1481 – 1496 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) derived the likelihood ratio test statistic for his scan method for both of these models. His scan statistic is widely used with freely available software, SaTScan^? (see Kulldorff, 2005 Kulldorff , M. ( 2005 ). SaTScan: software for the spatial, temporal, and space-time scan statistics , version 5.1 [computer program]. Information Management Services 2005; Available: http://www.satscan.org/. Accessed July 2007 . [Google Scholar]). We provide an alternative derivation of the likelihood ratio test statistic in the Poisson case. Our derivation is simpler and more general in the sense that it applies when the incidences are not aggregated into subregional counts.     

A Distribution-Free Control Chart Based on Order Statistics

In this article, we introduce a new distribution-free Shewhart-type control chart that takes into account the location of a single order statistic of the test sample (such as the median) as well as the number of observations in that test sample that lie between the control limits. Exact formulae for the alarm rate, the run length distribution, and the average run length (ARL) are all derived. A key advantage of the chart is that, due to its nonparametric nature, the false alarm rate and in-control run length distribution are the same for all continuous process distributions, and so will be naturally robust. Tables are provided for the implementation of the chart for some typical ARL values and false alarm rates. The empirical study carried out reveals that the new chart is preferable from a robustness point of view in comparison to a classical Shewhart-type chart and also the nonparametric chart of Chakraborti et al. (2004 Chakraborti , S. , van der Laan , P. , van de Wiel , M. A. ( 2004 ). A class of distribution-free control charts . J. Roy. Statist. Soc. Ser. C-Appl. Statist. 53 ( 3 ): 443 – 462 .[Web of Science ®] , [Google Scholar]).     

Goodness-of-Fit Test and Parameter Estimation for a Proportional Odds Model of Random Censorship

Censored data arise naturally in a number of fields, particularly in problems of reliability and survival analysis. There are several types of censoring, in this article, we will confine ourselves to the right randomly censoring type. Recently, Ahmadi et al. (2010 Ahmadi , J. , Doostparast , M. , Parsian , A. ( 2010 ). Bayes estimation based on random censored data for some life time models under symmetric and asymmetric loss functions . Communcations in Statistics-Theory and Methods , 39 : 3058 – 3071 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) considered the problem of estimating unknown parameters in a general framework based on the right randomly censored data. They assumed that the survival function of the censoring time is free of the unknown parameter. This assumption is sometimes inappropriate. In such cases, a proportional odds (PO) model may be more appropriate (Lam and Leung, 2001 Lam , K. F. , Leung , T. L. ( 2001 ). Marginal likelihood estimation for proportional odds models with right censored data . Lifetime Data Analysis 7 : 39 – 54 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]). Under this model, in this article, point and interval estimations for the unknown parameters are obtained. Since it is important to check the adequacy of models upon which inferences are based (Lawless, 2003 Lawless , J. F. (2003). Statistical Models and Methods for Lifetime Data. , 2nd ed. New York : John Wiley & Sons. [Google Scholar], p. 465), two new goodness-of-fit tests for PO model based on right randomly censored data are proposed. The proposed procedures are applied to two real data sets due to Smith (2002 Smith , P. J. ( 2002 ). Analysis of Failure and Survival Data . London : Chapman & Hall, CRC . [Google Scholar]). A Monte Carlo simulation study is conducted to carry out the behavior of the estimators obtained.     

Likelihood Ratio Order of m-Spacings in Multiple-Outlier Models

This article deals with the multiple-outlier exponential model. The likelihood ratio order between m-spacings of the combined sample is developed, some results extend the conclusions on simple spacings in Wen et al. (2007 Wen , S. , Lu , Q. , Hu , T. ( 2007 ). Likelihood ratio order of spacings of heterogeneous exponential random variables . J. Multivariate Anal. 98 : 743 – 756 .[Crossref], [Web of Science ®] , [Google Scholar]).     

Analyzing Supersaturated Designs by Means of an Information Based Criterion

The cost and time consumption of many industrial experimentations can be reduced using the class of supersaturated designs since this can be used for screening out the important factors from a large set of potentially active variables. A supersaturated design is a design for which there are fewer runs than effects to be estimated. Although there exists a wide study of construction methods for supersaturated designs, their analysis methods are yet in an early research stage. In this article, we propose a method for analyzing data using a correlation-based measure, named as symmetrical uncertainty. This method combines measures from the information theory field and is used as the main idea of variable selection algorithms developed in data mining. In this work, the symmetrical uncertainty is used from another viewpoint in order to determine more directly the important factors. The specific method enables us to use supersaturated designs for analyzing data of generalized linear models for a Bernoulli response. We evaluate our method by using some of the existing supersaturated designs, obtained according to methods proposed by Tang and Wu (1997 Tang , B. , Wu , C. F. J. (1997). A method for constructing supersaturated designs and its E(s ²)-optimality. Canadian Journal of Statistics 25:191–201.[Crossref], [Web of Science ®] , [Google Scholar]) as well as by Koukouvinos et al. (2008 Koukouvinos , C. , Mylona , K. , Simos , D. E. ( 2008 ). E(s ²)-optimal and minimax-optimal cyclic supersaturated designs via multi-objective simulated annealing . Journal of Statistical Planning and Inference 138 : 1639 – 1646 .[Crossref], [Web of Science ®] , [Google Scholar]). The comparison is performed by some simulating experiments and the Type I and Type II error rates are calculated. Additionally, Receiver Operating Characteristics (ROC) curves methodology is applied as an additional statistical tool for performance evaluation.     

On Testing Sample Selection Bias Under the Multicollinearity Problem

ABSTRACT

This paper reviews and extends the literature on the finite sample behavior of tests for sample selection bias. Monte Carlo results show that, when the “multicollinearity problem” identified by Nawata (1993 Nawata , K. ( 1993 ). A note on the estimation of models with sample-selection biases . Economics Letters 42 : 15 – 24 . [CSA] [CROSSREF] [Crossref], [Web of Science ®] , [Google Scholar]) is severe, (i) the t-test based on the Heckman–Greene variance estimator can be unreliable, (ii) the Likelihood Ratio test remains powerful, and (iii) nonnormality can be interpreted as severe sample selection bias by Maximum Likelihood methods, leading to negative Wald statistics. We also confirm previous findings (Leung and Yu, 1996 Leung , S. F. , Yu , S. ( 1996 ). On the choice between sample selection and two-part models . Journal of Econometrics 72 : 197 – 229 . [CSA] [CROSSREF] [Crossref], [Web of Science ®] , [Google Scholar]) that the standard regression-based t-test (Heckman, 1979 Heckman , J. J. ( 1979 ). Sample selection bias as a specification error . Econometrica 47 : 153 – 161 . [CSA] [Crossref], [Web of Science ®] , [Google Scholar]) and the asymptotically efficient Lagrange Multiplier test (Melino, 1982 Melino , A. ( 1982 ). Testing for sample selection bias . Review of Economic Studies 49 : 151 – 153 . [CSA] [Crossref], [Web of Science ®] , [Google Scholar]), are robust to nonnormality but have very little power.