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.     

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.     

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

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.     

On Joint Determination of the Number of States and the Number of Variables in Markov-Switching Models: A Monte Carlo Study

In this article, we examine the performance of two newly developed procedures that jointly select the number of states and variables in Markov-switching models by means of Monte Carlo simulations. They are Smith et al. (2006 Smith , A. , Naik , P. A. , Tsai , C. ( 2006 ). Markov-switching model selection using Kullback–Leibler divergence . Journal of Econometrics 134 ( 2 ): 553 – 577 .[Crossref], [Web of Science ®] , [Google Scholar]) and Psaradakis and Spagnolo (2006 Psaradakis , Z. , Spagnolo , N. ( 2006 ). Joint determination of the state dimension and autoregressive order for models with Markov regime switching . Journal of Time Series Analysis 27 ( 2 ): 753 – 766 .[Crossref], [Web of Science ®] , [Google Scholar]), respectively. The former develops Markov switching criterion (MSC) designed specifically for Markov-switching models, while the latter recommends the use of standard complexity-penalised information criteria (BIC, HQC, and AIC) in joint determination of the state dimension and the autoregressive order of Markov-switching models. The Monte Carlo evidence shows that BIC outperforms MSC while MSC and HQC are preferable over AIC.     

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

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.     

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.     

Bias-Corrected Maximum Likelihood Estimators in Nonlinear Heteroscedastic Models

Nonlinear heteroscedastic models are widely used in econometrics and statistical applications. We derive matrix formulae for the second-order biases of the maximum likelihood estimators of the parameters in the mean and variance response which generalize previous results by Cook et al. (1986 Cook , D. R. , Tsai , C. L. , Wei , B. C. ( 1986 ). Bias in nonlinear regression . Biometrika 73 : 615 – 623 .[Crossref], [Web of Science ®] , [Google Scholar]) and Cordeiro (1993 Cordeiro , G. M. ( 1993 ). Bartlett corrections and bias correction for two heteroscedastic regression models . Commun. Statist. Theor. Meth. 22 : 169 – 188 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]). The biases of the estimators are easily obtained as vectors of regression coefficients from suitable weighted linear regressions. The practical use of such biases is illustrated in a simulation study and in an application to a real data set.     

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.     

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.     

The Performance of Panel Cointegration Methods: Results from a Large Scale Simulation Study

This article presents results concerning the performance of both single equation and system panel cointegration tests and estimators. The study considers the tests developed in Pedroni (1999 Pedroni , P. ( 1999 ). Critical values for cointegration tests in heterogeneous panels with multiple regressors . Oxford Bulletin of Economics and Statistics 61 : 653 – 670 .[Crossref], [Web of Science ®] , [Google Scholar], 2004 Pedroni , P. ( 2004 ). Panel cointegration. Asymptotic and finite sample properties of pooled time series tests with an application to the PPP hypothesis . Econometric Theory 20 : 597 – 625 .[Crossref], [Web of Science ®] , [Google Scholar]), Westerlund (2005 Westerlund , J. ( 2005 ). New simple tests for panel cointegration . Econometric Reviews 24 : 297 – 316 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), Larsson et al. (2001 Larsson , R. , Lyhagen , J. , Löthgren , M. ( 2001 ). Likelihood-based cointegration tests in heterogeneous panels . Econometrics Journal 4 : 109 – 142 .[Crossref] , [Google Scholar]), and Breitung (2005 Breitung , J. ( 2005 ). A parametric approach to the estimation of cointegration vectors in panel data . Econometric Reviews 24 : 151 – 173 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) and the estimators developed in Phillips and Moon (1999 Phillips , P. C. B. , Moon , H. R. ( 1999 ). Linear regression limit theory for nonstationary panel data . Econometrica 67 : 1057 – 1111 .[Crossref], [Web of Science ®] , [Google Scholar]), Pedroni (2000 Pedroni , P. ( 2000 ). Fully modified OLS for heterogeneous cointegrated panels . In: Baltagi , B. H. , ed. Nonstationary Panels, Panel Cointegration, and Dynamic Panels . Amsterdam : Elsevier , pp. 93 – 130 .[Crossref] , [Google Scholar]), Kao and Chiang (2000 Kao , C. , Chiang , M.-H. ( 2000 ). On the estimation and inference of a cointegrated regression in panel data . In: Baltagi , B. H. , ed. Nonstationary Panels, Panel Cointegration, and Dynamic Panels . Amsterdam : Elsevier , pp. 179 – 222 .[Crossref] , [Google Scholar]), Mark and Sul (2003 Mark , N. C. , Sul , D. ( 2003 ). Cointegration vector estimation by panel dynamic OLS and long-run money demand . Oxford Bulletin of Economics and Statistics 65 : 655 – 680 .[Crossref], [Web of Science ®] , [Google Scholar]), Pedroni (2001 Pedroni , P. ( 2001 ). Purchasing power parity tests in cointegrated panels . Review of Economics and Statistics 83 : 1371 – 1375 . [Google Scholar]), and Breitung (2005 Breitung , J. ( 2005 ). A parametric approach to the estimation of cointegration vectors in panel data . Econometric Reviews 24 : 151 – 173 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]). We study the impact of stable autoregressive roots approaching the unit circle, of I(2) components, of short-run cross-sectional correlation and of cross-unit cointegration on the performance of the tests and estimators. The data are simulated from three-dimensional individual specific VAR systems with cointegrating ranks varying from zero to two for fourteen different panel dimensions. The usual specifications of deterministic components are considered.     

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.     

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.     

Testing for sphericity in a two-way error components panel data model

This article is concerned with sphericity test for the two-way error components panel data model. It is found that the John statistic and the bias-corrected LM statistic recently developed by Baltagi et al. (2011 Baltagi, B. H., Feng, Q., Kao, C. (2011). Testing for sphericity in a fixed effects panel data model. Econometrics Journal 14:25–47.[Crossref], [Web of Science ®] , [Google Scholar])Baltagi et al. (2012 Baltagi, B. H., Feng, Q., Kao, C. (2012). A Lagrange multiplier test for cross-sectional dependence in a fixed effects panel data model. Journal of Econometrics 170:164–177.[Crossref], [Web of Science ®] , [Google Scholar], which are based on the within residuals, are not helpful under the present circumstances even though they are in the one-way fixed effects model. However, we prove that when the within residuals are properly transformed, the resulting residuals can serve to construct useful statistics that are similar to those of Baltagi et al. (2011 Baltagi, B. H., Feng, Q., Kao, C. (2011). Testing for sphericity in a fixed effects panel data model. Econometrics Journal 14:25–47.[Crossref], [Web of Science ®] , [Google Scholar])Baltagi et al. (2012 Baltagi, B. H., Feng, Q., Kao, C. (2012). A Lagrange multiplier test for cross-sectional dependence in a fixed effects panel data model. Journal of Econometrics 170:164–177.[Crossref], [Web of Science ®] , [Google Scholar]). Simulation results show that the newly proposed statistics perform well under the null hypothesis and several typical alternatives.     

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.     

Skew-Normal ARMA Models with Nonlinear Heteroscedastic Predictors

Multivariate skew-normal (SN) distributions (Azzalini and Dalla Valle, 1996 Azzalini , A. , Dalla Valle , A. ( 1996 ). The multivariate skew-normal distribution . Biometrika 83 : 715 – 726 .[Crossref], [Web of Science ®] , [Google Scholar]) enjoy some of the useful properties of normal distributions, have nonlinear heteroscedastic predictors but lack the closure property of normal distributions (the sum of independent SN random variables is not SN). Recently, there has been a proliferation of classes of SN distributions with certain closure properties, one of the most promising being the closed skew-normal (CSN) distributions of González-Farías et al. (2004 González-Farías , G. , Dominguez-Molina , J. A. , Gupta , A. K. ( 2004 ). Additive properties of skew-normal random vectors . J. Statist. Plann. Infer. 126 : 521 – 534 .[Crossref], [Web of Science ®] , [Google Scholar]). We study the construction of stationary SN ARMA models for colored SN noise and show that their finite-dimensional distributions are skew-normal, seldom strictly stationary and their covariance functions differ from their normal ARMA counterparts in that they do not converge to zero for large lags. The situation is better for ARMA models driven by CSN noise, but at the additional cost of considerable computational complexity and a less explicit skewness parameter. In view of these results, the widespread use of such classes of SN distributions in the framework of ARMA models seem doubtful.     

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