Modeling censored data using modified lambda family

In the present article we propose the modified lambda family (MLF) which is the Freimer, Mudholkar, Kollia, and Lin (FMKL) parametrization of generalized lambda distribution (GLD) as a model for censored data. The expressions for probability weighted moments of MLF are derived and used to estimate the parameters of the distribution. We modified the estimation technique using probability weighted moments. It is shown that the distribution provides reasonable fit to a real censored data.     

The generalized transmuted-G family of distributions

We introduce a new class of continuous distributions called the generalized transmuted-G family which extends the transmuted-G class. We provide six special models of the new family. Some of its mathematical properties including explicit expressions for the ordinary and incomplete moments, generating function, Rényi and Shannon entropies, order statistics and probability weighted moments are derived. The estimation of the model parameters is performed by maximum likelihood. The flexibility of the proposed family is illustrated by means of three applications to real data sets.     

Estimation of the Lognormal-Pareto Distribution Using Probability Weighted Moments and Maximum Likelihood

This article deals with the estimation of the lognormal-Pareto and the lognormal-generalized Pareto distributions, for which a general result concerning asymptotic optimality of maximum likelihood estimation cannot be proved. We develop a method based on probability weighted moments, showing that it can be applied straightforwardly to the first distribution only. In the lognormal-generalized Pareto case, we propose a mixed approach combining maximum likelihood and probability weighted moments. Extensive simulations analyze the relative efficiencies of the methods in various setups. Finally, the techniques are applied to two real datasets in the actuarial and operational risk management fields.     

Theoretical L-moments and TL-moments using combinatorial identities and finite operators

Moments have been traditionally used to characterize a probability distribution. Recently, linear moments (L-moments) and trimmed L-moments (TL-moments) are appealing alternatives to the conventional moments. This paper focuses on the computation of theoretical L-moments and TL-moments and emphasizes the use of combinatorial identities. We are able to derive new closed-form formulas of L-moments and TL-moments for continuous probability distributions. Finally, closed-form formulas for the L-moments for the exponential distribution and the uniform distribution are also obtained.     

Testing the Stochastic Structure of Production: A Flexible Moment-Based Approach

Conventional production function specifications are shown to impose restrictions on the probability distribution of output that cannot be tested with the conventional models. These restrictions have important implications for firm behavior under uncertainty. A flexible representation of a firm's stochastic technology is developed based on the moments of the probability distribution of output. These moments are a unique representation of the technology and are functions of inputs. Large-sample estimators are developed for a linear moment model that is sufficiently flexible to test the restrictions implied by conventional production function specifications. The flexible moment-based approach is applied to milk production data. The first three moments of output are statistically significant functions of inputs. The cross-moment restrictions implied by conventional models are rejected.     

Fractional Moments and Maximum Entropy: Geometric Meaning

We consider the problem of recovering a probability density on a bounded or unbounded subset D of [0, ∞), from the knowledge of its sequence of fractional moments within a maximum entropy (MaxEnt) setup. Based upon entropy convergence results previously formulated, the fractional moments are selected so that the entropy of the MaxEnt approximation be minimum. A geometric interpretation of the reconstruction procedure is formulated as follows: the two moment curves generated by the unknown density and its MaxEnt approximation are interpolating in Hermite-Birkoff sense; that is, they are both interpolating and tangent at the selected nodes.     

GMM nonparametric correction methods for logistic regression with error‐contaminated covariates and partially observed instrumental variables

We consider logistic regression with covariate measurement error. Most existing approaches require certain replicates of the error‐contaminated covariates, which may not be available in the data. We propose generalized method of moments (GMM) nonparametric correction approaches that use instrumental variables observed in a calibration subsample. The instrumental variable is related to the underlying true covariates through a general nonparametric model, and the probability of being in the calibration subsample may depend on the observed variables. We first take a simple approach adopting the inverse selection probability weighting technique using the calibration subsample. We then improve the approach based on the GMM using the whole sample. The asymptotic properties are derived, and the finite sample performance is evaluated through simulation studies and an application to a real data set.     

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.     

Explicit expressions for gamma spacings and their moments

This work is concerned with the distributions of spacings from a two-parameter gamma distribution, when the shape parameter is a positive integer (or Erlang Distribution). We express the probability density functions of spacings and their moments in closed forms that are easy to implement using computer algebra systems like Mathematica.     

A note on the scale parameter of the dirichlet process

This paper gives an interpretation for the scale parameter of a Dirichlet process when the aim is to estimate a linear functional of an unknown probability distribution. We provide exact first and second posterior moments for such functionals under both informative and noninformative prior specifications. The noninformative case provides a normal approximation to the Bayesian bootstrap.     

The four-parameter Burr XII distribution: Properties,regression model,and applications

This paper introduces a new four-parameter lifetime model called the Weibull Burr XII distribution. The new model has the advantage of being capable of modeling various shapes of aging and failure criteria. We derive some of its structural properties including ordinary and incomplete moments, quantile and generating functions, probability weighted moments, and order statistics. The new density function can be expressed as a linear mixture of Burr XII densities. We propose a log-linear regression model using a new distribution so-called the log-Weibull Burr XII distribution. The maximum likelihood method is used to estimate the model parameters. Simulation results to assess the performance of the maximum likelihood estimation are discussed. We prove empirically the importance and flexibility of the new model in modeling various types of data.     

Distributional properties of continuous time processes: from CIR to bates

In this paper, we compute closed-form expressions of moments and comoments for the CIR process which allows us to provide a new construction of the transition probability density based on a moment argument that differs from the historic approach. For Bates’ model with stochastic volatility and jumps, we show that finite difference approximations of higher moments such as the skewness and the kurtosis are unstable and, as a remedy, provide exact analytic formulas for log-returns. Our approach does not assume a constant mean for log-price differentials but correctly incorporates volatility resulting from Ito’s lemma. We also provide R, MATLAB, and Mathematica modules with exact implementations of the theoretical conditional and unconditional moments. These modules should prove useful for empirical research.

     

New class of location-parameter discrete probability distributions and their characterizations

A new class of location-parameter discrete probability distributions (LDPD) has been defined where the population mean is the location parameter. It has been shown that some single parameter discrete distributions do not belong to this class and all discrete probability distributions belonging to this class can be characterized by their variances only. Expressions are given for the first four central moments and a recurrence formula for higher central moments has been obtained. Eight theorems are given to characterize the various distributions in the LDPD class.     

Pareto Poisson–Lindley distribution with applications

A new lifetime distribution is introduced based on compounding Pareto and Poisson–Lindley distributions. Several statistical properties of the distribution are established, including behavior of the probability density function and the failure rate function, heavy- and long-right tailedness, moments, the Laplace transform, quantiles, order statistics, moments of residual lifetime, conditional moments, conditional moment generating function, stress–strength parameter, Rényi entropy and Song's measure. We get maximum-likelihood estimators of the distribution parameters and investigate the asymptotic distribution of the estimators via Fisher's information matrix. Applications of the distribution using three real data sets are presented and it is shown that the distribution fits better than other related distributions in practical uses.     

Estimating the probability of obtaining nonfeasible parameter estimates of the generalized pareto distribution

In this paper we consider the problem of estimating the parameters of the generalized Pareto distribution. Both the method of moments and probability-weighted moments do not guarantee that their respective estimates will be consistent with the observed data. We present simple programs to predict the probability of obtaining such nonfeasible estimates. Our estimation techniques are based on results from intensive simulations and the successful modelling of the lower tail of the distribution of the upper bound of the support. More simulations are performed to validate the new procedure.     

Recurrence formula on the joint distribution of sample skewness and kurtosis under normality

Abstract

Two recurrence relations with respect to sample size are given concerning the joint distribution of skewness and kurtosis of random observations from a normal population: one between the probability density functions and the other between the product moments. As a consequence, the latter yields a recurrence formula for the moments of sample kurtosis. The exact moments of Jarque-Bera statistic is also given.     

Useful relationship between the log-survival function and truncated moments,with applications

We present a new characterization technique extracted from a well known idea in statistical inference. We use the partial derivative of the logarithm of the survival function in connection with truncated moments to characterize several probability distributions. Our methods introduce a unified technique to obtain several well known results in a unified way.     

Symbolic solutions of some linear recurrences

A symbolic method for solving linear recurrences of combinatorial and statistical interest is introduced. This method essentially relies on a representation of polynomial sequences as moments of a symbol that looks as the framework of a random variable with no reference to any probability space. We give several examples of applications and state an explicit form for the class of linear recurrences involving Sheffer sequences satisfying a special initial condition. The results here presented can be easily implemented in a symbolic software.     

Moment density estimation for positive random variables

An unknown moment-determinate cumulative distribution function or its density function can be recovered from corresponding moments and estimated from the empirical moments. This method of estimating an unknown density is natural in certain inverse estimation models like multiplicative censoring or biased sampling when the moments of unobserved distribution can be estimated via the transformed moments of the observed distribution. In this paper, we introduce a new nonparametric estimator of a probability density function defined on the positive real line, motivated by the above. Some fundamental properties of proposed estimator are studied. The comparison with traditional kernel density estimator is discussed.     

Evaluating the Calibration of Multi-Step-Ahead Density Forecasts Using Raw Moments

The evaluation of multi-step-ahead density forecasts is complicated by the serial correlation of the corresponding probability integral transforms. In the literature, three testing approaches can be found that take this problem into account. However, these approaches rely on data-dependent critical values, ignore important information and, therefore lack power, or suffer from size distortions even asymptotically. This article proposes a new testing approach based on raw moments. It is extremely easy to implement, uses standard critical values, can include all moments regarded as important, and has correct asymptotic size. It is found to have good size and power properties in finite samples if it is based on the (standardized) probability integral transforms.