Estimation and diagnostic for skew-normal partially linear models

Partially linear models (PLMs) are an important tool in modelling economic and biometric data and are considered as a flexible generalization of the linear model by including a nonparametric component of some covariate into the linear predictor. Usually, the error component is assumed to follow a normal distribution. However, the theory and application (through simulation or experimentation) often generate a great amount of data sets that are skewed. The objective of this paper is to extend the PLMs allowing the errors to follow a skew-normal distribution [A. Azzalini, A class of distributions which includes the normal ones, Scand. J. Statist. 12 (1985), pp. 171–178], increasing the flexibility of the model. In particular, we develop the expectation-maximization (EM) algorithm for linear regression models and diagnostic analysis via local influence as well as generalized leverage, following [H. Zhu and S. Lee, Local influence for incomplete-data models, J. R. Stat. Soc. Ser. B 63 (2001), pp. 111–126]. A simulation study is also conducted to evaluate the efficiency of the EM algorithm. Finally, a suitable transformation is applied in a data set on ragweed pollen concentration in order to fit PLMs under asymmetric distributions. An illustrative comparison is performed between normal and skew-normal errors.     

New goodness of fit tests for the Cauchy distribution

Some goodness-of-fit procedures for the Cauchy distribution are presented. The power comparisons indicate that the new tests possess good performances among the competitors, especially against symmetric alternatives. A financial data set is analyzed for illustration.     

An approximation to the convolution of gamma distributions

In general, the exact distribution of a convolution of independent gamma random variables is quite complicated and does not admit a closed form. Of all the distributions proposed, the gamma-series representation of Moschopoulos (1985 Moschopoulos, P. G. (1985). The distribution of the sum of independent gamma random variables. Annals of the Institute of Statistical Mathematics 37Part A:541–544. [Google Scholar]) is relatively simple to implement but for particular combinations of scale and/or shape parameters the computation of the weights of the series can result in complications with too much time consuming to allow a large-scale application. Recently, a compact random parameter representation of the convolution has been proposed by Vellaisamy and Upadhye (2009 Vellaisamy, P., Upadhye, N. S. (2009). On the sums of compound negative binomial and gamma random variables. Journal of Applied Probability 46:272–283.[Crossref], [Web of Science ®] , [Google Scholar]) and it allows to give an exact interpretation to the weights of the series. They describe an infinite discrete probability distribution. This result suggested to approximate Moschopoulos’s expression looking for an approximating theoretical discrete distribution for the weights of the series. More precisely, we propose a general negative binomial distribution. The result is an “excellent” approximation, fast and simple to implement for any parameter combination.     

A new generalized two-sided class of distributions with an emphasis on two-sided generalized normal distribution

The ordinary-G class of distributions is defined to have the cumulative distribution function (cdf) as the value of the cdf of the ordinary distribution F whose range is the unit interval at G, that is, F(G), and it generalizes the ordinary distribution. In this work, we consider the standard two-sided power distribution to define other classes like the beta-G and the Kumaraswamy-G classes. We extend the idea of two-sidedness to other ordinary distributions like normal. After studying the basic properties of the new class in general setting, we consider the two-sided generalized normal distribution with maximum likelihood estimation procedure.     

Comparison of the effectiveness of forecasts obtained by means of selected probability functions with respect to forecast error distributions

The forecasting of sales in a company is one of the crucial challenges that must be faced. Nowadays, there is a large spectrum of methods that enable making reliable forecasts. However, sometimes the nature of time series excludes many well-known and widely used forecasting methods (e.g., econometric models). Therefore, the authors decided to forecast on the basis of a seasonally adjusted median of selected probability distributions. The obtained forecasts were verified by means of distributions of the Theil U² coefficient and unbiasedness coefficient.     

Randomly weighted sums of linearly wide quadrant-dependent random variables with heavy tails

This paper investigates tail behavior of the randomly weighted sum ∑ⁿ_{k = 1}θ_kX_k and reaches an asymptotic formula, where X_k, 1 ? k ? n, are real-valued linearly wide quadrant-dependent (LWQD) random variables with a common heavy-tailed distribution, and θ_k, 1 ? k ? n, independent of X_k, 1 ? k ? n, are n non-negative random variables without any dependence assumptions. The LWQD structure includes the linearly negative quadrant-dependent structure, the negatively associated structure, and hence the independence structure. On the other hand, it also includes some positively dependent random variables and some other random variables. The obtained result coincides with the existing ones.     

The bivariate alpha-skew-normal distribution

In this paper, we propose a new bivariate distribution, namely bivariate alpha-skew-normal distribution. The proposed distribution is very flexible and capable of generalizing the univariate alpha-skew-normal distribution as its marginal component distributions; it features a probability density function with up to two modes and has the bivariate normal distribution as a special case. The joint moment generating function as well as the main moments are provided. Inference is based on a usual maximum-likelihood estimation approach. The asymptotic properties of the maximum-likelihood estimates are verified in light of a simulation study. The usefulness of the new model is illustrated in a real benchmark data.     

On fiducial generators

Fiducial inference has been gaining presence recently and it is the intention of the present article to look at the notion of fiducial generators; meaning procedures to simulate parameter values that in some sense correspond to simulations from some implicit fiducial distribution. It is well known that when the distribution has group structure, stemming from the natural pivotal associated, a fiducial may be obtained. It is in the non group distributions that there appears to be still room for finding a fiducial distribution. Recently some general procedures have been proposed for dealing with generalized fiducials, but these depend on certain choices for a structural equation or a fiducial equation, as in Hannig (2009 Hannig, J. (2009). On generalized fiducial inference. Stat. Sin. 19:491–544.[Web of Science ®] , [Google Scholar]) or Taraldsen and Lindqvist (2013 Taraldsen, G., Lindqvist, B.H. (2013). Fiducial theory and optimal inference. Ann. Stat. 41(1):323–341.[Crossref], [Web of Science ®] , [Google Scholar]), respectively. A brief presentation is made of an earlier approach to fiducial inference for multivariate parameters, as in Brillinger (1962 Brillinger, D.R. (1962). Examples bearing on the definition of fiducial probability with a bibliography. Ann. Math. Stat. 33(4):1349–1355.[Crossref] , [Google Scholar]), and the implied fiducial generator introduced in Engen and Lillegård (1997 Engen, S., Lillegård, M. (1997). Stochastic simulation conditioned on sufficient statistics. Biometrika 84(1):235–240.[Crossref], [Web of Science ®] , [Google Scholar]), trying to connect them. Three interesting non group distributions are seen; two of them, the truncated exponential and the two-parameter gamma, already reported in literature. A third non group distribution is analyzed; the inverse Gaussian, connecting the fiducial that results following Brillinger (1962 Brillinger, D.R. (1962). Examples bearing on the definition of fiducial probability with a bibliography. Ann. Math. Stat. 33(4):1349–1355.[Crossref] , [Google Scholar]), with a result pertaining confidence limits for the shape parameter in Hsieh (1990 Hsieh, H.K. (1990). Inferences on the coefficient of variation of an inverse-Gaussian distribution. Commun. Stat. - Theory Methods 19(5):1589–1605.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]). In the three cases, comparisons are made with the Bayesian posteriors that have been known to be close numerically. Some discussion is made on the issue of singularities of the fiducial density and its connection with densities that do not integrate to unity. As to the case of discrete observables, some comments are made for the Bernoulli distribution, only.