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.     

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.     

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.     

Regression M-estimators with One Modified Form of Doubly Censored Data

In this article, we consider the M-estimators for the linear regression model when both response and covariate variables are subject to double censoring. The proposed estimators are constructed as some functional of three types of estimators for a bivariate survival distribution. The first two estimators are the generalizations of the Campbell and Földes (1982 Campbell, G. and Földes, A. 1982. “Large sample properties of nonparametric statistical inference”. In Nonparametric Statistical Inference., Edited by: Gnredenko, B. V., Puri, M. L. and Vineze, I. 103–122. Amsterdam: North-Holland.  [Google Scholar]) and Dabrowska (1988 Dabrowska, D. M. 1988. Kaplan-Meier estimate on the plane. Annals of Statistics, 18: 1475–1489. [Crossref], [Web of Science ®] , [Google Scholar]) estimators proposed by Shen (2009 Shen, P. S. 2009. Nonparametric estimation of the bivariate survival function one modified form of doubly censored data. Computational Statistics, 25: 203–313. [Crossref], [Web of Science ®] , [Google Scholar]). The third estimator is the generalization of the Prentice and Cai (1992 Prentice, R. L. and Cai, J. 1992. Covariance and survivor function estimation using censored multivariate failure time data. Biometrika, 79: 495–512. [Crossref], [Web of Science ®] , [Google Scholar]) estimator. The consistency of the proposed M-estimators is established. A simulation study is conducted to investigate the performance of the proposed estimators. Furthermore, the simple bootstrap methods are used to estimate standard deviations and construct interval estimators.     

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.     

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.     

Additive Mixed Effect Model for Clustered Doubly Censored Data

Cai and Zeng (2011 Cai, J. and Zeng, D. 2011. Additive mixed effect model for clustered failure time data. Biometrics, 67(4): 1340–1351. [Crossref], [PubMed] , [Google Scholar]) proposed an additive mixed effect model to analyze clustered right-censored data. In this article, we demonstrate that the approach of Cai and Zeng (2011 Cai, J. and Zeng, D. 2011. Additive mixed effect model for clustered failure time data. Biometrics, 67(4): 1340–1351. [Crossref], [PubMed] , [Google Scholar]) can be extended to clustered doubly censored data. Furthermore, when both left- and right-censoring variables are always observed, we propose alternative estimators using the approach of Cai and Cheng (2004 Cai, T. and Cheng, S. C. 2004. Semiparametric regression analysis for doubly censored data. Biometrika, 91: 277–290. [Crossref], [Web of Science ®] , [Google Scholar]). A simulation study is conducted to investigate the performance of the proposed estimators.     

The performance of the truncated mixed control chart

Abstract

When the mixed chart proposed by Aslam et al. (2015 Aslam, M., M. Azam, N. Khan, and C.-H. Jun. 2015. A mixed control chart to monitor the process. International Journal of Production Research 53 (15):4684–93. doi:10.1080/00207543.2015.1031354.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) is in use, the sample items are classified as defective or not defective and, depending on the number of defectives, the quality characteristic X of the sample items are also measured. In this case, an Xbar chart decides the state of the process. The previous conforming/non-conforming classification truncates the X distribution and, because of that, the mathematical development to obtain the ARLs is complex. Aslam et al. (2015 Aslam, M., M. Azam, N. Khan, and C.-H. Jun. 2015. A mixed control chart to monitor the process. International Journal of Production Research 53 (15):4684–93. doi:10.1080/00207543.2015.1031354.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) didn’t pay attention to the fact that the X distribution is truncated and, due to that, they obtained incorrect ARLs.     

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.     

Semiparametric analysis of 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. In this note, using Martingale arguments of Chen et al. [3 Chen, K., Jin, Z. and Ying, Z. 2002. Semiparametric analysis of transformation models with censored data. Biometrika, 89: 659–668. [Crossref], [Web of Science ®] [Google Scholar]], we propose an estimator (denoted by ?β) for estimating regression coefficients of transformation model when L is always observed. Under Cox proportional hazards model, the proposed estimator is equivalent to the partial likelihood estimator for left-truncated and right-censored data if the left-censoring variables L were regarded as left-truncated variables. In this case, the estimator ?β can be obtained by the standard software. A simulation study is conducted to investigate the performance of ?β. For the purpose of comparison, the simulation study also includes the estimator proposed by Cai and Cheng [2 Cai, T. and Cheng, S. 2004. Semiparametric regression analysis for doubly censored data. Biometrika, 91: 277–290. [Crossref], [Web of Science ®] [Google Scholar]] for the case when L and U are always observed.     

Is Each NPMLE of a Continuous Bivariate Distribution Function with Singly Right-Censored Data Really Inconsistent?

We consider non-parametric estimation of a continuous cdf of a random vector (X ₁, X ₂). With bivariate RC data, it is stated in van der Laan (1996 Van der Laan , M. J. ( 1996 ) Efficient estimation in the bivariate censoring model and repairing NPMLE . Ann. Statist. 24 : 596 – 627 .[Crossref], [Web of Science ®] , [Google Scholar], p. 598¹⁰, Ann. Statist.), Quale et al. (2006 Quale , C. M. , van der Laan , M. J. , Robins , J. R. ( 2006 ). Locally efficient estimation with bivariate right-censored data . JASA. 101 : 1076 – 1084 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar], JASA) etc. that “it is well known that the NPMLE for continuous data is inconsistent (Tsai et al. (1986 Tsai , W. Y. , Leurgans , S. , Crowley , J. ( 1986 ). Nonparametric estimation of a bivariate survival function in the presence of censoring . Ann. Statist. 14 : 1351 – 1365 .[Crossref], [Web of Science ®] , [Google Scholar])).” The claim is based on a result in Tsai et al. (1986 Tsai , W. Y. , Leurgans , S. , Crowley , J. ( 1986 ). Nonparametric estimation of a bivariate survival function in the presence of censoring . Ann. Statist. 14 : 1351 – 1365 .[Crossref], [Web of Science ®] , [Google Scholar], p.1352, Ann. Statist.) that if X ₁ is right censored but not X ₂, then common ways for defining one NPMLE lead to inconsistency. If X ₁ is right censored and X ₂ is type I right-censored (which includes the case in Tsai et al.), we present a consistent NPMLE. The result corrects a common misinterpretation of Tsai's example (Tsai et al., 1986 Tsai , W. Y. , Leurgans , S. , Crowley , J. ( 1986 ). Nonparametric estimation of a bivariate survival function in the presence of censoring . Ann. Statist. 14 : 1351 – 1365 .[Crossref], [Web of Science ®] , [Google Scholar], Ann. Statist.).     

Aspects of Dependence in Generalized Farlie-Gumbel-Morgenstern Distributions

In this article, we obtained a dependence measure for generalized Farlie-Gumbel-Morgenstern (FGM) family in view of Kochar and Gupta (1987 Kochar , S. G. , Gupta , R. P. ( 1987 ). Competitors of Kendall-tau test for testing independence against PQD . Biometrika 74 ( 3 ): 664 – 669 .[Crossref], [Web of Science ®] , [Google Scholar]) and then compared this measure with Spearman's rho and Kendall's tau in FGM family. Moreover, we evaluated the empirical power of the class of distribution-free tests proposed by Kochar and Gupta (1987 Kochar , S. G. , Gupta , R. P. ( 1987 ). Competitors of Kendall-tau test for testing independence against PQD . Biometrika 74 ( 3 ): 664 – 669 .[Crossref], [Web of Science ®] , [Google Scholar], 1990 Kochar , S. G. , Gupta , R. P. ( 1990 ). Distribution-free tests based on sub-sample extrema for testing against positive dependence . Australian Journal of Statistics 32 : 45 – 51 .[Crossref] , [Google Scholar]) based on exact distribution of a U-statistics. This is derived via a simulation study for sample of sizes n = 6, 8, 10, 12, 16, and 20. Also, we compared our simulation results with those achieved by Amini et al. (2010 Amini , M. , Jabbari , H. , Mohtashami Borzadaran , G. R. , Azadbakhsh , M. ( 2010 ). Power comparison of independence test for the Farlie-Gumbel-Moregenstern family . Communications of the Korean Statistical Society 17 ( 4 ): 493 – 505 .[Crossref] , [Google Scholar]) and Güven and Kotz (2008 Güven , B. , Kotz , S. ( 2008 ). Test of independence for generalized Farlie-Gumbel-Morgenstern distributions . Journal of Computational and Applied Mathematics 212 : 102 – 111 .[Crossref], [Web of Science ®] , [Google Scholar]).     

On Complete Convergence in Marcinkiewicz-Zygmund Type SLLN for END Random Variables and Its Applications

Abstract

In this paper, the complete convergence for maximal weighted sums of extended negatively dependent (END, for short) random variables is investigated. Some sufficient conditions for the complete convergence and some applications to a nonparametric model are provided. The results obtained in the paper generalize and improve the corresponding ones of Wang et al. (2014 Wang, X. J., X. Deng, L. L. Zheng, and S. H. Hu. 2014. Complete convergence for arrays of rowwise negatively superadditive-dependent random variables and its applications. A Journal of Theoretical and Applied Statistics 48(4):834–50. [Google Scholar]b) and Shen, Xue, and Wang (2017 Shen, A., M. Xue, and W. Wang. 2017. Complete convergence for weighted sums of extended negatively dependent random variables. Communications in Statistics – Theory and Methods 46(3):1433–44.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]).     

MODIFICATION OF KELLY AND RICE METHOD FOR SYNERGISM TESTS

ABSTRACT

In this paper, we present a modified Kelly and Rice method for testing synergism. This approach is consistent with Berenbaum's1-3 Berenbaum, M.C. 1977. Synergism, Additivism, and Antagonism in Immuno-Suppression: A Critical Review. Clinical and Experimental Immunology, 28: 1–18. Berenbaum, M.C. 1985. The Expected Effect of a Combination of Agents: The General Solution. Journal of Theoretical Biology, 114: 413–431. Berenbaum, M.C. 1989. What is Synergism? Pharmacological Reviews. 41: 93–141.   framework for additivity. The delta method[4] Bishop, Y.M., Fienberg, S.E. and Holland, P.W. 1975. Discrete Multivariate Analysis: Theory and Practice Massachusetts: MIT Press.  [Google Scholar] is applied to obtain the estimated variance for the predicted additivity proportion. A Monte Carlo simulation study for the evaluation of the method's performance, i.e., global overall tests for synergism, is also discussed. Kelly and Rice[5] Kelly, C. and Rice, J. 1990. Monotone Smoothing with Application to Dose-Response Curves and the Assessment of Synergism. Biometrics, 46: 1071–1085. [Crossref], [PubMed], [Web of Science ®] , [Google Scholar] do not provide a correct test statistic because the variance is underestimated. Hence, the performance of the Kelly–Rice[5] Kelly, C. and Rice, J. 1990. Monotone Smoothing with Application to Dose-Response Curves and the Assessment of Synergism. Biometrics, 46: 1071–1085. [Crossref], [PubMed], [Web of Science ®] , [Google Scholar] method is generally anti-conservative, based on the simulation findings. In addition, the overall test of synergism with χ²(r) from the modified Kelly and Rice method for larger sample sizes is better than that with χ²(1) from the modified Kelly and Rice method.     

Optimal Halton Sequence via Inversive Scrambling

In this article, we develop the Halton sequence of generating quasi-random numbers to an optimal sequence which has the inversive property. The new constructed quasi-random number generator satisfies the extra uniformity condition on [0, 1]. We finally present the performances of this generator in contrast to the former optimal Halton sequence in Chi et al. (2005 Chi, H., Mascagni, M. and Warnock, T. 2005. On the Optimal Halton Sequences. Math. Comput. in Simul., 70(1): 9–21. [Crossref], [Web of Science ®] , [Google Scholar]) and modified optimal Halton sequence in Fathi et al. (2009 Fathi, B., Samimi, H. and Eskandari, A. 2009. A New Efficient Algorithm for Linear Scrambling Halton Sequence. International Journal of Applied Mathematics, 22(7): 1059–1065.  [Google Scholar]).     

Marginal correlation from an extended random-effects model for repeated and overdispersed counts

Vangeneugden et al. [15 Vangeneugden, T., Molenberghs, G., Laenen, A., Geys, H., Beunckens, C. and Sotto, C. 2007. Marginal correlation in longitudinal binary data based on generalized linear mixed models, Tech. Rep., Hasselt University. submitted for publication [Google Scholar]] derived approximate correlation functions for longitudinal sequences of general data type, Gaussian and non-Gaussian, based on generalized linear mixed-effects models (GLMM). Their focus was on binary sequences, as well as on a combination of binary and Gaussian sequences. Here, we focus on the specific case of repeated count data, important in two respects. First, we employ the model proposed by Molenberghs et al. [13 Molenberghs, G., Verbeke, G. and Demétrio, C. G.B. 2007. An extended random-effects approach to modeling repeated, overdispersed count data. Lifetime Data Anal., 13: 513–531. [Crossref], [PubMed], [Web of Science ®] [Google Scholar]], which generalizes at the same time the Poisson-normal GLMM and the conventional overdispersion models, in particular the negative-binomial model. The model flexibly accommodates data hierarchies, intra-sequence correlation, and overdispersion. Second, means, variances, and joint probabilities can be expressed in closed form, allowing for exact intra-sequence correlation expressions. Next to the general situation, some important special cases such as exchangeable clustered outcomes are considered, producing insightful expressions. The closed-form expressions are contrasted with the generic approximate expressions of Vangeneugden et al. [15 Vangeneugden, T., Molenberghs, G., Laenen, A., Geys, H., Beunckens, C. and Sotto, C. 2007. Marginal correlation in longitudinal binary data based on generalized linear mixed models, Tech. Rep., Hasselt University. submitted for publication [Google Scholar]]. Data from an epileptic-seizures trial are analyzed and correlation functions derived. It is shown that the proposed extension strongly outperforms the classical GLMM.     

Robust Bayesian inference via γ-divergence

Abstract

This paper presents the robust Bayesian inference based on the γ-divergence which is the same divergence as “type 0 divergence” in Jones et al. (2001 Jones, M. C., N. L. Hjort, I. R. Harris, and A. Basu. 2001. A comparison of related density-based minimum divergence estimators. Biometrika 88(3):865–73.[Crossref], [Web of Science ®] , [Google Scholar]) on the basis of Windham (1995 Windham, M. P. 1995. Robustifying model fitting. Journal of the Royal Statistical Society B57:599–609. [Google Scholar]). It is known that the minimum γ-divergence estimator works well to estimate the probability density for heavily contaminated data, and to estimate the variance parameters. In this paper, we propose a robust posterior distribution against outliers based on the γ-divergence and show the asymptotic properties of the proposed estimator. We also discuss some robustness properties of the proposed estimator and illustrate its performances in some simulation studies.     

Improved Bar-Lev,Bobovitch, and Boukai Randomized Response Models

In this article, two new improved randomized response models have been proposed. The proposed models are found to be more efficient than the recent randomized response model studied by Bar-Lev et al. (2004 Bar-Lev , S. K. , Bobovitch , E. , Boukai , B. ( 2004 ). A note on randomized response models for quantitative data . Metrika 60 : 255 – 260 .[Crossref], [Web of Science ®] , [Google Scholar]). The relative efficiency of the proposed models has been studied with respect to the Bar-Lev et al. (2004 Bar-Lev , S. K. , Bobovitch , E. , Boukai , B. ( 2004 ). A note on randomized response models for quantitative data . Metrika 60 : 255 – 260 .[Crossref], [Web of Science ®] , [Google Scholar]) model under different situations.     

An Asymmetric Kernel Estimator of Density Function for Stationary Associated Sequences

Here, we apply the smoothing technique proposed by Chaubey et al. (2007 Chaubey , Y. P. , Sen , A. , Sen , P. K. ( 2007 ). A new smooth density estimator for non-negative random variables. Technical Report No. 1/07. Department of Mathematics and Statistics, Concordia University, Montreal, Canada . [Google Scholar]) for the empirical survival function studied in Bagai and Prakasa Rao (1991 Bagai , I. , Prakasa Rao , B. L. S. ( 1991 ). Estimation of the survival function for stationary associated processes . Statist. Probab. Lett. 12 : 385 – 391 .[Crossref], [Web of Science ®] , [Google Scholar]) for a sequence of stationary non-negative associated random variables.The derivative of this estimator in turn is used to propose a nonparametric density estimator. The asymptotic properties of the resulting estimators are studied and contrasted with some other competing estimators. A simulation study is carried out comparing the recent estimator based on the Poisson weights (Chaubey et al., 2011 Chaubey , Y. P. , Dewan , I. , Li , J. ( 2011 ). Smooth estimation of survival and density functions for a stationary associated process using poisson weights . Statist. Probab. Lett. 81 : 267 – 276 .[Crossref], [Web of Science ®] , [Google Scholar]) showing that the two estimators have comparable finite sample global as well as local behavior.     

On a Generalization of Bivariate Cauchy Distribution

This paper addresses a generalization of the bivariate Cauchy distribution discussed by Fang et al. (1990 Fang , K. T. , Kotz , S. , Ng , K. W. ( 1990 ). Symmetric Multivariate and Related Distributions . London : Chapman and Hall .[Crossref] , [Google Scholar]), derived from a trivariate normal distribution with a general correlation matrix. We obtain explicit expressions for the joint distribution function and joint density function, and show that they reduce in a special case to the corresponding expressions of Fang et al. (1990 Fang , K. T. , Kotz , S. , Ng , K. W. ( 1990 ). Symmetric Multivariate and Related Distributions . London : Chapman and Hall .[Crossref] , [Google Scholar]). Finally, we show that this generalized distribution is useful in determining the orthant probability of a bivariate skew-normal distribution of 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]).