Likelihood procedure for testing changes in skew normal model with applications to stock returns

The skew normal distribution family is an attractive distribution family due to its mathematical tractability and inclusion of the normal distribution as the special case. It has wide applications in many applied fields such as finance, economics, and medical research. Such a distribution family has been studied extensively since it was introduced by Azzalini in 1985 Azzalini, A. (1985). A class of distributions which includes the normal ones. Scandinavian Journal of Statistics 12:171–178. [Google Scholar] for the first time. Yet, few work has been done on the study of change point problem related to this distribution family. In this article, we propose the likelihood ratio test (LRT) to detect changes in the parameters of the skew normal distribution associated with some asymptotic results of the test statistic. Simulations have been conducted under different scenarios to investigate the performance of the proposed method. Comparisons to some other existing method indicate the comparable power of the method in detecting changes in parameters of the skew normal distribution model. Applications on two real data: Brazilian and Tanzanian stock returns illustrate the detection procedure.     

A Note On Regions of Given Probability of the Extended Skew-normal Distribution

The construction of regions with assigned probability p and minimum geometric measure has theoretical and practical interests, such as the construction of tolerance regions. Following Azzalini (2001 Azzalini, A. (2001). A note on regions of given probability of the skew-normal distribution. Metron 59(3-4): 27–34. [Google Scholar]) and exploiting the normal approximation of the extended skew-normal distribution when some of its parameters go to infinity, we discuss an approach for the construction of regions with assigned probability p for the bivariate extended skew-normal distribution.     

A Generalization of the Skew-Normal Distribution: The Beta Skew-Normal

We consider a new generalization of the skew-normal distribution introduced by Azzalini (1985 Azzalini , A. ( 1985 ). A class of distributions which includes the normal ones . Scand. J. Statis. 12 ( 2 ): 171 – 178 .[Web of Science ®] , [Google Scholar]). We denote this distribution Beta skew-normal (BSN) since it is a special case of the Beta generated distribution (Jones, 2004 Jones , M. C. ( 2004 ). Families of distributions of order statistics . Test 13 ( 1 ): 1 – 43 .[Crossref], [Web of Science ®] , [Google Scholar]). Some properties of the BSN are studied. We pay attention to some generalizations of the skew-normal distribution (Bahrami et al., 2009 Bahrami , W. , Agahi , H. , Rangin , H. ( 2009 ). A two-parameter Balakrishnan skew-normal distribution . J. Statist. Res. Iran 6 : 231 – 242 . [Google Scholar]; Sharafi and Behboodian, 2008 Sharafi , M. , Behboodian , J. ( 2008 ). The Balakrishnan skew-normal density . Statist. Pap. 49 : 769 – 778 .[Crossref], [Web of Science ®] , [Google Scholar]; Yadegari et al., 2008 Yadegari , I. , Gerami , A. , Khaledi , M. J. ( 2008 ). A generalization of the Balakrishnan skew-normal distribution . Statist. Probab. Lett. 78 : 1165 – 1167 .[Crossref], [Web of Science ®] , [Google Scholar]) and to their relations with the BSN.     

New Influence Measures in Polytomous Logistic Regression Models Based on Phi-Divergence Measures

We consider the asymptotic distribution of divergence-based influence measures which are an extension for polytomous logistic regression of an influence measure proposed in Johnson (1985 Johnson, W.O. (1985). Influence measures for logistic regression: Another point of view. Biometrika 72: 5965.[Crossref], [Web of Science ®] , [Google Scholar]), for binary logistic regression. A numerical example compares the classical Cook’s distance with the divergence based influence measures.     

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.     

An objective Bayesian estimation in a two-period crossover design

Two-period crossover design is one of the commonly used designs in clinical trials. But, the estimation of treatment effect is complicated by the possible presence of carryover effect. It is known that ignoring the carryover effect when it exists can lead to poor estimates of the treatment effect. The classical approach by Grizzle (1965 Grizzle, J.E. (1965). The two-period change-over design and its use in clinical trials. Biometrics 21:467–480. See Grizzle (1974) for corrections.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) consists of two stages. First, a preliminary test is conducted on carryover effect. If the carryover effect is significant, analysis is based only on data from period one; otherwise, analysis is based on data from both periods. A Bayesian approach with improper priors was proposed by Grieve (1985 Grieve, A.P. (1985). A Bayesian analysis of the two-period crossover design for clinical trials. Biometrics 41:979–990.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) which uses a mixture of two models: a model with carryover effect and another without. The indeterminacy of the Bayes factor due to the arbitrary constant in the improper prior was addressed by assigning a minimally discriminatory value to the constant. In this article, we present an objective Bayesian estimation approach to the two-period crossover design which is also based on a mixture model, but using the commonly recommended Zellner–Siow g-prior. We provide simulation studies and a real data example and compare the numerical results with Grizzle (1965 Grizzle, J.E. (1965). The two-period change-over design and its use in clinical trials. Biometrics 21:467–480. See Grizzle (1974) for corrections.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar])’s and Grieve (1985 Grieve, A.P. (1985). A Bayesian analysis of the two-period crossover design for clinical trials. Biometrics 41:979–990.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar])’s approaches.     

Test for periodicity in restrictive EXPAR models

Abstract

This article is devoted to study the problem of test of periodicity in the restricted exponential autoregressive (EXPAR) model. The local asymptotic normality property, of this model, is shown via the adapted sufficient conditions due to Swensen (1985 Swensen, A.R. (1985). The asymptotic distribution of the likelihood ratio for autoregressive time series with a regression trend. J. Multivariate Anal. 16:54–70.[Crossref], [Web of Science ®] , [Google Scholar]). Using this result, in the case where the innovation density is specified, we obtain a parametric local asymptotic “most stringent” test.     

A new lifetime model with different types of failure rate

A new class of lifetime distributions, which can exhibit with upside-down bathtub-shaped, bathtub-shaped, decreasing, and increasing failure rates, is introduced. The new distribution is constructed by compounding generalized Weibull and logarithmic distributions, leading to improvement on the lifetime distribution considered in Dimitrakopoulou et al. (2007 Dimitrakopoulou, T., K. Adamidis, and S. Loukas. 2007. A lifetime distribution with an upside-down bathtub-shaped hazard function. IEEE Transactions on Reliability 56:308–11.[Crossref], [Web of Science ®] , [Google Scholar]) by having no restriction on the shape parameter and extending the result studied by Tahmasbi and Rezaei (2008 Tahmasbi, R., and S. Rezaei. 2008. A two-parameter lifetime distribution with decreasing failure rate. Computational Statistics and Data Analysis 52:3889–901.[Crossref], [Web of Science ®] , [Google Scholar]) in the general form. The proposed model includes the exponential–logarithmic and Weibull–logarithmic distributions as special cases. Various statistical properties of the proposed class are discussed. Furthermore, estimation via the maximum likelihood method and the Fisher information matrix are discussed. Applications to real data demonstrate that the new class of distributions is more flexible than other recently proposed classes.     

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.     

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.     

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.     

Skew-Normal ARMA Models with Nonlinear Heteroscedastic Predictors

Multivariate skew-normal (SN) distributions (Azzalini and Dalla Valle, 1996 Azzalini , A. , Dalla Valle , A. ( 1996 ). The multivariate skew-normal distribution . Biometrika 83 : 715 – 726 .[Crossref], [Web of Science ®] , [Google Scholar]) enjoy some of the useful properties of normal distributions, have nonlinear heteroscedastic predictors but lack the closure property of normal distributions (the sum of independent SN random variables is not SN). Recently, there has been a proliferation of classes of SN distributions with certain closure properties, one of the most promising being the closed skew-normal (CSN) distributions of González-Farías et al. (2004 González-Farías , G. , Dominguez-Molina , J. A. , Gupta , A. K. ( 2004 ). Additive properties of skew-normal random vectors . J. Statist. Plann. Infer. 126 : 521 – 534 .[Crossref], [Web of Science ®] , [Google Scholar]). We study the construction of stationary SN ARMA models for colored SN noise and show that their finite-dimensional distributions are skew-normal, seldom strictly stationary and their covariance functions differ from their normal ARMA counterparts in that they do not converge to zero for large lags. The situation is better for ARMA models driven by CSN noise, but at the additional cost of considerable computational complexity and a less explicit skewness parameter. In view of these results, the widespread use of such classes of SN distributions in the framework of ARMA models seem doubtful.     

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.     

On General Continuous Triangular and Two-sided Power Distributions

In this article, we directly introduce the continuous version of the general discrete triangular distributions (Kokonendji and Zocchi, 2010 Kokonendji, C.C., Zocchi, S.S. (2010). Extensions of discrete triangular distribution and boundary bias in kernel estimation for discrete functions. Statist. Probab. Lett. 80:1655–1662.[Crossref], [Web of Science ®] , [Google Scholar]). It is bounded and, in general, unimodal with pike. It contains thus a very useful class of two-sided power distributions (van Dorp and Kotz, 2002a Van Dorp, J.R., Kotz, S. (2002a). A novel extension of the triangular distribution and its parameter estimation. Statistician 51:1–17. [Google Scholar],b Van Dorp, J.R., Kotz, S. (2002b). The standard two-sided power distribution and its properties; with applications in financial engineering. Amer. Statistician 56:90–99.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar], 2003 Van Dorp, J.R., Kotz, S. (2003). Generalization of two-sided power distributions and their convolution. Commun. Statist. Theor. Meth. 32:1703–1723.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]). Moments, particular cases, limit distributions, and relations between parameters are straightforwardly derived.     

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.     

Fisher Information Matrix for the Generalized Feller-Pareto Distribution

In this article, the exact form of Fisher information matrix for the generalized Feller-Pareto (GFP) distribution is determined. The GFP family is a general distribution which includes a variety of distributions as special cases. For example:

??generalized Singh-Maddala distribution which in turn includes Burr, Fisk, and Lomax distribution (see Kleiber and Kotz, 2003 Kleiber, C., Kotz, S. (2003). Statistical Size Distributions in Economics and Actuarial Sciences. Hoboken, NJ: John Wiley.[Crossref] , [Google Scholar]);

??a Pareto IV distribution which includes a hierarchy of Pareto models, omitted an additional location parameter (see Arnold, 1983 Arnold, B.C. (1983). Pareto Distributions. Fairland, MD: International Cooperative Publishing House. [Google Scholar], 2008 Arnold, B.C. (2008). Pareto and generalized pareto distributions. In: Modeling Income Distributions and Lorenz Curves, Economic Studies in Equality, Social Exclusion and Well-Being, Chotikapanich, D. (Ed.), New York: Springer. pp. 119–145.[Crossref] , [Google Scholar]); and

??beta Lomax distribution which includes, for example, beta II and Lomax distributions.

Application of these distributions covers a wide spectrum of areas ranging from actuarial science, economics, finance to bioscience, telecommunications, and medicine.     

Analysis of Order Statistics from Distributions with Regularly Varying Tails

In this article, necessary conditions for comparing order statistics from distributions with regularly varying tails are discussed in terms of various stochastic orders. A necessary and sufficient condition for stochastically comparing tail behaviors of order statistics is derived. The main results generalize and recover some results in Kleiber (2002 Kleiber, C. (2002). Variability ordering of heavy-tailed distributions with applications to order statistics. Statist. Probab. Lett. 58:381388.[Crossref], [Web of Science ®] , [Google Scholar], 2004 Kleiber, C. (2004). Lorenz ordering of order statistics from log-logistic and related distributions. J. Statist. Plann. Infer. 120:2004.[Crossref], [Web of Science ®] , [Google Scholar]). Extensions to coherent systems are mentioned as well.     

Subjective Bayesian Analysis of the Elliptical Model

For the multivariate elliptical model subjective Bayesian estimators of the location vector and some functions of the characteristic matrix with the normal-inverse Wishart and the normal-Wishart as prior, respectively, are derived. Fang and Li (1999 Fang, K.T., Li, R.Z. (1999). Bayesian statistical inference on elliptical matrix distributions. J. Multivariate Anal. 70: 66–85.[Crossref], [Web of Science ®] , [Google Scholar]) considered the elliptical model for Bayesian analysis for an objective prior structure. In addition, the newly developed results are applied to the multivariate normal- and t-distribution. A performance study is done to evaluate the normal-gamma and normal-inverse gamma distributions as suitable priors. A practical application for the posterior distributions of the multivariate t-distribution is included by means of Gibbs sampling and a Metropolis-Hastings algorithm.