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.     

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.     

Note on posterior inference for the Bingham distribution

The properties of high-dimensional Bingham distributions have been studied by Kume and Walker (2014 Kume, A., and S. G. Walker. 2014. On the Bingham distribution with large dimension. Journal of Multivariate Analysis 124:345–52.[Crossref], [Web of Science ®] , [Google Scholar]). Fallaize and Kypraios (2016 Fallaize, C. J., and T. Kypraios. 2016. Exact Bayesian inference for the Bingham distribution. Statistics and Computing 26:349–60.[Crossref], [Web of Science ®] , [Google Scholar]) propose the Bayesian inference for the Bingham distribution and they use developments in Bayesian computation for distributions with doubly intractable normalizing constants (Møller et al. 2006 Møller, J., A. N. Pettitt, R. Reeves, and K. K. Berthelsen. 2006. An efficient Markov chain Monte Carlo method for distributions with intractable normalising constants. Biometrika 93 (2):451–458.[Crossref], [Web of Science ®] , [Google Scholar]; Murray, Ghahramani, and MacKay 2006 Murray, I., Z. Ghahramani, and D. J. C. MacKay. 2006. MCMC for doubly intractable distributions. In Proceedings of the 22nd annual conference on uncertainty in artificial intelligence (UAI-06), 359–66. AUAI Press. [Google Scholar]). However, they rely heavily on two Metropolis updates that they need to tune. In this article, we propose instead a model selection with the marginal likelihood.     

A quantile-based study on ageing intensity function

Jiang, Ji, and Xiao (2003 Jiang, R., P. Ji, and X. Xiao. 2003. Aging property of unimodal failure rate models. Reliability Engineering and System Safety 79(1):113–6.[Crossref], [Web of Science ®] , [Google Scholar]) has introduced a quantitative measure known as the ageing intensity function for evaluating the ageing properties of a component/system. In recent years, there has been a great interest on the study of quantile function, an equivalent alternative to the distribution function approach. Unlike the distribution function approach, the quantile method possess some unique properties (see Gilchrist 2000 Gilchrist, W. 2000. Statistical modelling with quantile functions. Boca Raton, Florida: Chapman and Hall/CRC.[Crossref] , [Google Scholar], Nair, Sankaran, and Balakrishnan 2013 Nair N. U., P. G. Sankaran, N. Balakrishnan. 2013. Quantile-based reliability concepts. In: Quantile-Based Reliability Analysis. Statistics for Industry and Technology. New York, NY: Birkhäuser.[Crossref] , [Google Scholar]). Motivated with this, in the present paper we introduce a quantile-based ageing intensity function and study its various ageing properties. We also study some stochastic comparison of random variables based on the proposed measure.     

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

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.     

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

ESTIMATION OF PARAMETERS IN ARFIMA PROCESSES: A SIMULATION STUDY

It is known that, in the presence of short memory components, the estimation of the fractional parameter d in an Autoregressive Fractionally Integrated Moving Average, ARFIMA(p, d, q), process has some difficulties (see [1] Smith, J., Taylor, N. and Yadav, S. 1997. Comparing the bias and misspecification in ARFIMA models. Journal of Time Series Analysis, 18(5): 507–527. [Crossref] , [Google Scholar]). In this paper, we continue the efforts made by Smith et al. [1] Smith, J., Taylor, N. and Yadav, S. 1997. Comparing the bias and misspecification in ARFIMA models. Journal of Time Series Analysis, 18(5): 507–527. [Crossref] , [Google Scholar] and Beveridge and Oickle [2] Beveridge, S. and Oickle, C. 1993. Estimating fractionally integrated time series models. Economics Letters, 43: 137–142.  [Google Scholar] by conducting a simulation study to evaluate the convergence properties of the iterative estimation procedure suggested by Hosking [3] Hosking, J. 1981. Fractional differencing. Biometrika, 68(1): 165–176. [Crossref], [Web of Science ®] , [Google Scholar]. In this context we consider some semiparametric approaches and a parametric method proposed by Fox-Taqqu[4] Fox, R. and Taqqu, M. S. 1986. Large-sample properties of parameter estimates for strongly dependent stationary gaussian time series. The Annals of Statistics, 14(2): 517–532. [Crossref], [Web of Science ®] , [Google Scholar]. We also investigate the method proposed by Robinson [5] Robinson, P. M. 1995a. Log-periodogram regression of time series with long range dependence. The Annals of Statistics, 23(3): 1048–1072. [Crossref], [Web of Science ®] , [Google Scholar] and a modification using the smoothed periodogram function.     

Bayesian estimation of sensitivity level and population proportion of a sensitive characteristic in a binary optional unrelated question RRT model

Sihm et al. (2016 Sihm, J. S., A. Chhabra, and S. N. Gupta. 2016. An optional unrelated question RRT model. Involve: A Journal of Mathematics 9 (2):195–209.[Crossref] , [Google Scholar]) proposed an unrelated question binary optional randomized response technique (RRT) model for estimating the proportion of population that possess a sensitive characteristic and the sensitivity level of the question. In our work, decision theoretic approach has been followed to obtain Bayes estimates of the two parameters along with their corresponding minimal Bayes posterior expected losses (BPEL) using beta prior and squared error loss function (SELF). Relative losses are also examined to compare the performances of the Bayes estimates with those of the classical estimates obtained by Sihm et al. (2016 Sihm, J. S., A. Chhabra, and S. N. Gupta. 2016. An optional unrelated question RRT model. Involve: A Journal of Mathematics 9 (2):195–209.[Crossref] , [Google Scholar]). The results obtained are illustrated with the help of real survey data using non informative prior.     

On generalized exponential chain ratio estimators under stratified two-phase random sampling

Sanaullah et al. (2014 Sanaullah, A., Ali, H.M., Noor ul Amin, M., Hanif, M. (2014). Generalized exponential chain ratio estimators under stratified two-phase random sampling. Appl. Math. Comput. 226:541–547.[Crossref], [Web of Science ®] , [Google Scholar]) have suggested generalized exponential chain ratio estimators under stratified two-phase sampling scheme for estimating the finite population mean. However, the bias and mean square error (MSE) expressions presented in that work need some corrections, and consequently the study based on efficiency comparison also requires corrections. In this article, we revisit Sanaullah et al. (2014 Sanaullah, A., Ali, H.M., Noor ul Amin, M., Hanif, M. (2014). Generalized exponential chain ratio estimators under stratified two-phase random sampling. Appl. Math. Comput. 226:541–547.[Crossref], [Web of Science ®] , [Google Scholar]) estimator and provide the correct bias and MSE expressions of their estimator. We also propose an estimator which is more efficient than several competing estimators including the classes of estimators in Sanaullah et al. (2014 Sanaullah, A., Ali, H.M., Noor ul Amin, M., Hanif, M. (2014). Generalized exponential chain ratio estimators under stratified two-phase random sampling. Appl. Math. Comput. 226:541–547.[Crossref], [Web of Science ®] , [Google Scholar]). Three real datasets are used for efficiency comparisons.     

Quantile-based Chernoff distance for truncated random variables

Several probability distributions such as power-Pareto distribution (see Gilchrist 2000 Gilchrist, W. 2000. Statistical modelling with quantile functions. Boca Raton, FL: Chapman and Hall/CRC.[Crossref] , [Google Scholar] and Hankin and Lee 2006 Hankin, R. K. S., and A. Lee. 2006. A new family of non-negative distributions. Australian and New Zealand Journal of Statistics 48:67–78.[Crossref], [Web of Science ®] , [Google Scholar]), various forms of lambda distributions (see Ramberg and Schmeiser 1974 Ramberg, J. S., and B. W. Schmeiser. 1974. An appropriate method for generating asymmetric random variables. Communications of the ACM 17:78–82.[Crossref], [Web of Science ®] , [Google Scholar] and Freimer et al. 1988 Freimer, M., S. Mudholkar, G. Kollia, and C. T. Lin. 1988. A study of the generalized lambda family. Communications in Statistics - Theory and Methods 17:3547–67.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), Govindarajulu distribution (see Nair, Sankaran, and Vineshkumar 2012 Nair, U. N., P. G. Sankaran, and B. Vineshkumar. 2012. The Govindarajulu distribution: some properties and applications. Communications in Statistics—Theory and Methods 41:4391–406.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), etc., do not have manageable distribution functions, though they have tractable quantile functions. Hence, analytical study of the properties of Chernoff distance of two random variables associated with these distributions via traditional distribution function-based tool becomes difficult. To make this simple, in this paper, we introduce quantile-based Chernoff distance for (left or right) truncated random variables and study its various properties. Some useful bounds as well as characterization results are obtained.     

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.     

A Bottom-Up Dynamic Model of Portfolio Credit Risk with Stochastic Intensities and Random Recoveries

In Bielecki et al. (2014a Bielecki , T. R. , Cousin , A. , Crépey , S. , Herbertsson , A. ( 2014a ). Dynamic hedging of portfolio credit risk in a markov copula model . J. Optimiz. Theor. Applic . doi: DOI 10.1007/s10957-013-0318-4 (forthcoming) .[Crossref] , [Google Scholar]), the authors introduced a Markov copula model of portfolio credit risk where pricing and hedging can be done in a sound theoretical and practical way. Further theoretical backgrounds and practical details are developed in Bielecki et al. (2014b Bielecki , T. R. , Cousin , A. , Crépey , S. , Herbertsson , A. ( 2014b ). A bottom-up dynamic model of portfolio credit risk - Part I: Markov copula perspective . In: Recent Adv. Fin. Eng. 2012 , World Scientific (preprint version available at http://dx.doi.org/10.2139/ssrn.1844574) . [Google Scholar],c) where numerical illustrations assumed deterministic intensities and constant recoveries. In the present paper, we show how to incorporate stochastic default intensities and random recoveries in the bottom-up modeling framework of Bielecki et al. (2014a Bielecki , T. R. , Cousin , A. , Crépey , S. , Herbertsson , A. ( 2014a ). Dynamic hedging of portfolio credit risk in a markov copula model . J. Optimiz. Theor. Applic . doi: DOI 10.1007/s10957-013-0318-4 (forthcoming) .[Crossref] , [Google Scholar]) while preserving numerical tractability. These two features are of primary importance for applications like CVA computations on credit derivatives (Assefa et al., 2011 Assefa , S. , Bielecki , T. R. , Crépey , S. , Jeanblanc , M. ( 2011 ). CVA computation for counterparty risk assessment in credit portfolios . In: Bielecki , T.R. , Brigo , D. , Patras , F. , Eds., Credit Risk Frontiers . Hoboken : Wiley/Bloomberg-Press . [Google Scholar]; Bielecki et al., 2012 Bielecki , T. R. , Crépey , S. , Jeanblanc , M. , Zargari , B. ( 2012 ). Valuation and Hedging of CDS counterparty exposure in a markov copula model . Int. J. Theoret. Appl. Fin. 15 ( 1 ): 1250004 .[Crossref] , [Google Scholar]), as CVA is sensitive to the stochastic nature of credit spreads and random recoveries allow to achieve satisfactory calibration even for “badly behaved” data sets. This article is thus a complement to Bielecki et al. (2014a Bielecki , T. R. , Cousin , A. , Crépey , S. , Herbertsson , A. ( 2014a ). Dynamic hedging of portfolio credit risk in a markov copula model . J. Optimiz. Theor. Applic . doi: DOI 10.1007/s10957-013-0318-4 (forthcoming) .[Crossref] , [Google Scholar]), Bielecki et al. (2014b Bielecki , T. R. , Cousin , A. , Crépey , S. , Herbertsson , A. ( 2014b ). A bottom-up dynamic model of portfolio credit risk - Part I: Markov copula perspective . In: Recent Adv. Fin. Eng. 2012 , World Scientific (preprint version available at http://dx.doi.org/10.2139/ssrn.1844574) . [Google Scholar]) and Bielecki et al. (2014c Bielecki , T. R. , Cousin , A. , Crépey , S. , Herbertsson , A. ( 2014c ). A bottom-up dynamic model of portfolio credit risk - Part II: Common-shock interpretation, calibration and hedging issues . Recent Adv. Fin. Eng. 2012 , World Scientific (preprint version available at http://dx.doi.org/10.2139/ssrn.2245130) . [Google Scholar]).     

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.     

MULTI-LEVEL FACTORIAL DESIGNS WITH MINIMUM NUMBERS OF LEVEL CHANGES

The order of experimental runs in a fractional factorial experiment is essential when the cost of level changes in factors is considered. The generalized foldover scheme given by [1] Coster, D. C. and Cheng, C. S. 1988. Minimum cost trend free run orders of fractional factorial designs. The Annals of Statistics, 16: 1188–1205. [Crossref], [Web of Science ®] , [Google Scholar]gives an optimal order to experimental runs in an experiment with specified defining contrasts. An experiment can be specified by a design requirement such as resolution or estimation of some interactions. To meet such a requirement, we can find several sets of defining contrasts. Applying the generalized foldover scheme to these sets of defining contrasts, we obtain designs with different numbers of level changes and then the design with minimum number of level changes. The difficulty is to find all the sets of defining contrasts. An alternative approach is investigated by [2] Cheng, C. S., Martin, R. J. and Tang, B. 1998. Two-level factorial designs with extreme numbers of level changes. The Annals of Statistics, 26: 1522–1539. [Crossref], [Web of Science ®] , [Google Scholar]for two-level fractional factorial experiments. In this paper, we investigate experiments with all factors in slevels.     

An adaptive multivariate CUSUM control chart for signaling a range of location shifts

ABSTRACT

In this work, we proposed an adaptive multivariate cumulative sum (CUSUM) statistical process control chart for signaling a range of location shifts. This method was based on the multivariate CUSUM control chart proposed by Pignatiello and Runger (1990 Pignatiello, J.J., Runger, G.C. (1990). Comparisons of multivariate CUSUM charts. J. Qual. Technol. 22(3):173–186.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), but we adopted the adaptive approach similar to that discussed by Dai et al. (2011 Dai, Y., Luo, Y., Li, Z., Wang, Z. (2011). A new adaptive CUSUM control chart for detecting the multivariate process mean. Qual. Reliab. Eng. Int. 27(7):877–884.[Crossref], [Web of Science ®] , [Google Scholar]), which was based on a different CUSUM method introduced by Crosier (1988 Crosier, R.B. (1988). Multivariate generalizations of cumulative sum quality-control schemes. Technometrics 30(3):291–303.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]). The reference value in this proposed procedure was changed adaptively in each run, with the current mean shift estimated by exponentially weighted moving average (EWMA) statistic. By specifying the minimal magnitude of the mean shift, our proposed control chart achieved a good overall performance for detecting a range of shifts rather than a single value. We compared our adaptive multivariate CUSUM method with that of Dai et al. (2001 Dai, Y., Luo, Y., Li, Z., Wang, Z. (2011). A new adaptive CUSUM control chart for detecting the multivariate process mean. Qual. Reliab. Eng. Int. 27(7):877–884.[Crossref], [Web of Science ®] , [Google Scholar]) and the non adaptive versions of these two methods, by evaluating both the steady state and zero state average run length (ARL) values. The detection efficiency of our method showed improvements over the comparative methods when the location shift is unknown but falls within an expected range.     

An improved estimation in stratified random sampling

ABSTRACT

The article suggests a class of estimators of population mean in stratified random sampling using auxiliary information with its properties. In addition, various known estimators/classes of estimators are identified as members of the suggested class. It has been shown that the suggested class of estimators under optimum condition performs better than the usual unbiased, usual combined ratio, usual combined regression, Kadilar and Cingi (2005 Kadilar, C., Cingi, H. (2005). A new ratio estimator in stratified sampling. Commun. Stat. Theory Methods 34:597–602.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), Singh and Vishwakarma (2006 Singh, H.P., Vishwakarma, G.K. (2006). Combined ratio-product estimator of finite population mean in stratified sampling. Metodologia de Encuestas Monografico: Incidencias en el trabjo de Campo 7(1):32–40. [Google Scholar]) estimators and the members belonging to the classes of estimators envisaged by Kadilar and Cingi (2003 Kadilar, C., Cingi, H. (2003). Ratio estimator in stratified sampling. Biomet. J. 45:218–225.[Crossref], [Web of Science ®] , [Google Scholar]), Singh, Tailor et al. (2008 Singh, H.P., Agnihotri, N. (2008). A general procedure of estimating population mean using auxiliary information in sample surveys. Stat. Trans. 9(1):71–87. [Google Scholar]), Singh et al. (2009 Singh, R., Kumar, M., Chaudhary, M.K., Kadilar, C. (2009). Improved exponential estimator in stratified random sampling. Pak. J. Stat. Oper. Res. 5(2):67–82.[Crossref] , [Google Scholar]), Singh and Vishwakarma (2010 Singh, H.P., Vishwakarma, G.K. (2010). A general procedure for estimating the population mean in stratified sampling using auxiliary information. METRON 67(1):47–65.[Crossref] , [Google Scholar]) and Koyuncu and Kadilar (2010) Koyuncu, N., Kadilar, C. (2010). On improvement in estimating population mean in stratified random sampling. J. Appl. Stat. 37(6):999–1013.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar].     

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.     

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.     

Weighting variables in Kohonen competitive learning algorithms

This paper presents a new variable weight method, called the singular value decomposition (SVD) approach, for Kohonen competitive learning (KCL) algorithms based on the concept of Varshavsky et al. [18 R. Varshavsky, A. Gottlieb, M. Linial, and D. Horn, Novel unsupervised feature filtering of bilogical data, Bioinformatics 22 (2006), pp. 507–513.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]]. Integrating the weighted fuzzy c-means (FCM) algorithm with KCL, in this paper, we propose a weighted fuzzy KCL (WFKCL) algorithm. The goal of the proposed WFKCL algorithm is to reduce the clustering error rate when data contain some noise variables. Compared with the k-means, FCM and KCL with existing variable-weight methods, the proposed WFKCL algorithm with the proposed SVD's weight method provides a better clustering performance based on the error rate criterion. Furthermore, the complexity of the proposed SVD's approach is less than Pal et al. [17 S.K. Pal, R.K. De, and J. Basak, Unsupervised feature evaluation: a neuro-fuzzy approach, IEEE. Trans. Neural Netw. 11 (2000), pp. 366–376.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]], Wang et al. [19 X.Z. Wang, Y.D. Wang, and L.J. Wang, Improving fuzzy c-means clustering based on feature-weight learning, Pattern Recognit. Lett. 25 (2004), pp. 1123–1132.[Crossref], [Web of Science ®] , [Google Scholar]] and Hung et al. [9 W. -L. Hung, M. -S. Yang, and D. -H. Chen, Bootstrapping approach to feature-weight selection in fuzzy c-means algorithms with an application in color image segmentation, Pattern Recognit. Lett. 29 (2008), pp. 1317–1325.[Crossref], [Web of Science ®] , [Google Scholar]].