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.     

The Finite Time Ruin Probability of a New Risk Model Based on Entrance Process

In this article, we consider a new insurance risk model based on the entrance process proposed in Li et al. (2005 Li , Z. , Zhu , J. , Chen , F. ( 2005 ). Study of a risk model based on the entrance process . Statist. Probab. Lett. 72 : 1 – 10 .[Crossref], [Web of Science ®] , [Google Scholar]), and investigate the finite time ruin probabilities of this model. It is showed that an exponential upper bound for the finite time ruin probability exists, when the distributions of the claim size are light tailed. Furthermore, when the distributions of the claim size are heavy tailed, an asymptotic formula for the finite time ruin probability is obtained.     

Performance Evaluation of Different Tests for Location Parameters

In this article, our objective is to evaluate the performance of different tests which are used to compare the equality of more than two location parameters. We have considered six tests (including some commonly used) in this study, one of which is parametric and the others are nonparametric. These tests include the usual F test (Fisher and Mackenzie, 1923 Fisher , R. A. , Mackenzie , M. A. ( 1923 ). Studies in crop variation. II. The manurial response of different potato . Journal of Agricultural Science 13 : 311 – 320 .[Crossref], [Web of Science ®] , [Google Scholar]), Kruskal–Wallis test (Kruskall and Wallis, 1952 Kruskall , W. H. , Wallis , W. A. ( 1952 ). Use of ranks in one-criterion variance analysis . Journal of American Statistical Association 47 : 583 – 621 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), Kolmogorov–Smirnov test (David, 1958 David , H. T. ( 1958 ). Three-sample Kolnogorov–Smirnov test . The Annals of Mathematical Statistics 29 : 842 – 851 .[Crossref] , [Google Scholar]), the g test (Stekler, 1987 Stekler , H. O. ( 1987 ). Who forecasters better? Journal of Business and Economic Statistics 5 : 155 – 158 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), f test (Batchelor, 1990 Batchelor , R. A. ( 1990 ). All forecasters are equal . Journal of Business and Economic Statistics 8 : 143 – 144 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), and Extension of Median test (as given in Daniel, 1990 Daniel , W. W. ( 1990 ). Applied Nonparametric Statistics. , 2nd ed. Duxbury Classic Series , Boston . [Google Scholar]). Performance of these tests are compared under different symmetric, skewed and contaminated probability distributions that include Normal, Cauchy, Uniform, Laplace, Lognormal, Exponential, Weibull, Gamma, t, Chi-square, Half Normal, Mixed Weibull, and Mixed Normal. Performances of these tests are measured in terms of power. We have suggested appropriate tests which may perform better under different situations. It is expected that researchers will find these results useful in decision making.     

Influence diagnostics in Birnbaum–Saunders nonlinear regression models

We consider the issue of assessing influence of observations in the class of Birnbaum–Saunders nonlinear regression models, which is useful in lifetime data analysis. Our results generalize those in Galea et al. [8 Galea, M., Leiva, V. and Paula, G. A. 2004. Influence diagnostics in log-Birnbaum–Saunders regression models. J. Appl. Stat., 31: 1049–1064. [Taylor & Francis Online], [Web of Science ®] [Google Scholar]] which are confined to Birnbaum–Saunders linear regression models. Some influence methods, such as the local influence, total local influence of an individual and generalized leverage are discussed. Additionally, the normal curvatures for studying local influence are derived under some perturbation schemes. We also give an application to a real fatigue data set.     

Tests of Multivariate Independence for Ordinal Data

Population and sample versions of Kendall and Spearman measures of association suitable for multivariate ordinal data are defined. The latter generalize the indices of dependence of Ruymgaart and van Zuijlen (1978 Ruymgaart , F. H. , van Zuijlen , M. C. A. ( 1978 ). Asymptotic normality of multivariate linear rank statistics in the non-i.i.d. case . Ann. Statist. 6 : 588 – 602 .[Crossref], [Web of Science ®] , [Google Scholar]), Joe (1990 Joe , H. ( 1990 ). Multivariate concordance . J. Multivariate Anal. 35 : 12 – 30 .[Crossref], [Web of Science ®] , [Google Scholar]), and Schmid and Schmidt (2007 Schmid , F. , Schmidt , R. ( 2007 ). Multivariate extensions of Spearman's rho and related statistics . Statist. Probab. Lett. 77 : 407 – 416 .[Crossref], [Web of Science ®] , [Google Scholar]) by allowing atoms in the underlying distribution. The representation of the proposed empirical measures as U-statistics enables to establish their asymptotic normality under general distributions. A special attention is given to tests of independence for multivariate ordinal data, where the power of the new methodologies are investigated under fixed and contiguous alternatives.     

A Simple Method for Constructing Multidimensional Distributions of Correlated Categorical Data

We suggest a method for constructing a multidimensional distribution of correlated categorical data with fixed marginal distributions and specified degrees of association based on the log-linear models. A convex combination approach by Lee (1997 Lee , A. J. ( 1997 ). Some simple methods for generating correlated categorical variates . Computational Statistics & Data Analysis 26 : 133 – 148 .[Crossref], [Web of Science ®] , [Google Scholar]) is applied to get a joint distribution with fixed Pearson chi-square coefficient. By using the suggested method, we can generate three-dimensional distributions which have a fixed association among three variables. Therefore, the suggested method could be extended to higher dimensions.     

Empirical Likelihood Method for Censored Median Regression Models

Based on the semiparametric median regression analysis for the right-censored data developed by Ying et al. (1995 Ying , Z. , Jung , S. H. , Wei , L. J. ( 1995 ). Survival analysis with median regression models . J. Amer. Statist. Assoc. 90 : 178 – 184 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), an empirical likelihood based inferential procedure for the regression coefficients is proposed. The limiting distribution of the proposed log-empirical likelihood ratio test statistic follows a chi-squared distribution, which corresponds to the standard asymptotic results of the empirical likelihood method. The inference about the subsets of the entire regression coefficients vector is discussed. The proposed method is illustrated by some simulation studies.     

Combining Unbiased Ridge and Principal Component Regression Estimators

In the presence of multicollinearity problem, ordinary least squares (OLS) estimation is inadequate. To circumvent this problem, two well-known estimation procedures often suggested are the unbiased ridge regression (URR) estimator given by Crouse et al. (1995 Crouse , R. , Jin , C. , Hanumara , R. ( 1995 ). Unbiased ridge estimation with prior information and ridge trace . Commun. Statist. Theor. Meth. 24 : 2341 – 2354 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) and the (r, k) class estimator given by Baye and Parker (1984 Baye , M. , Parker , D. ( 1984 ). Combining ridge and principal component regression: a money demand illustration . Commun. Statist. Theor. Meth. 13 : 197 – 205 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]). In this article, we proposed a new class of estimators, namely modified (r, k) class ridge regression (MCRR) which includes the OLS, the URR, the (r, k) class, and the principal components regression (PCR) estimators. It is based on a criterion that combines the ideas underlying the URR and the PCR estimators. The standard properties of this new class estimator have been investigated and a numerical illustration is done. The conditions under which the MCRR estimator is better than the other two estimators have been investigated.     

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.     

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.     

Saddlepoint Approximations to the Limiting Distribution of the Modified Anderson–Darling Test Statistic

The approximation for the distribution function of test statistic is extremely important in statistics. The standard and higher-order saddlepoint approximations are considered in tails of the limiting distribution for the modified Anderson–Darling test. The saddlepoint approximations are compared with the approximation of Sinclair et al. (1990 Sinclair , C. D. , Spurr , B. D. , Ahmad , M. I. ( 1990 ). Modified Anderson Darling test . Communication Statistics—Theory and Methods 19 : 3677 – 3686 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) for upper tail area. An empirical function is derived to estimate the critical values of a saddlepoint approximation.     

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.     

A Monte Carlo Study of Recent Ridge Parameters

A number of procedures have been developed for finding biased estimators of regression parameters. One of these procedures is the ridge regression. In this article, a new approach to obtain the ridge parameter K is suggested and then evaluated by Monte Carlo simulations. A number of different models are investigated for different number of observations, the strength of correlation between the explanatory variables, and distribution of the error terms. The mean squared error (MSE) criterion is used to examine the performance of the proposed estimators when compared with other well-known estimators. Under certain conditions, it is shown that at least one of the proposed estimators have a smaller MSE than the ordinary least squared estimator (OLS) and Hoerl and Kennard (1970a Hoerl , A. E. , Kennard , R. W. ( 1970a ). Ridge regression: biased estimation for non-orthogonal problems . Technometrics 12 : 55 – 67 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) estimator (HK).     

Strong Limit Theorems for Pairwise NQD Random Variables

A number of strong laws of large numbers for sequences of pairwise negative quadrant dependent (NQD) random variables have been established by using the generalized three series theorem. In this article, we obtain a strong law of large numbers by using the truncation technique and method of subsequences instead of the generalized three series theorem. Our result generalizes and improves on the corresponding one in Li and Yang (2008 Li , R. , Yang , W. ( 2008 ). Strong convergence of pairwise NQD random sequences . J. Math. Anal. Appl. 344 : 741 – 747 .[Crossref], [Web of Science ®] , [Google Scholar]). We also obtain a complete convergence result for an array of rowwise pairwise NQD random variables.     

Nonparametric Knn estimation with monotone constraints

ABSTRACT

The K-nearest-neighbor (Knn) method is known to be more suitable in fitting nonparametrically specified curves than the kernel method (with a globally fixed smoothing parameter) when data sets are highly unevenly distributed. In this paper, we propose to estimate a nonparametric regression function subject to a monotonicity restriction using the Knn method. We also propose using a new convergence criterion to measure the closeness between an unconstrained and the (monotone) constrained Knn-estimated curves. This method is an alternative to the monotone kernel methods proposed by Hall and Huang (2001 Hall, P., Huang, L.-S. (2001). Nonparametric kernel regression subject to monotonicity constraints. The Annals of Statistics 29:624–647.[Crossref], [Web of Science ®] , [Google Scholar]), and Du et al. (2013 Du, P., Parmeter, C. F., Racine, J. S. (2013). Nonparametric kernel regression with multiple predictors and multiple shape constraints. Statistic Sinica 23:1347–1371.[Web of Science ®] , [Google Scholar]). We use a bootstrap procedure for testing the validity of the monotone restriction. We apply our method to the “Job Market Matching” data taken from Gan and Li (2016 Gan, L., Li, Q. (2016). Efficiency of thin and thick market. Journal of Econometrics 192(1):40–54.[Crossref], [Web of Science ®] , [Google Scholar]) and find that the unconstrained/constrained Knn estimators work better than kernel estimators for this type of highly unevenly distributed data.     

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.     

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.     

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.”     

Generalized Symmetric Divergence Measures and the Probability of Error

There are three classical divergence measures exist in the literature on information theory and statistics. These are namely, Jeffryes-Kullback-Leiber (Jeffreys, 1946 Jeffreys , H. ( 1946 ). An invariant form for the prior probability in estimation problems . Proc. Roy. Soc. Lon. Ser. A 186 : 453 – 461 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]; Kullback and Leibler, 1951 Kullback , S. , Leibler , R. A. ( 1951 ). On information and sufficiency . Ann. Math. Statist. 22 : 79 – 86 .[Crossref] , [Google Scholar]) J-divergence. Sibson-Burbea-Rao (Sibson, 1969 Sibson , R. ( 1969 ). Information radius . Z. Wahrs. und verw Geb. 14 : 149 – 160 .[Crossref], [Web of Science ®] , [Google Scholar]), Jensen-Shannon divegernce, (Burbea and Rao, 1982 Burbea , J. , Rao , C. R. ( 1982 ). On the convexity of some divergence measures based on entropy functions . IEEE Trans. Inform. Theor. IT-28 : 489 – 495 .[Crossref], [Web of Science ®] , [Google Scholar]), and Taneja (1995 Taneja , I.J. ( 1995 ). New developments in generalized information measures , In: Hawkes , P. W. , Ed., Advances in Imaging and Electron Physics . 91:37–136 . [Google Scholar]). Arithmetic-Geometric divergence. These three measures bear an interesting relationship among each other. The divergence measures like Hellinger (1909 Hellinger , H. ( 1909 ). Neue Begründung der theorie der quadratischen formen von unendlichen vielen Veränderlichen . J. Reine Aug. Math. 136 : 210 – 271 .[Crossref] , [Google Scholar]) discrimination, symmetric χ²-divergence, and triangular discrimination are also known in the literature. In this article, we have considered generalized symmetric divergence measures having the measures given above as particular cases. Bounds on the probability of error are obtained in terms of generalized symmetric divergence measures. Study of bounds on probability of error is extended for the difference of divergence measures.     

Simultaneous Confidence Intervals for the Parameters of Pareto Distribution under Progressive Censoring

In this article, we consider the progressive Type II right censored sample from Pareto distribution. We introduce a new approach for constructing the simultaneous confidence interval of the unknown parameters of this distribution under progressive censoring. A Monte Carlo study is also presented for illustration. It is shown that this confidence region has a smaller area than that introduced by Ku? and Kaya (2007 Ku? , C. , Kaya , M. F. ( 2007 ). Estimation for the parameters of the Pareto distribution under progressive censoring . Commun. Statist. Theor. Meth. 36 : 1359 – 1365 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]).