Styjun

Let $(X_i)_i$ be a sequence of iid positive variables of mean 1 and variance $sigma^2$. Let $barX_n = fracsum_i=1^n X_in$.

My question is: Can we can bound $mathbbE(1/barX_n)$ as a function of $sigma$ and $n$?

There seems to be some strategy that may work based on the taylor extension, but

• I'm not sure about the hypothesis that need to be met;


• if it works in this case; and


• if we can say something definite on $barX_n$ or if we need to use the central limit theorem and can only say this for the normal approximation?

More details about the Taylor expansion. According to this wikipedia article,
$$mathbbE(f(X)) approx f(mu_X) +fracf''(mu_X)2sigma_X^2$$

So in my case it would give something like:
$$mathbbE(1/barX_n) approx 1 +fracsigma^24 n$$
I'm trying to find maybe a formal proof of a similar result, or hypothesis so that it works. Maybe references?
Thanks

EDIT: if needed, we can consider that the $(X_i)_i$ are discrete, there exists $v_1<cdots<v_K$ such that $mathbbP(X=v_k)=p_k$ and $sum p_k = 1$. In this case we know that $barX_n geq v_1$. Although I believe something can be said in the general case.

PS: this is almost a cross-post of this on Math.SE.

You cannot bound that expectation in $sigma, n$. That's because there is the distinct possibility that the expectation do not exist at all (or, is $infty$.) See I've heard that ratios or inverses of random variables often are problematic, in not having expectations. Why is that?. If the conditions given there is fulfilled for the density of $X_1$, it will so be for the density of $barX_n$. If densities do not exist, but probability mass functions do, it is simpler, since your assumptions prohibit a probability atom at zero, but a probability density can still be positive at zero even if $P(X >0)=1$.

For a useful bound you will at least need to restrict the common distribution of $X_1, dotsc, X_n$ much more.

EDIT

After your new information, and with $v_1>0$, the expectation of $1/barX_n$ certainly will exist (irrespective if $K$ is finite or not.) And, since the function $xmapsto 1/x$ is convex for $x>0$, we can use the Jensen Inequality to conclude that $DeclareMathOperatorEmathbbEE 1/barX_n ge 1/E barX_n$.

I think I have the gist of it.
Given that $f(x)=1/x$ is infinitely differentiable in 1. Taylor's theorem tells us:

There exists $varepsilon>0$ such that $f(x) = f(1) + f'(1) (x-1)+ fracf''(1)(x-1)^22 + fracf'''(varepsilon) (x-1)^22$.

In our case, if $X_i$ belongs in the domaine $[v_1;+infty[$, then $barX_n$ has the same domain and we have $varepsilon geq v_1$.

Hence
$mathbbE(1/barX_n) = mathbbEleft (1 - (barX_n-1) + frac(x-1)^24+ fracf'''(varepsilon) (barX_n-1)^22 right)$, and
beginalign*
mathbbE(1/barX_n) &= 1 + fracf'''(varepsilon) mathbbEleft ((barX_n-1)^2right )2 = 1 +frac V(barX_n)4 - frac V(barX_n)12 varepsilon^4\
endalign*
and hence
$$1 + frac sigma^24 n- fracsigma^212 v_1^4 n leq mathbbE(1/barX_n) leq 1 + frac sigma^24 n.$$

For the case where $X_i$ do not admit a minimum but has an unlimited number of moments, one can do a similar transformation using the full taylor expansion:

beginalign*
mathbbE(1/barX_n) &= sum_i=0^+infty fracf^(i)(1)i!mathbbEleft((barX_n-1)^iright)\
&= sum_i=0^+infty frac(-1)^ii!i!mathbbEleft((barX_n-1)^iright)
endalign*

Now if we can say something about the $k^th$ moment of $tildeX_n$ being $O(1/n^k/2)$ this validates that $mathbbE(1/barX_n) approx 1 + frac sigma^24 n$.