First-Order Integer-Valued Autoregressive (INAR (1)) Process: Distributional and Regression Properties

Some properties of a first-order integer-valued autoregressive process (INAR)) are investigated. The approach begins with discussing the self-decomposability and unimodality of the 1-dimensional marginals of the process {X_n} generated according to the scheme X_n=α° X _n-i +e_n, where α° X _n-1 denotes a sum of X_{n - 1}, independent 0 - 1 random variables Y^(n-1), independent of X _n-1 with Pr -( y ^{(n - 1)}= 1) = 1 - Pr ( y ^(n-i)= 0) =α. The distribution of the innovation process ( e _n) is obtained when the marginal distribution of the process ( X _n) is geometric. Regression behavior of the INAR(1) process shows that the linear regression property in the backward direction is true only for the Poisson INAR(1) process.     

Uitbijternegerende schatters bij duplobepalingen

Outlyer-ignoring estimators for measurement in duplo.
By hypothesis a measurement u is the sum of two independent random variables, the normal random variable with expectation μ, and standard error σ, and a random error φ:

Basically two independent measurements u₁ and u₂ over u are to give the estimate x=1/2(u₁+ u₂) over μ.
However, to reduce the effect of the error φ on a final estimate of μ, one adds, according to a common practice, a third or even a fourth measurement u₃, u₄, in the case that the basic pair differs by more than a number A. For this extended set of measurements two outlyer-ignoring estimator y and z of μ are defined, and investigated against three specifications fo the error φ. Also an outlyer-ignoring estimate of σ is considered, and its application is illustrated by an example.     

An alternate development of conditioning

Abstract  The "classical" development of conditioning, due to K olmogorov , does not agree with the "practical" (more intuitive, but unrigorous) way in which probabilists and statisticians actually think about conditioning. This paper describes an alternative to the classical development. It is shown that standard concepts and results can be developed, rigorously, along lines, which correspond to the "practical" approach, and so as to include the classical material as a special case. More specifically, let Xand Y be random variables (r.v.'s) from (Ω, f, P) to ( x, f_x ) and (y. f^y.), respectively. In this paper, the fundamental concept is the conditional probability P(AX = x ), a function of xε x which satisfies a "natural" defining condition. This is used to define a conditional distribution P_y/x, as a mapping x × f^y-R such that, as a function of B, P_ylx=x,(B ) is a probability measure on f^y. Then, for a numerical r.v. Y , conditional expectation E(Y/X) is defined as a mapping x → whose value at x isE(Y/X = x) = ydP_Y/x=i(y ). Basic properties of conditional probabilities, distributions, and expectations, are derived and their existence and uniqueness are discussed. Finally, for a sub-o-algebra and a numerical r.v. Y , the classical conditional expectation E(Y) is obtained as E(Y/X) with X = i , the identity mapping from (Ω, f) to (Ω).     

A central limit theorem for sums of correlated products

Consider a sequence of random points placed on the nonnegative integers with i.i.d. geometric (1/2) interpoint spacings y _i . Let x _i denote the numbers of points placed at integer i . We prove a central limit theorem for the partial sums of the sequence x ₀ y ₀, x ₁ y ₁, . . . The problem is connected with a question concerning different bootstrap procedures.     

PSEUDO-R2 MEASURES FOR SOME COMMON LIMITED DEPENDENT VARIABLE MODELS

Abstract. A large number of different Pseudo- R ² measures for some common limited dependent variable models are surveyed. Measures include those based solely on the maximized likelihoods with and without the restriction that slope coefficients are zero, those which require further calculations based on parameter estimates of the coefficients and variances and those that are based solely on whether the qualitative predictions of the model are correct or not. The theme of the survey is that while there is no obvious criterion for choosing which Pseudo- R ² to use, if the estimation is in the context of an underlying latent dependent variable model, a case can be made for basing the choice on the strength of the numerical relationship to the OLS- R ² in the latent dependent variable. As such an OLS- R ² can be known in a Monte Carlo simulation, we summarize Monte Carlo results for some important latent dependent variable models (binary probit, ordinal probit and Tobit) and find that a Pseudo- R ² measure due to McKelvey and Zavoina scores consistently well under our criterion. We also very briefly discuss Pseudo- R ² measures for count data, for duration models and for prediction-realization tables.     

On divergence and convergence of sums of nonnegative random variables

Abstract  If X ₁, X ₂,… are exponentially distributed random variables thenσ^∞_{k= 1} X_k=∞ with probability 1 iff σ^∞_{k= 1} EX_k=∞. This result, which is basic for a criterion in the theory of Markov jump processes for ruling out explosions (infinitely many transitions within a finite time) is usually proved under the assumption of independence (see FREEDMAN (1971), p. 153–154 or BREI-MAN (1968), p. 337–338), but is shown in this note to hold without any assumption on the joint distribution. More generally, it is investigated when sums of nonnegative random variables with given marginal distributions converge or diverge whatever are their joint distributions.     

Generalized densities of order statistics

Let X ₁, ... , X _n be independent identically distributed random variables with distribution F . We derive expressions for generalized joint 'densities' of order statistics of X ₁, ... , X _n , for arbitrary distributions F , in terms of Radon–Nikodym derivatives with respect to product measures based on F . We then give formulae for conditional distributions of order statistics and use them to derive results concerning Markov properties of order statistics, formulae for distributions of trimmed sums, and other useful representations. Our approach leads to simple and natural expressions which appear not to have been given before.     

Approximation for the cumulative and inverse gamma distribution

The gamma distribution function can be expressed in terms of the Normal distribution and density functions with sufficient accuracy for most practical purposes.
The distribution function for the density x^Λ-1e^-x/μ/μ^ΛΓ(A) on 0 -R(Λ){(1 + 1/1 2Λ) φ(z) + 11 -z/4Λ^1/2+2(z²+ 2)/45Λ] φ(z) /3 Λ^1/2} where φ(z)≅1/[1 +_e^{-2z(√2/π+z2 /28)}] and φ(z) = e^{-z2 /2}/√2π are the Normal distribution and density functions, y is the appropriate root of y-y²/6+y³/36-y⁴/270= In (x/Λμ), z= Λ^1/2 y, and R( Λ) is the remainder term in Stirling's approximation for In Γ(Λ).     

Time Series Analysis of Repeated Surveys: The State–space Approach

Cross sectional estimates from repeated surveys form a time series { y_t }. These estimates can be viewed as the sum y _t = Y _t + e _t of two processes, { Y _t }, the population process and { e _t }, the survey error process. Serial correlations in the latter series are usually present, mainly due to sample overlap. Other sources of data such as censuses, administrative records and demographic population counts are also available. The state–space modelling approach to the analysis of repeated surveys allows combining information from different sources, incorporating benchmarking constraints in a natural way. Results from these methods seem to compare favourably with those from X-11-ARIMA in filtering out survey errors.     

On a crossroad of resampling plans: bootstrapping elementary symmetric polynomials

We investigate the validity of the bootstrap method for the elementary symmetric polynomials S ^{( k )} _n =( ⁿ _k )⁻¹Σ_{1≤ i 1}< ... < i _k ≤ n X _{i 1} ... X _{i k} of i.i.d. random variables X ₁, ..., X _n . For both fixed and increasing order k , as n→∞ the cases where μ=E X ₁[moe2]0, the nondegenerate case, and where μ=E X ₁=0, the degenerate case, are considered.     

Inequalities relating maximal moments to other measures of dispersion

Let X , X ₁, ..., X_k be i.i.d. random variables, and for k ∈ N let D_k ( X ) = E ( X ₁ V ... V X _{k +1}) − EX be the k th centralized maximal moment. A sharp lower bound is given for D ₁( X ) in terms of the Lévy concentration Q_l ( X ) = sup _{x ∈ R} P ( X ∈[ x , x + l ]). This inequality, which is analogous to P. Levy's concentration-variance inequality, illustrates the fact that maximal moments are a gauge of how much spread out the underlying distribution is. It is also shown that the centralized maximal moments are increased under convolution.     

The asymptotic distribution of the sum of weighted squared residuals in binary choice models

It is claimed by some authors that the distribution of the sum of weighted squared residuals, used as a goodness of fit measure in binary choice models, behaves for large n as a x²_n– k–1 distribution. This claim seems to be based on a false analogy with the well–known Pearson x² statistic for frequency tables with a fixed number of cells and cell sizes tending to infinity. We derive the asymptotic (normal) distribution and show that the approximation by the x² distribution in general will not be valid. A new x² test is proposed based on the asymptotic normality of the measure.     

Minimax character of the two-sample χ2-test

Abstract  We consider the problem of testing equality of the mean vectors of two multivariate normal populations when covariates are present and the covariance matrix is known. As an application of the H unt -S tein theorem it is shown that the χ²-^-test of level a maximizes, among all level a tests, the minimum power on each of the contours where the χ²-^-test has constant power. A corollary is that the χ²-^-test is most stringent level a .     

Substituut-variabelen in correlatieberekingen

The recently repeated assertion that in correlation analysis it makes little difference whether one variable (x₂) is used instead of another one (x₃), provided the coefficient of correlation (r₂₃) between x₂ and x₃ is high, is scrutinized.
To that purpose the ranges of coefficients of correlation with respect to the substitute variable are expressed in formula 3. Moreover, by way of example, extreme values of coefficients of simple correlation (r₁₃ and r₃₄), of multiple correlation (R_1.34 and R_3.14) and of regression (α₁₃ and α₁₄, α₃₁ and α₃₄) relating to the substitute variable, are calculated on the basis of empirical values of coefficients of simple correlation relating to the substituted and the remaining variables.
The outcome of those calculations are summarized in the tables 1 and 3, and in the graph.
Table 1 presents ranges of r₁₃ for given values of r₁₂ and r₂₃, table 3 shows extreme values of coefficients of single and multiple correlation and regression in case an additional variable x₄ is introduced and r₁₂, r₁₄, r₂₄ and r₂₃ are given. The graph shows an ellipse as the boundary of the inner closed domain of compatible values of r₁₃ and r₃₄.
Those results clearly indicate the need for caution in substituting one variable by another.     

Selection from Logistic populations

Assume k ( k ≥ 2) independent populations π₁, π₂μ_k are given. The associated independent random variables X_i,( i = 1,2,… k ) are Logistically distributed with unknown means μ₁, μ₂, μ_k and equal variances. The goal is to select that population which has the largest mean. The procedure is to select that population which yielded the maximal sample value. Let μ₍₁₎≤μ₍₂₎≤…≤μ_(k) denote the ordered means. The probability of correct selection has been determined for the Least Favourable Configuration μ₍₁₎=μ₍₂₎==μ_{(k – 1)}=μ_(k)–δ where δ > 0. An exact formula for the probability of correct selection is given.     

On the asymptotically uniform distribution modulo 1 of extreme order statistics

Let (X_m)∞₁ be a sequence of independent and identically distributed random variables. We give sufficient conditions for the fractional part of rnax (X₁., X_n) to converge in distribution, as n ←∞ to a random variable with a uniform distribution on [0, 1).     

On Bayesian selection of the best normal population using theKullback–Leibler divergence measure

In this paper, we use the Bayesian approach to study the problem of selecting the best population among k different populations π₁, ..., π_k (k≥2) relative to some standard (or control) population π₀. Here, π₀ is considered to be the population with the desired characteristics. The best population is defined to be the one which is closest to the ideal population π₀ . The procedure uses the idea of minimizing the posterior expected value of the Kullback–Leibler (KL) divergence measure of π _i from π₀. The populations under consideration are assumed to be multivariate normal. An application to regression problems is also presented. Finally, a numerical example using real data set is provided to illustrate the implementation of the selection procedure.     

A simple measure of skewness

Groeneveld (1986) in discussing the skewness for the Weibull family has pointed out the shortcomings of the classical measures of asymmetry—the standardized third moment and the Pearson measure of skewness. He has shown that a modified form of the Pearson measure b₃= (μ-m)/E|X-m| portrays the skewness of Weibull family quite well. We give another competitive measure of skewness T that is easy to interpret and is based on conditional expectations. The proposed measure satisfies the desirable properties of a skewness measure.     

Characteristic properties of order statistics based on random sample size from an exponential distribution

Suppose X₁, X₂, X_m is a random sample of size m from a population with probability density function f (x), x > 0), and let X_{1, m}< × 2, m <… < X_{m, m} be the corresponding order statistics.
We assume m is an integer-valued random variable with P( m = k ) = p (1- p )^k-1, k = 1,2,… and 0 < p < 1. Two characterizations of the exponential distribution are given based on the distributional properties of X_{l, m}.     

Stability and self-decomposability of semi-group valued random variables

A random variable X on IR₊ is said to be self-decomposable, dif for all c∈ (0, 1) there exists a random variable X_c on IR₊ such that X=^dcX+X_c . It is said to be stable if it is self-decomposable and X_c=^d (1 - c)X' , where X and X' are identically and independently distributed. The notions of stability and self-decomposability for infinitely divisible random variables are generalised to abelian semi-groups ( S, + ) with S having an identical involution, by using characteristic functions. The generalised definitions involve semi-groups of scaling operators T . There operators can be interpreted in a slightly different context as generalised continuous-time branching processes (with immigration). The underlying importance of the generator of the semi-groups T in the characterisation of stability and self-decomposability is stressed.