¿Cómo se calcula la expectativa de


12

Si se distribuye exponencialmente con el parámetro y son mutuamente independientes, ¿cuál es la expectativa de ( i = 1 , . . . , N ) λ X iXi(i=1,...,n)λXi

(i=1nXi)2

en términos de y y posiblemente otras constantes?λnλ

Nota: Esta pregunta obtuvo una respuesta matemática en /math//q/12068/4051 . Los lectores también lo echarían un vistazo.


55
Las dos copias de esta pregunta se refieren entre sí y, apropiadamente, el sitio de estadísticas (aquí) tiene una respuesta estadística y el sitio de matemáticas tiene una respuesta matemática. Parece una buena división: ¡déjelo en pie!
Whuber

Respuestas:


31

Si , entonces (bajo independencia), , entonces se distribuye gamma (ver wikipedia ). Entonces, solo necesitamos . Como , sabemos que . Por lo tanto, (consulte wikipedia para conocer las expectativas y la variación de la distribución gamma).y = x iG a m m a ( n , 1 / λ ) y E [ y 2 ] V a r [ y ] = E [ y 2 ] - E [ y ] 2 E [ y 2 ] = V a r [xiExp(λ)y=xiGamma(n,1/λ)yE[y2]Var[y]=E[y2]E[y]2 E [ y 2 ] = n / λ 2 + n 2 / λ 2 = n ( 1 + n ) / λ 2E[y2]=Var[y]+E[y]2E[y2]=n/λ2+n2/λ2=n(1+n)/λ2


Gracias. También se proporcionó una forma muy ordenada de responder la pregunta (que conduce a la misma respuesta) en math.stackexchange (enlace anterior en la pregunta) hace unos minutos.
Wolfgang el

2
La respuesta matemática calcula las integrales utilizando la linealidad de la expectativa. De alguna manera es más simple. Pero me gusta su solución porque explota el conocimiento estadístico : porque sabe que una suma de variables exponenciales independientes tiene una distribución Gamma, ya está.
whuber

1
Lo disfruté bastante y de ninguna manera soy estadístico o matemático.
Kortuk

Respuesta muy elegante.
Cyrus S

1
@Dilip El matemático tiende a ver esta pregunta como una pregunta integral y procede directamente a integrarla. El estadístico lo reexpresa en términos de cantidades estadísticas familiares, como la varianza, y relaciones estadísticas familiares, como que el Exponencial es Gamma y la familia Gamma está cerrada por convolución. Las respuestas son las mismas, pero los enfoques son completamente diferentes. Luego está la cuestión de lo que realmente significa "hacer una integración". Por ejemplo, esta integral complicada se realiza puramente algebraicamente.
whuber

9

La respuesta anterior es muy agradable y responde completamente a la pregunta, pero, en cambio, proporcionaré una fórmula general para el cuadrado esperado de una suma y la aplicaré al ejemplo específico mencionado aquí.

Para cualquier conjunto de constantes es un hecho quea1,...,an

(i=1nai)2=i=1nj=1naiaj

Esto es cierto por la propiedad distributiva y queda claro cuando considera lo que está haciendo cuando calcula a mano.(a1+...+an)(a1+...+an)

Por lo tanto, para una muestra de variables aleatorias , independientemente de las distribuciones,X1,...,Xn

E([i=1nXi]2)=E(i=1nj=1nXiXj)=i=1nj=1nE(XiXj)

siempre que existan estas expectativas.

En el ejemplo del problema, son variables aleatorias iid , lo que nos dice que y para cada . Por independencia, para , tenemosX1,...,Xnexponential(λ)E(Xi)=1/λvar(Xi)=1/λ2iij

E(XiXj)=E(Xi)E(Xj)=1λ2

Hay de estos términos en la suma. Cuando , tenemosn2ni=j

E(XiXj)=E(Xi2)=var(Xi)+E(Xi)2=2λ2

y hay de estos términos en la suma. Por lo tanto, usando la fórmula anterior,n

E(i=1nXi)2=i=1nj=1nE(XiXj)=(n2n)1λ2+n2λ2=n2+nλ2

es tu respuesta


3

Este problema es solo un caso especial del problema mucho más general de 'momentos de momentos' que generalmente se definen en términos de notación de suma de potencia. En particular, en notación de suma de potencia:

s1=i=1nXi

Luego, independientemente de la distribución , el póster original busca (siempre que existan los momentos). Dado que el operador de expectativas es solo el primer Momento sin procesar, la solución se da en el software mathStatica por:E[s12]

ingrese la descripción de la imagen aquí

[El '___ToRaw' significa que queremos que la solución se presente en términos de momentos crudos de la población (en lugar de decir momentos centrales o acumulativos). ]

Finalmente, si ~ Exponencial ( ) con pdf :λ f ( x )Xλf(x)

f = Exp[-x/λ]/λ;      domain[f] = {x, 0, ∞} &&  > 0};

entonces podemos reemplazar los momentos en la solución general con los valores reales para una variable aleatoria exponencial, así:μisol

ingrese la descripción de la imagen aquí

Todo listo.


PD La razón por la que las otras soluciones publicadas aquí producen una respuesta con en el denominador en lugar del numerador es, por supuesto, porque están usando una parametrización diferente de la distribución exponencial. Como el OP no indicó qué versión estaba usando, decidí usar la definición de libro de texto de teoría de distribución estándar Johnson Kotz et al ... solo para equilibrar las cosas :)λ2

Al usar nuestro sitio, usted reconoce que ha leído y comprende nuestra Política de Cookies y Política de Privacidad.
Licensed under cc by-sa 3.0 with attribution required.