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.     

DSGE Models with Student-t Errors

This paper deals with Dynamic Stochastic General Equilibrium (DSGE) models under a multivariate student-t distribution for the structural shocks. Based on the solution algorithm of Klein (2000) and the gamma-normal representation of the t-distribution, the TaRB-MH algorithm of Chib and Ramamurthy (2010 Chib , S. , Ramamurthy , S. ( 2010 ). Tailored randomized block MCMC methods with application to DSGE models . Journal of Econometrics 108 : 19 – 38 .[Crossref], [Web of Science ®] , [Google Scholar]) is used to estimate the model. A technique for estimating the marginal likelihood of the DSGE student-t model is also provided. The methodologies are illustrated first with simulated data and then with the DSGE model of Ireland (2004 Ireland , P. N. ( 2004 ). Technology shocks in the new keynesian model . Review of Economics and Statistics 86 ( 4 ): 923 – 936 .[Crossref], [Web of Science ®] , [Google Scholar]) where the results support the t-error model in relation to the Gaussian model.     

An Improved Solution for Multiple Comparisons with a Control in Response Surface Methodology with Two Predictor Variables

Sa and Edwards (1993 Sa , P. , Edwards , D. ( 1993 ). Multiple comparisons with a control in response surface methodology . Technometrics 35 ( 4 ): 436 – 445 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) first proposed the Multiple Comparisons with a Control problem in Response Surface Methodology. They provided an exact solution for one predictor variable and a conservative solution when number of predictor variables is more than one. Merchant et al. (1998 Merchant , A. , McCann , M. , Edwards , D. ( 1998 ). Improved multiple comparisons with a control in response surface analysis . Technometrics 40 ( 4 ): 291 – 303 .[Crossref], [Web of Science ®] , [Google Scholar]) improved the solution for the latter case. This article improves Merchant et al.'s solution for the case of rotatable designs in two predictor variables.     

Empirical Likelihood Ratio for Linear Transformation Models with Doubly Censored Data

Double censoring arises when T represents an outcome variable that can only be accurately measured within a certain range, [L, U], where L and U are the left- and right-censoring variables, respectively. When L is always observed, we consider the empirical likelihood inference for linear transformation models, based on the martingale-type estimating equation proposed by Chen et al. (2002 Chen , K. , Jin , Z. , Ying , Z. ( 2002 ). Semiparametric analysis of transformation models with censored data . Biometrika 89 : 659 – 668 .[Crossref], [Web of Science ®] , [Google Scholar]). It is demonstrated that both the approach of Lu and Liang (2006 Lu , W. , Liang , Y. ( 2006 ). Empirical likelihood inference for linear transformation models . Journal of Multivariate Analysis 97 : 1586 – 1599 .[Crossref], [Web of Science ®] , [Google Scholar]) and that of Yu et al. (2011 Yu , W. , Sun , Y. , Zheng , M. ( 2011 ). Empirical likelihood method for linear transformation models . Annals of the Institute of Statistical Mathematics 63 : 331 – 346 .[Crossref], [Web of Science ®] , [Google Scholar]) can be extended to doubly censored data. Simulation studies are conducted to investigate the performance of the empirical likelihood ratio methods.     

Generalized Scrambling in Quantitative Optional Randomized Response Models

Huang (2010 Huang , K. C. ( 2010 ). Unbiased estimators of mean, variance and sensitivity level for quantitative characteristics in finite population sampling . Metrika 71 : 341 – 352 .[Crossref], [Web of Science ®] , [Google Scholar]) proposed an optional randomized response model using a linear combination scrambling which is a generalization of the multiplicative scrambling of 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]) and the additive scrambling of Gupta et al. (2006, 2010). In this article, we discuss two main issues. (1) Can the Huang (2010 Huang , K. C. ( 2010 ). Unbiased estimators of mean, variance and sensitivity level for quantitative characteristics in finite population sampling . Metrika 71 : 341 – 352 .[Crossref], [Web of Science ®] , [Google Scholar]) model be improved further by using a two-stage approach?; (2) Does the linear combination scrambling provide any benefit over the additive scrambling of Gupta et al. (2010 Gupta , S. N. , Shabbir , J. , Sehra , S. ( 2010 ). Mean and sensitivity estimation in optional randomized response models . J. Statist. Plann. Infer. 140 : 2870 – 2874 .[Crossref], [Web of Science ®] , [Google Scholar])? We will note that the answer to the first question is “yes” but the answer to the second question is “no.”     

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.     

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.     

Bayesian Inference in Dynamic Disequilibrium Models: An Application to the Polish Credit Market

We propose a Bayesian approach for inference in a dynamic disequilibrium model. To circumvent the difficulties raised by the Maddala and Nelson (1974 Maddala , G. , Nelson , F. ( 1974 ). Maximum likelihood methods for models of markets in disequilibrium . Econometrica 42 ( 6 ): 1013 – 1030 .[Crossref], [Web of Science ®] , [Google Scholar]) specification in the dynamic case, we analyze a dynamic extended version of the disequilibrium model of Ginsburgh et al. (1980 Ginsburgh , V. , Tishler , A. , Zang , I. ( 1980 ). Alternative estimation methods for two regime models . European Economic Review 13 : 207 – 228 .[Crossref], [Web of Science ®] , [Google Scholar]). We develop a Gibbs sampler based on the simulation of the missing observations. The feasibility of the approach is illustrated by an empirical analysis of the Polish credit market, for which we conduct a specification search using the posterior deviance criterion of Spiegelhalter et al. (2002 Spiegelhalter , D. J. , Best , N. G. , Carlin , B. P. , Van Der Linde , A. ( 2002 ). Bayesian measures of complexity and fit . Journal of the Royal Statistical Society, Ser. B 64 : 583 – 639 .[Crossref] , [Google Scholar]).     

Identification-robust moment-based tests for Markov switching in autoregressive models

ABSTRACT

This paper develops tests of the null hypothesis of linearity in the context of autoregressive models with Markov-switching means and variances. These tests are robust to the identification failures that plague conventional likelihood-based inference methods. The approach exploits the moments of normal mixtures implied by the regime-switching process and uses Monte Carlo test techniques to deal with the presence of an autoregressive component in the model specification. The proposed tests have very respectable power in comparison with the optimal tests for Markov-switching parameters of Carrasco et al. (2014 Carrasco, M., Hu, L., Ploberger, W. (2014). Optimal test for Markov switching parameters. Econometrica 82(2):765–784.[Crossref], [Web of Science ®] , [Google Scholar]), and they are also quite attractive owing to their computational simplicity. The new tests are illustrated with an empirical application to an autoregressive model of USA output growth.     

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.     

Designs Robust Against the Presence of an Aberration in a Diallel Cross Design Setup

In this article, we find designs insensitive to the presence of an outlier in a diallel cross design setup for estimating a complete set of orthonormal contrasts among the effects of the general combining abilities of a set of parental lines. The criterion of robustness, suggested by Mandal (1989 Mandal , N. K. ( 1989 ). On robust designs . Cal. Stat. Assoc. Bull. 38 : 115 – 119 . [Google Scholar]) in block design setup and used by Biswas (2012 Biswas , A. ( 2012 ). Block designs robust against the presence of an aberration in a treatment-control setup . Commun. Statist. Theor. Meth. 41 : 920 – 933 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) in treatment-control setup, is adapted here. Complete diallel cross designs, suggested by Gupta and Kageyama (1994 Gupta , S. , Kageyama , S. ( 1994 ). Optimal complete diallel crosses . Biometrika 81 : 420 – 424 .[Crossref], [Web of Science ®] , [Google Scholar]), and partial diallel cross designs, suggested by Gupta et al. (1995 Gupta , S. , Das , A. , Kageyama , S. ( 1995 ). Single replicate orthogonal block designs for circulant partial diallel crosses . Commun. Statist. A 24 : 2601 – 2607 .[Taylor & Francis Online] , [Google Scholar]) and Mukerjee (1997 Mukerjee , R. ( 1997 ). Optimal partial diallel crosses . Biometrika 84 : 939 – 948 .[Crossref], [Web of Science ®] , [Google Scholar]), are found to be robust under certain conditions.     

A Unified Approach to Ordering Comparison of GPS Distributions with their Mixtures

Recently, several authors have been concerned with ordering comparison of known distributions of the family of generalized power series (GPS) distributions with their mixtures in various senses. In this article, we shall employ a unified approach and obtain similar results, more generally, for all members of the class of the GPS distributions. Some of the previous findings of Misra et al. (2003 Misra , N. , Singh , H. , Harner , E. J. ( 2003 ). Stochastic comparisons of Poisson and binomial random variables with their mixtures . Statist. Probab. Lett. 65 : 279 – 290 .[Crossref], [Web of Science ®] , [Google Scholar]), Alamatsaz and Abbasi (2008 Alamatsaz , M. H. , Abbasi , S. ( 2008 ). Ordering comparison of negative binomial random variables with their mixtures . Statist. Probab. Lett. 78 : 2234 – 2239 . [Google Scholar]), and Aghababaei Jazi and Alamatsaz (2010 Aghababaei Jazi , M. , Alamatsaz , M. H. ( 2010 ). Ordering comparison of logarithmic series random variables with their mixtures . Commun. Statist. Theor. Meth. 39 : 3252 – 3263 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) in this connection, then, follow as corollaries. Further, we have derived some more ordering comparison results.     

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.     

Robustness of the EWMA control chart for individual observations

The traditional exponentially weighted moving average (EWMA) chart is one of the most popular control charts used in practice today. The in-control robustness is the key to the proper design and implementation of any control chart, lack of which can render its out-of-control shift detection capability almost meaningless. To this end, Borror et al. [5 Borror, C. M., Montgomery, D. C. and Runger, G. C. 1999. Robustness of the EWMA control chart to non-normality. J. Qual. Technol., 31(3): 309–316. [Taylor & Francis Online], [Web of Science ®] [Google Scholar]] studied the performance of the traditional EWMA chart for the mean for i.i.d. data. We use a more extensive simulation study to further investigate the in-control robustness (to non-normality) of the three different EWMA designs studied by Borror et al. [5 Borror, C. M., Montgomery, D. C. and Runger, G. C. 1999. Robustness of the EWMA control chart to non-normality. J. Qual. Technol., 31(3): 309–316. [Taylor & Francis Online], [Web of Science ®] [Google Scholar]]. Our study includes a much wider collection of non-normal distributions including light- and heavy-tailed and symmetric and asymmetric bi-modal as well as the contaminated normal, which is particularly useful to study the effects of outliers. Also, we consider two separate cases: (i) when the process mean and standard deviation are both known and (ii) when they are both unknown and estimated from an in-control Phase I sample. In addition, unlike in the study done by Borror et al. [5 Borror, C. M., Montgomery, D. C. and Runger, G. C. 1999. Robustness of the EWMA control chart to non-normality. J. Qual. Technol., 31(3): 309–316. [Taylor & Francis Online], [Web of Science ®] [Google Scholar]], the average run-length (ARL) is not used as the sole performance measure in our study, we consider the standard deviation of the run-length (SDRL), the median run-length (MDRL), and the first and the third quartiles as well as the first and the 99th percentiles of the in-control run-length distribution for a better overall assessment of the traditional EWMA chart's in-control performance. Our findings sound a cautionary note to the (over) use of the EWMA chart in practice, at least with some types of non-normal data. A summary and recommendations are provided.     

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

A Benchmarking Approach to Temporal Disaggregation of Economic Time Series by Related Series

Several methods have been devised to deal with the problem of temporal disaggregation of economic time series (a) either when related series are available or (b) when only aggregate figures exist. In this article, we propose a statistical model-based approach to temporal disaggregation of economic time series by related series. The proposed approach is performed in two stages. In the first stage, we evaluate a preliminary estimate of the disaggregated series using a regression model for the disaggregated series and related series observed in the same frequency. The preliminary estimate of disaggregated series obtained in the first step is not consistent with aggregate figures. To ensure consistency we propose in the second stage, the use of a modified benchmarking approach based on signal extraction (Hillmer and Trabelsi, 1987 Hillmer , S. C. , Trabelsi , A. ( 1987 ). Benchmarking of economic time series . J. Amer. Statist. Assoc. 82 ( 400 ): 1064 – 1071 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]; Trabelsi and Hillmer, 1990 Trabelsi , A. , Hillmer , S. C. ( 1990 ). Benchmarking time series with reliable benchmarks . Appl. Statist. 39 ( 3 ): 367 – 379 .[Crossref], [Web of Science ®] , [Google Scholar]) to adjust the preliminary estimate of disaggregate series. The approach developed here is used for Seasonally Adjusted (SA) and Not Seasonally Adjusted (NSA) data. A comparison with previous temporal disaggregation methods has been done.     

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

Improved Estimation of Finite Population Variance Using Auxiliary Information

This article addresses the problem of estimating of finite population variance using auxiliary information in simple random sampling. A ratio-cum-difference type class of estimators for population variance has been suggested with its properties under large sample approximation. It has been shown that the suggested class of estimators is more efficient than usual unbiased, difference, Das and Tripathi (1978 Das , A. K. , Tripathi , T. P. ( 1978 ). Use of auxiliary information in estimating the finite population variance . Sankhya C 40 : 139 – 148 . [Google Scholar]), Isaki (1983 Isaki , C. T. ( 1983 ). Variance estimation using auxiliary information . J. Amer. Statist. Assoc. 78 : 117 – 123 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), Singh et al. (1988 Singh , H. P. , Upadhyaya , L. N. , Namjoshi , U. D. ( 1988 ). Estimation of finite population variance . Curr. Sci. 57 : 1331 – 1334 .[Web of Science ®] , [Google Scholar]), Kadilar and Cingi (2006 Kadilar , C. , Cingi , H. ( 2006 ). Ratio estimators for the population variance in simple and stratified random sampling . Appl. Math. Comp. 173 ( 2 ): 1047 – 1059 .[Crossref], [Web of Science ®] , [Google Scholar]), and other estimators/classes of estimators. In addition, we support this theoretical result with the aid of a empirical study.     

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.     

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.