Time series models associated with Mittag-Leffler type distributions and its properties

ABSTRACT

We have provided a fractional generalization of the Poisson renewal processes by replacing the first time derivative in the relaxation equation of the survival probability by a fractional derivative of order α(0 < α ? 1). A generalized Laplacian model associated with the Mittag-Leffler distribution is examined. We also discuss some properties of this new model and its relevance to time series. Distribution of gliding sums, regression behaviors, and sample path properties are studied. Finally we introduce the q-Mittag-Leffler process associated with the q-Mittag-Leffler distribution.     

THE q-BETA-GEOMETRIC DISTRIBUTION AS A MODEL FOR FECUNDABILITY

The number of sterile couples in a retrospective study of the number of cycles to conception is necessarily zero; this is not so for a prospective study. The paper puts forward a modification of Weinberg and Gladen's beta geometric model for cycles to conception that is suitable for both types of investigation. The probability that a couple achieves conception at the xth cycle, but not earlier, is assumed to take the form R_x = (1 ? ρ)/(1 ? m ^x?1 ρ/u), instead of μ/(1 ? θ + θx). The set of parameter restraints (0 < m < 1, 0< ρ < 1, 1 < u) is appropriate for retrospective data, whilst the alternative set of restraints (1 < m, 1 < ρ, 0 < u < 1) is appropriate for prospective data. The decrease in R_x over time can be interpreted not only as a time effect, but also as a heterogeneity effect by replacing Weinberg and Gladen's beta mixture of geometric distributions by a q-beta mixture.     

Tests for properties of residual life

The lOOα -percentile (0 < α < 1) residual life function at time x is defined to be the lOOα -percentile of the remaining life given survival up to time x. Joe and Proschan (1982b) develop tests for testing the alternatives representing decreasing 100α-percentile residual life (DPRL-α ) and the property ‘new better than used with respect to the lOOα -percentile’ (NBUP-α ). In this paper, tests are developed for DPRL[α, l) and NBUP[α, l) alternatives, where DPRL[α, l) is the class of life distributions which are DPRL-β distributions for all a ≤ β < 1 if 0 ≤ α < 1 and for all 0 < β < 1 if α = 0, and NBUP[α, l) is similarly defined. When α = 0, the DPRL[α, l) class of life distributions corresponds to the increasing failure rate class and the NBUP[α, l) class of life distributions corresponds to the new better than used class, and the test statistics are respectively asymptotically equivalent to the Hollander and Proschan (1975) test statistics for decreasing mean residual life and new better than used in expectation alternatives.     

z Test for the significance of the mean of a stable probability distribution with 1<α≤2

The analysis of data using a stable probability distribution with tail parameter α<2 (sometimes called a Pareto–Levy distribution) seems to have been avoided in the past in part because of the lack of a significance test for the mean, even though it appears to be the correct distribution to use for describing returns in the financial markets. A z test for the significance of the mean of a stable distribution with tail parameter 1<α≤2 is defined. Tables are calculated and displayed for the 5% and 1% significance levels for a range of tail and skew parameters α and β. Through the use of maximum likelihood estimates, the test becomes a practical tool even when α and β are not that accurately determined. As an example, the z test is applied to the daily closing prices for the Dow Jones Industrial average from 2 January 1940 to 19 March 2010.     

ASYMPTOTIC AND SMALL SAMPLE STATISTICAL PROPERTIES OF RANDOM FRAILTY VARIANCE ESTIMATES FOR SHARED GAMMA FRAILTY MODELS

This paper concerns maximum likelihood estimation for the semiparametric shared gamma frailty model; that is the Cox proportional hazards model with the hazard function multiplied by a gamma random variable with mean 1 and variance θ. A hybrid ML-EM algorithm is applied to 26 400 simulated samples of 400 to 8000 observations with Weibull hazards. The hybrid algorithm is much faster than the standard EM algorithm, faster than standard direct maximum likelihood (ML, Newton Raphson) for large samples, and gives almost identical results to the penalised likelihood method in S-PLUS 2000. When the true value θ₀ of θ is zero, the estimates of θ are asymptotically distributed as a 50–50 mixture between a point mass at zero and a normal random variable on the positive axis. When θ₀ > 0, the asymptotic distribution is normal. However, for small samples, simulations suggest that the estimates of θ are approximately distributed as an x ? (100 ? x)% mixture, 0 ≤ x ≤ 50, between a point mass at zero and a normal random variable on the positive axis even for θ₀ > 0. In light of this, p-values and confidence intervals need to be adjusted accordingly. We indicate an approximate method for carrying out the adjustment.     

A generalization of the binomial distribution

This paper presents a new departure in the generalization of the binomial distribution by adopting the assumption that the underlying Bernoulli trials take on the values α or β where α < β, rather than the conventional values 0 or 1. The adoption of this more general assumption renders the binomial distribution a four-parameter distribution of the form B(n,p,α,β), and requires the generalization of Romanovsky's (1923) reduction formula for central moments. This paper assesses the usefulness of B(n,p,α,β), and its reduction formula, in the numerical analysis of two problems of interest to decision theorists.     

Analysis of mixed correlated bivariate zero-inflated count and (k,l)-inflated beta responses with application to social network datasets

This paper presents a new model that monitors the basic network formation mechanisms via the attributes through time. It considers the issue of joint modeling of longitudinal inflated (0, 1)-support continuous and inflated count response variables. For joint model of mentioned response variables, a correlated generalized linear mixed model is studied. The fraction response is inflated in two points k and l (k < l) and a k and l inflated beta distribution is introduced to use as its distribution. Also, the count response is inflated in zero and we use some members of zero-inflated power series distributions, hurdle-at-zero, members of zero-inflated double power series distributions and zero-inflated generalized Poisson distribution as our count response distribution. A full likelihood-based approach is used to yield maximum likelihood estimates of the model parameters and the model is applied to a real social network obtained from an observational study where the rate of the ith node’s responsiveness to the jth node and the number of arrows or edges with some specific characteristics from the ith node to the jth node are the correlated inflated (0, 1)-support continuous and inflated count response variables, respectively. The effect of the sender and receiver positions in an office environment on the responses are investigated simultaneously.     

A discrete distribution including the Poisson

Abstract

This article presents a new generalization of the Poisson distribution, with the parameters α > 0 and θ > 0, using the Marshall and Olkin (1997 Marshall, A.W., Olkin, I. (1997). A new method for adding a parameter to a family of distributions with application to the exponential and Weibull families. Biometrika 84(3):641–652.[Crossref], [Web of Science ®] , [Google Scholar]) scheme and adding a parameter to the classical Poisson distribution. The particular case of α = 1 gives the Poisson distribution. The new distribution is unimodal and has a failure rate that monotonically increases or decreases depending on the value of the parameter α. After reviewing some of the properties of this distribution, we investigated the question of parameter estimation. Expected frequencies were calculated for two data sets, one with an index of dispersion larger than one and the other with an index of dispersion smaller than one. In both cases the distribution provided a very satisfactory fit.     

On computer generation of gamma random variables by rejection and composition procedures 2

The paper presents various algorithms for generating gamma random variables, by combining rejection and composition procedures. Two efficient algorithms are given for the case when the parameter of the gamma distribution is 0<v<l. For the case vl, several algorithms are given but they (except one), work reliably only for small values of v. Results of some computer tests together with FORTRAN subroutines are also presented.     

A note on a strong renewal theorem

Let {S_n, n ≥ 1} be a sequence of partial sums of independent and identically distributed non-negative random variables with a common distribution function F. Let F belong to the domain of attraction of a stable law with exponent α, 0 < α < 1. Suppose H(t) = ? N(t), t ? 0, where N(t) = max(n : S_n ≥ t). Under some additional assumptions on F, the difference between H(t) and its asymptotic value as t → ∞ is estimated.     

Peculiar bias properties of the ols estimator when applied to a dynamic model with autocorrelated disturbances

This study reveals that contrary to the conventional wisdom among econometricians, the bias of the OLS estimator can be quite small when the estimator is applied to a geometrically distributed lag model, y_t<ce:glyph name="dbnd6"/> α + βx _t+ λy _t-₁. + u_t, with autocorrelated disturbances, be they AR(1), MA(1), MA(2), AR(2), and ARMA(1,1). This happens when λ is large and x_tis smoothly trended (e.g., a real GNP series). In fact, the bias of the OLS estimator becomes zero at one parameter combination, and the OLS estimator performs well over a wide range around this parameter combination. By decomposing the disturbance term into two parts, the paper also explains why OLS shows such an unexpected property. These findings have both pedagogical and practical significance.     

Mean square error efficient estimation of an exponential mean under an exchangeable single outlier model

This paper is concerned with estimating θ, the mean of an exponential distribution under a single outlier exchangeable model. It is a.ssumed that the single outlying observation is also exponential with mean θ/α, where 0 < α < 1. The estirnators proposed are weighted averages of the order statistics. The formulas for the weights minimizing the mean square error are presented. These weights are calculated for certain combinations of the sample size n and of α. It is found that the optimal weights very nearly have a certain form. The mean square errors of a simplified estitnator are compared lo those of Joshi (1972, 1988) and of Clhikkagoudar and Kunchur (1980). A nlodification of Joshi's iterative procedure is suggested.     

CONDITIONING VS. STANDARDIZATION FOR CONTRASTS ON ERRORS ATTRACTED TO STABLE LAWS

In multiple regression and other settings one encounters the problem of estimating sampling distributions for contrast operations applied to i.i.d. errors. Permutation bootstrap applied to least squares residuals has been proven to consistently estimate conditionalsampling distributions of contrasts, conditional upon order statistics of errors, even for long-tailed error distributions. How does this compare with the unconditional sampling distribution of the contrast when standardizing by the sample s.d. of the errors (or the residuals)? For errors belonging to the domain of attraction of a normal we present a limit theorem proving that these distributions are far closer to one another than they are to the limiting standard normal distribution. For errors attracted to α-stable laws with α ≤ 2 we construct random variables possessing these conditional and unconditional sampling distributions and develop a Poisson representation for their a.s. limit correlation ρ_α. We prove that ρ₂= 1, ρ_α→ 1 for α → 0 + or 2 ?, and ρ_α< 1 a.s. for α < 2.     

Statistics of Shape,Direction and Cylindrical Variables

In statistical shape analysis, the shape of an object is understood to be what remains after the effects of location, scale and rotation are removed. We consider the distributional problem of triangular shape and an associated direction; motivated by a data set of microscopic fossils. We begin by constructing a parallel transport system such that the data transform onto the space 𝒮²×𝒮². A joint shape distribution on 𝒮²×𝒮¹ is proposed based on Jupp & Mardia's bivariate distribution on 𝒮²×𝒮¹. For concentrated data, an approximation to the distribution on 𝒮²×𝒮¹ is given by a distribution on ?¹×𝒮¹, and we explore a distribution on this space by extending Mardia & Sutton's distribution on ?²×𝒮¹. In this distribution, the expected edgel direction varies linearly in the shape coordinates. This is found to be a useful model for the microfossil data.     

The Distribution of Sums of Certain I.I.D. Pareto Variates

ABSTRACT

Though the Pareto distribution is important to actuaries and economists, an exact expression for the distribution of the sum of n i.i.d. Pareto variates has been difficult to obtain in general. This article considers Pareto random variables with common probability density function (pdf) f(x) = (α/β) (1 + x/β)^α+1 for x > 0, where α = 1,2,… and β > 0 is a scale parameter. To date, explicit expressions are known only for a few special cases: (i) α = 1 and n = 1,2,3; (ii) 0 < α < 1 and n = 1,2,…; and (iii) 1 < α < 2 and n = 1,2,…. New expressions are provided for the more general case where β > 0, and α and n are positive integers. Laplace transforms and generalized exponential integrals are used to derive these expressions, which involve integrals of real valued functions on the positive real line. An important attribute of these expressions is that the integrands involved are non oscillating.     

Utilizing the Flexibility of the Epsilon-Skew-Normal Distribution for Common Regression Problems

In this paper we illustrate the properties of the epsilon-skew-normal (ESN) distribution with respect to developing more flexible regression models. The ESN model is a simple one-parameter extension of the standard normal model. The additional parameter ~ corresponds to the degree of skewness in the model. In the fitting process we take advantage of relatively new powerful routines that are now available in standard software packages such as SAS. It is illustrated that even if the true underlying error distribution is exactly normal there is no practical loss n power with respect to testing for non-zero regression coefficients. If the true underlying error distribution is slightly skewed, the ESN model is superior in terms of statistical power for tests about the regression coefficient. This model has good asymptotic properties for samples of size n>50.     

Empirical bayes interval estimates involving uniform distributions

Bayesian and empirical Bayesian decision rules are exhibited for the interval estimation of the parameter 0 of a Uniform (0,θ) distribution. The estimate ?,δ>resulting in the interval [?,?+δ]suffers loss given by L(?,δ>,θ)=1-[?≦e≦?+δ]+c₁((?-θ)²+(?+δ?θ)²))+c₂δ. The solution is presented for prior distributions G which have bounded support, no point masses,∫θ^?mdG(θ)<∞ and for some integer m. An example is presented involving a particular parametric form for G and rates of risk convergence in the empirical Bayes problem for this example are calculated.     

Rate of consistency of one sample tests of location

Let X₁,…,X_n be a sample from a population with continuous distribution function F(x?θ) such that F(x)+F(-x)=1 and 0<F(x)<1, x?R¹. It is shown that the power- function of a monotone test of H: θ=θ₀ against K: θ>θ₀ cannot tend to 1 as θ?θ₀ → ∞ more than n times faster than the tails of F tend to 0. Some standard as well as robust tests are considered with respect to this rate of convergence.     

A Note on an Approximation of the Mobile Communications Expenditures Distribution Function Using a Mixture Model

Approximating the distribution of mobile communications expenditures (MCE) is complicated by zero observations in the sample. To deal with the zero observations by allowing a point mass at zero, a mixture model of MCE distributions is proposed and applied. The MCE distribution is specified as a mixture of two distributions, one with a point mass at zero and the other with full support on the positive half of the real line. The model is empirically verified for individual MCE survey data collected in Seoul, Korea. The mixture model can easily capture the common bimodality feature of the MCE distribution. In addition, when covariates were added to the model, it was found that the probability that an individual has non-expenditure significantly varies with some variables. Finally, the goodness-of-fit test suggests that the data are well represented by the mixture model.     

Inference for Random Coefficient INAR(1) Process Based on Frequency Domain Analysis

In this article, we present the explicit expressions for the higher-order moments and cumulants of the first-order random coefficient integer-valued autoregressive (RCINAR(1)) process. The spectral and bispectral density functions are also obtained, which can characterize the RCINAR(1) process in the frequency domain. We use a frequency domain approach which is named Whittle criterion to estimate the parameters of the process. We propose a test statistic which is based on the frequency domain approach for the hypothesis test, H₀: α = 0?H₁: 0 < α < 1, where α is the mean of the random coefficient in the process. The asymptotic distribution of the test statistic is obtained. We compare the proposed test statistic with other statistics that can test serial dependence in time series of count via a typically numerical simulation, which indicates that our proposed test statistic has a good power.