Bayesian Robustness of the Positive False Discovery Rate

The positive false discovery rate was introduced by Storey (2003 Storey , J. D. (2003). The positive false discovery rate: a Bayesian interpretation and the q-value. Ann. Statist. 31:2013–2035.[Crossref], [Web of Science ®] , [Google Scholar]) as an alternative to the family wise error rate for the case in which we are simultaneously testing a large amount of hypotheses. The positive false discovery rate has a very nice Bayesian interpretation (as it was shown by Storey, 2003 Storey , J. D. (2003). The positive false discovery rate: a Bayesian interpretation and the q-value. Ann. Statist. 31:2013–2035.[Crossref], [Web of Science ®] , [Google Scholar]) and its robustness is analyzed. The emphasis is on the ε-contamination class (one of the most used classes of priors for Bayesian robustness) and it is shown that robustness is not obtained when the basic prior concentrates the probability on the null hypothesis.     

Symmetric quantiles and their applications

To develop estimators with stronger efficiencies than the trimmed means which use the empirical quantile, Kim (1992) Kim, S. J. 1992. The metrically trimmed means as a robust estimator of location. Annals of Statistics, 20: 1534–1547. [Crossref], [Web of Science ®] [Google Scholar] and Chen & Chiang (1996) Chen, L. A. and Chiang, Y. C. 1996. Symmetric type quantile and trimmed means for location and linear regression model. Journal of Nonparametric Statistics, 7: 171–185. [Taylor & Francis Online] [Google Scholar], implicitly or explicitly used the symmetric quantile, and thus introduced new trimmed means for location and linear regression models, respectively. This study further investigates the properties of the symmetric quantile and extends its application in several aspects. (a) The symmetric quantile is more efficient than the empirical quantiles in asymptotic variances when quantile percentage α is either small or large. This reveals that for any proposal involving the α th quantile of small or large α s, the symmetric quantile is the right choice; (b) a trimmed mean based on it has asymptotic variance achieving a Cramer-Rao lower bound in one heavy tail distribution; (c) an improvement of the quantiles-based control chart by Grimshaw & Alt (1997) Grimshaw, S. D. and Alt, F. B. 1997. Control charts for quantile function values. Journal of Quality Technology, 29: 1–7. [Taylor & Francis Online] [Google Scholar] is discussed; (d) Monte Carlo simulations of two new scale estimators based on symmetric quantiles also support this new quantile.     

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.     

Estimators for the Population Mean in the Case of Missing Data

We propose a class of estimators for the population mean when there are missing data in the data set. Obtaining the mean square error equations of the proposed estimators, we show the conditions where the proposed estimators are more efficient than the sample mean, ratio-type estimators, and the estimators in Singh and Horn (2000 Singh , S. , Horn , S. ( 2000 ). Compromised imputation in survey sampling . Metrika 51 : 267 – 276 .[Crossref], [Web of Science ®] , [Google Scholar]) and Singh and Deo (2003 Singh , S. , Deo , B. (2003). Imputation by power transformation. Statist. Pap. 44:555–579.[Crossref], [Web of Science ®] , [Google Scholar]) in the case of missing data. These conditions are also supported by a numerical example.     

The Alpha-power Tobit Model

The main object of this article is to propose an extension of the tobit model for which the error distribution follows the power-normal distribution (Gupta and Gupta, 2008 Gupta , D. , Gupta , R. C. ( 2008 ). Analyzing skewed data by power normal model . Test 17 : 197 – 210 .[Crossref], [Web of Science ®] , [Google Scholar]). Inference is dealt with by using the likelihood approach. Simulation studies and application to a real data set are used to demonstrate the usefulness of the extension.     

Computing the Exact Fisher Information Matrix of Periodic State-Space Models

The purpose of this article is to develop algorithms for computing the exact Fisher information matrix of periodic time-varying state-space models. We first present a relatively simple recursive algorithm which computes the elements of the exact information matrix without involving numerical differentiation, since all required derivatives are analytically evaluated. The proposed algorithm extends the procedure due to Cavanaugh and Shumway (1996 Cavanaugh , J. E. , Shumway , R. H. ( 1996 ). On computing the expected Fisher information matrix for state-space model parameters . Statist. Probab. Lett. 26 : 347 – 355 .[Crossref], [Web of Science ®] , [Google Scholar]) to the periodic state-space framework. Exploiting the approach used in Klein et al. (2000 Klein , A. , Mélard , G. , Zahaf , T. ( 2000 ). Construction of the exact Fisher information matrix of Gaussian time series models by means of matrix differential rules . Linear Alg. Applic. 321 : 209 – 232 .[Crossref], [Web of Science ®] , [Google Scholar]), a second algorithm is proposed in order to obtain the exact information matrix as a whole instead of element by element. The algorithms are first developed in a general framework and then specialized to the case of a periodic Gaussian vector autoregressive moving-average (PVARMA) model.     

Modeling of Inverse Gaussian Frailty Model for Bivariate Survival Data

Shared frailty models are often used to model heterogeneity in survival analysis. The most common shared frailty model is a model in which hazard function is a product of random factor (frailty) and baseline hazard function which is common to all individuals. There are certain assumptions about the baseline distribution and distribution of frailty. In this article, we consider inverse Gaussian distribution as frailty distribution and three different baseline distributions namely, Weibull, generalized exponential, and exponential power distribution. With these three baseline distributions, we propose three different inverse Gaussian shared frailty models. To estimate the parameters involved in these models we adopt Markov Chain Monte Carlo (MCMC) approach. We present a simulation study to compare the true values of the parameters with the estimated values. Also, we apply these three models to a real life bivariate survival data set of McGilchrist and Aisbett (1991 McGilchrist , C. A. , Aisbett , C. W. ( 1991 ). Regression with frailty in survival analysis . Biometrics 47 : 461 – 466 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) related to kidney infection and a better model is suggested for the data.     

Curve Fitting Under Jump and Peak Irregularities Using Local Linear Regression

In Desmet and Gijbels (2009 Desmet , L. , Gijbels , I. ( 2009 ). Local linear fitting and improved estimation near peaks . Canad. J. Statist. 37 ( 3 ): 453 – 475 .[Crossref] , [Google Scholar]), the problem of curve fitting on functions with peaks was addressed and a method was proposed that was building further on the one used in Gijbels et al. (2007 Gijbels , I. , Lambert , A. , Qiu , P. ( 2007 ). Jump-preserving regression and smoothing using local linear fitting: a compromise . Ann. Instit. Statist. Math. 59 : 235 – 272 .[Crossref], [Web of Science ®] , [Google Scholar]) when dealing with functions that have jump discontinuities. In this article, we propose a common framework for both estimation problems and an integrated procedure where an appropriate diagnostic quantity is a key ingredient. In particular we discuss practical procedure parameter selection considering target functions with an unknown number of irregularities of unknown type. The method is illustrated both on simulated and real life data sets.     

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.     

Relative Errors of Difference-Based Variance Estimators in Nonparametric Regression

Difference-based estimators for the error variance are popular since they do not require the estimation of the mean function. Unlike most existing difference-based estimators, new estimators proposed by Müller et al. (2003 Müller , U. , Schick , A. , Wefelmeyer , W. ( 2003 ). Estimating the error variance in nonparametric regression by a covariate-matched U-statistic . Statistics 37 : 179 – 188 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) and Tong and Wang (2005 Tong , T. , Wang , Y. ( 2005 ). Estimating residual variance in nonparametric regression using least squares . Biometrika 92 : 821 – 830 .[Crossref], [Web of Science ®] , [Google Scholar]) achieved the asymptotic optimal rate as residual-based estimators. In this article, we study the relative errors of these difference-based estimators which lead to better understanding of the differences between them and residual-based estimators. To compute the relative error of the covariate-matched U-statistic estimator proposed by Müller et al. (2003 Müller , U. , Schick , A. , Wefelmeyer , W. ( 2003 ). Estimating the error variance in nonparametric regression by a covariate-matched U-statistic . Statistics 37 : 179 – 188 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), we develop a modified version by using simpler weights. We further investigate its asymptotic property for both equidistant and random designs and show that our modified estimator is asymptotically efficient.     

Efficient Ratio and Product Estimators in Stratified Random Sampling

This paper suggests an efficient class of ratio and product estimators for estimating the population mean in stratified random sampling using auxiliary information. It is interesting to mention that, in addition to many, Koyuncu and Kadilar (2009 Koyuncu , N. , Kadilar , C. ( 2009 ). Ratio and product estimators in stratified random sampling . J. Statist. Plann. Infer. 139 : 2552 – 2558 .[Crossref], [Web of Science ®] , [Google Scholar]), Kadilar and Cingi (2003 Kadilar , C. , Cingi , H. ( 2003 ). Ratio estimator in stratified sampling . Biometr. J. 45 : 218 – 225 .[Crossref], [Web of Science ®] , [Google Scholar], 2005 Kadilar , C. , Cingi , H. ( 2005 ). A new estimator in stratified random sampling . Commun. Statist. Theor. Meth. 34 : 597 – 602 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), and Singh and Vishwakarma (2007 Singh , H. P. , Vishwakarma , G. K. ( 2007 ). Modified exponential ratio and product estimators for finite population mean in double sampling . Austr. J. Statist. 36 ( 3 ): 217 – 225 . [Google Scholar]) estimators are identified as members of the proposed class of estimators. The expressions of bias and mean square error (MSE) of the proposed estimators are derived under large sample approximation in general form. Asymptotically optimum estimator (AOE) in the class is identified alongwith its MSE formula. It has been shown that the proposed class of estimators is more efficient than combined regression estimator and Koyuncu and Kadilar (2009 Koyuncu , N. , Kadilar , C. ( 2009 ). Ratio and product estimators in stratified random sampling . J. Statist. Plann. Infer. 139 : 2552 – 2558 .[Crossref], [Web of Science ®] , [Google Scholar]) estimator. Moreover, theoretical findings are supported through a numerical example.     

The Asymptotic Approximation of EPMC for Linear Discriminant Rules Using a Moore-Penrose Inverse Matrix in High Dimension

We consider the discriminant rule in a high-dimensional setting, i.e., when the number of feature variables p is comparable to or larger than the number of observations N. The discriminant rule must be modified in order to cope with singular sample covariance matrix in high-dimension. One way to do so is by considering the Moor-Penrose inverse matrix. Recently, Srivastava (2006 Srivastava , M. S. ( 2006 ). Minimum distance classification rules for high dimensional data . J. Multivariate Anal. 97 : 2057 – 2070 .[Crossref], [Web of Science ®] , [Google Scholar]) proposed maximum likelihood ratio rule by using Moor-Penrose inverse matrix of sample covariance matrix. In this article, we consider the linear discriminant rule by using Moor-Penrose inverse matrix of sample covariance matrix (LDRMP). With the discriminant rule, the expected probability of misclassification (EPMC) is commonly used as measure of the classification accuracy. We investigate properties of EPMC for LDRMP in high-dimension as well as the one of the maximum likelihood rule given by Srivastava (2006 Srivastava , M. S. ( 2006 ). Minimum distance classification rules for high dimensional data . J. Multivariate Anal. 97 : 2057 – 2070 .[Crossref], [Web of Science ®] , [Google Scholar]). From our asymptotic results, we show that the classification accuracy of LDRMP depends on new distance. Additionally, our asymptotic result is verified by using the Monte Carlo simulation.     

Bayesian Capture-Recapture Analysis: An Application in Modeling Semelparity of a Neotropical Didelphid Marsupial

We develop a series of Bayesian statistical models for estimating survival of a neotropic didelphid marsupial, the Brazilian gracile mouse opossum (Gracilinanus microtarsus). These models are based on the Cormack–Jolly–Seber model (Cormack, 1964 Cormack , R. M. ( 1964 ). Estimates of survival from the sighting of marked animals . Biometrika 51 : 429 – 438 .[Crossref], [Web of Science ®] , [Google Scholar]; Jolly 1965 Jolly , G. M. ( 1965 ). Explicit estimates from capture-recapture data with both death and immigration stochastic model . Biometrika 52 : 225 – 247 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]; Seber 1965 Seber , G. A. F. ( 1965 ). A note on the multiple recapture census . Biometrika 52 : 249 – 259 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) with both survival and recapture rates expressed as a function of covariates using a logit link. The proposed models allow taking into account heterogeneity in capture probability caused by the existence of different groups of individuals in the population. The models were applied to two cohorts (Cohort, 2000, 2001) with the first one including 14 and the second one 15 sampling occasions. The best models for each of the cohorts indicate that G. microtarsus is best described as partially semelparous, a condition in which mortality after the first mating is high but graded over time, with a fraction of males surviving for a second breeding season (Boonstra, 2005 Boonstra , R. ( 2005 ). Equipped for life: the adaptive role of the stress axis in male mammals . Journal of Mammalogy 86 : 236 – 247 .[Crossref], [Web of Science ®] , [Google Scholar]).     

Data-Driven Tests for Trend

This article proposes a new nonparametric test for the ordered alternatives problem in the k-sample setting for null hypothesis of lack of trend. This article further elaborates upon and extends the results of Ledwina and Wy?upek (2012a Ledwina , T. , Wy?upek , G. ( 2012a ). Two-sample test against one-sided alternatives . Scand. J. Statist. 39 : 358 – 381 .[Crossref], [Web of Science ®] , [Google Scholar]) obtained for k = 2. Simulations show that the new test has high and stable power and is able to control the Type I error to satisfactory extent, thus solving the problem posed in Terpstra and Magel (2003 Terpstra , J. T. , Magel , R. C. ( 2003 ). A new nonparametric test for the ordered alternative problem . J. Nonparametr. Statist. 15 : 289 – 301 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]). Our theoretical results say that asymptotic errors of both kinds do not exceed significance level, thus implying that the test is asymptotically unbiased.     

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]).     

Asymptotic Properties of Conditional Quantile Estimator Under Left-Truncated and α-Mixing Conditions

In this article, we establish strong uniform convergence and asymptotic normality of estimators of conditional quantile and conditional distribution function for a left truncated model when the data exhibit some kind of dependence. It is assumed that the observations form a stationary α-mixing sequence. The results of Lemdani et al. (2009 Lemdani , M. , Ould-Saïd , E. , Poulin , N. ( 2009 ). Asymptotic properties of a conditional quantile estimator with randomly truncated data . J. Multivariate Anal. 100 : 546 – 559 .[Crossref], [Web of Science ®] , [Google Scholar]) are relaxed from the i.i.d. assumption to α-mixing setting. Finite sample behavior of the estimators is investigated via simulations as well.     

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.     

An Efficiency and Protection Degree-Based Comparison Among the Quantitative Randomized Response Strategies

There have been many alternative strategies for implementing sampling survey on quantitative characteristic of sensitive issues by using randomized response (RR) technique. The efficiency of most of those strategies has been improved by choosing the suitable design parameters of model. However, the two different procedures with pre-assigned design parameter values cannot ensure that they possess the same protection degree to the respondents. Some earlier comparisons of those strategies are inadequate (as in Eichhorn and Hayre, 1983 Eichhorn , B. H. , Hayre , L. S. ( 1983 ). Scrambled randomized response methods for obtaining sensitive quantitative data . J. Statist. Plann. Infer. 7 : 307 – 316 .[Crossref], [Web of Science ®] , [Google Scholar]; Gupta et al., 2002 Gupta , S. , Gupta , B. , Singh , S. ( 2002 ). Estimation of sensitivity level of personal interview survey questions . J. Statist. Plann. Infer. 100 ( 2 ): 239 – 247 .[Crossref], [Web of Science ®] , [Google Scholar]). Some literature contains a more comprehensive comparison based on efficiency and protection degree to the respondents among the qualitative characteristic RR techniques (see Bhargava and Singh, 2002 Bhargava , M. , Singh , R. ( 2002 ). On the efficiency comparison of certain randomized response strategies . Metrika 55 : 191 – 197 .[Crossref], [Web of Science ®] , [Google Scholar]; Nayak, 1994 Nayak , T. K. ( 1994 ). On randomized response surveys for estimating a proportion . Commun. Statist. Theor. Meth. 23 ( 11 ): 3303 – 3321 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]; Zaizai and Zankan, 2004 Zaizai , Y. , Zankan , N. ( 2004 ). A fair comparison of the randomized response strategies . Acta Mathematica Scientia 24A ( 3 ): 362 – 368 . [Google Scholar]). As far as the comparisons are concerned that are based on efficiency and protection degree to the respondents among the quantitative characteristic RR techniques, very few related studies have been found so far. The purpose of this article is to give a more adequate comparison among those earlier quantitative characteristic RR strategies. It is found that several important differences between the results obtained in this article and some known results exist. Therefore, these earlier RR strategies should be reevaluated.     

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.     

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.