Articulo de referencia

Media aritmética ponderada

La media aritmética ponderada es similar a la media aritmética ordinaria (el tipo de promedio más común ), con la diferencia de que, en lugar de que cada dato contribuya por igu...

La media aritmética ponderada es similar a la media aritmética ordinaria (el tipo de promedio más común ), con la diferencia de que, en lugar de que cada dato contribuya por igual al promedio final, algunos datos contribuyen más que otros. El concepto de media ponderada se utiliza en estadística descriptiva y también aparece de forma más general en otras áreas de las matemáticas .

Si todos los pesos son iguales, entonces la media ponderada es igual a la media aritmética . Si bien las medias ponderadas generalmente se comportan de manera similar a las medias aritméticas, tienen algunas propiedades contraintuitivas, como se refleja, por ejemplo, en la paradoja de Simpson .

Ejemplos

Ejemplo básico

Dadas dos clases escolares —una con 20 alumnos y otra con 30— y las calificaciones de las pruebas en cada clase son las siguientes:

Clase matutina = {62, 67, 71, 74, 76, 77, 78, 79, 79, 80, 80, 81, 81, 82, 83, 84, 86, 89, 93, 98}

Clase de la tarde = {81, 82, 83, 84, 85, 86, 87, 87, 88, 88, 89, 89, 89, 90, 90, 90, 90, 91, 91, 91, 92, 92, 93, 93, 94, 95, 96, 97, 98, 99}

La media de la clase de la mañana es 80 y la de la clase de la tarde es 90. La media simple de ambas es 85. Sin embargo, esto no tiene en cuenta la diferencia en el número de alumnos por clase (20 frente a 30); por lo tanto, el valor de 85 no refleja la nota media de los alumnos (independientemente de la clase). La nota media de los alumnos se puede obtener promediando todas las notas, sin importar la clase (sumar todas las notas y dividirlas por el número total de alumnos). incógnita¯=430050=86.{\displaystyle {\bar {x}}={\frac {4300}{50}}=86.}

O bien, esto se puede lograr ponderando las medias de las clases según el número de alumnos en cada una. A la clase más numerosa se le da mayor "peso":

incógnita¯=(20×80)+(30×90)20+30=86.{\displaystyle {\bar {x}}={\frac {(20\times 80)+(30\times 90)}{20+30}}=86.}

De este modo, la media ponderada permite calcular la media de las calificaciones de los estudiantes sin necesidad de conocer la nota individual de cada uno. Solo se requieren las medias de las clases y el número de alumnos en cada una.

Ejemplo de combinación convexa

Dado que solo los pesos relativos son relevantes, cualquier media ponderada puede expresarse mediante coeficientes que suman uno. Dicha combinación lineal se denomina combinación convexa .

Utilizando el ejemplo anterior, obtendríamos los siguientes pesos:

2020+30=0,4{\displaystyle {\frac {20}{20+30}}=0,4}
3020+30=0,6{\displaystyle {\frac {30}{20+30}}=0,6}

Luego, aplica los pesos de esta manera:

incógnita¯=(0,4×80)+(0,6×90)=86.{\displaystyle {\bar {x}}=(0.4\times 80)+(0.6\times 90)=86.}

Definición matemática

Formalmente, la media ponderada de una tupla finita no vacía de datos(incógnita1,incógnita2,,incógnitanorte){\displaystyle \left(x_{1},x_{2},\dots ,x_{n}\right)}, con pesos no negativos correspondientes(w1,w2,,wnorte){\displaystyle \left(w_{1},w_{2},\dots ,w_{n}\right)}es

incógnita¯=i=1nortewiincógnitaii=1nortewi,{\displaystyle {\bar {x}}={\frac {\sum \limits _{i=1}^{n}w_{i}x_{i}}{\sum \limits _{i=1}^{n}w_{i}}},}

que se expande a:

incógnita¯=w1incógnita1+w2incógnita2++wnorteincógnitanortew1+w2++wnorte.{\displaystyle {\bar {x}}={\frac {w_{1}x_{1}+w_{2}x_{2}+\cdots +w_{n}x_{n}}{w_{1}+w_{2}+\cdots +w_{n}}}.}

Por lo tanto, los elementos de datos con un peso alto contribuyen más a la media ponderada que los elementos con un peso bajo. Los pesos no pueden ser negativos para que la ecuación funcione. [ a ] ​​Algunos pueden ser cero, pero no todos (ya que no se permite la división por cero).

Las fórmulas se simplifican cuando los pesos se normalizan de manera que sumen 1, es decir,i=1nortewi=1{\textstyle \sum \limits _{i=1}^{n}{w_{i}'}=1}Para tales pesos normalizados, la media ponderada es equivalentemente:

incógnita¯=i=1nortewiincógnitai{\displaystyle {\bar {x}}=\sum \limits _{i=1}^{n}{w_{i}'x_{i}}}.

Siempre se pueden normalizar los pesos realizando la siguiente transformación sobre los pesos originales:

wi=wij=1nortewj{\displaystyle w_{i}'={\frac {w_{i}}{\sum \limits _{j=1}^{n}{w_{j}}}}}.

El significado ordinario1nortei=1norteincógnitai{\textstyle {\frac {1}{n}}\sum \limits _{i=1}^{n}{x_{i}}}es un caso especial de la media ponderada donde todos los datos tienen el mismo peso.

Si los elementos de datos son variables aleatorias independientes e idénticamente distribuidas con varianzaσ2{\displaystyle \sigma ^{2}}, el error estándar de la media ponderada ,σincógnita¯{\displaystyle \sigma _{\bar {x}}}, se puede demostrar mediante la propagación de la incertidumbre que es:

σincógnita¯=σi=1nortewi2{\textstyle \sigma _{\bar {x}}=\sigma {\sqrt {\sum \limits _{i=1}^{n}w_{i}'^{2}}}}

Pesos definidos por la varianza

Para la media ponderada de una lista de datos para la cual cada elementoincógnitai{\displaystyle x_{i}}potencialmente proviene de una distribución de probabilidad diferente con varianza conocidaσi2{\displaystyle \sigma _{i}^{2}}Dado que todos tienen la misma media, una posible elección para los pesos viene dada por el recíproco de la varianza:

wi=1σi2.{\displaystyle w_{i}={\frac {1}{\sigma _{i}^{2}}}.}

La media ponderada en este caso es:

incógnita¯=i=1norte(incógnitaiσi2)i=1norte1σi2=i=1norte(incógnitaiwi)i=1nortewi,{\displaystyle {\bar {x}}={\frac {\sum _{i=1}^{n}\left({\dfrac {x_{i}}{\sigma _{i}^{2}}}\right)}{\sum _{i=1}^{n}{\dfrac {1}{\sigma _{i}^{2}}}}}={\frac {\sum _{i=1}^{n}\left(x_{i}\cdot w_{i}\right)}{\sum _{i=1}^{n}w_{i}}},}

y el error estándar de la media ponderada (con ponderaciones de varianza inversa) es:

σincógnita¯=1i=1norteσi2=1i=1nortewi,{\displaystyle \sigma _{\bar {x}}={\sqrt {\frac {1}{\sum _{i=1}^{n}\sigma _{i}^{-2}}}}={\sqrt {\frac {1}{\sum _{i=1}^{n}w_{i}}}},}

Tenga en cuenta que esto se reduce aσincógnita¯2=σ02/norte{\displaystyle \sigma _{\bar {x}}^{2}=\sigma _{0}^{2}/n}cuando todosσi=σ0{\displaystyle \sigma _{i}=\sigma _{0}}Es un caso especial de la fórmula general de la sección anterior.

σincógnita¯2=i=1nortewi2σi2=i=1norteσi4σi2(i=1norteσi2)2.{\displaystyle \sigma _{\bar {x}}^{2}=\sum _{i=1}^{n}{w_{i}'^{2}\sigma _{i}^{2}}={\frac {\sum _{i=1}^{n}{\sigma _{i}^{-4}\sigma _{i}^{2}}}{\left(\sum _{i=1}^{n}\sigma _{i}^{-2}\right)^{2}}}.}

Las ecuaciones anteriores se pueden combinar para obtener:

incógnita¯=σincógnita¯2i=1norteincógnitaiσi2.{\displaystyle {\bar {x}}=\sigma _{\bar {x}}^{2}\sum _{i=1}^{n}{\frac {x_{i}}{\sigma _{i}^{2}}}.}

La importancia de esta elección radica en que esta media ponderada es el estimador de máxima verosimilitud de la media de las distribuciones de probabilidad, bajo el supuesto de que son independientes y siguen una distribución normal con la misma media.

Propiedades estadísticas

Expectativa

La media ponderada de la muestra,incógnita¯{\displaystyle {\bar {x}}}, es en sí misma una variable aleatoria. Su valor esperado y desviación estándar están relacionados con los valores esperados y las desviaciones estándar de las observaciones, como se indica a continuación. Para simplificar, asumimos ponderaciones normalizadas (ponderaciones que suman uno).

Si las observaciones tienen valores esperados mi(incógnitai)=μi,{\displaystyle E(x_{i})={\mu _{i}},} entonces la media muestral ponderada tiene expectativa mi(incógnita¯)=i=1nortewiμi.{\displaystyle E({\bar {x}})=\sum _{i=1}^{n}{w_{i}'\mu _{i}}.} En particular, si las medias son iguales,μi=μ{\displaystyle \mu _{i}=\mu }, entonces la esperanza de la media muestral ponderada será ese valor, mi(incógnita¯)=μ.{\displaystyle E({\bar {x}})=\mu .}

Diferencia

Caso simple iid

Cuando se tratan los pesos como constantes y se tiene una muestra de n observaciones de variables aleatorias no correlacionadas , todas con la misma varianza y esperanza (como es el caso de las variables aleatorias i.i.d. ), entonces la varianza de la media ponderada se puede estimar como la multiplicación de la varianza no ponderada por el efecto de diseño de Kish (ver demostración ):

Var(y¯w)=σ^y2w2¯w¯2{\displaystyle \operatorname {Var} ({\bar {y}}_{w})={\hat {\sigma }}_{y}^{2}{\frac {\overline {w^{2}}}{{\bar {w}}^{2}}}}

Conσ^y2=i=1norte(yiy¯)2norte1{\displaystyle {\hat {\sigma }}_{y}^{2}={\frac {\sum _{i=1}^{n}(y_{i}-{\bar {y}})^{2}}{n-1}}},w¯=i=1nortewinorte{\displaystyle {\bar {w}}={\frac {\sum _{i=1}^{n}w_{i}}{n}}}, yw2¯=i=1nortewi2norte{\displaystyle {\overline {w^{2}}}={\frac {\sum _{i=1}^{n}w_{i}^{2}}{n}}}

Sin embargo, esta estimación es bastante limitada debido a la fuerte suposición sobre las observaciones de y . Esto ha llevado al desarrollo de estimadores alternativos y más generales.

Perspectiva del muestreo de encuestas

Desde una perspectiva basada en modelos , nos interesa estimar la varianza de la media ponderada cuando los diferentesyi{\displaystyle y_{i}}no son variables aleatorias i.i.d. Una perspectiva alternativa para este problema es la de un diseño de muestreo arbitrario de los datos en el que las unidades se seleccionan con probabilidades desiguales (con reemplazo). [ 1 ] : 306

En la metodología de encuestas , la media poblacional de alguna cantidad de interés y se calcula tomando una estimación del total de y sobre todos los elementos de la población ( Y o a veces T ) y dividiéndola por el tamaño de la población, ya sea conocido (norte{\displaystyle N}) o estimado (norte^{\displaystyle {\hat {N}}}). En este contexto, cada valor de y se considera constante, y la variabilidad proviene del procedimiento de selección. Esto contrasta con los enfoques "basados ​​en modelos" en los que la aleatoriedad a menudo se describe en los valores de y. El procedimiento de muestreo de la encuesta produce una serie de valores indicadores de Bernoulli (Ii{\displaystyle I_{i}}) que obtienen 1 si alguna observación i está en la muestra y 0 si no fue seleccionada. Esto puede ocurrir con un tamaño de muestra fijo o con un tamaño de muestra variable (por ejemplo, muestreo de Poisson ). La probabilidad de que un elemento sea elegido, dada una muestra, se denota comoPAG(Ii=1Algunos tamaños de muestra norte)=πi{\displaystyle P(I_{i}=1\mid {\text{Some sample of size }}n)=\pi _{i}}y la probabilidad de selección en una sola extracción esPAG(Ii=1|una muestra extraída)=pagiπinorte{\displaystyle P(I_{i}=1|{\text{one sample draw}})=p_{i}\approx {\frac {\pi _{i}}{n}}}(Si N es muy grande y cadapagi{\displaystyle p_{i}}es muy pequeño). Para la siguiente derivación asumiremos que la probabilidad de seleccionar cada elemento está totalmente representada por estas probabilidades. [ 2 ] : 42, 43, 51 Es decir: seleccionar algún elemento no influirá en la probabilidad de extraer otro elemento (esto no se aplica a cosas como el diseño de muestreo por conglomerados ).

Dado que cada elemento (yi{\displaystyle y_{i}}) es fijo, y la aleatoriedad proviene de si se incluye o no en la muestra (Ii{\displaystyle I_{i}}), a menudo hablamos de la multiplicación de los dos, que es una variable aleatoria. Para evitar confusiones en la siguiente sección, llamaremos a este término:yi=yiIi{\displaystyle y'_{i}=y_{i}I_{i}}Con la siguiente expectativa:mi[yi]=yimi[Ii]=yiπi{\displaystyle E[y'_{i}]=y_{i}E[I_{i}]=y_{i}\pi _{i}}; y varianza:V[yi]=yi2V[Ii]=yi2πi(1πi){\displaystyle V[y'_{i}]=y_{i}^{2}V[I_{i}]=y_{i}^{2}\pi _{i}(1-\pi _{i})}.

Cuando cada elemento de la muestra se infla por el inverso de su probabilidad de selección, se denominaπ{\displaystyle \pi }-valores y expandidos , es decir:yˇi=yiπi{\displaystyle {\check {y}}_{i}={\frac {y_{i}}{\pi _{i}}}}. Una cantidad relacionada espag{\displaystyle p}-valores y expandidos :yipagi=norteyˇi{\displaystyle {\frac {y_{i}}{p_{i}}}=n{\check {y}}_{i}}. [ 2 ] : 42, 43, 51, 52 Como se indicó anteriormente, podemos agregar una marca de graduación si multiplicamos por la función indicadora . Es decir:yˇi=Iiyˇi=Iiyiπi{\displaystyle {\check {y}}'_{i}=I_{i}{\check {y}}_{i}={\frac {I_{i}y_{i}}{\pi _{i}}}}

En esta perspectiva basada en el diseño , los pesos, utilizados en el numerador de la media ponderada, se obtienen tomando el inverso de la probabilidad de selección (es decir, el factor de inflación). Es decir:wi=1πi1norte×pagi{\displaystyle w_{i}={\frac {1}{\pi _{i}}}\approx {\frac {1}{n\times p_{i}}}}.

Varianza de la suma ponderada ( estimador pwr para totales)

Si se conoce el tamaño de la población N, podemos estimar la media de la población utilizandoY¯^conocido norte=Y^pagwrnortei=1nortewiyinorte{\displaystyle {\hat {\bar {Y}}}_{{\text{known }}N}={\frac {{\hat {Y}}_{pwr}}{N}}\approx {\frac {\sum _{i=1}^{n}w_{i}y'_{i}}{N}}}.

Si el diseño de muestreo es uno que da como resultado un tamaño de muestra fijo n (como en el muestreo pps ), entonces la varianza de este estimador es:

Var(Y¯^conocido norte)=1norte2nortenorte1i=1norte(wiyiwy¯)2{\displaystyle \operatorname {Var} \left({\hat {\bar {Y}}}_{{\text{known }}N}\right)={\frac {1}{N^{2}}}{\frac {n}{n-1}}\sum _{i=1}^{n}\left(w_{i}y_{i}-{\overline {wy}}\right)^{2}}
Prueba

La fórmula general se puede desarrollar de la siguiente manera:

Y¯^conocido norte=Y^pagwrnorte=1nortei=1norteyipaginortei=1norteyiπinorte=i=1nortewiyinorte.{\displaystyle {\hat {\bar {Y}}}_{{\text{known }}N}={\frac {{\hat {Y}}_{pwr}}{N}}={\frac {{\frac {1}{n}}\sum _{i=1}^{n}{\frac {y'_{i}}{p_{i}}}}{N}}\approx {\frac {\sum _{i=1}^{n}{\frac {y'_{i}}{\pi _{i}}}}{N}}={\frac {\sum _{i=1}^{n}w_{i}y'_{i}}{N}}.}

La población total se denota comoY=i=1norteyi{\displaystyle Y=\sum _{i=1}^{N}y_{i}}y puede estimarse mediante el estimador de Horvitz-Thompson (insesgado) , también llamadoπ{\displaystyle \pi }-estimador. Este estimador puede ser estimado a su vez utilizando el estimador pwr (es decir:pag{\displaystyle p}-estimador expandido con reemplazo, o estimador de "probabilidad con reemplazo"). Con la notación anterior, es:Y^pagwr=1nortei=1norteyipagi=i=1norteyinortepagii=1norteyiπi=i=1nortewiyi{\displaystyle {\hat {Y}}_{pwr}={\frac {1}{n}}\sum _{i=1}^{n}{\frac {y'_{i}}{p_{i}}}=\sum _{i=1}^{n}{\frac {y'_{i}}{np_{i}}}\approx \sum _{i=1}^{n}{\frac {y'_{i}}{\pi _{i}}}=\sum _{i=1}^{n}w_{i}y'_{i}}. [ 2 ] : 51

La varianza estimada del estimador pwr viene dada por: [ 2 ] : 52Var(Y^pagwr)=nortenorte1i=1norte(wiyiwy¯)2{\displaystyle \operatorname {Var} ({\hat {Y}}_{pwr})={\frac {n}{n-1}}\sum _{i=1}^{n}\left(w_{i}y_{i}-{\overline {wy}}\right)^{2}} dóndewy¯=i=1nortewiyinorte{\displaystyle {\overline {wy}}=\sum _{i=1}^{n}{\frac {w_{i}y_{i}}{n}}}.

La fórmula anterior se tomó de Sarndal et al. (1992) (también presentada en Cochran 1977), pero se escribió de manera diferente. [ 2 ] : 52 [ 1 ] : 307 (11.35) El lado izquierdo muestra cómo se escribió la varianza y el lado derecho muestra cómo hemos desarrollado la versión ponderada:

Var(Y^potencia)=1norte1norte1i=1norte(yipagiY^pagwr)2=1norte1norte1i=1norte(nortenorteyipaginortenortei=1nortewiyi)2=1norte1norte1i=1norte(norteyiπinortei=1nortewiyinorte)2=norte2norte1norte1i=1norte(wiyiwy¯)2=nortenorte1i=1norte(wiyiwy¯)2{\displaystyle {\begin{aligned}\operatorname {Var} ({\hat {Y}}_{\text{pwr}})&={\frac {1}{n}}{\frac {1}{n-1}}\sum _{i=1}^{n}\left({\frac {y_{i}}{p_{i}}}-{\hat {Y}}_{pwr}\right)^{2}\\&={\frac {1}{n}}{\frac {1}{n-1}}\sum _{i=1}^{n}\left({\frac {n}{n}}{\frac {y_{i}}{p_{i}}}-{\frac {n}{n}}\sum _{i=1}^{n}w_{i}y_{i}\right)^{2}={\frac {1}{n}}{\frac {1}{n-1}}\sum _{i=1}^{n}\left(n{\frac {y_{i}}{\pi _{i}}}-n{\frac {\sum _{i=1}^{n}w_{i}y_{i}}{n}}\right)^{2}\\&={\frac {n^{2}}{n}}{\frac {1}{n-1}}\sum _{i=1}^{n}\left(w_{i}y_{i}-{\overline {wy}}\right)^{2}\\&={\frac {n}{n-1}}\sum _{i=1}^{n}\left(w_{i}y_{i}-{\overline {wy}}\right)^{2}\end{aligned}}}

Y así llegamos a la fórmula de arriba.

Un término alternativo, para cuando el muestreo tiene un tamaño de muestra aleatorio (como en el muestreo de Poisson ), se presenta en Sarndal et al. (1992) como: [ 2 ] : 182

Var(Y¯^potencia (conocida) norte))=1norte2i=1nortej=1norte(Δˇijyˇiyˇj){\displaystyle \operatorname {Var} ({\hat {\bar {Y}}}_{{\text{pwr (known }}N{\text{)}}})={\frac {1}{N^{2}}}\sum _{i=1}^{n}\sum _{j=1}^{n}\left({\check {\Delta }}_{ij}{\check {y}}_{i}{\check {y}}_{j}\right)}

Conyˇi=yiπi{\displaystyle {\check {y}}_{i}={\frac {y_{i}}{\pi _{i}}}}. También,do(Ii,Ij)=πijπiπj=Δij{\displaystyle C(I_{i},I_{j})=\pi _{ij}-\pi _{i}\pi _{j}=\Delta _{ij}}dóndeπij{\displaystyle \pi _{ij}}es la probabilidad de seleccionar tanto i como j. [ 2 ] : 36 YΔˇij=1πiπjπij{\displaystyle {\check {\Delta }}_{ij}=1-{\frac {\pi _{i}\pi _{j}}{\pi _{ij}}}}y para i=j:Δˇii=1πiπiπi=1πi{\displaystyle {\check {\Delta }}_{ii}=1-{\frac {\pi _{i}\pi _{i}}{\pi _{i}}}=1-\pi _{i}}. [ 2 ] : 43

Si las probabilidades de selección no están correlacionadas (es decir:ij:do(Ii,Ij)=0{\displaystyle \forall i\neq j:C(I_{i},I_{j})=0}), y suponiendo que la probabilidad de cada elemento es muy pequeña, entonces:

Var(Y¯^potencia (conocida) norte))=1norte2i=1norte(wiyi)2{\displaystyle \operatorname {Var} ({\hat {\bar {Y}}}_{{\text{pwr (known }}N{\text{)}}})={\frac {1}{N^{2}}}\sum _{i=1}^{n}\left(w_{i}y_{i}\right)^{2}}
Prueba

Suponemos que(1πi)1{\displaystyle (1-\pi _{i})\approx 1}y eso Var(Y^potencia (conocida) norte))=1norte2i=1nortej=1norte(Δˇijyˇiyˇj)=1norte2i=1norte(Δˇiiyˇiyˇi)=1norte2i=1norte((1πi)yiπiyiπi)=1norte2i=1norte(wiyi)2{\displaystyle {\begin{aligned}\operatorname {Var} ({\hat {Y}}_{{\text{pwr (known }}N{\text{)}}})&={\frac {1}{N^{2}}}\sum _{i=1}^{n}\sum _{j=1}^{n}\left({\check {\Delta }}_{ij}{\check {y}}_{i}{\check {y}}_{j}\right)\\&={\frac {1}{N^{2}}}\sum _{i=1}^{n}\left({\check {\Delta }}_{ii}{\check {y}}_{i}{\check {y}}_{i}\right)\\&={\frac {1}{N^{2}}}\sum _{i=1}^{n}\left((1-\pi _{i}){\frac {y_{i}}{\pi _{i}}}{\frac {y_{i}}{\pi _{i}}}\right)\\&={\frac {1}{N^{2}}}\sum _{i=1}^{n}\left(w_{i}y_{i}\right)^{2}\end{aligned}}}

Varianza de la media ponderada ( estimador π para la media de razón)

La sección anterior trató sobre la estimación de la media poblacional como una razón de un total poblacional estimado (Y^{\displaystyle {\hat {Y}}}) con un tamaño de población conocido (norte{\displaystyle N}), y la varianza se estimó en ese contexto. Otro caso común es que el tamaño de la población en sí (norte{\displaystyle N}) es desconocido y se estima utilizando la muestra (es decir:norte^{\displaystyle {\hat {N}}}). La estimación denorte{\displaystyle N}se puede describir como la suma de pesos. Entonces, cuandowi=1πi{\displaystyle w_{i}={\frac {1}{\pi _{i}}}}obtenemosnorte^=i=1nortewiIi=i=1norteIiπi=i=1norte1ˇi{\displaystyle {\hat {N}}=\sum _{i=1}^{n}w_{i}I_{i}=\sum _{i=1}^{n}{\frac {I_{i}}{\pi _{i}}}=\sum _{i=1}^{n}{\check {1}}'_{i}}Con la notación anterior, el parámetro que nos interesa es la razón de las sumas deyi{\displaystyle y_{i}}s y 1s. Es decir:R=Y¯=i=1norteyiπii=1norte1πi=i=1norteyˇii=1norte1ˇi=i=1nortewiyii=1nortewi{\displaystyle R={\bar {Y}}={\frac {\sum _{i=1}^{N}{\frac {y_{i}}{\pi _{i}}}}{\sum _{i=1}^{N}{\frac {1}{\pi _{i}}}}}={\frac {\sum _{i=1}^{N}{\check {y}}_{i}}{\sum _{i=1}^{N}{\check {1}}_{i}}}={\frac {\sum _{i=1}^{N}w_{i}y_{i}}{\sum _{i=1}^{N}w_{i}}}}Podemos estimarlo utilizando nuestra muestra con:R^=Y¯^=i=1norteIiyiπii=1norteIi1πi=i=1norteyˇii=1norte1ˇi=i=1nortewiyii=1nortewi1i=i=1nortewiyii=1nortewi1i=y¯w{\displaystyle {\hat {R}}={\hat {\bar {Y}}}={\frac {\sum _{i=1}^{N}I_{i}{\frac {y_{i}}{\pi _{i}}}}{\sum _{i=1}^{N}I_{i}{\frac {1}{\pi _{i}}}}}={\frac {\sum _{i=1}^{N}{\check {y}}'_{i}}{\sum _{i=1}^{N}{\check {1}}'_{i}}}={\frac {\sum _{i=1}^{N}w_{i}y'_{i}}{\sum _{i=1}^{N}w_{i}1'_{i}}}={\frac {\sum _{i=1}^{n}w_{i}y'_{i}}{\sum _{i=1}^{n}w_{i}1'_{i}}}={\bar {y}}_{w}}. Al pasar de usar N a usar n, sabemos que todas las variables indicadoras obtienen 1, por lo que podríamos simplemente escribir:y¯w=i=1nortewiyii=1nortewi{\displaystyle {\bar {y}}_{w}={\frac {\sum _{i=1}^{n}w_{i}y_{i}}{\sum _{i=1}^{n}w_{i}}}}Este será el valor estimado para valores específicos de y y w, pero las propiedades estadísticas se obtienen al incluir la variable indicadora.y¯w=i=1nortewiyii=1nortewi1i{\displaystyle {\bar {y}}_{w}={\frac {\sum _{i=1}^{n}w_{i}y'_{i}}{\sum _{i=1}^{n}w_{i}1'_{i}}}}. [ 2 ] : 162, 163, 176

Esto se denomina estimador de razón y es aproximadamente insesgado para R. [ 2 ] : 182

En este caso, la variabilidad de la razón depende de la variabilidad de las variables aleatorias tanto en el numerador como en el denominador, así como de su correlación. Dado que no existe una forma analítica cerrada para calcular esta varianza, se utilizan varios métodos para la estimación aproximada. Principalmente, la linealización de primer orden de la serie de Taylor , los métodos asintóticos y el bootstrap/jackknife. [ 2 ] : 172 El método de linealización de Taylor podría conducir a una subestimación de la varianza para tamaños de muestra pequeños en general, pero eso depende de la complejidad del estadístico. Para la media ponderada, se supone que la varianza aproximada es relativamente precisa incluso para tamaños de muestra medianos. [ 2 ] : 176 Cuando el muestreo tiene un tamaño de muestra aleatorio (como en el muestreo de Poisson ), es como sigue: [ 2 ] : 182

V(y¯w)^=1(i=1nortewi)2i=1nortewi2(yiy¯w)2{\displaystyle {\widehat {V({\bar {y}}_{w})}}={\frac {1}{(\sum _{i=1}^{n}w_{i})^{2}}}\sum _{i=1}^{n}w_{i}^{2}(y_{i}-{\bar {y}}_{w})^{2}}.

Siπipaginorte{\displaystyle \pi _{i}\approx p_{i}n}, entonces ya sea usandowi=1πi{\displaystyle w_{i}={\frac {1}{\pi _{i}}}}owi=1pagi{\displaystyle w_{i}={\frac {1}{p_{i}}}}daría el mismo estimador, ya que multiplicarwi{\displaystyle w_{i}}Por algún factor, se obtendría el mismo estimador. Esto también significa que si escalamos la suma de los pesos para que sea igual a un tamaño de población N conocido de antemano , el cálculo de la varianza sería el mismo. Cuando todos los pesos son iguales entre sí, esta fórmula se reduce al estimador de varianza insesgado estándar.

Prueba

La linealización de Taylor establece que para un estimador de razón general de dos sumas (R^=Y^Z^{\displaystyle {\hat {R}}={\frac {\hat {Y}}{\hat {Z}}}}), se pueden expandir alrededor del valor verdadero R, y dan: [ 2 ] : 178

R^=Y^Z^=i=1nortewiyii=1nortewiziR+1Zi=1norte(yiπiRziπi){\displaystyle {\hat {R}}={\frac {\hat {Y}}{\hat {Z}}}={\frac {\sum _{i=1}^{n}w_{i}y'_{i}}{\sum _{i=1}^{n}w_{i}z'_{i}}}\approx R+{\frac {1}{Z}}\sum _{i=1}^{n}\left({\frac {y'_{i}}{\pi _{i}}}-R{\frac {z'_{i}}{\pi _{i}}}\right)}

Y la varianza se puede aproximar mediante: [ 2 ] : 178, 179

V(R^)^=1Z^2i=1nortej=1norte(ΔˇijyiR^ziπiyjR^zjπj)=1Z^2[V(Y^)^+R^V(Z^)^2R^do^(Y^,Z^)]{\displaystyle {\widehat {V({\hat {R}})}}={\frac {1}{{\hat {Z}}^{2}}}\sum _{i=1}^{n}\sum _{j=1}^{n}\left({\check {\Delta }}_{ij}{\frac {y_{i}-{\hat {R}}z_{i}}{\pi _{i}}}{\frac {y_{j}-{\hat {R}}z_{j}}{\pi _{j}}}\right)={\frac {1}{{\hat {Z}}^{2}}}\left[{\widehat {V({\hat {Y}})}}+{\hat {R}}{\widehat {V({\hat {Z}})}}-2{\hat {R}}{\hat {C}}({\hat {Y}},{\hat {Z}})\right]}.

El términodo^(Y^,Z^){\displaystyle {\hat {C}}({\hat {Y}},{\hat {Z}})}es la covarianza estimada entre la suma estimada de Y y la suma estimada de Z. Dado que esta es la covarianza de dos sumas de variables aleatorias , incluiría muchas combinaciones de covarianzas que dependerán de las variables indicadoras. Si la probabilidad de selección no está correlacionada (es decir:ij:Δij=do(Ii,Ij)=0{\displaystyle \forall i\neq j:\Delta _{ij}=C(I_{i},I_{j})=0}), este término aún incluiría una suma de n covarianzas para cada elemento i entreyi=Iiyi{\displaystyle y'_{i}=I_{i}y_{i}}yzi=Iizi{\displaystyle z'_{i}=I_{i}z_{i}}Esto ayuda a ilustrar que esta fórmula incorpora el efecto de la correlación entre y y z sobre la varianza de los estimadores de razón.

Al definirzi=1{\displaystyle z_{i}=1}lo anterior se convierte en: [ 2 ] : 182

V(R^)^=V(y¯w)^=1norte^2i=1nortej=1norte(Δˇijyiy¯wπiyjy¯wπj).{\displaystyle {\widehat {V({\hat {R}})}}={\widehat {V({\bar {y}}_{w})}}={\frac {1}{{\hat {N}}^{2}}}\sum _{i=1}^{n}\sum _{j=1}^{n}\left({\check {\Delta }}_{ij}{\frac {y_{i}-{\bar {y}}_{w}}{\pi _{i}}}{\frac {y_{j}-{\bar {y}}_{w}}{\pi _{j}}}\right).}

Si las probabilidades de selección no están correlacionadas (es decir:ij:Δij=do(Ii,Ij)=0{\displaystyle \forall i\neq j:\Delta _{ij}=C(I_{i},I_{j})=0}), y cuando se supone que la probabilidad de cada elemento es muy pequeña (es decir:(1πi)1{\displaystyle (1-\pi _{i})\approx 1}), entonces lo anterior se reduce a lo siguiente: V(y¯w)^=1norte^2i=1norte((1πi)yiy¯wπi)2=1(i=1nortewi)2i=1nortewi2(yiy¯w)2.{\displaystyle {\widehat {V({\bar {y}}_{w})}}={\frac {1}{{\hat {N}}^{2}}}\sum _{i=1}^{n}\left((1-\pi _{i}){\frac {y_{i}-{\bar {y}}_{w}}{\pi _{i}}}\right)^{2}={\frac {1}{(\sum _{i=1}^{n}w_{i})^{2}}}\sum _{i=1}^{n}w_{i}^{2}(y_{i}-{\bar {y}}_{w})^{2}.}

Thomas Lumley proporcionó una recreación similar de la prueba (con algunos errores al final) en crossvalidated. [ 3 ]

Tenemos (al menos) dos versiones de la varianza para la media ponderada: una con estimación del tamaño de la población conocida y otra con estimación del tamaño de la población desconocida. No hay un enfoque uniformemente mejor, pero la literatura presenta varios argumentos para preferir el uso de la versión de estimación de la población (incluso cuando se conoce el tamaño de la población). [ 2 ] : 188 Por ejemplo: si todos los valores de y son constantes, el estimador con tamaño de población desconocido dará el resultado correcto, mientras que el que tiene tamaño de población conocido tendrá cierta variabilidad. Además, cuando el tamaño de la muestra en sí es aleatorio (por ejemplo, en el muestreo de Poisson ), la versión con media de población desconocida se considera más estable. Por último, si la proporción de muestreo está correlacionada negativamente con los valores (es decir, menor probabilidad de muestrear una observación que es grande), entonces la versión con tamaño de población desconocido compensa ligeramente eso.

Para el caso trivial en el que todos los pesos son iguales a 1, la fórmula anterior es igual que la fórmula habitual para la varianza de la media (pero observe que utiliza el estimador de máxima verosimilitud para la varianza en lugar de la varianza insesgada. Es decir: la divide por n en lugar de (n-1)).

Validación por arranque

Gatz et al. (1995) demostraron que, en comparación con los métodos de remuestreo (bootstrapping ) , la siguiente estimación de la varianza de la razón-media mediante la linealización de la serie de Taylor es razonable para el cuadrado del error estándar de la media (cuando se utiliza en el contexto de la medición de constituyentes químicos): [ 4 ] : 1186

σincógnita¯w2^=norte(norte1)(nortew¯)2[(wiincógnitaiw¯incógnita¯w)22incógnita¯w(wiw¯)(wiincógnitaiw¯incógnita¯w)+incógnita¯w2(wiw¯)2]{\displaystyle {\widehat {\sigma _{{\bar {x}}_{w}}^{2}}}={\frac {n}{(n-1)(n{\bar {w}})^{2}}}\left[\sum (w_{i}x_{i}-{\bar {w}}{\bar {x}}_{w})^{2}-2{\bar {x}}_{w}\sum (w_{i}-{\bar {w}})(w_{i}x_{i}-{\bar {w}}{\bar {x}}_{w})+{\bar {x}}_{w}^{2}\sum (w_{i}-{\bar {w}})^{2}\right]}

dóndew¯=winorte{\displaystyle {\bar {w}}={\frac {\sum w_{i}}{n}}}Una mayor simplificación conduce a

σincógnita¯2^=norte(norte1)(nortew¯)2wi2(incógnitaiincógnita¯w)2{\displaystyle {\widehat {\sigma _{\bar {x}}^{2}}}={\frac {n}{(n-1)(n{\bar {w}})^{2}}}\sum w_{i}^{2}(x_{i}-{\bar {x}}_{w})^{2}}

Gatz et al. mencionan que la formulación anterior fue publicada por Endlich et al. (1988) cuando trataron la media ponderada como una combinación de un estimador total ponderado dividido por un estimador del tamaño de la población, [ 5 ] basado en la formulación publicada por Cochran (1977), como una aproximación a la media de razón. Sin embargo, Endlich et al. no parecieron publicar esta derivación en su artículo (aunque mencionan que la usaron), y el libro de Cochran incluye una formulación ligeramente diferente. [ 1 ] : 155 Aun así, es casi idéntica a las formulaciones descritas en secciones anteriores.

Estimadores basados ​​en replicación

Debido a que no existe una forma analítica cerrada para la varianza de la media ponderada, se propuso en la literatura recurrir a métodos de replicación como el Jackknife y el Bootstrapping . [ 1 ] : 321

Otras notas

Para observaciones no correlacionadas con varianzasσi2{\displaystyle \sigma _{i}^{2}}, la varianza de la media muestral ponderada es

σincógnita¯2=i=1nortewi2σi2{\displaystyle \sigma _{\bar {x}}^{2}=\sum _{i=1}^{n}{w_{i}'^{2}\sigma _{i}^{2}}}

cuya raíz cuadradaσincógnita¯{\displaystyle \sigma _{\bar {x}}}se puede denominar error estándar de la media ponderada (caso general) .

En consecuencia, si todas las observaciones tienen la misma varianza,σi2=σ02{\displaystyle \sigma _{i}^{2}=\sigma _{0}^{2}}, la media muestral ponderada tendrá varianza

σincógnita¯2=σ02i=1nortewi2,{\displaystyle \sigma _{\bar {x}}^{2}=\sigma _{0}^{2}\sum _{i=1}^{n}{w_{i}'^{2}},}

dónde1/nortei=1nortewi21{\textstyle 1/n\leq \sum _{i=1}^{n}{w_{i}'^{2}}\leq 1}La varianza alcanza su valor máximo,σ02{\displaystyle \sigma _{0}^{2}}, cuando todos los pesos excepto uno son cero. Su valor mínimo se encuentra cuando todos los pesos son iguales (es decir, media no ponderada), en cuyo caso tenemosσincógnita¯=σ0/norte{\textstyle \sigma _{\bar {x}}=\sigma _{0}/{\sqrt {n}}}, es decir, degenera en el error estándar de la media , al cuadrado.

Dado que siempre se pueden transformar pesos no normalizados en pesos normalizados, todas las fórmulas de esta sección se pueden adaptar a pesos no normalizados reemplazando todoswi=wii=1nortewi{\displaystyle w_{i}'={\frac {w_{i}}{\sum _{i=1}^{n}{w_{i}}}}}.

Mediana

Cuando los pesos sondo{\displaystyle c}yd{\displaystyle d}, el promedio aritmético ponderado de dos fraccionesa/do{\displaystyle a/c}yb/d{\displaystyle b/d}es igual a su mediana .

Varianza de muestra ponderada

Normalmente, cuando se calcula una media, es importante conocer la varianza y la desviación estándar de esa media. Cuando se calcula una media ponderadaμ{\displaystyle \mu ^{*}}Si se utiliza, la varianza de la muestra ponderada es diferente de la varianza de la muestra no ponderada.

La varianza de muestra ponderada sesgadaσ^w2{\displaystyle {\hat {\sigma }}_{\mathrm {w} }^{2}}se define de forma similar a la varianza de muestra sesgada normalσ^2{\displaystyle {\hat {\sigma }}^{2}}:

σ^2 =i=1norte(incógnitaiμ)2norteσ^w2=i=1nortewi(incógnitaiμ)2i=1nortewi{\displaystyle {\begin{aligned}{\hat {\sigma }}^{2}\ &={\frac {\sum \limits _{i=1}^{N}\left(x_{i}-\mu \right)^{2}}{N}}\\{\hat {\sigma }}_{\mathrm {w} }^{2}&={\frac {\sum \limits _{i=1}^{N}w_{i}\left(x_{i}-\mu ^{*}\right)^{2}}{\sum _{i=1}^{N}w_{i}}}\end{aligned}}}

dóndei=1nortewi=1{\displaystyle \sum _{i=1}^{N}w_{i}=1}para pesos normalizados. Si los pesos son pesos de frecuencia (y por lo tanto son variables aleatorias), se puede demostrar queσ^w2{\displaystyle {\hat {\sigma }}_{\mathrm {w} }^{2}}es el estimador de máxima verosimilitud deσ2{\displaystyle \sigma ^{2}}para observaciones gaussianas i.i.d.

Para muestras pequeñas, es habitual utilizar un estimador insesgado para la varianza poblacional. En muestras normales no ponderadas, la N del denominador (que corresponde al tamaño de la muestra) se cambia a N  1 (véase la corrección de Bessel ). En el caso de ponderaciones, existen dos estimadores insesgados diferentes: uno para ponderaciones de frecuencia y otro para ponderaciones de fiabilidad .

ponderaciones de frecuencia

Si los pesos son pesos de frecuencia (donde un peso es igual al número de ocurrencias), entonces el estimador insesgado es:

s2 =i=1nortewi(incógnitaiμ)2i=1nortewi1{\displaystyle s^{2}\ ={\frac {\sum \limits _{i=1}^{N}w_{i}\left(x_{i}-\mu ^{*}\right)^{2}}{\sum _{i=1}^{N}w_{i}-1}}}

Esto aplica efectivamente la corrección de Bessel para ponderaciones de frecuencia. Por ejemplo, si los valores{2,2,4,5,5,5}{\displaystyle \{2,2,4,5,5,5\}}Si se extraen de la misma distribución, podemos tratar este conjunto como una muestra no ponderada o como una muestra ponderada.{2,4,5}{\displaystyle \{2,4,5\}}con pesos correspondientes{2,1,3}{\displaystyle \{2,1,3\}}y obtenemos el mismo resultado en ambos casos.

Si los pesos de frecuencia{wi}{\displaystyle \{w_{i}\}}se normalizan a 1, entonces la expresión correcta después de la corrección de Bessel se convierte en

s2 =i=1nortewii=1nortewi1i=1nortewi(incógnitaiμ)2{\displaystyle s^{2}\ ={\frac {\sum _{i=1}^{N}w_{i}}{\sum _{i=1}^{N}w_{i}-1}}\sum _{i=1}^{N}w_{i}\left(x_{i}-\mu ^{*}\right)^{2}}

donde el número total de muestras esi=1nortewi{\displaystyle \sum _{i=1}^{N}w_{i}}(nonorte{\displaystyle N}). En cualquier caso, la información sobre el número total de muestras es necesaria para obtener una corrección imparcial, incluso siwi{\displaystyle w_{i}}tiene un significado diferente al de peso de frecuencia.

El estimador solo puede ser insesgado si los pesos no están estandarizados ni normalizados , ya que estos procesos cambian la media y la varianza de los datos y, por lo tanto, conllevan una pérdida de la tasa base (el recuento de la población, que es un requisito para la corrección de Bessel).

Pesos de confiabilidad

Si en cambio los pesos son pesos de confiabilidad (valores no aleatorios que reflejan la confiabilidad relativa de la muestra, a menudo derivados de la varianza de la muestra), podemos determinar un factor de corrección para obtener un estimador insesgado. Suponiendo que cada variable aleatoria se muestrea de la misma distribución con mediaμ{\displaystyle \mu }y variación realσactual2{\displaystyle \sigma _{\text{actual}}^{2}}, teniendo en cuenta las expectativas que tenemos,

mi[σ^2]=i=1nortemi[(incógnitaiμ)2]norte=mi[(incógnitami[incógnita])2]1nortemi[(incógnitami[incógnita])2]=(norte1norte)σactual2mi[σ^w2]=i=1nortewimi[(incógnitaiμ)2]V1=mi[(incógnitami[incógnita])2]V2V12mi[(incógnitami[incógnita])2]=(1V2V12)σactual2{\displaystyle {\begin{aligned}\operatorname {E} [{\hat {\sigma }}^{2}]&={\frac {\sum \limits _{i=1}^{N}\operatorname {E} [(x_{i}-\mu )^{2}]}{N}}\\&=\operatorname {E} [(X-\operatorname {E} [X])^{2}]-{\frac {1}{N}}\operatorname {E} [(X-\operatorname {E} [X])^{2}]\\&=\left({\frac {N-1}{N}}\right)\sigma _{\text{actual}}^{2}\\\operatorname {E} [{\hat {\sigma }}_{\mathrm {w} }^{2}]&={\frac {\sum \limits _{i=1}^{N}w_{i}\operatorname {E} [(x_{i}-\mu ^{*})^{2}]}{V_{1}}}\\&=\operatorname {E} [(X-\operatorname {E} [X])^{2}]-{\frac {V_{2}}{V_{1}^{2}}}\operatorname {E} [(X-\operatorname {E} [X])^{2}]\\&=\left(1-{\frac {V_{2}}{V_{1}^{2}}}\right)\sigma _{\text{actual}}^{2}\end{aligned}}}

dóndeV1=i=1nortewi{\displaystyle V_{1}=\sum _{i=1}^{N}w_{i}}yV2=i=1nortewi2{\displaystyle V_{2}=\sum _{i=1}^{N}w_{i}^{2}}Por lo tanto, el sesgo en nuestro estimador es(1V2V12){\displaystyle \left(1-{\frac {V_{2}}{V_{1}^{2}}}\right)}, análogo al(norte1norte){\displaystyle \left({\frac {N-1}{N}}\right)}sesgo en el estimador no ponderado (nótese también que V12/V2=nortemiFF{\displaystyle \ V_{1}^{2}/V_{2}=N_{eff}}es el tamaño efectivo de la muestra ). Esto significa que para corregir el sesgo de nuestro estimador necesitamos dividir previamente por1(V2/V12){\displaystyle 1-\left(V_{2}/V_{1}^{2}\right)}, asegurando que el valor esperado de la varianza estimada sea igual a la varianza real de la distribución muestral . La estimación final insesgada de la varianza muestral es:

sw2 =σ^w21(V2/V12)=i=1nortewi(incógnitaiμ)2V1(V2/V1),{\displaystyle {\begin{aligned}s_{\mathrm {w} }^{2}\ &={\frac {{\hat {\sigma }}_{\mathrm {w} }^{2}}{1-(V_{2}/V_{1}^{2})}}\\[4pt]&={\frac {\sum \limits _{i=1}^{N}w_{i}(x_{i}-\mu ^{*})^{2}}{V_{1}-(V_{2}/V_{1})}},\end{aligned}}}[ 6 ]

dóndemi[sw2]=σactual2{\displaystyle \operatorname {E} [s_{\mathrm {w} }^{2}]=\sigma _{\text{actual}}^{2}}Los grados de libertad de esta varianza muestral ponderada e insesgada varían en consecuencia desde N  1 hasta  0. La desviación estándar es simplemente la raíz cuadrada de la varianza anterior.

Como nota al margen, se han descrito otros enfoques para calcular la varianza muestral ponderada. [ 7 ]

Covarianza de muestra ponderada

En una muestra ponderada, cada vector filaincógnitai{\displaystyle \mathbf {x} _{i}}(a cada conjunto de observaciones individuales en cada una de las K variables aleatorias) se le asigna un pesowi0{\displaystyle w_{i}\geq 0}.

Luego, el vector de media ponderadaμ{\displaystyle \mathbf {\mu ^{*}} }es dado por

μ=i=1nortewiincógnitaii=1nortewi.{\displaystyle \mathbf {\mu ^{*}} ={\frac {\sum _{i=1}^{N}w_{i}\mathbf {x} _{i}}{\sum _{i=1}^{N}w_{i}}}.}

Y la matriz de covarianza ponderada viene dada por: [ 8 ]

do=i=1nortewi(incógnitaiμ)T(incógnitaiμ)V1.{\displaystyle \mathbf {C} ={\frac {\sum _{i=1}^{N}w_{i}\left(\mathbf {x} _{i}-\mu ^{*}\right)^{T}\left(\mathbf {x} _{i}-\mu ^{*}\right)}{V_{1}}}.}

De forma similar a la varianza muestral ponderada, existen dos estimadores insesgados diferentes dependiendo del tipo de ponderaciones.

ponderaciones de frecuencia

Si los pesos son pesos de frecuencia , la estimación ponderada insesgada de la matriz de covarianzado{\displaystyle \textstyle \mathbf {C} }, con la corrección de Bessel, viene dada por: [ 8 ]

do=i=1nortewi(incógnitaiμ)T(incógnitaiμ)V11.{\displaystyle \mathbf {C} ={\frac {\sum _{i=1}^{N}w_{i}\left(\mathbf {x} _{i}-\mu ^{*}\right)^{T}\left(\mathbf {x} _{i}-\mu ^{*}\right)}{V_{1}-1}}.}

Este estimador solo puede ser insesgado si los pesos no están estandarizados ni normalizados , ya que estos procesos cambian la media y la varianza de los datos y, por lo tanto, conllevan una pérdida de la tasa base (el recuento de la población, que es un requisito para la corrección de Bessel).

Pesos de confiabilidad

En el caso de los pesos de confiabilidad , los pesos se normalizan :

V1=i=1nortewi=1.{\displaystyle V_{1}=\sum _{i=1}^{N}w_{i}=1.}

(Si no lo son, divida los pesos por su suma para normalizarlos antes de calcularlos)V1{\displaystyle V_{1}}:

wi=wii=1nortewi{\displaystyle w_{i}'={\frac {w_{i}}{\sum _{i=1}^{N}w_{i}}}}

Luego, el vector de media ponderadaμ{\displaystyle \mathbf {\mu ^{*}} }se puede simplificar a

μ=i=1nortewiincógnitai.{\displaystyle \mathbf {\mu ^{*}} =\sum _{i=1}^{N}w_{i}\mathbf {x} _{i}.}

y la estimación ponderada insesgada de la matriz de covarianzado{\displaystyle \mathbf {C} }es: [ 9 ]

do=i=1nortewi(i=1nortewi)2i=1nortewi2i=1nortewi(incógnitaiμ)T(incógnitaiμ)=i=1nortewi(incógnitaiμ)T(incógnitaiμ)V1(V2/V1).{\displaystyle {\begin{aligned}\mathbf {C} &={\frac {\sum _{i=1}^{N}w_{i}}{\left(\sum _{i=1}^{N}w_{i}\right)^{2}-\sum _{i=1}^{N}w_{i}^{2}}}\sum _{i=1}^{N}w_{i}\left(\mathbf {x} _{i}-\mu ^{*}\right)^{T}\left(\mathbf {x} _{i}-\mu ^{*}\right)\\&={\frac {\sum _{i=1}^{N}w_{i}\left(\mathbf {x} _{i}-\mu ^{*}\right)^{T}\left(\mathbf {x} _{i}-\mu ^{*}\right)}{V_{1}-(V_{2}/V_{1})}}.\end{aligned}}}

El razonamiento aquí es el mismo que en la sección anterior.

Dado que asumimos que los pesos están normalizados, entoncesV1=1{\displaystyle V_{1}=1}y esto se reduce a:

do=i=1nortewi(incógnitaiμ)T(incógnitaiμ)1V2.{\displaystyle \mathbf {C} ={\frac {\sum _{i=1}^{N}w_{i}\left(\mathbf {x} _{i}-\mu ^{*}\right)^{T}\left(\mathbf {x} _{i}-\mu ^{*}\right)}{1-V_{2}}}.}

Si todos los pesos son iguales, es decirwi/V1=1/norte{\displaystyle w_{i}/V_{1}=1/N}, entonces la media ponderada y la covarianza se reducen a la media muestral no ponderada y la covarianza mencionadas anteriormente.

Estimaciones con valores vectoriales

Lo anterior se generaliza fácilmente al caso de tomar la media de estimaciones con valores vectoriales. Por ejemplo, las estimaciones de posición en un plano pueden tener menos certeza en una dirección que en otra. Al igual que en el caso escalar, la media ponderada de múltiples estimaciones puede proporcionar una estimación de máxima verosimilitud . Simplemente reemplazamos la varianza.σ2{\displaystyle \sigma ^{2}}por la matriz de covarianzado{\displaystyle \mathbf {C} }y la inversa aritmética por la inversa de la matriz (ambas denotadas de la misma manera, mediante superíndices); la matriz de pesos queda entonces: [ 10 ]

Wi=doi1.{\displaystyle \mathbf {W} _{i}=\mathbf {C} _{i}^{-1}.}

La media ponderada en este caso es: incógnita¯=doincógnita¯(i=1norteWiincógnitai),{\displaystyle {\bar {\mathbf {x} }}=\mathbf {C} _{\bar {\mathbf {x} }}\left(\sum _{i=1}^{n}\mathbf {W} _{i}\mathbf {x} _{i}\right),} (donde el orden del producto matriz-vector no es conmutativo ), en términos de la covarianza de la media ponderada: doincógnita¯=(i=1norteWi)1,{\displaystyle \mathbf {C} _{\bar {\mathbf {x} }}=\left(\sum _{i=1}^{n}\mathbf {W} _{i}\right)^{-1},}

Por ejemplo, consideremos la media ponderada del punto [1 0] con alta varianza en el segundo componente y [0 1] con alta varianza en el primer componente. Entonces

incógnita1:=[10],do1:=[100100]{\displaystyle \mathbf {x} _{1}:={\begin{bmatrix}1&0\end{bmatrix}}^{\top },\qquad \mathbf {C} _{1}:={\begin{bmatrix}1&0\\0&100\end{bmatrix}}}
incógnita2:=[01],do2:=[100001]{\displaystyle \mathbf {x} _{2}:={\begin{bmatrix}0&1\end{bmatrix}}^{\top },\qquad \mathbf {C} _{2}:={\begin{bmatrix}100&0\\0&1\end{bmatrix}}}

Entonces, la media ponderada es:

incógnita¯=(do11+do21)1(do11incógnita1+do21incógnita2)=[0,9901000,9901][11]=[0,99010,9901]{\displaystyle {\begin{aligned}{\bar {\mathbf {x} }}&=\left(\mathbf {C} _{1}^{-1}+\mathbf {C} _{2}^{-1}\right)^{-1}\left(\mathbf {C} _{1}^{-1}\mathbf {x} _{1}+\mathbf {C} _{2}^{-1}\mathbf {x} _{2}\right)\\[5pt]&={\begin{bmatrix}0.9901&0\\0&0.9901\end{bmatrix}}{\begin{bmatrix}1\\1\end{bmatrix}}={\begin{bmatrix}0.9901\\0.9901\end{bmatrix}}\end{aligned}}}

lo cual tiene sentido: la estimación [1 0] es "conforme" en el segundo componente y la estimación [0 1] es conforme en el primer componente, por lo que la media ponderada es casi [1 1].

Consideración de las correlaciones

En el caso general, supongamos queincógnita=[incógnita1,,incógnitanorte]T{\displaystyle \mathbf {X} =[x_{1},\dots ,x_{n}]^{T}},do{\displaystyle \mathbf {C} }es la matriz de covarianza que relaciona las cantidadesincógnitai{\displaystyle x_{i}},incógnita¯{\displaystyle {\bar {x}}}es la media común que se debe estimar, yJ{\displaystyle \mathbf {J} }es una matriz de diseño igual a un vector de unos[1,,1]T{\displaystyle [1,\dots ,1]^{T}}(de longitudnorte{\displaystyle n}El teorema de Gauss-Markov establece que la estimación de la media con varianza mínima viene dada por:

σincógnita¯2=(JTWJ)1,{\displaystyle \sigma _{\bar {x}}^{2}=(\mathbf {J} ^{T}\mathbf {W} \mathbf {J} )^{-1},}

y

incógnita¯=σincógnita¯2(JTWincógnita),{\displaystyle {\bar {x}}=\sigma _{\bar {x}}^{2}(\mathbf {J} ^{T}\mathbf {W} \mathbf {X} ),}

dónde:

W=do1.{\displaystyle \mathbf {W} =\mathbf {C} ^{-1}.}

Disminución de la fuerza de las interacciones

Consideremos la serie temporal de una variable independiente.incógnita{\displaystyle x}y una variable dependientey{\displaystyle y}, connorte{\displaystyle n}observaciones muestreadas en momentos discretosti{\displaystyle t_{i}}. En muchas situaciones comunes, el valor dey{\displaystyle y}en ese momentoti{\displaystyle t_{i}}depende no solo deincógnitai{\displaystyle x_{i}}pero también en sus valores pasados. Comúnmente, la fuerza de esta dependencia disminuye a medida que aumenta la separación de las observaciones en el tiempo. Para modelar esta situación, se puede reemplazar la variable independiente por su media móvil.z{\displaystyle z}para un tamaño de ventanametro{\displaystyle m}.

zk=i=1metrowiincógnitak+1i.{\displaystyle z_{k}=\sum _{i=1}^{m}w_{i}x_{k+1-i}.}

Pesos que disminuyen exponencialmente

En el escenario descrito en la sección anterior, con mayor frecuencia la disminución de la fuerza de interacción obedece a una ley exponencial negativa. Si las observaciones se muestrean en tiempos equidistantes, entonces la disminución exponencial es equivalente a una disminución en una fracción constante.0<Δ<1{\displaystyle 0<\Delta <1}en cada paso de tiempo. Configuraciónw=1Δ{\displaystyle w=1-\Delta }podemos definirmetro{\displaystyle m}pesos normalizados por

wi=wi1V1,{\displaystyle w_{i}={\frac {w^{i-1}}{V_{1}}},}

dóndeV1{\displaystyle V_{1}}es la suma de los pesos no normalizados. En este casoV1{\displaystyle V_{1}}es simplemente

V1=i=1metrowi1=1wmetro1w,{\displaystyle V_{1}=\sum _{i=1}^{m}{w^{i-1}}={\frac {1-w^{m}}{1-w}},}

que se acercaV1=1/(1w){\displaystyle V_{1}=1/(1-w)}para valores grandes demetro{\displaystyle m}.

La constante de amortiguaciónw{\displaystyle w}debe corresponder a la disminución real de la fuerza de interacción. Si esto no se puede determinar a partir de consideraciones teóricas, entonces las siguientes propiedades de pesos que disminuyen exponencialmente son útiles para hacer una elección adecuada: en el paso(1w)1{\displaystyle (1-w)^{-1}}, el peso es aproximadamente igual ami1(1w)=0,39(1w){\displaystyle {e^{-1}}(1-w)=0.39(1-w)}, el área de la cola el valormi1{\displaystyle e^{-1}}la zona de la cabeza1mi1=0,61{\displaystyle {1-e^{-1}}=0.61}. El área de la cola en el pasonorte{\displaystyle n}esminorte(1w){\displaystyle \leq {e^{-n(1-w)}}}. Donde principalmente el más cercanonorte{\displaystyle n}Las observaciones importan y el efecto de las observaciones restantes puede ignorarse con seguridad, entonces elijaw{\displaystyle w}de manera que la zona de la cola sea suficientemente pequeña.

Promedios ponderados de funciones

El concepto de promedio ponderado puede extenderse a funciones. [ 11 ] Los promedios ponderados de funciones desempeñan un papel importante en los sistemas de cálculo diferencial e integral ponderado. [ 12 ]

Corrección de la sobredispersión o subdispersión

Las medias ponderadas se utilizan normalmente para encontrar la media ponderada de datos históricos, en lugar de datos generados teóricamente. En este caso, habrá algún error en la varianza de cada punto de datos. Normalmente, los errores experimentales pueden subestimarse debido a que el experimentador no tiene en cuenta todas las fuentes de error al calcular la varianza de cada punto de datos. En este caso, la varianza en la media ponderada debe corregirse para tener en cuenta el hecho de queχ2{\displaystyle \chi ^{2}}es demasiado grande. La corrección que debe hacerse es

σ^incógnita¯2=σincógnita¯2χν2{\displaystyle {\hat {\sigma }}_{\bar {x}}^{2}=\sigma _{\bar {x}}^{2}\chi _{\nu }^{2}}

dóndeχν2{\displaystyle \chi _{\nu }^{2}}es el chi-cuadrado reducido :

χν2=1(norte1)i=1norte(incógnitaiincógnita¯)2σi2;{\displaystyle \chi _{\nu }^{2}={\frac {1}{(n-1)}}\sum _{i=1}^{n}{\frac {(x_{i}-{\bar {x}})^{2}}{\sigma _{i}^{2}}};}

La raíz cuadradaσ^incógnita¯{\displaystyle {\hat {\sigma }}_{\bar {x}}} can be called the standard error of the weighted mean (variance weights, scale corrected).

When all data variances are equal, σi=σ0{\displaystyle \sigma _{i}=\sigma _{0}}, they cancel out in the weighted mean variance, σx¯2{\displaystyle \sigma _{\bar {x}}^{2}}, which again reduces to the standard error of the mean (squared), σx¯2=σ2/n{\displaystyle \sigma _{\bar {x}}^{2}=\sigma ^{2}/n}, formulated in terms of the sample standard deviation (squared),

σ2=i=1n(xix¯)2n1.{\displaystyle \sigma ^{2}={\frac {\sum _{i=1}^{n}(x_{i}-{\bar {x}})^{2}}{n-1}}.}

See also

Notes

  1. Technically, negatives may be used if all the values are either zero or negatives. This fills no function however as the weights work as absolute values.

References

  1. 1234Cochran, W. G. (1977). Sampling Techniques (3rd ed.). Nashville, TN: John Wiley & Sons. ISBN 978-0-471-16240-7
  2. 1234567891011121314151617Carl-Erik Sarndal; Bengt Swensson; Jan Wretman (1992). Model Assisted Survey Sampling. Springer. ISBN 978-0-387-97528-3.
  3. Thomas Lumley (https://stats.stackexchange.com/users/249135/thomas-lumley), How to estimate the (approximate) variance of the weighted mean?, URL (version: 2021-06-08): https://stats.stackexchange.com/q/525770
  4. Gatz, Donald F.; Smith, Luther (June 1995). "The standard error of a weighted mean concentration—I. Bootstrapping vs other methods". Atmospheric Environment. 29 (11): 1185–1193. Bibcode:1995AtmEn..29.1185G. doi:10.1016/1352-2310(94)00210-C. - pdf link
  5. Endlich, R. M.; Eymon, B. P.; Ferek, R. J.; Valdes, A. D.; Maxwell, C. (1988-12-01). "Statistical Analysis of Precipitation Chemistry Measurements over the Eastern United States. Part I: Seasonal and Regional Patterns and Correlations". Journal of Applied Meteorology and Climatology. 27 (12): 1322–1333. Bibcode:1988JApMe..27.1322E. doi:10.1175/1520-0450(1988)027<1322:SAOPCM>2.0.CO;2.
  6. "Biblioteca Científica GNU – Manual de Referencia: Muestras Ponderadas" . Gnu.org . Consultado el 22 de diciembre de 2017 .
  7. "Error estándar ponderado y su impacto en las pruebas de significancia (WinCross vs. Quantum y SPSS), Dr. Albert Madansky" (PDF) . Analyticalgroup.com . Consultado el 22 de diciembre de 2017 .
  8. 1 2 Price, George R. (abril de 1972). "Extensión de las matemáticas de selección de covarianza" ( PDF) . Annals of Human Genetics . 35 (4): 485– 490. doi : 10.1111/j.1469-1809.1957.tb01874.x . PMID 5073694. S2CID 37828617 .  
  9. Mark Galassi, Jim Davies, James Theiler, Brian Gough, Gerard Jungman, Michael Booth y Fabrice Rossi. Biblioteca Científica GNU - Manual de referencia, Versión 1.15 , 2011. Sec. 21.7 Muestras ponderadas
  10. James, Frederick (2006). Métodos estadísticos en física experimental (2.ª ed.). Singapur: World Scientific. pág. 324. ISBN   981-270-527-9.
  11. GH Hardy , JE Littlewood y G. Pólya . Desigualdades (2.ª ed.), Cambridge University Press, ISBN 978-0-521-35880-4, 1988.
  12. Jane Grossman, Michael Grossman, Robert Katz. Los primeros sistemas de cálculo diferencial e integral ponderado , ISBN 0-9771170-1-4, 1980.

Lecturas adicionales

  • Bevington, Philip R (1969). Reducción de datos y análisis de errores para las ciencias físicas . Nueva York, NY: McGraw-Hill. OCLC 300283069 . 
  • Strutz, T. (2010). Ajuste de datos e incertidumbre (Una introducción práctica a los mínimos cuadrados ponderados y más allá) . Vieweg+Teubner. ISBN 978-3-8348-1022-9.
  • David Terr. "Media ponderada" . MathWorld .
  • Herramienta para calcular el promedio ponderado