Suma limitante de iid Gamma variantes

Sea una secuencia de variables aleatorias distribuidas de forma independiente e idéntica con la función de densidad de probabilidad; Muestre que $X_1,X_2,\ldots$

f (x) = {\begin{cases} \frac{1}{2} x^{2} e^{- x} & if x > 0; \\ 0 & otherwise . \end{cases}

$f(x) = \left\{ \begin{array}{ll} \frac{1}{2}x^2 e^{-x} & \mbox{if $x>0$};\\ 0 & \mbox{otherwise}.\end{array} \right.$

lim_{n \to \infty} P [X_{1} + X_{2} + \dots + X_{n} \geq 3 (n - \sqrt{n})] \geq \frac{1}{2}

$\lim_{n\to \infty} P[X_1+X_2+\ldots+X_n\ge 3(n-\sqrt{n})] \ge \frac{1}{2}$

Lo que intenté

A primera vista, pensé que debería usar la desigualdad de Chebyshev, ya que la pregunta es mostrar un límite inferior $X_1+X_2+\ldots +X_n$ . Sin embargo, pensé en el signo de límite que indica claramente que el problema puede estar relacionado de alguna manera con el Teorema del límite central (CLT)

Deje $S_n=X_1+X_2+\ldots +X_n$

E (S_{n}) = \sum_{i = 0}^{n} E (X_{i}) = 3 n (since E (X_{i}) = 3) V (S_{n}) = \sum_{i = 0}^{n} V (X_{i}) = 3 n (since V (X_{i}) = 3 and X_{i} are i.i.d)

$E(S_n)=\sum_{i=0}^{n} E(X_i)=3n \ (\text{since } E(X_i)=3) \\ V(S_n)=\sum_{i=0}^{n} V(X_i)=3n \ (\text{since } V(X_i)=3 \text{ and } X_i \text{ are i.i.d})$

Ahora, usando CLT, para grande $n$ , $X_1+X_2+........+X_n \sim N(3n,3n)$
O,

z = \frac{S_{n} - 3 n}{\sqrt{3 n}} \sim N (0, 1) as n \to \infty

$z=\frac{S_n-3n}{\sqrt{3n}} \sim N(0,1) \text{ as } n\to \infty$

Ahora,

lim_{n \to \infty} P [X_{1} + X_{2} + . . . . . . . . + X_{n} \geq 3 (n - \sqrt{n})] = lim_{n \to \infty} P (S_{n} - 3 n \geq - 3 \sqrt{n}) = lim_{n \to \infty} P (\frac{S_{n} - 3 n}{\sqrt{3 n}} \geq - \sqrt{3}) = P (z \geq - \sqrt{3}) = P (- \sqrt{3} \leq z < 0) + P (z \geq 0) = P (- \sqrt{3} \leq z < 0) + \frac{1}{2} \dots (1)

$\lim_\limits{n\to \infty} P[X_1+X_2+........+X_n\ge 3(n-\sqrt{n})] = \lim_\limits{n\to \infty}P(S_n-3n \ge -3\sqrt{n}) = \lim_\limits{n\to \infty} P\left(\frac{S_n-3n}{\sqrt{3n}} \ge -\sqrt{3}\right) =P(z\ge -\sqrt{3}) =P\left(-\sqrt{3}\le z<0\right)+P(z\ge 0 ) =P\left(-\sqrt{3}\le z<0\right)+\frac{1}{2}\cdots(1)$

Desde $P(-\sqrt{3}\le z<0) \ge 0$ , por lo tanto desde $(1)$ ,

lim_{n \to \infty} P [X_{1} + X_{2} + . . . . . . . . + X_{n} \geq 3 (n - \sqrt{n})] \geq \frac{1}{2}

$\lim_\limits{n\to \infty} P[X_1+X_2+........+X_n\ge 3(n-\sqrt{n})]\ge \frac{1}{2}$

¿Estoy en lo correcto?

probability self-study

CLT parece un enfoque razonable pero " "no tiene sentido ..

lim_{n \to \infty} P [X_{1} + X_{2} + . . . . . . . . + X_{n} \geq 3 (n - \sqrt{n})] = P (S_{n} - 3 n \geq - 3 \sqrt{n})

$\lim_\limits{n\to \infty} P[X_1+X_2+........+X_n\ge 3(n-\sqrt{n})] =P(S_n-3n \ge -3\sqrt{n})$

— P.Windridge

Creo que debería ser

lim_{n \to \infty} P [X_{1} + X_{2} + . . . . . . . . + X_{n} \geq 3 (n - \sqrt{n})] = lim_{n \to \infty} P (S_{n} - 3 n \geq - 3 \sqrt{n}) = lim_{n \to \infty} P (\frac{S_{n} - 3 n}{\sqrt{3 n}} \geq - \sqrt{3}) = P (z \geq - \sqrt{3})

Como alternativa, considere que iid y así . La mediana de una variable aleatoria Gamma no se conoce en forma cerrada pero se conoce (cf. Wikipedia ) que para grandes , la mediana de un mentiras variable aleatoria entre y . Dado que , debe ser que al menos la mitad de la masa de probabilidad se encuentra a la derecha de .

X_{i} \sim Γ (3, 1)

$X_i \sim\Gamma(3,1)$

X_{1} + X_{2} + \cdot + X_{n} \sim Γ (3 n, 1)

$X_1+X_2+\cdot+X_n \sim \Gamma(3n,1)$

n

$n$

Γ (3 n, 1)

$\Gamma(3n,1)$

3 n - \frac{1}{3}

$3n-\frac 13$

3 n

$3n$

3 (n - \sqrt{n}) < 3 n - \frac{1}{3}

$3(n-\sqrt{n}) < 3n-\frac 13$

3 (n - \sqrt{n})

$3(n-\sqrt{n})$

— Dilip Sarwate

Respuestas:

Tenías razón en que la desigualdad de Chebyshev funcionará. Proporciona un límite algo crudo pero efectivo que se aplica a muchas de esas secuencias, revelando que la característica crucial de esta secuencia es que la varianza de las sumas parciales crece como máximo linealmente con . $n$

Consideremos, entonces, el caso extremadamente general de cualquier secuencia de variables no correlacionadas con medias y variaciones finitas Deje ser la suma de los primeros de ellos, $X_i$ $\mu_i$ $\sigma_i^2.$ $Y_n$ $n$

Y_{n} = \sum_{i = 1}^{n} X_{i} .

$Y_n = \sum_{i=1}^n X_i.$

En consecuencia, la media de es $Y_n$

m_{n} = \sum_{i = 1}^{n} μ_{n}

$m_n = \sum_{i=1}^n \mu_n$

y su varianza es

s_{n}^{2} = Var (Y_{n}) = \sum_{i = 1}^{n} Var (X_{i}) + 2 \sum_{j > i} Cov (X_{i}, X_{j}) = \sum_{i = 1}^{n} σ_{i}^{2} .

$s_n^2 = \operatorname{Var}(Y_n) = \sum_{i=1}^n \operatorname{Var}(X_i) + 2\sum_{j > i}\operatorname{Cov}(X_i,X_j) = \sum_{i=1}^n \sigma_i^2.$

Supongamos que crece como máximo linealmente con : $s_n^2$ $n$ es decir, existe un número tal que para todos los suficientemente grandes Deje (aún por determinar), observe que $\lambda\gt 0$ $n,$ $s_n^2 \le \lambda^2 n.$ $k\gt 0$

m - k \sqrt{n} \leq m - \frac{k}{λ} s_{n},

$m - k \sqrt{n} \le m - \frac{k}{\lambda}s_n,$

y aplicar la desigualdad de Chebyshev a para obtener $Y_n$

\begin{aligned} Pr (Y_{n} \geq m_{n} - k \sqrt{n}) & \geq Pr (Y_{n} \geq m_{n} - \frac{k}{λ} s_{n}) \\ \geq Pr (| Y_{n} - m_{n} | \leq \frac{k}{λ} s_{n}) \\ \geq 1 - \frac{λ^{2}}{k^{2}} . \end{aligned}

$\eqalign{ \Pr(Y_n \ge m_n - k\sqrt{n}) &\ge \Pr(Y_n \ge m_n - \frac{k}{\lambda}s_n) \\ &\ge \Pr(|Y_n - m_n| \le \frac{k}{\lambda} s_n) \\ &\ge 1 - \frac{\lambda^2}{k^2}. }$

Las dos primeras desigualdades son básicas: siguen porque cada evento sucesivo es un subconjunto del anterior.

En el caso en cuestión, donde es independiente (y por lo tanto no correlacionado) con los medios y las varianzas tenemos y $X_i$ $\mu_i=3$ $\sigma_i^2=3,$ $m_n=3n$

s_{n} = \sqrt{3} \sqrt{n},

$s_n = \sqrt{3}\sqrt{n},$

de donde podemos tomar tan pequeño como El evento en la pregunta corresponde a donde $\lambda$ $\sqrt{3}.$ $3(n-\sqrt{n}) = \mu_n - 3\sqrt{n}$ $k=3,$

Pr (Y_{n} \geq 3 n - 3 \sqrt{n}) \geq 1 - \frac{{\sqrt{3}}^{2}}{3^{2}} = \frac{2}{3} > \frac{1}{2},

$\Pr(Y_n \ge 3n - 3\sqrt{n}) \ge 1 - \frac{\sqrt{3}^{\ 2}}{3^2} = \frac{2}{3}\gt \frac{1}{2},$

QED

— whuber
fuente

Como alternativa a la excelente respuesta de Whuber, intentaré derivar el límite exacto de la probabilidad en cuestión. Una de las propiedades de la distribución gamma es que las sumas de variables aleatorias gamma independientes con el mismo parámetro de velocidad / escala también son variables aleatorias gamma con forma igual a la suma de las formas de esas variables. (Eso se puede probar fácilmente utilizando las funciones generadoras de la distribución). En el presente caso tenemos , por lo que obtenemos la suma: $X_1,...X_n \sim \text{IID Gamma}(3,1)$

S_{n} \equiv X_{1} + \dots + X_{n} \sim Gamma (3 n, 1) .

$S_n \equiv X_1 + \cdots + X_n \sim \text{Gamma}(3n, 1).$

Por lo tanto, podemos escribir la probabilidad exacta de interés utilizando el CDF de la distribución gamma. Si denota el parámetro de forma y denota el argumento de interés, tenemos: $a = 3n$ $x = 3(n-\sqrt{n})$

\begin{aligned} H (n) & \equiv P (S_{n} \geq 3 (n - \sqrt{n})) \\ = \frac{Γ (a, x)}{Γ (a)} \\ = \frac{a Γ (a)}{a Γ (a) + x^{a} e^{- x}} \cdot \frac{Γ (a + 1, x)}{Γ (a + 1)} . \end{aligned}

$\begin{equation} \begin{aligned} H(n) &\equiv \mathbb{P}(S_n \geq 3(n-\sqrt{n})) \\[12pt] &= \frac{\Gamma(a, x)}{\Gamma(a)} \\[6pt] &= \frac{a \Gamma(a)}{a \Gamma(a) + x^a e^{-x}} \cdot \frac{\Gamma(a+1, x)}{\Gamma(a+1)}. \\[6pt] \end{aligned} \end{equation}$

Para encontrar el límite de esta probabilidad, primero notamos que podemos escribir el segundo parámetro en términos del primero como donde . Usando un resultado que se muestra en Temme (1975) (Ec. 1.4, p. 1109) tenemos la equivalencia asintótica: $x = a + \sqrt{2a} \cdot y$ $y = -\sqrt{3/2}$

\begin{aligned} \frac{Γ (a + 1, x)}{Γ (a + 1)} & \sim \frac{1}{2} + \frac{1}{2} \cdot erf (- y) + \sqrt{\frac{2}{9 a π}} (1 + y^{2}) \exp (- y^{2}) . \end{aligned}

$\begin{aligned} \frac{\Gamma(a+1, x)}{\Gamma(a+1)} &\sim \frac{1}{2} + \frac{1}{2} \cdot \text{erf}(-y) + \sqrt{\frac{2}{9a \pi}} (1+y^2) \exp( - y^2). \end{aligned}$

Usando la aproximación de Stirling y la definición limitante del número exponencial, también se puede demostrar que:

\begin{aligned} \frac{a Γ (a)}{a Γ (a) + x^{a} e^{- x}} & \sim \frac{\sqrt{2 π} \cdot a \cdot (a - 1)^{a - 1 / 2}}{\sqrt{2 π} \cdot a \cdot (a - 1)^{a - 1 / 2} + x^{a} \cdot e^{a - x - 1}} \\ = \frac{\sqrt{2 π} \cdot a \cdot (1 - \frac{1}{a})^{a - 1 / 2}}{\sqrt{2 π} \cdot a \cdot (1 - \frac{1}{a})^{a - 1 / 2} + \sqrt{x} \cdot (\frac{x}{a})^{a - 1 / 2} \cdot e^{a - x - 1}} \\ = \frac{\sqrt{2 π} \cdot a \cdot e^{- 1}}{\sqrt{2 π} \cdot a \cdot e^{- 1} + \sqrt{x} \cdot e^{x - a} \cdot e^{a - x - 1}} \\ = \frac{\sqrt{2 π} \cdot a}{\sqrt{2 π} \cdot a + \sqrt{x}} \\ \sim \frac{\sqrt{2 π a}}{\sqrt{2 π a} + 1} . \end{aligned}

$\begin{aligned} \frac{a \Gamma(a)}{a \Gamma(a) + x^a e^{-x}} &\sim \frac{\sqrt{2 \pi} \cdot a \cdot (a-1)^{a-1/2}}{\sqrt{2 \pi} \cdot a \cdot (a-1)^{a-1/2} + x^a \cdot e^{a-x-1}} \\[6pt] &= \frac{\sqrt{2 \pi} \cdot a \cdot (1-\tfrac{1}{a})^{a-1/2}}{\sqrt{2 \pi} \cdot a \cdot (1-\tfrac{1}{a})^{a-1/2} + \sqrt{x} \cdot (\tfrac{x}{a})^{a-1/2} \cdot e^{a-x-1}} \\[6pt] &= \frac{\sqrt{2 \pi} \cdot a \cdot e^{-1}}{\sqrt{2 \pi} \cdot a \cdot e^{-1} + \sqrt{x} \cdot e^{x-a} \cdot e^{a-x-1}} \\[6pt] &= \frac{\sqrt{2 \pi} \cdot a}{\sqrt{2 \pi} \cdot a + \sqrt{x}} \\[6pt] &\sim \frac{\sqrt{2 \pi a}}{\sqrt{2 \pi a} + 1}. \\[6pt] \end{aligned}$

Sustituyendo los valores relevantes, obtenemos:

\begin{aligned} H (n) & = \frac{a Γ (a)}{a Γ (a) + x^{a} e^{- x}} \cdot \frac{Γ (a + 1, x)}{Γ (a + 1)} \\ \sim \frac{\sqrt{2 π a}}{\sqrt{2 π a} + 1} \cdot [\frac{1}{2} + \frac{1}{2} \cdot erf (\sqrt{\frac{3}{2}}) + \sqrt{\frac{2}{9 a π}} \cdot \frac{5}{2} \cdot \exp (\frac{3}{2})] . \end{aligned}

$\begin{equation} \begin{aligned} H(n) &= \frac{a \Gamma(a)}{a \Gamma(a) + x^a e^{-x}} \cdot \frac{\Gamma(a+1, x)}{\Gamma(a+1)} \\[6pt] &\sim \frac{\sqrt{2 \pi a}}{\sqrt{2 \pi a} + 1} \cdot \Bigg[ \frac{1}{2} + \frac{1}{2} \cdot \text{erf} \Big( \sqrt{\frac{3}{2}} \Big) + \sqrt{\frac{2}{9a \pi}} \cdot \frac{5}{2} \cdot \exp \Big( \frac{3}{2} \Big) \Bigg]. \\[6pt] \end{aligned} \end{equation}$

Esto nos da el límite:

lim_{n \to \infty} H (n) = \frac{1}{2} + \frac{1}{2} \cdot erf (\sqrt{\frac{3}{2}}) = 0.9583677.

$\lim_{n \rightarrow \infty} H(n) = \frac{1}{2} + \frac{1}{2} \cdot \text{erf} \Big( \sqrt{\frac{3}{2}} \Big) = 0.9583677.$

Esto nos da el límite exacto de la probabilidad de interés, que es mayor que la mitad.

— Ben - Restablece a Monica
fuente