A Comparison of Procedures for Controlling the False Discovery Rate in the Presence of Small Variance Genes: A Simulation Study

The Significance Analysis of Microarrays (SAM; Tusher et al., 2001 Tusher , V. G. , Tibshirani , R. , Chu , G. ( 2001 ). Significance analysis of microarrys applied to the ionizing radiation response . Proceedings of the National Academy of Sciences 98 : 5116 – 5121 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) method is widely used in analyzing gene expression data while controlling the FDR by using resampling-based procedure in the microarray setting. One of the main components of the SAM procedure is the adjustment of the test statistic. The introduction of the fudge factor to the test statistic aims at deflating the large value of test statistics due to the small standard error of gene-expression. Lin et al. (2008 Lin , D. , Shkedy , Z. , Burzykowski , T. , Göhlmann , H. W. H. , De Bondt , A. , Perera , T. , Geerts , T. , Bijnens , L. ( 2008 ). Significance analysis of microarray (SAM) for comparisons of several treatments with one control . Biometric Journal, MCP 50 ( 5 ): 801 – 823 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) pointed out that the fudge factor does not effectively improve the power and the control of the FDR as compared to the SAM procedure without the fudge factor in the presence of small variance genes. Motivated by the simulation results presented in Lin et al. (2008 Lin , D. , Shkedy , Z. , Burzykowski , T. , Göhlmann , H. W. H. , De Bondt , A. , Perera , T. , Geerts , T. , Bijnens , L. ( 2008 ). Significance analysis of microarray (SAM) for comparisons of several treatments with one control . Biometric Journal, MCP 50 ( 5 ): 801 – 823 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]), in this article, we extend our study to compare several methods for choosing the fudge factor in the modified t-type test statistics and use simulation studies to investigate the power and the control of the FDR of the considered methods.     

The Identification of Potential Longitudinal Biomarkers and Measurements of Effectiveness for Biomarkers as Surrogates in Multivariate Survival Data

A proposed method based on frailty models is used to identify longitudinal biomarkers or surrogates for a multivariate survival. This method is an extention of earlier models by Wulfsohn and Tsiatis (1997 Wulfsohn , M. S. , Tsiatis , A. A. ( 1997 ). A joint model for survival and longitudinal data measured with error . Biometrics 53 ( 1 ): 330 – 339 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) and Song et al. (2002 Song , X. , Davidian , M. , Tsiatis , A. A. ( 2002 ). A Semiparametric likelihood approach to joint modeling of longitudinal and time-to-event data . Biometrics 58 ( 4 ): 742 – 753 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]). In this article, similar to Henderson et al. (2002 Henderson , R. , Diggle , P. J. , Dobson , A. ( 2002 ). Identification and efficacy of longitudinal markers for survival . Biostatistics 3 ( 1 ): 33 – 50 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]), a joint likelihood function combines the likelihood functions of the longitudinal biomarkers and the multivariate survival times. We use simulations to explore how the number of individuals, the number of time points per individual and the functional form of the random effects from the longitudianl biomarkers influence the power to detect the association of a longitudinal biomarker and the multivariate survival time. The proposed method is illustrate by using the gastric cancer data.     

Estimating Incidence Rate on Current Status Data with Application to a Phase IV Cancer Trial

New drug discovery in the pediatrics has dramatically improved survival, but with long- term adverse events. This motivates the examination of adverse outcomes such as long-term toxicity in a phase IV trial. An ideal approach to monitor long-term toxicity is to systematically follow the survivors, which is generally not feasible. Instead, cross-sectional surveys are conducted in Hudson et al. (2007 Hudson , M. M. , Rai , S. N. , Nunez , C. , Merchant , T. E. , Marina , N. M. , Zalamea , N. , Cox , C. , Phipps , S. , Pompeu , R. , Rosenthal , D. ( 2007 ). Noninvasive evaluation of late anthracycline cardiac toxicity in childhood cancer survivors . J. Clin. Oncol. 25 : 3635 – 3643 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]), with one of the objectives to estimate the cumulative incidence rates along with specific interest in fixed-term (5 or 10 year) rates. We present inference procedures based on current status data to our motivating example with very interesting findings.     

Is Each NPMLE of a Continuous Bivariate Distribution Function with Singly Right-Censored Data Really Inconsistent?

We consider non-parametric estimation of a continuous cdf of a random vector (X ₁, X ₂). With bivariate RC data, it is stated in van der Laan (1996 Van der Laan , M. J. ( 1996 ) Efficient estimation in the bivariate censoring model and repairing NPMLE . Ann. Statist. 24 : 596 – 627 .[Crossref], [Web of Science ®] , [Google Scholar], p. 598¹⁰, Ann. Statist.), Quale et al. (2006 Quale , C. M. , van der Laan , M. J. , Robins , J. R. ( 2006 ). Locally efficient estimation with bivariate right-censored data . JASA. 101 : 1076 – 1084 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar], JASA) etc. that “it is well known that the NPMLE for continuous data is inconsistent (Tsai et al. (1986 Tsai , W. Y. , Leurgans , S. , Crowley , J. ( 1986 ). Nonparametric estimation of a bivariate survival function in the presence of censoring . Ann. Statist. 14 : 1351 – 1365 .[Crossref], [Web of Science ®] , [Google Scholar])).” The claim is based on a result in Tsai et al. (1986 Tsai , W. Y. , Leurgans , S. , Crowley , J. ( 1986 ). Nonparametric estimation of a bivariate survival function in the presence of censoring . Ann. Statist. 14 : 1351 – 1365 .[Crossref], [Web of Science ®] , [Google Scholar], p.1352, Ann. Statist.) that if X ₁ is right censored but not X ₂, then common ways for defining one NPMLE lead to inconsistency. If X ₁ is right censored and X ₂ is type I right-censored (which includes the case in Tsai et al.), we present a consistent NPMLE. The result corrects a common misinterpretation of Tsai's example (Tsai et al., 1986 Tsai , W. Y. , Leurgans , S. , Crowley , J. ( 1986 ). Nonparametric estimation of a bivariate survival function in the presence of censoring . Ann. Statist. 14 : 1351 – 1365 .[Crossref], [Web of Science ®] , [Google Scholar], Ann. Statist.).     

On Complete Convergence for Arrays of Dependent Random Variables

A complete convergence theorem for an array of rowwise independent random variables was established by Sung et al. (2005 Sung , S. H. , Volodin , A. I. , Hu , T.-C. ( 2005 ). More on complete convergence for arrays . Statist. Probab. Lett. 71 : 303 – 311 .[Crossref], [Web of Science ®] , [Google Scholar]). This result has been generalized and extended by Kruglov et al. (2006 Kruglov , V. M. , Volodin , A. I. , Hu , T.-C. ( 2006 ). On complete convergence for arrays . Statist. Probab. Lett. 76 : 1631 – 1640 .[Crossref], [Web of Science ®] , [Google Scholar]) and Chen et al. (2007 Chen , P. , Hu , T.-C. , Liu , X. , Volodin , A. ( 2007 ). On complete convergence for arrays of rowwise negatively associated random variables . Theor. Probab. Appl. 52 : 393 – 397 . [Google Scholar]). In this article, we extend the results of Sung et al. (2005 Sung , S. H. , Volodin , A. I. , Hu , T.-C. ( 2005 ). More on complete convergence for arrays . Statist. Probab. Lett. 71 : 303 – 311 .[Crossref], [Web of Science ®] , [Google Scholar]), Kruglov et al. (2006 Kruglov , V. M. , Volodin , A. I. , Hu , T.-C. ( 2006 ). On complete convergence for arrays . Statist. Probab. Lett. 76 : 1631 – 1640 .[Crossref], [Web of Science ®] , [Google Scholar]), and Chen et al. (2007 Chen , P. , Hu , T.-C. , Liu , X. , Volodin , A. ( 2007 ). On complete convergence for arrays of rowwise negatively associated random variables . Theor. Probab. Appl. 52 : 393 – 397 . [Google Scholar]) to an array of dependent random variables satisfying Hoffmann-Jørgensen type inequalities.     

Aspects of Dependence in Generalized Farlie-Gumbel-Morgenstern Distributions

In this article, we obtained a dependence measure for generalized Farlie-Gumbel-Morgenstern (FGM) family in view of Kochar and Gupta (1987 Kochar , S. G. , Gupta , R. P. ( 1987 ). Competitors of Kendall-tau test for testing independence against PQD . Biometrika 74 ( 3 ): 664 – 669 .[Crossref], [Web of Science ®] , [Google Scholar]) and then compared this measure with Spearman's rho and Kendall's tau in FGM family. Moreover, we evaluated the empirical power of the class of distribution-free tests proposed by Kochar and Gupta (1987 Kochar , S. G. , Gupta , R. P. ( 1987 ). Competitors of Kendall-tau test for testing independence against PQD . Biometrika 74 ( 3 ): 664 – 669 .[Crossref], [Web of Science ®] , [Google Scholar], 1990 Kochar , S. G. , Gupta , R. P. ( 1990 ). Distribution-free tests based on sub-sample extrema for testing against positive dependence . Australian Journal of Statistics 32 : 45 – 51 .[Crossref] , [Google Scholar]) based on exact distribution of a U-statistics. This is derived via a simulation study for sample of sizes n = 6, 8, 10, 12, 16, and 20. Also, we compared our simulation results with those achieved by Amini et al. (2010 Amini , M. , Jabbari , H. , Mohtashami Borzadaran , G. R. , Azadbakhsh , M. ( 2010 ). Power comparison of independence test for the Farlie-Gumbel-Moregenstern family . Communications of the Korean Statistical Society 17 ( 4 ): 493 – 505 .[Crossref] , [Google Scholar]) and Güven and Kotz (2008 Güven , B. , Kotz , S. ( 2008 ). Test of independence for generalized Farlie-Gumbel-Morgenstern distributions . Journal of Computational and Applied Mathematics 212 : 102 – 111 .[Crossref], [Web of Science ®] , [Google Scholar]).     

Improved strategy to collect sensitive data by using geometric distribution as a randomization device

Abstract

In this article, we improvise Singh and Grewal (2013 Singh, S., and I. S. Grewal. 2013. Geometric distribution as a randomization device implemented in the Kuk’s model. International Journal of Contemporary Mathematical Sciences 8:243–8.[Crossref] , [Google Scholar]) and Hussain et al. (2016 Hussain, Z., J. Shabbir, Z. Pervez, S. F. Shah, and M. Khan. 2016. Generalized geometric distribution of order k: A flexible choice to randomize the response. Communications in Statistics: Simulation and Computation 46:4708–21.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) techniques by introducing a new two-stage randomization response process. Using the proposed new technique, we achieve better efficiency and increasing protection of privacy of respondents than the Kuk (1990 Kuk, A. Y. C. 1990. Asking sensitive questions indirectly. Biometrika 77 (2):436–8.[Crossref], [Web of Science ®] , [Google Scholar]), Singh and Grewal (2013 Singh, S., and I. S. Grewal. 2013. Geometric distribution as a randomization device implemented in the Kuk’s model. International Journal of Contemporary Mathematical Sciences 8:243–8.[Crossref] , [Google Scholar]) and Hussain et al. (2016 Hussain, Z., J. Shabbir, Z. Pervez, S. F. Shah, and M. Khan. 2016. Generalized geometric distribution of order k: A flexible choice to randomize the response. Communications in Statistics: Simulation and Computation 46:4708–21.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) models. The relative efficiency and protection of the respondents of the proposed two-stage randomization device have been investigated through simulation study, and the situations are reported where the proposed estimator performs better than its competitors. The SAS code used to investigate the performance of the proposed strategy are also provided.     

A Generalization of the Skew-Normal Distribution: The Beta Skew-Normal

We consider a new generalization of the skew-normal distribution introduced by Azzalini (1985 Azzalini , A. ( 1985 ). A class of distributions which includes the normal ones . Scand. J. Statis. 12 ( 2 ): 171 – 178 .[Web of Science ®] , [Google Scholar]). We denote this distribution Beta skew-normal (BSN) since it is a special case of the Beta generated distribution (Jones, 2004 Jones , M. C. ( 2004 ). Families of distributions of order statistics . Test 13 ( 1 ): 1 – 43 .[Crossref], [Web of Science ®] , [Google Scholar]). Some properties of the BSN are studied. We pay attention to some generalizations of the skew-normal distribution (Bahrami et al., 2009 Bahrami , W. , Agahi , H. , Rangin , H. ( 2009 ). A two-parameter Balakrishnan skew-normal distribution . J. Statist. Res. Iran 6 : 231 – 242 . [Google Scholar]; Sharafi and Behboodian, 2008 Sharafi , M. , Behboodian , J. ( 2008 ). The Balakrishnan skew-normal density . Statist. Pap. 49 : 769 – 778 .[Crossref], [Web of Science ®] , [Google Scholar]; Yadegari et al., 2008 Yadegari , I. , Gerami , A. , Khaledi , M. J. ( 2008 ). A generalization of the Balakrishnan skew-normal distribution . Statist. Probab. Lett. 78 : 1165 – 1167 .[Crossref], [Web of Science ®] , [Google Scholar]) and to their relations with the BSN.     

The performance of the truncated mixed control chart

Abstract

When the mixed chart proposed by Aslam et al. (2015 Aslam, M., M. Azam, N. Khan, and C.-H. Jun. 2015. A mixed control chart to monitor the process. International Journal of Production Research 53 (15):4684–93. doi:10.1080/00207543.2015.1031354.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) is in use, the sample items are classified as defective or not defective and, depending on the number of defectives, the quality characteristic X of the sample items are also measured. In this case, an Xbar chart decides the state of the process. The previous conforming/non-conforming classification truncates the X distribution and, because of that, the mathematical development to obtain the ARLs is complex. Aslam et al. (2015 Aslam, M., M. Azam, N. Khan, and C.-H. Jun. 2015. A mixed control chart to monitor the process. International Journal of Production Research 53 (15):4684–93. doi:10.1080/00207543.2015.1031354.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) didn’t pay attention to the fact that the X distribution is truncated and, due to that, they obtained incorrect ARLs.     

A Broad Class of Partially Specified Autoregressions on Multi-Casting Data

Extending the bifurcating autoregressive (BAR) process (cf. Cowan and Staudte, 1986 Cowan , R. , Staudte , R. G. ( 1986 ). The bifurcating autoregression model in cell lineage studies . Biometrics 42 : 769 – 783 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) to multi-casting (multi-splitting) data, Hwang and Choi (2009 Hwang , S. Y. , Choi , M. S. ( 2009 ). Modeling and large sample estimation for multi-casting autoregression . Statist. Prob. Lett. 79 : 1943 – 1950 .[Crossref], [Web of Science ®] , [Google Scholar]) introduced multi-casting autoregression (MCAR, for short) defined on multi-casting tree structured data. This article is concerned with the case when the MCAR model is partially specified only through conditional mean and variance without directly imposing autoregressive (AR) structure. The resulting class of models will be referred to as P-MCAR (partially specified MCAR). The P-MCAR considerably enlarges the class of multi-casting models including (as special cases) MCAR, random coefficient MCAR, conditionally heteroscedastic multi-casting models and binomial-thinning processes. Moment structures for this broad P-MCAR class are investigated. Least squares (LS) estimation method is discussed and asymptotic relative efficiency (ARE) of the generalized-LS over ordinary-LS is obtained in a closed form. A simulation study is conducted to illustrate results.     

A Variance Estimator for Systematic Sampling from a Deliberately Ordered Population

ABSTRACT

The systematic sampling (SYS) design (Madow and Madow, 1944 Madow , L. H. , Madow , W. G. ( 1944 ). On the theory of systematic sampling . Ann. Math. Statist. 15 : 1 – 24 .[Crossref] , [Google Scholar]) is widely used by statistical offices due to its simplicity and efficiency (e.g., Iachan, 1982 Iachan , R. ( 1982 ). Systematic sampling a critical review . Int. Statist. Rev. 50 : 293 – 303 .[Crossref], [Web of Science ®] , [Google Scholar]). But it suffers from a serious defect, namely, that it is impossible to unbiasedly estimate the sampling variance (Iachan, 1982 Iachan , R. ( 1982 ). Systematic sampling a critical review . Int. Statist. Rev. 50 : 293 – 303 .[Crossref], [Web of Science ®] , [Google Scholar]) and usual variance estimators (Yates and Grundy, 1953 Yates , F. , Grundy , P. M. ( 1953 ). Selection without replacement from within strata with probability proportional to size . J. Roy. Statist. Soc. Series B 1 : 253 – 261 . [Google Scholar]) are inadequate and can overestimate the variance significantly (Särndal et al., 1992 Särndal , C. E. , Swenson , B. , Wretman , J. H. ( 1992 ). Model Assisted Survey Sampling . New York : Springer-Verlag , Ch. 3 .[Crossref] , [Google Scholar]). We propose a novel variance estimator which is less biased and that can be implemented with any given population order. We will justify this estimator theoretically and with a Monte Carlo simulation study.     

Methods for Bivariate Survival Data with Mismeasured Covariates Under an Accelerated Failure Time Model

Accelerated failure time models are useful in survival data analysis, but such models have received little attention in the context of measurement error. In this paper we discuss an accelerated failure time model for bivariate survival data with covariates subject to measurement error. In particular, methods based on the marginal and joint models are considered. Consistency and efficiency of the resultant estimators are investigated. Simulation studies are carried out to evaluate the performance of the estimators as well as the impact of ignoring the measurement error of covariates. As an illustration we apply the proposed methods to analyze a data set arising from the Busselton Health Study (Knuiman et al., 1994 Knuiman , M. W. , Cullent , K. J. , Bulsara , M. K. , Welborn , T. A. , Hobbs , M. S. T. ( 1994 ). Mortality trends, 1965 to 1989, in Busselton, the site of repeated health surveys and interventions . Austral. J. Public Health 18 : 129 – 135 . [CSA] [Crossref], [PubMed] , [Google Scholar]).     

Improved Predictive Inference for Responses Based on Location Parameter Median from the Rayleigh Model Given a Doubly Censored Sample

The purpose of this article is to investigate the predictive inference for responses from the location parameter mean as well as from the median given a doubly censored sample from the two-parameter Rayleigh model. The predictive results by Khan et al. (2010 Khan , H. M. R. , Provost , S. B. , Ashima , S. ( 2010 ). Predictive Inference from a Two-Parameter Rayleigh Life Model Given a Doubly Censored Sample . Commun. Statist. Theor. Meth. 39 ( 7 ): 1237 – 1246 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) are used to obtain the predictive inference for responses from the median, where Khan et al. (2010 Khan , H. M. R. , Provost , S. B. , Ashima , S. ( 2010 ). Predictive Inference from a Two-Parameter Rayleigh Life Model Given a Doubly Censored Sample . Commun. Statist. Theor. Meth. 39 ( 7 ): 1237 – 1246 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) obtained the future estimates from the mean. A numerical example representing 66 liver cancer patients is used for predictive analysis. It is concluded that the predictive inference from the median gives precise results as compared with the location parameter mean.     

On Simultaneous Estimation of the Number of Signals and Frequencies of Multiple Sinusoids When Some Observations Are Missing

In an earlier article (Bai et al., 1999 Bai , Z. D. , Rao , C. R. , Wu , Y. H. , Zen , M. M. , Zhao , L. C. ( 1999 ). The simultaneous estimation of the number of signals and frequencies of multiple sinusods when some observations are missing: I. Asymptotics . Proc. Natl. Acad. Sci. 96 : 11,106 – 11,110 . [CSA] [Crossref], [Web of Science ®] , [Google Scholar]), the problem of simultaneous estimation of the number of signals and frequencies of multiple sinusoids is considered in the case that some observations are missing. The number of signals is estimated with an information theoretic criterion and the frequencies are estimated with eigenvariation linear prediction. Asymptotic properties of the procedure are investigated but the Monte Carlo simulation is not performed. In this article, a slightly different but scale invariant criterion for detection is proposed and the estimation of frequencies remains the same. Asymptotic properties of this new procedure are provided. Monte Carlo Simulation for both procedures is carried out. Furthermore, comparison on the real signals is also given.     

Using Frailty Models to Identify the Longitudinal Biomarkers in Survival Analysis

We introduce a score test to identify longitudinal biomarkers or surrogates for a time to event outcome. This method is an extension of Henderson et al. (2000 Henderson , R. , Diggle , P. , Dobson , A. ( 2000 ). Joint modelling of longitudinal measurements and event time data . Biostatistics 1 ( 4 ): 465 – 480 .[Crossref], [PubMed] , [Google Scholar], 2002 Henderson , R. , Diggle , P. , Dobson , A. ( 2002 ). Identification and efficacy of longitudinal markers for survival . Biostatistics 3 ( 1 ): 33 – 50 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]). In this article, a score test is based on a joint likelihood function which combines the likelihood functions of the longitudinal biomarkers and the survival times. Henderson et al. (2000 Henderson , R. , Diggle , P. , Dobson , A. ( 2000 ). Joint modelling of longitudinal measurements and event time data . Biostatistics 1 ( 4 ): 465 – 480 .[Crossref], [PubMed] , [Google Scholar], 2002 Henderson , R. , Diggle , P. , Dobson , A. ( 2002 ). Identification and efficacy of longitudinal markers for survival . Biostatistics 3 ( 1 ): 33 – 50 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) assumed that the same random effect exists in the longitudinal component and in the Cox model and then they can derive a score test to determine if a longitudinal biomarker is associated with time to an event. We extend this work and our score test is based on a joint likelihood function which allows other random effects to be present in the survival function.

Considering heterogeneous baseline hazards in individuals, we use simulations to explore how the factors can influence the power of a score test to detect the association of a longitudinal biomarker and the survival time. These factors include the functional form of the random effects from the longitudinal biomarkers, in the different number of individuals, and time points per individual. We illustrate our method using a prothrombin index as a predictor of survival in liver cirrhosis patients.     

On the Feasibility of Basing Horvitz and Thompson's Estimator on a Sample by Rao,Hartley, and Cochran's Scheme

Formulae for the first and second order inclusion probabilities for the Rao et al. (1962 Rao , J. N. K. , Hartley , H. O. , Cochran , W. G. ( 1962 ). On a simple procedure of unequal probability sampling without replacement . J. Roy. Statist. Soc. B 24 : 482 – 491 . [CSA]  [Google Scholar]) (RHC) scheme of sampling are derived. They enable one to evaluate, for a sample drawn according to the RHC scheme, the Horvitz and Thompson's (1952 Horvitz , D. G. , Thompson , D. J. ( 1952 ). A generalization of sampling without replacement from a finite universe . J. Amer. Statist. Assoc. 47 : 663 – 685 . [CSA] [Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) estimator (HTE) along with its unbiased variance estimator given by Yates and Grundy (1953 Yates , F. , Grundy , P. M. ( 1953 ). Selection without replacement from within strata with probability proportional to size . J. Amer. Statist. Assoc. 75 : 206 – 211 . [CSA]  [Google Scholar]). So, for a sample at hand thus drawn one may choose between the RHCE and the HTE for use on finding which one has the smaller coefficient of variation.     

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.     

Second-Order Covariance Matrix Formula for Heteroskedastic Generalized Linear Models

For the first time, we provide a matrix formula for second-order covariances of maximum likelihood estimates in heteroskedastic generalized linear models, thus generalizing the results of Cordeiro (2004 Cordeiro , G. M. ( 2004 ). Second-order covariance matrix of maximum likelihood estimates in generalized linear models . Statist. Probab. Lett. 66 : 153 – 160 .[Crossref], [Web of Science ®] , [Google Scholar]) and Cordeiro et al. (2006 Cordeiro , G. M. , Barroso , L. P. , Botter , D. A. (2006). Covariance matrix formula for generalized linear models with unknown dispersion. Commun. Statist. Theor. Meth. 35:113–120.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) related to the generalized linear models with known and unknown dispersion parameter, respectively. The covariance matrix formula does not involve cumulants of log-likelihood derivatives and can be easily obtained using simple matrix operations. We apply our main result to a simple model. Some simulations show that the second-order covariances can be quite pronounced in small to moderate samples. The usual covariances of the maximum likelihood estimates can be corrected by these second-order covariances.     

A side sensitive modified group runs control chart to detect shifts in the process mean

Gadre and Rattihalli [5 Gadre, M. P. and Rattihalli, R. N. 2006. Modified group runs control charts to detect increases in fraction non-conforming and shifts in the process mean. Commun. Stat. Simul. Comput., 35: 225–240. [Taylor & Francis Online], [Web of Science ®] [Google Scholar]] have introduced the Modified Group Runs (MGR) control chart to identify the increases in fraction non-conforming and to detect shifts in the process mean. The MGR chart reduces the out-of-control average time-to-signal (ATS), as compared with most of the well-known control charts. In this article, we develop the Side Sensitive Modified Group Runs (SSMGR) chart to detect shifts in the process mean. With the help of numerical examples, it is illustrated that the SSMGR chart performs better than the Shewhart's X¯ chart, the synthetic chart [12 Wu, Z. and Spedding, T. A. 2000. A synthetic control chart for detecting small shifts in the process mean. J. Qual. Technol., 32: 32–38. [Taylor & Francis Online], [Web of Science ®] [Google Scholar]], the Group Runs chart [4 Gadre, M. P. and Rattihalli, R. N. 2004. A group runs control chart for detecting shifts in the process mean. Econ. Qual. Control, 19: 29–43. [Crossref] [Google Scholar]], the Side Sensitive Group Runs chart [6 Gadre, M. P. and Rattihalli, R. N. 2007. A side sensitive group runs control chart for detecting shifts in the process mean. Stat. Methods Appl., 16: 27–37. [Crossref], [Web of Science ®] [Google Scholar]], as well as the MGR chart [5 Gadre, M. P. and Rattihalli, R. N. 2006. Modified group runs control charts to detect increases in fraction non-conforming and shifts in the process mean. Commun. Stat. Simul. Comput., 35: 225–240. [Taylor & Francis Online], [Web of Science ®] [Google Scholar]]. In some situations it is also superior to the Cumulative Sum chart p9 Page, E. S. 1954. Continuous inspection schemes. Biometrika, 41: 100–114. [Crossref], [Web of Science ®] [Google Scholar]] and the exponentially weighed moving average chart [10 Roberts, S. W. 1959. Control chart tests based on geometric moving averages. Technometrics, 1: 239–250. [Taylor & Francis Online] [Google Scholar]]. In the steady state also, its performance is better than the above charts.     

A Note on Importance Resampling for Multi-Dimensional Statistics

Johns (1988 Johns , M. V. (1988). Importance sampling for bootstrap confidence intervals. Journal of the American Statistical Association 83:709–714.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), Davison (1988 Davison , A. C. ( 1988 ). Discussion of paper by D. V. Hinkley . Journal of the Royal Statistical Society Series B 50 : 356 – 357 . [Google Scholar]), and Do and Hall (1991 Do , K. A. , Hall , P. ( 1991 ). On importance sampling for the bootstrap . Biometrika 78 : 161 – 167 .[Crossref], [Web of Science ®] , [Google Scholar]) used importance sampling for calculating bootstrap distributions of one-dimensional statistics. Realizing that their methods can not be extended easily to multi-dimensional statistics, Fuh and Hu (2004 Fuh , C. D. , Hu , I. ( 2004 ). Efficient importance sampling for events of moderate deviations with applications . Biometrika 91 : 471 – 490 .[Crossref], [Web of Science ®] , [Google Scholar]) proposed an exponential tilting formula for statistics of multi-dimension, which is optimal in the sense that the asymptotic variance is minimized for estimating tail probabilities of asymptotically normal statistics. For one-dimensional statistics, Hu and Su (2008 Hu , J. , Su , Z. ( 2008 ). Adaptive resampling algorithms for estimating bootstrap distributions . Journal of Statistical Planning and Inference 138 ( 6 ): 1763 – 1777 .[Crossref], [Web of Science ®] , [Google Scholar]) proposed a multi-step variance minimization approach that can be viewed as a generalization of the two-step variance minimization approach proposed by Do and Hall (1991 Do , K. A. , Hall , P. ( 1991 ). On importance sampling for the bootstrap . Biometrika 78 : 161 – 167 .[Crossref], [Web of Science ®] , [Google Scholar]). In this article, we generalize the approach of Hu and Su (2008 Hu , J. , Su , Z. ( 2008 ). Adaptive resampling algorithms for estimating bootstrap distributions . Journal of Statistical Planning and Inference 138 ( 6 ): 1763 – 1777 .[Crossref], [Web of Science ®] , [Google Scholar]) to multi-dimensional statistics, which applies to general statistics and does not resort to asymptotics. Empirical results on a real survival data set show that the proposed algorithm provides significant computational efficiency gains.