THE GENERALIZED SECANT HYPERBOLIC DISTRIBUTION AND ITS PROPERTIES

ABSTRACT

The properties of a family of distributions generalizing the secant hyperbolic are developed. This family consists of symmetric distributions, with kurtosis ranging from 1.8 to infinity, and includes the logistic as a special case, the uniform as a limiting case, and closely approximates the normal and Student's t-distributions with corresponding kurtosis. A significant difference between this family and Student's t is that for any member of the generalized secant hyperbolic family, all moments are finite. Further, technical difficulties associated with evaluating moments of Student's t (especially for fractional degrees of freedom) are not present with this family. The properties of the maximum likelihood and modified maximum likelihood estimates of the location and scale parameters for complete samples are considered. Examples illustrate the methods developed in this work.     

Tail-Adaptive Location Rank Test for the Generalized Secant Hyperbolic Distribution

The generalized secant hyperbolic distribution (GSHD) was recently introduced as a modeling tool in data analysis. The GSHD is a unimodal distribution that is completely specified by location, scale, and shape parameters. It has also been shown elsewhere that the rank procedures of location are regular, robust, and asymptotically fully efficient. In this article, we study certain tail weight measures for the GSHD and introduce a tail-adaptive rank procedure of location based on those tail weight measures. We investigate the properties of the new adaptive rank procedure and compare it to some conventional estimators.     

ODD MAN OUT—THE MEIXNER HYPERGEOMETRIC DISTRIBUTION

Two open problems are described for the hyperbolic secant distribution, as a special case of the more general Meixner hypergeomet-ric distribution. The first concerns the completeness of the family of functions sech θ x , θ≥ 0. The second concerns the characterization of the class of canonical correlation sequences in bivariate distributions with marginals sech θ x , sech θ y. In both cases some partial results are put forward.     

GSH Dependence Modeling with an Application to Risk Management

The generalized secant hyperbolic distribution (GSH) can be used to represent financial data with heavy tails as an alternative to the Student-t, because it guarantees the existence of all moments, also with a high kurtosis value. In order to obtain a multivariate extension of the GSH distribution, in this article we present two approaches to model the dependence, the copula approach and independent component analysis. Since the methodologies considered allow to simulate the GSH dependence, we show also the empirical results obtained in the estimation of risk of a financial portfolio by the Monte Carlo method.     

Asymmetric generalizations of symmetric univariate probability distributions obtained through quantile splicing

Abstract

Balakrishnan et al. proposed a two-piece skew logistic distribution by making use of the cumulative distribution function (CDF) of half distributions as the building block, to give rise to an asymmetric family of two-piece distributions, through the inclusion of a single shape parameter. This paper proposes the construction of asymmetric families of two-piece distributions by making use of quantile functions of symmetric distributions as building blocks. This proposition will enable the derivation of a general formula for the L-moments of two-piece distributions. Examples will be presented, where the logistic, normal, Student’s t(2) and hyperbolic secant distributions are considered.     

Choosing the link function and accounting for link uncertainty in generalized linear models using Bayes factors

One important component of model selection using generalized linear models (GLM) is the choice of a link function. We propose using approximate Bayes factors to assess the improvement in fit over a GLM with canonical link when a parametric link family is used. The approximate Bayes factors are calculated using the Laplace approximations given in [32], together with a reference set of prior distributions. This methodology can be used to differentiate between different parametric link families, as well as allowing one to jointly select the link family and the independent variables. This involves comparing nonnested models and so standard significance tests cannot be used. The approach also accounts explicitly for uncertainty about the link function. The methods are illustrated using parametric link families studied in [12] for two data sets involving binomial responses. The first author was supported by Sonderforschungsbereich 386 Statistische Analyse Diskreter Strukturen, and the second author by NIH Grant 1R01CA094212-01 and ONR Grant N00014-01-10745.     

Hodges-Lehmann Scale Estimator for Cauchy Distribution

We draw here on the relation between the Cauchy and hyperbolic secant distributions to prove that the MLE of the scale parameter of the Cauchy distribution is log-normally distributed and to study the properties of a Hodges-Lehmann type estimator for the scale parameter. This scale estimator is slightly biased but performs well even on small samples regardless of the location parameter. The asymptotic efficiency of the estimator is 98%.     

On the matrix-variate generalized hyperbolic distribution and its Bayesian applications

In the first part of the paper, we introduce the matrix-variate generalized hyperbolic distribution by mixing the matrix normal distribution with the matrix generalized inverse Gaussian density. The p-dimensional generalized hyperbolic distribution of [Barndorff-Nielsen, O. (1978). Hyperbolic distributions and distributions on hyperbolae. Scand. J. Stat., 5, 151–157], the matrix-T distribution and many well-known distributions are shown to be special cases of the new distribution. Some properties of the distribution are also studied. The second part of the paper deals with the application of the distribution in the Bayesian analysis of the normal multivariate linear model.     

A new alternative quantile regression model for the bounded response with educational measurements applications of OECD countries

This article introduces a new distribution with two tuning parameters specified on the unit interval. It follows from a ‘hyperbolic secant transformation’ of a random variable following the Weibull distribution. The lack of research on the prospect of hyperbolic transformations providing flexible distributions over the unit interval is a motivation for the study. The main distributional structural properties of the new distribution are established. The different estimation methods and two simulation works have been derived for model parameters. Subsequently, we develop a related quantile regression model for further statistical perspectives. We consider two real data applications based on the educational measurements of both OECD and some non-members of OECD countries. Our regression model aims to relate the desire to get top grades on certain young students in the OECD countries with some of their Education and School Life Index such as reading performance, work environment at home, and paid work experience. It is shown that the elaborated quantile regression model has a better fitting power than famous regression models when the unit response variable possesses skewed distribution as well as two independent variables are significant in the statistical sense at any standard significance level for the median response.     

Sampling nested Archimedean copulas

We give algorithms for sampling from non-exchangeable Archimedean copulas created by the nesting of Archimedean copula generators, where in the most general algorithm the generators may be nested to an arbitrary depth. These algorithms are based on mixture representations of these copulas using Laplace transforms. While in principle the approach applies to all nested Archimedean copulas, in practice the approach is restricted to certain cases where we are able to sample distributions with given Laplace transforms. Precise instructions are given for the case when all generators are taken from the Gumbel parametric family or the Clayton family; the Gumbel case in particular proves very easy to simulate.     

On the property of multivariate generalized hyperbolic distribution and the Stein-type inequality

This note consists of two parts . In the first part, we provide a pedagogic review on the multivariate generalized hyperbolic (MGH) distribution. We show that this probability family is close under margining, conditioning, and linear transforms; however, such property does not hold for its subclasses. In the second part, we obtain the Stein-type inequality in the context of MGH distribution. Moreover, we apply the Stein-type inequality to prove a lower bound for Var[h(X)]. Particularly, we present examples when X belongs to some well-known subclasses in MGH family.     

Heavy or semi-heavy tail,that is the question

While there has been considerable research on the analysis of extreme values and outliers by using heavy-tailed distributions, little is known about the semi-heavy-tailed behaviors of data when there are a few suspicious outliers. To address the situation where data are skewed possessing semi-heavy tails, we introduce two new skewed distribution families of the hyperbolic secant with exciting properties. We extend the semi-heavy-tailedness property of data to a linear regression model. In particular, we investigate the asymptotic properties of the ML estimators of the regression parameters when the error term has a semi-heavy-tailed distribution. We conduct simulation studies comparing the ML estimators of the regression parameters under various assumptions for the distribution of the error term. We also provide three real examples to show the priority of the semi-heavy-tailedness of the error term comparing to heavy-tailedness. Online supplementary materials for this article are available. All the new proposed models in this work are implemented by the shs R package, which can be found on the GitHub webpage.     

On an exponential family

This paper considers a generalization of the exponential type distributions in the class of exponential families. A characterization and a method of generating an exponential family from a given family are given. In particular the generalized gamma, the generalized Poisson, the inverse Gaussian distributions belonging to this family are discussed. The approximations of the cumulative sums for the generalized gamma and the generalized Poisson by the Chi-square are considered. Some of the results are extended to the bivariate case.     

A multiple stream generator based on de Bruijn digraph homomorphisms

We propose a method to obtain several streams of pseudorandom numbers based on a backbone generator of the generalized shift register type. The method is based on inverting one cycle in a de Bruijn digraph into many sequences in a higher-order de Bruijn graph via an appropriate graph homomorphism. We apply this technique to twisted generalized feedback shift register generators and to the Mersenne Twister MT19937. Positive results of statistical testing are reported.     

Robust finite mixture modeling of multivariate unrestricted skew-normal generalized hyperbolic distributions

In this paper, we introduce an unrestricted skew-normal generalized hyperbolic (SUNGH) distribution for use in finite mixture modeling or clustering problems. The SUNGH is a broad class of flexible distributions that includes various other well-known asymmetric and symmetric families such as the scale mixtures of skew-normal, the skew-normal generalized hyperbolic and its corresponding symmetric versions. The class of distributions provides a much needed unified framework where the choice of the best fitting distribution can proceed quite naturally through either parameter estimation or by placing constraints on specific parameters and assessing through model choice criteria. The class has several desirable properties, including an analytically tractable density and ease of computation for simulation and estimation of parameters. We illustrate the flexibility of the proposed class of distributions in a mixture modeling context using a Bayesian framework and assess the performance using simulated and real data.

     

Regression Properties for Asymmetric Generalized Scale Mixtures of Multivariate Gaussian Variables

Analytical properties of regression and the variance–covariance matrix of asymmetric generalized scale mixture of multivariate Gaussian variables are presented. The analysis includes an in-depth analytical investigation of the first two conditional moments of the mixing variable. Exact computable expressions for the prediction and the conditional variance are presented for the generalized hyperbolic distribution using the inversion theorem for Fourier transforms. An application to financial log returns is demonstrated via the classical Euler approximation. The methodology is illustrated by analyzing the regression of intraday log returns for CISCO against the corresponding data from S&P 500.     

上海银行间同业拆放利率ES风险度量研究

随着基准利率地位的不断变化,上海银行间同行业拆放利率（SHIBOR）市场风险管理对金融机构将会越来越重要。然而同正态分布相比而言,SHIBOR收益率变量具有偏态等特征。提出采用广义双曲线分布来拟合收益率序列。为了解决参数估计难的问题,提出利用强有力的EM算法对于解决像包含Bessel函数这样复杂、具有大量局部最优解的优化问题,具有很现实的意义,同时利用蒙特卡罗模拟方法来计算广义双曲线分布下的VaR值、ES值,最后讨论广义双曲线分布在SHIBOR市场风险度量中的应用。     

Random number generators for the generalized Birnbaum–Saunders distribution

The generalized Birnbaum–Saunders distribution pertains to a class of lifetime models including both lighter and heavier tailed distributions. This model adapts well to lifetime data, even when outliers exist, and has other good theoretical properties and application perspectives. However, statistical inference tools may not exist in closed form for this model. Hence, simulation and numerical studies are needed, which require a random number generator. Three different ways to generate observations from this model are considered here. These generators are compared by utilizing a goodness-of-fit procedure as well as their effectiveness in predicting the true parameter values by using Monte Carlo simulations. This goodness-of-fit procedure may also be used as an estimation method. The quality of this estimation method is studied here. Finally, through a real data set, the generalized and classical Birnbaum–Saunders models are compared by using this estimation method.     

Generalized elliptical distributions

A family of distributions generated by an operator acting on generalized normal density is introduced. This family contains as particular cases many known distributions, including the generalized normal, generalized t, and generalized gamma distributions. Several mathematical properties of the family (including expansions, characteristic function, moments, cumulants, and order statistics properties) are derived. Estimation procedures are derived too by the method of moments, method of maximum likelihood, and the method of empirical characteristic function. A real data application is presented. Finally, extensions to the multivariate case are outlined.     

Minimum phi-divergence estimators for loglinear models with linear constraints and multinomial sampling

In this paper the family ofφ-divergence estimators for loglinear models with linear constraints and multinomial sampling is studied. This family is an extension of the maximum likelihood estimator studied by Haber and Brown (1986). A simulation study is presented and some alternative estimators to the maximum likelihood are obtained. This work was parcially supported by Grant DGES PB2003-892