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

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

On Testing Sample Selection Bias Under the Multicollinearity Problem

ABSTRACT

This paper reviews and extends the literature on the finite sample behavior of tests for sample selection bias. Monte Carlo results show that, when the “multicollinearity problem” identified by Nawata (1993 Nawata , K. ( 1993 ). A note on the estimation of models with sample-selection biases . Economics Letters 42 : 15 – 24 . [CSA] [CROSSREF] [Crossref], [Web of Science ®] , [Google Scholar]) is severe, (i) the t-test based on the Heckman–Greene variance estimator can be unreliable, (ii) the Likelihood Ratio test remains powerful, and (iii) nonnormality can be interpreted as severe sample selection bias by Maximum Likelihood methods, leading to negative Wald statistics. We also confirm previous findings (Leung and Yu, 1996 Leung , S. F. , Yu , S. ( 1996 ). On the choice between sample selection and two-part models . Journal of Econometrics 72 : 197 – 229 . [CSA] [CROSSREF] [Crossref], [Web of Science ®] , [Google Scholar]) that the standard regression-based t-test (Heckman, 1979 Heckman , J. J. ( 1979 ). Sample selection bias as a specification error . Econometrica 47 : 153 – 161 . [CSA] [Crossref], [Web of Science ®] , [Google Scholar]) and the asymptotically efficient Lagrange Multiplier test (Melino, 1982 Melino , A. ( 1982 ). Testing for sample selection bias . Review of Economic Studies 49 : 151 – 153 . [CSA] [Crossref], [Web of Science ®] , [Google Scholar]), are robust to nonnormality but have very little power.     

The Asymptotic Estimate of Ruin Probability Under a Class of Risk Model in the Presence of Heavy Tails

In contrast with the classical Cramér–Lundberg model where the premium process is a linear function of time, we consider the ruin probability under the risk model where the aggregate premium consists of both a compound Poisson process and a linear process of time. Moreover, a constant interest force is also taken into account in our model. We restrict ourselves to the case where the claim size is heavy-tailed, i.e., the equilibrium distribution function of the claim size belongs to a wide subclass of the subexponential distribution. An asymptotic formula for the ruin probability is obtained by using the similar method of Kalashnikov and Konstantinides (2000 Kalashnikov , V. , Konstantinides , D. ( 2000 ). Ruin under interest force and subexponential claims: a simple treament . Insur. Math. Econ. 27 : 145 – 149 .[Crossref], [Web of Science ®] , [Google Scholar]). The asymptotic formula we get here is the same as the one in Asmussen (1998 Asmussen , S. ( 1998 ). Subexponential asymptotics for stochastic processes: extremal behaviour, stationary distribution and first passage probabilities . Ann. Appl. Probab. 8 : 354 – 374 .[Crossref], [Web of Science ®] , [Google Scholar]), Klüppelberg and Stadtmüller (1998 Klüppelberg , C. , Stadtmüller , U. ( 1998 ). Ruin probabilities in the presence of heavy-tails and interest rates . Scand. Actuarial J. 1 : 49 – 58 .[Taylor & Francis Online] , [Google Scholar]), and Kalashnikov and Konstantinides (2000 Kalashnikov , V. , Konstantinides , D. ( 2000 ). Ruin under interest force and subexponential claims: a simple treament . Insur. Math. Econ. 27 : 145 – 149 .[Crossref], [Web of Science ®] , [Google Scholar]) which did not consider the stochastic premium.     

A Simulation Study on the Probability of Correct Selection for Large k Populations

An increasing number of contemporary datasets are high dimensional. Applications require these datasets be screened (or filtered) to select a subset for further study. Multiple testing is the standard tool in such applications, although alternatives have begun to be explored. In order to assess the quality of selection in these high-dimensional contexts, Cui and Wilson (2008b Cui , X. , Wilson , J. ( 2008b ). On the probability of correct selection for large k populations with application to microarray data . Biometrical Journal 50 ( 5 ): 870 – 883 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) proposed two viable methods of calculating the probability that any such selection is correct (PCS). PCS thereby serves as a measure of the quality of competing statistics used for selection. The first simulation study of this article investigates the two PCS statistics of the above article. It shows that in the high-dimensional case PCS can be accurately estimated and is robust under certain conditions. The second simulation study investigates a nonparametric estimator of PCS.     

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.     

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

D s -Optimal Designs for Quadratic Log Contrast Model for Experiments with Mixtures

The approximate D _s -optimal design for discriminating between linear and quadratic log contrast models for experiments with mixtures suggested by Aitchison and Bacon-Shone (1984 Aitchison , J. , Bacon-Shone , J. ( 1984 ). Log contrast models for experiments with mixtures . Biometrika 71 : 323 – 330 .[Crossref], [Web of Science ®] , [Google Scholar]) is investigated, where the experimental domain is restricted further as in Chan (1992 Chan , L. Y. ( 1992 ). D-Optimal design for a quadratic log contrast model for experiments with mixtures . Commun. Statist. Theor. Meth. 21 ( 10 ): 2909 – 2930 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]). It is found that for a symmetric subspace of the finite dimensional simplex, there is a D _s -optimal design with the nice structure that puts a weight 1/2 ^k?1 on the centroid of this subspace and the remaining weight is uniformly distributed on the vertices of the experimental domain. Finally, the D _s -efficiency of the D-optimal design for quadratic model and the design given by Aitchison and Bacon-Shone (1984 Aitchison , J. , Bacon-Shone , J. ( 1984 ). Log contrast models for experiments with mixtures . Biometrika 71 : 323 – 330 .[Crossref], [Web of Science ®] , [Google Scholar]) are also discussed.     

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.     

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

A Nonparametric Test for Symmetry Based on Freeman and Halton's Ideas on Contingency Tables

In this article, we propose a nonparametric method to test for symmetry in bivariate data. By using the extension of Fisher's exact treatment for 2 × 2 contingency tables proposed by Freeman and Halton (1951 Freeman , G. H. , Halton , J. H. ( 1951 ). Note on an exact treatment of contingency tables, goodness of fit and other problems of significance . Biometrika 38 : 141 – 149 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]), we can test the hypothesis of equal distribution for two samples of integer valued variables. Then, by counting the number of observations belonging to each cell of a symmetric, appropriately built grid, we can produce the two samples of integers required to use this test for equal distribution. The resulting test for symmetry is potentially extendible to higher dimensions. A simulation study is performed to compare with some known tests (Bowker, 1948 Bowker , A. H. ( 1948 ). A test for symmetry in contingency tables . Journal of the American Statistical Association 43 : 572 – 574 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]; Hollander, 1971 Hollander , M. ( 1971 ). A nonparametric test for bivariate symmetry . Biometrika 58-1 : 203 – 212 .[Crossref], [Web of Science ®] , [Google Scholar]; and its improvement given in Krampe and Kuhnt, 2007 Krampe , A. , Kuhnt , S. ( 2007 ). Bowker's test for symmetry and modifications within the algebraic framework . Computational Statistics and Data Analysis 51 : 4124 – 4142 .[Crossref], [Web of Science ®] , [Google Scholar]). Our proposal represents a competitive option as a test for symmetry.     

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.     

Covariate-Adjusted Partially Linear Regression Models

Motivated by covariate-adjusted regression (CAR) proposed by Sentürk and Müller (2005 Sentürk , D. , Müller , H. G. ( 2005 ). Covariate-adjusted regression . Biometrika 92 : 75 – 89 .[Crossref], [Web of Science ®] , [Google Scholar]) and an application problem, in this article we introduce and investigate a covariate-adjusted partially linear regression model (CAPLM), in which both response and predictor vector can only be observed after being distorted by some multiplicative factors, and an additional variable such as age or period is taken into account. Although our model seems to be a special case of covariate-adjusted varying coefficient model (CAVCM) given by Sentürk (2006 Sentürk , D. ( 2006 ). Covariate-adjusted varying coefficient models . Biostatistics 7 : 235 – 251 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]), the data types of CAPLM and CAVCM are basically different and then the methods for inferring the two models are different. In this article, the estimate method motivated by Cui et al. (2008 Cui , X. , Guo , W. S. , Lin , L. , Zhu , L. X. ( 2008 ). Covariate-adjusted nonlinear regression . Ann. Statist. 37 : 1839 – 1870 . [Google Scholar]) is employed to infer the new model. Furthermore, under some mild conditions, the asymptotic normality of estimator for the parametric component is obtained. Combined with the consistent estimate of asymptotic covariance, we obtain confidence intervals for the regression coefficients. Also, some simulations and a real data analysis are made to illustrate the new model and methods.     

Hazards Regression for Cancer Studies

The complication in analyzing tumor data is that the tumors detected in a screening program tend to be slowly progressive tumors, which is the so-called length-biased sampling that is inherent in screening studies. Under the assumption that all subjects have the same tumor growth function, Ghosh (2008 Ghosh , D. ( 2008 ). Proportional hazards regression for cancer studies . Biometrics 64 : 141 – 148 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) developed estimation procedures for proportional hazards model. In this article, by modeling growth function as a function of covariates, we demonstrate that Ghosh (2008 Ghosh , D. ( 2008 ). Proportional hazards regression for cancer studies . Biometrics 64 : 141 – 148 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar])'s approach can be extended to the case when each subject has a specific growth function. A simulation study is conducted to demonstrate the potential usefulness of the proposed estimators for the regression parameters in the proportional and additive hazards model.     

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.     

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.     

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.