Asymptotic properties for the moment estimators in Rayleigh distribution with two parameters

ABSTRACT

We consider the asymptotic properties for the moment estimators in Rayleigh distribution with two parameters. The law of the iterated logarithm for the estimators can be obtained. Moreover, we can give a simple proof of the asymptotic normality which has been obtained by Li and Li (2012) Li, Y.W., Li, M.H. (2012). Moment estimation of the parameters in Rayleigh distribution with two parameters. Commun. Stat.-Theor. Methods 41:2643–2660.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar].     

On Fitting Transformation Model to Survey Data

In this article, we consider fitting a semiparametric linear model to survey data with censored observations. The specific goal of the paper is to extend the methods of Cheng et al. (1995 Cheng, S.C., Wei, L.J., Ying, Z. (1995). Analysis of transformation models with censored data. Biometrika 82(4):835–845.[Crossref], [Web of Science ®] , [Google Scholar]) and Chen et al. (2002 Chen, K., Jin, Z. Ying, Z. (2002). Semiparametric analysis of transformation models with censored data. Biometrika 89:659–668.[Crossref], [Web of Science ®] , [Google Scholar]) to the case when the sample has been drawn from a population using a complex sampling design. Similar to the approach of Lin (2000 Lin, D.Y. (2000). On fitting Cox’s proportional hazards models to survey data. Biometrika 87:37–47.[Crossref], [Web of Science ®] , [Google Scholar]), we regard the survey population as a random sample from an infinite universe and accounts for this randomness in the statistical inference. A simulation study is conducted to investigate the performance of the proposed estimators.     

Group Variable Selection with Oracle Property by Weight-Fused Adaptive Elastic Net Model for Strongly Correlated Data

This paper is the generalization of weight-fused elastic net (Fu and Xu, 2012 Fu, G., Xu, Q. (2012). Grouping variable selection by weight fused elastic net for multi-collinear data. Communications in Statistics-Simulation and Computation 41(2):205–221.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), which performs group variable selection by combining weight-fused LASSO(wfLasso) and elastic net (Zou and Hastie, 2005 Zou, H., Hastie, T. (2005). Regularization and variable selection via the elastic net. Journal of the Royal Statistical Society: Series B (Statistical Methodology) 67(2):301–320.[Crossref], [Web of Science ®] , [Google Scholar]) penalties. In this study, the elastic net penalty is replaced by adaptive elastic net penalty (AdaEnet) (Zou and Zhang, 2009 Zou, H., Zhang, H. (2009). On the adaptive elastic-net with a diverging number of parameters. Annals of Statistics 37(4):1733–1751.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]), and a new group variable selection algorithm with oracle property (Fan and Li, 2001 Fan, J., Li, R. (2001). Variable selection via nonconcave penalized likelihood and its oracle properties. Journal of the American Statistical Association 96(456):1348–1360.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]; Zou, 2006 Zou, H. (2006). The adaptive lasso and its oracle properties. Journal of the American Statistical Association 101(476):1418–1429.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) is obtained.     

Combined estimating function for random coefficient models with correlated errors

ABSTRACT

Estimating functionshave been shown to be convenient to study inference for non linear time series models. Recently, Thavaneswaran et al. (2012 Thavaneswaran, A., Liang, Y., Frank, J. (2012). Inference for random coefficient volatility models. Stat. Probab. Lett. 82(12):2086–2090.[Crossref], [Web of Science ®] , [Google Scholar]) used combined estimating functions to study inference for random coefficient autoregressive (RCA) models with generalized autoregressive heteroscedasticity errors. While most RCA modeling assumes that the random term and the error are independent, Chandra and Taniguchi (2001 Chandra, S.A., Taniguchi, M. (2001). Estimating functions for nonlinear time series models. Ann. Inst. Stat. Math 53(1):125–141.[Crossref], [Web of Science ®] , [Google Scholar]) studied inference for RCA models with correlated errors using linear estimating functions. In this paper, we derive the quadratic estimating functions for the joint estimation of the conditional mean, variance, and correlation parameters of the RCA models with correlated errors.     

Generalized linear failure rate power series distribution

Abstract

In this article, we introduce a new class of lifetime distributions. This new class includes several previously known distributions such as those of Chahkandi and Ganjali (2009 Chahkandi, M., Ganjali, M. (2009). On some lifetime distributions with decreasing failure rate. Computat. Statist. Data Anal. 53:4433–4440.[Crossref], [Web of Science ®] , [Google Scholar]), Mahmoudi and Jafari (2012 Mahmoudi, E., Jafari, A.A. (2012). Generalized exponential power series distributions. Comput. Statist. Data Anal. 56(12):4047–4066.[Crossref], [Web of Science ®] , [Google Scholar]), and Nadarajah et al. (2012 Nadarajah, S., Shahsanaei, F., Rezaei, S. (2012). A new four-parameter lifetime distribution. J. Statist. Computat. Simul.. ifirst, 1–16. [Google Scholar]). This new class of four-parameter distributions allows for flexible failure rate behavior. Indeed, the failure rate function here can be increasing, decreasing, bathtub-shaped or upside-down bathtub-shaped. Several distributional properties of the new class including moments, quantiles and order statistics are studied. An EM algorithm for computing the estimates of the parameters involved is proposed and some maximum entropy characterizations are discussed. Finally, to show the flexibility and potential of the new class of distributions, applications to two real data sets are provided.     

A Strong Limit Theorem for Weighted Sums of Negatively Dependent Random Variables

In this article, we establish a new complete convergence theorem for weighted sums of negatively dependent random variables. As corollaries, many results on the almost sure convergence and complete convergence for weighted sums of negatively dependent random variables are obtained. In particular, the results of Jing and Liang (2008 Jing, B.Y., Liang, H.Y. (2008). Strong limit theorems for weighted sums of negatively associated random variables. J. Theor. Probab. 21:890–909.[Crossref], [Web of Science ®] , [Google Scholar]), Sung (2012 Sung, S.H. (2012). Complete convergence for weighted sums of negatively dependent random variables. Stat. Pap. 53:73–82.[Crossref], [Web of Science ®] , [Google Scholar]), and Wu (2010) can be obtained.     

Berry-Esseen Bounds for Random Index Non Linear Statistics via Stein's Method

Under the second moment condition, we obtain Berry-Esseen bounds for random index non linear statistics by using a technique discussed in Chen and Shao (2007 Chen, L. H.Y., Shao, Q.-M. (2007). Normal approximation for nonlinear statistics using a concentration inequality approach. Bernoulli 13(2):581–599.[Crossref], [Web of Science ®] , [Google Scholar]). A concept in this article is to approximate any random index non-linear statistic by a random index linear statistic. The bounds for random sums of independent random variables are also provided. Applications are the bounds for random U-statistics and random sums of the present values in investment analysis.     

Generalized linear regression models incorporating original outcome distributions

ABSTRACT

In this article, we propose an approach for incorporating continuous and discrete original outcome distributions into the usual exponential family regression models. The new approach is an extension of the works of Suissa (1991 Suissa, S. (1991). Binary methods for continuous outcomes: A parametric alternative. J. Clin. Epidemiol. 44:241–248.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) and Suissa and Blais (1995 Suissa, S., Blais, L. (1995). Binary regression with continuous outcomes. Stat. Med. 14:247–255.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]), which present methods to estimate the risk of an event defined in a sample subspace of an original continuous outcome variable. Simulation studies are presented in order to illustrate the performance of the developed methodology. Real data sets are analyzed by using the proposed models.     

On Augmented Franses Tests for Seasonal Unit Roots

This article extends the results reported in del Barrio Castro, Osborn and Taylor (2012 del Barrio Castro, T., Osborn, D.R., Taylor, A. M.R. (2012). On augmented HEGY tests for seasonal unit roots. Econometric Theor. 18:1121–1143.[Crossref], [Web of Science ®] , [Google Scholar]) to the approach followed by Franses (1991a Franses, P. H. (1991a). Model selection and seasonality in time series. Tibergen Institute Series, 18. [Google Scholar],b Franses, P.H. (1991b). Seasonality, non-stationarity and the forecasting of monthly time series. Int. J. Forecast. 7:199–208.[Crossref], [Web of Science ®] , [Google Scholar]) to test for seasonal unit roots, providing the asymptotic representation to the seasonal unit roots tests proposed by Franses for a general number of seasons S.     

Designs Robust against Violation of Normality Assumption in the Standard F-test

In this article, designs are found for which the F-test of analysis of variance is insensitive to violation of normality assumption. Atiqullah (1962 Atiquallah, M. (1962). The estimation of residual variance in quadratically balanced least-squares problems and the robustness of the F-test. Biometrika 49: 83–91.[Crossref], [Web of Science ®] , [Google Scholar]) proved that the F-test for treatments adjusting for blocks in the intra-block analysis of a balanced incomplete block design is robust against non-normality in the observations. Here an attempt has been made to identify other designs robust in this sense. In particular, it is observed that for testing relevant hypothesis, a partially balanced incomplete block design in block design setup, under certain conditions, is robust. Robustness of a balanced treatment incomplete block design and a partially balanced treatment incomplete block design (Biswas, 2012 Biswas, A. (2012). Block designs robust against the presence of an aberration in a treatment-control setup. Commun Statist.Theor Methods 41: 920–933.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), in treatment-control design setup, is also studied. Moreover, a new measure of robustness is introduced for further study. The performance of the F-test in presence of non-normality in the observations for a quadratically balanced design is also examined.     

Interval Estimation of the Stress-Strength Reliability with Independent Normal Random Variables

This article develops a procedure to obtain highly accurate confidence interval estimates for the stress-strength reliability R = P(X > Y) where X and Y are data from independent normal distributions of unknown means and variances. Our method is based on third-order likelihood analysis and is compared to the conventional first-order likelihood ratio procedure as well as the approximate methods of Reiser and Guttman (1986 Reiser, B., Guttman, I. (1986). Statistical inference for Pr(Y < X): the normal case. Technometrics 28: 253257.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) and Guo and Krishnamoorthy (2004 Guo, H., Krishnamoorthy, K. (2004). New approximate inferential methods for the reliability parameter in a stress-strength model: the normal case. Commun. Statist. Theor. Meth. 33: 17151731.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]). The use of our proposed method is illustrated by an empirical example and its superior accuracy in terms of coverage probability and error rate are examined through Monte Carlo simulation studies.     

Optimal Simultaneous Confidence Bands in Multiple Linear Regression with Predictor Variables Constrained in an Ellipsoidal Region

A simultaneous confidence band provides useful information on the plausible range of an unknown regression model function, just as a confidence interval gives the plausible range of an unknown parameter. For a multiple linear regression model, confidence bands of different shapes, such as the hyperbolic band and the constant width band, can be constructed and the predictor variable region over which a confidence band is constructed can take various forms. One interesting but unsolved problem is to find the optimal (shape) confidence band over an ellipsoidal region χ_E under the Minimum Volume Confidence Set (MVCS) criterion of Liu and Hayter (2007 Liu, W., Hayter, A.J. (2007). Minimum area confidence set optimality for confidence bands in simple linear regression. J. Amer. Statist. Assoc. 102:181–190.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) and Liu et al. (2009 Liu, W., Bretz, F., Hayter, A.J., Wynn, H.P. (2009). Assessing non-superiority, non-inferiority or equivalence when comparing two regression models over a restricted covariate region. Biometrics 65:1279–1287.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]). This problem is challenging as it involves optimization over an unknown function that determines the shape of the confidence band over χ_E. As a step towards solving this difficult problem, in this paper, we introduce a family of confidence bands over χ_E, called the inner-hyperbolic bands, which includes the hyperbolic and constant-width bands as special cases. We then search for the optimal confidence band within this family under the MVCS criterion. The conclusion from this study is that the hyperbolic band is not optimal even within this family of inner-hyperbolic bands and so cannot be the overall optimal band. On the other hand, the constant width band can be optimal within the family of inner-hyperbolic bands when the region χ_E is small and so might be the overall optimal band.     

A proportional hazards model for the analysis of doubly censored competing risks data

ABSTRACT

Competing risks data are common in medical research in which lifetime of individuals can be classified in terms of causes of failure. In survival or reliability studies, it is common that the patients (objects) are subjected to both left censoring and right censoring, which is refereed as double censoring. The analysis of doubly censored competing risks data in presence of covariates is the objective of this study. We propose a proportional hazards model for the analysis of doubly censored competing risks data, using the hazard rate functions of Gray (1988 Gray, R.J. (1988). A class of k-sample tests for comparing the cumulative incidence of a competing risk. Ann. Statist. 16:1141–1154.[Crossref], [Web of Science ®] , [Google Scholar]), while focusing upon one major cause of failure. We derive estimators for regression parameter vector and cumulative baseline cause specific hazard rate function. Asymptotic properties of the estimators are discussed. A simulation study is conducted to assess the finite sample behavior of the proposed estimators. We illustrate the method using a real life doubly censored competing risks data.     

First- and Second-order Asymptotics for the Tail Distortion Risk Measure of Extreme Risks

The tail distortion risk measure at level p was first introduced in Zhu and Li (2012 Zhu, L., Li, H. (2012). Tail distortion risk and its asymptotic analysis. Insur. Math. Econ. 51(1):115–121.[Crossref], [Web of Science ®] , [Google Scholar]), where the parameter p ∈ (0, 1) indicates the confidence level. They established first-order asymptotics for this risk measure, as p↑1, for the Fréchet case. In this article, we extend their work by establishing both first-order and second-order asymptotics for the Fréchet, Weibull, and Gumbel cases. Numerical studies are also carried out to examine the accuracy of both asymptotics.     

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.     

Inference in log-alpha-power and log-skew-normal multivariate models

ABSTRACT

Random vectors with positive components are common in many applied fields, for example, in meteorology, when daily precipitation is measured through a region Marchenko and Genton (2010 Marchenko, Y., Genton, M. (2010). Multivariate log-skew-elliptical distributions with applications to precipitation data. Environmetrics 21:318–340.[Crossref], [Web of Science ®] , [Google Scholar]). Frequently, the log-normal multivariate distribution is used for modeling this type of data. This modeling approach is not appropriate for data with high asymmetry or kurtosis. Consequently, more flexible multivariate distributions than the log-normal multivariate are required. As an alternative to this distribution, we propose the log-alpha-power multivariate and log-skew-normal multivariate models. The first model is an extension for positive data of the fractional order statistics model Durrans (1992 Durrans, S. (1992). Distributions of fractional order statistics in hydrology. Water Resour. Res. 28:1649–1655.[Crossref], [Web of Science ®] , [Google Scholar]). The second one is an extension of the log-skew-normal model studied by Mateu-Figueras and Pawlowsky-Glahn (2007 Mateu-Figueras, G., Pawlowsky-Glahn, V. (2007). The skew-normal distribution on the simplex. Commun. Stat.-Theory Methods 36:1787–1802.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]). We study parameter estimation for these models by means of pseudo-likelihood and maximum likelihood methods. We illustrate the proposal analyzing a real dataset.     

Baker- Lin-Huang Type Bivariate Distributions Based on Order Statistics

Baker (2008 Baker, R. (2008). An order-statistics-based method for constructing multivariate distributions with fixed marginals. J. Multivariate Anal. 99: 23122327.[Crossref], [Web of Science ®] , [Google Scholar]) introduced a new class of bivariate distributions based on distributions of order statistics from two independent samples of size n. Lin and Huang (2010 Lin, G.D., Huang, J.S. (2010). A note on the maximum correlation for Baker’s bivariate distributions with fixed marginals. J. Multivariate Anal. 101: 22272233.[Crossref], [Web of Science ®] , [Google Scholar]) discovered an important property of Baker’s distribution and showed that the Pearson’s correlation coefficient for this distribution converges to maximum attainable value, i.e., the correlation coefficient of the Fréchet upper bound, as n increases to infinity. Bairamov and Bayramoglu (2013 Bairamov, I., Bayramoglu, K. (2013). From Huang-Kotz distribution to Baker’s distribution. J. Multivariate Anal. 113: 106115.[Crossref], [Web of Science ®] , [Google Scholar]) investigated a new class of bivariate distributions constructed by using Baker’s model and distributions of order statistics from dependent random variables, allowing higher correlation than that of Baker’s distribution. In this article, a new class of Baker’s type bivariate distributions with high correlation are constructed based on distributions of order statistics by using an arbitrary continuous copula instead of the product copula.     

Alternative Approximations to Value-At-Risk: A Comparison

This article compares three value-at-risk (VaR) approximation methods suggested in the literature: Cornish and Fisher (1937 Cornish, E.A., Fisher, R.A. (1937). Moments and cumulants in the specification of distributions. Revue de l’Institut International de Statistique 5:307–320.[Crossref] , [Google Scholar]), Sillitto (1969 Sillitto, G.P. (1969). Derivation of approximants to the inverse distribution function of a continuous univariate population from the order statistics of a sample. Biometrika 56:641–650.[Crossref], [Web of Science ®] , [Google Scholar]), and Liu (2010 Liu, W.-H. (2010). Estimation and testing of portfolio value-at-risk based on L-comoment matrices. Journal of Futures Markets 30:897–908.[Crossref], [Web of Science ®] , [Google Scholar]). Simulation results are obtained for three families of distributions: student-t, skewed-normal, and skewed-t. We recommend the Sillitto approximation as the best method to evaluate the VaR when the financial return has an unknown, skewed, and heavy-tailed distribution.     

A Universal Prior Distribution for Bayesian Consistency of Non parametric Procedures

The introduction of the Hausdorff α-entropy in Xing (2008a Xing, Y. (2008a). Convergence rates of posterior distributions for observations without the iid structure, 38 pages. Available at: www.arxiv.org:0811.4677v1. [Google Scholar]), Xing (2008b Xing, Y. (2008b). On adaptive Bayesian inference. Electron. J. Stat. 2:848–862.[Crossref] , [Google Scholar]), Xing (2010 Xing, Y. (2010). Rates of posterior convergence for iid Observations. Commun. Stat. Theory Methods. 39(19):3389–3398.[Taylor & Francis Online] , [Google Scholar]), Xing (2011 Xing, Y. (2011). Convergence rates of nonparametric posterior distributions. J. Stat. Plann. Inference 141:3382–3390.[Crossref], [Web of Science ®] , [Google Scholar]), and Xing and Ranneby (2009 Xing, Y., Ranneby, B. (2009). Sufficient conditions for Bayesian consistency. J. Stat. Plann. Inference. 139:2479–2489.[Crossref], [Web of Science ®] , [Google Scholar]) has lead a series of improvements of well-known results on posterior consistency. In this paper we discuss an application of the Hausdorff α-entropy. We construct a universal prior distribution such that the corresponding posterior distribution is almost surely consistent. The approach of the construction of this type of prior distribution is natural, but it works very well for all separable models. We illustrate such prior distributions by examples. In particular, we obtain that if the true density function is known to be some normal probability density function with unknown mean and unknown variance then without any additional assumption one can construct a prior distribution which leads to posterior consistency.     

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.