Cantidad estadística
La media aritmética ponderada es similar a una media aritmética ordinaria (el tipo más común de promedio ), excepto que en lugar de que cada uno de los puntos de datos contribuya por igual al promedio final, algunos puntos de datos contribuyen más que otros. La noción de media ponderada desempeña un papel en la 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 la misma que 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
Se le dan dos clases escolares , una con 20 estudiantes y otra con 30 estudiantes , y las calificaciones de las pruebas en cada clase son las siguientes:
- Clase de la mañana = {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 media de la clase de la tarde es 90. La media no ponderada de las dos medias es 85. Sin embargo, esto no tiene en cuenta la diferencia en el número de estudiantes en cada clase (20 versus 30); por lo tanto, el valor de 85 no refleja la calificación promedio del estudiante (independientemente de la clase). La calificación promedio de un estudiante se puede obtener promediando todas las calificaciones, sin importar las clases (sume todas las calificaciones y divida por el número total de estudiantes):![{\displaystyle {\bar {x}}={\frac {4300}{50}}=86.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
O esto se puede lograr ponderando las medias de las clases por el número de estudiantes en cada clase. A la clase más grande se le da más "peso":
![{\displaystyle {\bar {x}}={\frac {(20\times 80)+(30\times 90)}{20+30}}=86.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Así, la media ponderada permite encontrar la nota media de los alumnos sin conocer la puntuación de cada alumno. Sólo se necesitan las medias de clase y el número de estudiantes en cada clase.
Ejemplo de combinación convexa
Dado que sólo los pesos relativos son relevantes, cualquier media ponderada puede expresarse utilizando coeficientes que sumen uno. Esta combinación lineal se llama combinación convexa .
Usando el ejemplo anterior, obtendríamos los siguientes pesos:
![{\displaystyle {\frac {20}{20+30}}=0,4}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\frac {30}{20+30}}=0,6}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Luego, aplica los pesos de esta manera:
![{\displaystyle {\bar {x}}=(0,4\times 80)+(0,6\times 90)=86.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Definición matemática
Formalmente, la media ponderada de una tupla finita de datos no vacía , con sus correspondientes pesos no negativos es
![{\displaystyle \left(w_{1},w_{2},\dots,w_{n}\right)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\bar {x}}={\frac {\sum \limits _{i=1}^{n}w_{i}x_{i}}{\sum \limits _{i=1}^ {n}w_ {i}}},}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
que se expande a:
![{\displaystyle {\bar {x}}={\frac {w_{1}x_{1}+w_{2}x_{2}+\cdots +w_{n}x_{n}}{w_{1} +w_{2}+\cdots +w_{n}}}.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
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 modo que sumen 1, es decir, . Para tales ponderaciones normalizadas, la media ponderada es equivalente:![{\textstyle \sum \limits _ {i=1}^{n}{w_{i}'}=1}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
.
Siempre se pueden normalizar los pesos realizando la siguiente transformación sobre los pesos originales:
.
La media ordinaria es un caso especial de la media ponderada donde todos los datos tienen pesos iguales.![{\textstyle {\frac {1}{n}}\sum \limits _{i=1}^{n}{x_{i}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Si los elementos de datos son variables aleatorias independientes y distribuidas idénticamente con varianza , se puede demostrar que el error estándar de la media ponderada , mediante la propagación de la incertidumbre , es:![{\displaystyle \sigma ^{2}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \sigma _{\bar {x}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\textstyle \sigma _{\bar {x}}=\sigma {\sqrt {\sum \limits _ {i=1}^{n}w_{i}'^{2}}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Pesos definidos por la varianza
Para la media ponderada de una lista de datos en la que cada elemento proviene potencialmente de una distribución de probabilidad diferente con varianza conocida , todos con la misma media, una posible elección para las ponderaciones viene dada por el recíproco de la varianza:
![{\displaystyle \sigma _ {i}^{2}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle w_{i}={\frac {1}{\sigma _{i}^{2}}}.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
La media ponderada en este caso es:
![{\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}}},}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
y el error estándar de la media ponderada (con ponderaciones de varianza inversa) es:
![{\displaystyle \sigma _{\bar {x}}={\sqrt {\frac {1}{\sum _{i=1}^{n}\sigma _{i}^{-2}}}} ={\sqrt {\frac {1}{\sum _{i=1}^{n}w_{i}}}},}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Tenga en cuenta que esto se reduce a cuando todos . Es un caso especial de la fórmula general de la sección anterior,![{\displaystyle \sigma _{\bar {x}}^{2}=\sigma _{0}^{2}/n}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \sigma _{i}=\sigma _{0}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\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}}}.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Las ecuaciones anteriores se pueden combinar para obtener:
![{\displaystyle {\bar {x}}=\sigma _{\bar {x}}^{2}\sum _{i=1}^{n}{\frac {x_{i}}{\sigma _ {i}^{2}}}.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
La importancia de esta elección es 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 están normalmente distribuidas con la misma media.
Propiedades estadísticas
Expectativa
La media muestral ponderada, , 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, de la siguiente manera. Para simplificar, asumimos ponderaciones normalizadas (ponderaciones que suman uno).![{\displaystyle {\bar {x}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Si las observaciones tienen valores esperados,
entonces la media muestral ponderada tiene expectativa.
En particular, si las medias son iguales, entonces la expectativa de la media muestral ponderada será ese valor,![{\displaystyle E(x_{i})={\mu _{i}},}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle E({\bar {x}})=\sum _{i=1}^{n}{w_{i}'\mu _{i}}.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \mu _{i}=\mu }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle E({\bar {x}})=\mu.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Diferencia
caso iid simple
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 expectativa (como es el caso de las variables aleatorias iid ), 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 prueba ):
![{\displaystyle \operatorname {Var} ({\bar {y}}_{w})={\hat {\sigma }}_{y}^{2}{\frac {\overline {w^{2} }}{{\bar {w}}^{2}}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Con , , y![{\displaystyle {\hat {\sigma }}_{y}^{2}={\frac {\sum _{i=1}^{n}(y_{i}-{\bar {y}}) ^{2}}{n-1}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\bar {w}}={\frac {\sum _{i=1}^{n}w_{i}}{n}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\overline {w^{2}}}={\frac {\sum _{i=1}^{n}w_{i}^{2}}{n}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Sin embargo, esta estimación es bastante limitada debido a los fuertes supuestos sobre las observaciones y . Esto ha llevado al desarrollo de estimadores alternativos, más generales.
Perspectiva del muestreo de la encuesta
Desde una perspectiva basada en modelos , nos interesa estimar la varianza de la media ponderada cuando las diferentes no son variables aleatorias iid . Una perspectiva alternativa para este problema es la de algún diseño de muestreo arbitrario de los datos en el que las unidades se seleccionan con probabilidades desiguales (con reemplazo). [1] : 306 ![{\ Displaystyle y_ {i}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
En la metodología de encuesta , 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 ( ). o estimado ( ). 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 del indicador Bernoulli ( ) 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 muestreo de tamaño de muestra variado (por ejemplo, muestreo de Poisson ). La probabilidad de que se elija algún elemento, dada una muestra, se denota como , y la probabilidad de selección de un sorteo es (si N es muy grande y cada uno es muy pequeño). Para la siguiente derivación asumiremos que la probabilidad de seleccionar cada elemento está completamente 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 ).![{\displaystyle N}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\sombrero {N}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\ Displaystyle I_ {i}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle P(I_{i}=1\mid {\text{Alguna muestra de tamaño }}n)=\pi _{i}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle P(I_{i}=1|{\text{un sorteo de muestra}})=p_{i}\approx {\frac {\pi _{i}}{n}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\ Displaystyle p_ {i}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Como cada elemento ( ) es fijo, y la aleatoriedad proviene de estar incluido o no en la muestra ( ), solemos hablar de la multiplicación de los dos, que es una variable aleatoria. Para evitar confusiones en la siguiente sección, llamemos a este término: . Con la siguiente expectativa: ; y variación: .![{\ Displaystyle y_ {i}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\ Displaystyle I_ {i}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\ Displaystyle y'_ {i} = y_ {i} I_ {i}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle E[y'_{i}]=y_{i}E[I_{i}]=y_{i}\pi _{i}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle V[y'_{i}]=y_{i}^{2}V[I_{i}]=y_{i}^{2}\pi _{i}(1-\pi _{ i})}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Cuando cada elemento de la muestra se infla por la inversa de su probabilidad de selección, se denomina valores de y expandidos , es decir: . Una cantidad relacionada tiene valores de y expandidos : . [2] : 42, 43, 51, 52 Como arriba, podemos agregar una marca si multiplicamos por la función del indicador. Es decir:![{\displaystyle \pi }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\check {y}}_{i}={\frac {y_{i}}{\pi _{i}}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle p}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\frac {y_{i}}{p_{i}}}=n{\check {y}}_{i}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\check {y}}'_{i}=I_{i}{\check {y}}_{i}={\frac {I_{i}y_{i}}{\pi _ i}}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
En esta perspectiva basada en el diseño , las ponderaciones, utilizadas en el numerador de la media ponderada, se obtienen tomando la inversa de la probabilidad de selección (es decir, el factor de inflación). Es decir: .![{\displaystyle w_{i}={\frac {1}{\pi _{i}}}\aprox {\frac {1}{n\times p_{i}}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Varianza de la suma ponderada (poder-estimador de totales)
Si se conoce el tamaño de la población N, podemos estimar la media poblacional usando .![{\displaystyle {\hat {\bar {Y}}}_{{\text{conocido }}N}={\frac {{\hat {Y}}_{pwr}}{N}}\approx {\ frac {\sum _{i=1}^{n}w_{i}y'_{i}}{N}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Si el diseño de muestreo da como resultado un tamaño de muestra fijo n (como en el muestreo pps ), entonces la varianza de este estimador es:
![{\displaystyle \operatorname {Var} \left({\hat {\bar {Y}}}_{{\text{conocido }}N}\right)={\frac {1}{N^{2}} }{\frac {n}{n-1}}\sum _{i=1}^{n}\left(w_{i}y_{i}-{\overline {wy}}\right)^{2 }}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
PruebaLa fórmula general se puede desarrollar así:
![{\displaystyle {\hat {\bar {Y}}}_{{\text{conocido }}N}={\frac {{\hat {Y}}_{pwr}}{N}}={\frac {{\frac {1}{n}}\sum _{i=1}^{n}{\frac {y'_{i}}{p_{i}}}}{N}}\aprox {\ frac {\sum _{i=1}^{n}{\frac {y'_{i}}{\pi _{i}}}}{N}}={\frac {\sum _{i= 1}^{n}w_{i}y'_{i}}{N}}.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
La población total se denota como y puede estimarse mediante el estimador (insesgado) de Horvitz-Thompson , también llamado estimador. Este estimador puede estimarse a su vez utilizando el estimador pwr (es decir: estimador ampliado con reemplazo o estimador de "probabilidad con reemplazo"). Con la notación anterior, es: . [2] : 51 ![{\displaystyle Y=\sum _ {i=1}^{N}y_ {i}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \pi }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle p}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\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}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
La varianza estimada del estimador pwr viene dada por: [2] : 52
donde .![{\displaystyle \operatorname {Var} ({\hat {Y}}_{pwr})={\frac {n}{n-1}}\sum _{i=1}^{n}\left(w_ {i}y_{i}-{\overline {wy}}\right)^{2}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\overline {wy}}=\sum _{i=1}^{n}{\frac {w_{i}y_{i}}{n}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
La fórmula anterior fue tomada de Sarndal et al. (1992) (también presentado en Cochran 1977), pero fue escrito de manera diferente. [2] : 52 [1] : 307 (11.35) El lado izquierdo es cómo se escribió la varianza y el lado derecho es cómo desarrollamos la versión ponderada:
![{\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}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Y llegamos a la fórmula desde 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
![{\displaystyle \operatorname {Var} ({\hat {\bar {Y}}}_{{\text{pwr (conocido }}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)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Con . Además, ¿dónde está la probabilidad de seleccionar tanto i como j? [2] : 36 Y , y para i=j: . [2] : 43 ![{\displaystyle {\check {y}}_{i}={\frac {y_{i}}{\pi _{i}}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle C(I_{i},I_{j})=\pi _{ij}-\pi _{i}\pi _{j}=\Delta _{ij}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \pi _{ij}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\check {\Delta }}_{ij}=1-{\frac {\pi _{i}\pi _{j}}{\pi _{ij}}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\check {\Delta }}_{ii}=1-{\frac {\pi _{i}\pi _{i}}{\pi _{i}}}=1-\pi _ {i}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Si las probabilidades de selección no están correlacionadas (es decir: ), y al asumir que la probabilidad de cada elemento es muy pequeña, entonces:![{\displaystyle \forall i\neq j:C(I_{i},I_{j})=0}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \operatorname {Var} ({\hat {\bar {Y}}}_{{\text{pwr (conocido }}N{\text{)}}})={\frac {1}{N ^{2}}}\sum _{i=1}^{n}\left(w_{i}y_{i}\right)^{2}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
PruebaAsumimos eso y eso![{\displaystyle (1-\pi _ {i})\aproximadamente 1}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\begin{aligned}\operatorname {Var} ({\hat {Y}}_{{\text{pwr (conocido }}N{\text{)}}})&={\frac {1 }{N^{2}}}\sum _{i=1}^{n}\sum _{j=1}^{n}\left({\check {\Delta }}_{ij}{\ comprobar {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}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Varianza de la media ponderada (π-estimador de relación-media)
La sección anterior trató sobre la estimación de la media poblacional como una proporción de una población total estimada ( ) con un tamaño de población conocido ( ), y la varianza se estimó en ese contexto. Otro caso común es que el tamaño de la población en sí ( ) se desconoce y se estima utilizando la muestra (es decir: ). La estimación de puede describirse como la suma de pesos. Así que cuando lleguemos . Con la notación anterior, el parámetro que nos interesa es la relación de las sumas de sy 1s. Es decir: . Podemos estimarlo usando nuestra muestra con: . A medida que pasamos del uso de N al uso de n, en realidad sabemos que todas las variables del indicador obtienen 1, por lo que simplemente podríamos escribir: . Esta será la estimación para valores específicos de y y w, pero las propiedades estadísticas surgen al incluir la variable indicadora . [2] : 162, 163, 176 ![{\displaystyle {\sombrero {Y}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle N}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle N}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\sombrero {N}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle N}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle w_{i}={\frac {1}{\pi _{i}}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\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}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\ Displaystyle y_ {i}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle R={\bar {Y}}={\frac {\sum _{i=1}^{N}{\frac {y_{i}}{\pi _{i}}}}{\ suma _{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}}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\sombrero {R}}={\sombrero {\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}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\bar {y}}_{w}={\frac {\sum _{i=1}^{n}w_{i}y_{i}}{\sum _{i=1}^ {n}w_{i}}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\bar {y}}_{w}={\frac {\sum _{i=1}^{n}w_{i}y'_{i}}{\sum _{i=1 }^{n}w_{i}1'_{i}}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Esto se llama estimador de razón y es aproximadamente insesgado para R. [2] : 182
En este caso, la variabilidad de la relació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 una estimación aproximada. Principalmente linealización de primer orden de series de Taylor , asintóticas y bootstrap/jackknife. [2] : 172 El método de linealización de Taylor podría llevar a una subestimación de la varianza para tamaños de muestra pequeños en general, pero eso depende de la complejidad de la estadística. Para la media ponderada, se supone que la varianza aproximada es relativamente precisa incluso para tamaños de muestra medianos. [2] : 176 Para cuando el muestreo tiene un tamaño de muestra aleatorio (como en el muestreo de Poisson ), queda como sigue: [2] : 182
.
Si , entonces usar o daría el mismo estimador, ya que multiplicar por algún factor conduciría al mismo estimador. También significa que si escalamos la suma de ponderaciones para que sea igual a un tamaño de población N conocido desde antes , 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.![{\displaystyle \pi _{i}\aproximadamente p_{i}n}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle w_{i}={\frac {1}{\pi _{i}}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle w_{i}={\frac {1}{p_{i}}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle w_{i}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
PruebaLa linealización de Taylor establece que para un estimador de razón general de dos sumas ( ), se pueden expandir alrededor del valor verdadero R y dar: [2] : 178 ![{\displaystyle {\sombrero {R}}={\frac {\sombrero {Y}}{\sombrero {Z}}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\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}}}\aprox R+{\frac {1}{Z}}\sum _{i =1}^{n}\left({\frac {y'_{i}}{\pi _{i}}}-R{\frac {z'_{i}}{\pi _{i} }}\bien)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Y la varianza se puede aproximar mediante: [2] : 178, 179
.
El término 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: ), este término aún incluiría una suma de n covarianzas para cada elemento i entre y . Esto ayuda a ilustrar que esta fórmula incorpora el efecto de la correlación entre y y z en la varianza de los estimadores de razones.![{\displaystyle {\sombrero {C}}({\sombrero {Y}},{\sombrero {Z}})}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \forall i\neq j:\Delta _{ij}=C(I_{i},I_{j})=0}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle y'_{i}=I_{i}y_{i}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle z'_{i}=I_{i}z_{i}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Al definir lo anterior queda: [2] : 182 ![{\displaystyle z_{i}=1}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\widehat {V({\hat {R}})}}={\widehat {V({\bar {y}}_{w})}}={\frac {1}{{\ sombrero {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}}}\derecha).}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Si las probabilidades de selección no están correlacionadas (es decir: ), y cuando se supone que la probabilidad de cada elemento es muy pequeña (es decir: ), entonces lo anterior se reduce a lo siguiente:![{\displaystyle \forall i\neq j:\Delta _{ij}=C(I_{i},I_{j})=0}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle (1-\pi _ {i})\aproximadamente 1}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\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}}}\ derecha)^{2}={\frac {1}{(\sum _{i=1}^{n}w_{i})^{2}}}\sum _{i=1}^{n} w_{i}^{2}(y_{i}-{\bar {y}}_{w})^{2}.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Thomas Lumley proporcionó una recreación similar de la prueba (con algunos errores al final) en validación cruzada. [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 desconocida. No existe un enfoque uniformemente mejor, pero la literatura presenta varios argumentos para preferir usar la versión de estimación de 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 poblacional desconocida se considera más estable. Por último, si la proporción del muestreo está correlacionada negativamente con los valores (es decir, menor probabilidad de muestrear una observación que es grande), entonces la versión de tamaño de población desconocido lo compensa ligeramente.
Para el caso trivial en el que todos los pesos son iguales a 1, la fórmula anterior es igual a la fórmula normal 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: dividiéndolo por n en lugar de (n-1)).
Validación de arranque
Se ha demostrado, por Gatz et al. (1995), que en comparación con los métodos de arranque , lo siguiente (estimación de la varianza de la relación media utilizando la linealización de la serie de Taylor ) es una estimación razonable para el cuadrado del error estándar de la media (cuando se utiliza en el contexto de la medición de componentes químicos) : [4] : 1186
![{\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{\ barra {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]}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
dónde . Una mayor simplificación conduce a![{\displaystyle {\bar {w}}={\frac {\sum w_{i}}{n}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\widehat {\sigma _{\bar {x}}^{2}}}={\frac {n}{(n-1)(n{\bar {w}})^{2} }}\sum w_{i}^{2}(x_{i}-{\bar {x}}_{w})^{2}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Gatz et al. mencionar que la formulación anterior fue publicada por Endlich et al. (1988) al tratar la media ponderada como una combinación de un estimador total ponderado dividido por un estimador del tamaño de la población, [5] basándose en la formulación publicada por Cochran (1977), como una aproximación a la media del ratio. Sin embargo, Endlich et al. no pareció 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 Aún así, es casi idéntico 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, en la literatura se propuso confiar en métodos de replicación como Jackknife y Bootstrapping . [1] : 321
Otras notas
Para observaciones no correlacionadas con varianzas , la varianza de la media muestral ponderada es [ cita necesaria ]![{\displaystyle \sigma _ {i}^{2}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \sigma _{\bar {x}}^{2}=\sum _{i=1}^{n}{w_{i}'^{2}\sigma _{i}^{2} }}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
cuya raíz cuadrada puede denominarse error estándar de la media ponderada (caso general) . [ cita necesaria ]![{\displaystyle \sigma _{\bar {x}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
En consecuencia, si todas las observaciones tienen la misma varianza, la media muestral ponderada tendrá varianza![{\displaystyle \sigma _ {i}^{2}=\sigma _ {0}^{2}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \sigma _{\bar {x}}^{2}=\sigma _{0}^{2}\sum _{i=1}^{n}{w_{i}'^{2} },}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
dónde . La varianza alcanza su valor máximo, 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 , es decir, degenera en el error estándar de la media , al cuadrado.![{\textstyle 1/n\leq \sum _ {i=1}^{n}{w_{i}'^{2}}\leq 1}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \sigma _{0}^{2}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\textstyle \sigma _{\bar {x}}=\sigma _{0}/{\sqrt {n}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Debido a 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 todo .![{\displaystyle w_{i}'={\frac {w_{i}}{\sum _{i=1}^{n}{w_{i}}}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Conceptos relacionados
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 utiliza una media ponderada , la varianza de la muestra ponderada es diferente de la varianza de la muestra no ponderada.![{\displaystyle \mu ^{*}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
La varianza de la muestra ponderada sesgada se define de manera similar a la varianza de la muestra sesgada normal : ![{\displaystyle {\hat {\sigma }}_{\mathrm {w} }^{2}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\sombrero {\sigma }}^{2}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\begin{alineado}{\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{alineado}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
donde para pesos normalizados. Si las ponderaciones son ponderaciones de frecuencia (y, por lo tanto, son variables aleatorias), se puede demostrar [ cita necesaria ] que es el estimador de máxima verosimilitud de las observaciones gaussianas iid .![{\displaystyle \sum _{i=1}^{N}w_{i}=1}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\hat {\sigma }}_{\mathrm {w} }^{2}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \sigma ^{2}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Para muestras pequeñas, se acostumbra utilizar un estimador insesgado de la varianza poblacional. En muestras normales no ponderadas, la N en el denominador (correspondiente al tamaño de la muestra) se cambia a N − 1 (consulte la corrección de Bessel ). En el entorno ponderado, en realidad existen dos estimadores insesgados diferentes, uno para el caso de las ponderaciones de frecuencia y otro para el caso de las ponderaciones de confiabilidad .
Pesos de frecuencia
Si las ponderaciones son ponderaciones de frecuencia (donde una ponderación es igual al número de ocurrencias), entonces el estimador insesgado es:
![{\displaystyle s^{2}\ ={\frac {\sum \limits _{i=1}^{N}w_{i}\left(x_{i}-\mu ^{*}\right)^ {2}}{\sum _{i=1}^{N}w_{i}-1}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Esto aplica efectivamente la corrección de Bessel para ponderaciones de frecuencia.
Por ejemplo, si los valores se extraen de la misma distribución, entonces podemos tratar este conjunto como una muestra no ponderada, o podemos tratarlo como una muestra ponderada con los pesos correspondientes , y obtenemos el mismo resultado de cualquier manera.![{\displaystyle \{2,2,4,5,5,5\}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \{2,4,5\}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \{2,1,3\}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Si las ponderaciones de frecuencia se normalizan a 1, entonces la expresión correcta después de la corrección de Bessel se convierte en![{\displaystyle \{w_{i}\}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\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}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
donde el número total de muestras es (no ). En cualquier caso, la información sobre el número total de muestras es necesaria para obtener una corrección insesgada, incluso si tiene un significado diferente al de ponderación frecuencial.![{\displaystyle \sum _{i=1}^{N}w_{i}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle N}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle w_{i}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
El estimador puede ser insesgado sólo si las ponderaciones no están estandarizadas ni normalizadas , estos procesos cambian la media y la varianza de los datos y, por lo tanto, conducen a una pérdida de la tasa base (el recuento de población, que es un requisito para la corrección de Bessel).
Pesos de confiabilidad
Si, en cambio, las ponderaciones no son aleatorias ( ponderaciones de confiabilidad [ definición necesaria ] ), podemos determinar un factor de corrección para producir un estimador insesgado. Suponiendo que cada variable aleatoria se muestrea de la misma distribución con media y varianza real , tomando las expectativas que tenemos,![{\displaystyle \mu}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \sigma _{\text{actual}}^{2}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\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}}}\\&=\nombreoperador {E} [(X-\nombreoperador {E} [X])^{2}]-{\frac {V_{2}}{V_{1}^{2 }}}\nombredeloperador {E} [(X-\nombredeloperador {E} [X])^{2}]\\&=\left(1-{\frac {V_{2}}{V_{1}^ {2}}}\right)\sigma _ {\text{actual}}^{2}\end{aligned}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
donde y . Por lo tanto, el sesgo en nuestro estimador es , análogo al sesgo en el estimador no ponderado (obsérvese también que es el tamaño de muestra efectivo ). Esto significa que para eliminar el sesgo de nuestro estimador debemos dividir previamente entre , asegurando que el valor esperado de la varianza estimada sea igual a la varianza real de la distribución muestral.![{\displaystyle V_{1}=\sum _{i=1}^{N}w_{i}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle V_{2}=\sum _{i=1}^{N}w_{i}^{2}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \left(1-{\frac {V_{2}}{V_{1}^{2}}}\right)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \left({\frac {N-1}{N}}\right)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \ V_{1}^{2}/V_{2}=N_{eff}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle 1-\left(V_{2}/V_{1}^{2}\right)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
La estimación final insesgada de la varianza muestral es:
[6]
dónde .![{\displaystyle \operatorname {E} [s_{\mathrm {w} }^{2}]=\sigma _ {\text{actual}}^{2}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Los grados de libertad de la 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, a cada vector de fila (cada conjunto de observaciones individuales en cada una de las K variables aleatorias) se le asigna un peso .![{\displaystyle \mathbf {x} _ {i}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle w_{i}\geq 0}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Entonces el vector medio ponderado viene dado por![{\displaystyle \mathbf {\mu ^{*}} }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \mathbf {\mu ^{*}} ={\frac {\sum _{i=1}^{N}w_{i}\mathbf {x} _{i}}{\sum _{i =1}^{N}w_{i}}}.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Y la matriz de covarianza ponderada viene dada por: [8]
![{\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}}}.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
De manera similar a la varianza de la muestra ponderada, existen dos estimadores insesgados diferentes según el tipo de ponderaciones.
Pesos de frecuencia
Si las ponderaciones son ponderaciones de frecuencia , la estimación ponderada insesgada de la matriz de covarianza , con la corrección de Bessel, viene dada por: [8]![{\displaystyle \textstyle \mathbf {C} }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\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}}.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Este estimador puede ser insesgado sólo si las ponderaciones no están estandarizadas ni normalizadas , estos procesos cambian la media y la varianza de los datos y, por lo tanto, conducen a una pérdida de la tasa base (el recuento de población, que es un requisito para la corrección de Bessel).
Pesos de confiabilidad
En el caso de las ponderaciones de confiabilidad , las ponderaciones están normalizadas :
![{\displaystyle V_{1}=\sum _{i=1}^{N}w_{i}=1.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
(Si no es así, divida los pesos por su suma para normalizarlos antes de calcular :![{\ Displaystyle V_ {1}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle w_{i}'={\frac {w_{i}}{\sum _{i=1}^{N}w_{i}}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Entonces el vector medio ponderado se puede simplificar a![{\displaystyle \mathbf {\mu ^{*}} }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \mathbf {\mu ^{*}} =\sum _{i=1}^{N}w_{i}\mathbf {x} _{i}.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
y la estimación ponderada insesgada de la matriz de covarianza es: [9]![{\displaystyle \mathbf {C} }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\begin{alineado}\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}} }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
El razonamiento aquí es el mismo que en el apartado anterior.
Como suponemos que los pesos están normalizados, esto se reduce a:![{\displaystyle V_{1}=1}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\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}}}.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Si todas las ponderaciones son iguales, es decir , entonces la media ponderada y la covarianza se reducen a la media muestral no ponderada y la covarianza anteriores.![{\displaystyle w_{i}/V_{1}=1/N}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Estimaciones valoradas por vectores
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 avión pueden tener menos certeza en una dirección que en otra. Como en el caso escalar, la media ponderada de estimaciones múltiples puede proporcionar una estimación de máxima verosimilitud . Simplemente reemplazamos la varianza por la matriz de covarianza y la inversa aritmética por la matriz inversa (ambas indicadas de la misma manera, mediante superíndices); la matriz de peso luego dice: [10]
![{\displaystyle \mathbf {C} }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \mathbf {W} _{i}=\mathbf {C} _{i}^{-1}.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
La media ponderada en este caso es:
(donde el orden del producto matriz-vector no es conmutativo ), en términos de la covarianza de la media ponderada:![{\displaystyle {\bar {\mathbf {x} }}=\mathbf {C} _ {\bar {\mathbf {x} }}\left(\sum _{i=1}^{n}\mathbf { W} _{i}\mathbf {x} _{i}\right),}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \mathbf {C} _{\bar {\mathbf {x} }}=\left(\sum _{i=1}^{n}\mathbf {W} _{i}\right)^{ -1},}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Por ejemplo, considere 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
![{\displaystyle \mathbf {x} _{1}:={\begin{bmatrix}1&0\end{bmatrix}}^{\top },\qquad \mathbf {C} _{1}:={\begin{ bmatrix}1&0\\0&100\end{bmatrix}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \mathbf {x} _{2}:={\begin{bmatrix}0&1\end{bmatrix}}^{\top },\qquad \mathbf {C} _{2}:={\begin{ bmatrix}100&0\\0&1\end{bmatrix}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
entonces la media ponderada es:
![{\displaystyle {\begin{alineado}{\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}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
lo cual tiene sentido: la estimación [1 0] es "cumple" en el segundo componente y la estimación [0 1] es compatible en el primer componente, por lo que la media ponderada es casi [1 1].
Contabilización de correlaciones
En el caso general, supongamos que , es la matriz de covarianza que relaciona las cantidades , es la media común a estimar y es una matriz de diseño igual a un vector de unos (de longitud ). El teorema de Gauss-Markov establece que la estimación de la media con varianza mínima viene dada por:![{\displaystyle \mathbf {X} =[x_{1},\dots,x_{n}]^{T}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \mathbf {C} }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle x_{i}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\bar {x}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle [1,\puntos,1]^{T}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle n}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \sigma _{\bar {x}}^{2}=(\mathbf {J} ^{T}\mathbf {W} \mathbf {J} )^{-1},}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
y
![{\displaystyle {\bar {x}}=\sigma _ {\bar {x}}^{2}(\mathbf {J} ^{T}\mathbf {W} \mathbf {X} ),}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
dónde:
![{\displaystyle \mathbf {W} =\mathbf {C} ^{-1}.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Disminución de la fuerza de las interacciones.
Considere la serie de tiempo de una variable independiente y una variable dependiente , con observaciones muestreadas en tiempos discretos . En muchas situaciones comunes, el valor de at time depende no sólo de sus valores pasados sino también de ellos. 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 para un tamaño de ventana .![{\displaystyle x}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle y}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle n}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\ Displaystyle t_ {i}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle y}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\ Displaystyle t_ {i}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle x_{i}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle z}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle m}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle z_{k}=\sum _{i=1}^{m}w_{i}x_{k+1-i}.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Pesos exponencialmente decrecientes
En el escenario descrito en la sección anterior, lo más frecuente es que la disminución de la fuerza de interacción obedezca una ley exponencial negativa. Si las observaciones se muestrean en tiempos equidistantes, entonces la disminución exponencial equivale a una disminución en una fracción constante en cada paso de tiempo. Configuración podemos definir pesos normalizados por![{\displaystyle 0<\Delta <1}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle w=1-\Delta}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle m}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle w_{i}={\frac {w^{i-1}}{V_{1}}},}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
donde es la suma de los pesos no normalizados. En este caso es simplemente![{\ Displaystyle V_ {1}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\ Displaystyle V_ {1}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle V_{1}=\sum _{i=1}^{m}{w^{i-1}}={\frac {1-w^{m}}{1-w}},}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
acercándose para valores grandes de .![{\displaystyle V_{1}=1/(1-w)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle m}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
La constante de amortiguación 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 los pesos exponencialmente decrecientes son útiles para hacer una elección adecuada: en el paso , el peso es aproximadamente igual a , el área de la cola al valor , el área de la cabeza . El área de la cola en el paso es . Cuando importan principalmente las observaciones más cercanas y el efecto de las observaciones restantes se puede ignorar con seguridad, entonces elija una opción tal que el área de la cola sea suficientemente pequeña.![{\displaystyle w}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle (1-w)^{-1}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {e^{-1}}(1-w)=0,39(1-w)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle e^{-1}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {1-e^{-1}}=0,61}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle n}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \leq {e^{-n(1-w)}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle n}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle w}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Promedios ponderados de funciones
El concepto de promedio ponderado se puede extender a funciones. [11] Los promedios ponderados de funciones juegan un papel importante en los sistemas de cálculo diferencial e integral ponderado. [12]
Corrección de dispersión excesiva o insuficiente
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. Por lo general, 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 es demasiado grande. La corrección que se debe hacer es![{\displaystyle \chi^{2}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\hat {\sigma }}_{\bar {x}}^{2}=\sigma _{\bar {x}}^{2}\chi _{\nu }^{2}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
¿Dónde está el chi-cuadrado reducido ?![{\displaystyle \chi _{\nu }^{2}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \chi _{\nu }^{2}={\frac {1}{(n-1)}}\sum _{i=1}^{n}{\frac {(x_{i} -{\bar {x}})^{2}}{\sigma _{i}^{2}}};}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
La raíz cuadrada puede denominarse error estándar de la media ponderada (ponderaciones de varianza, escala corregida) .![{\displaystyle {\sombrero {\sigma }}_{\bar {x}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Cuando todas las varianzas de los datos son iguales, se cancelan en la varianza media ponderada, que nuevamente se reduce al error estándar de la media (al cuadrado), formulado en términos de la desviación estándar de la muestra (al cuadrado),![{\displaystyle \sigma _{i}=\sigma _{0}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \sigma _{\bar {x}}^{2}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \sigma _{\bar {x}}^{2}=\sigma ^{2}/n}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \sigma ^{2}={\frac {\sum _{i=1}^{n}(x_{i}-{\bar {x}})^{2}}{n-1} }.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Ver también
Notas
- ^ Técnicamente, se pueden utilizar negativos si todos los valores son cero o negativos. Sin embargo, esto no cumple ninguna función ya que los pesos funcionan como valores absolutos .
Referencias
- ^ abcd Cochran, WG (1977). Técnicas de muestreo (3ª ed.). Nashville, Tennessee: John Wiley & Sons. ISBN 978-0-471-16240-7
- ^ abcdefghijklmnopq Carl-Erik Sarndal; Bengt Swensson; Jan Wretman (1992). Muestreo de encuestas asistido por modelos . ISBN 978-0-387-97528-3.
- ^ Thomas Lumley (https://stats.stackexchange.com/users/249135/thomas-lumley), ¿Cómo estimar la varianza (aproximada) de la media ponderada?, URL (versión: 2021-06-08): https: //stats.stackexchange.com/q/525770
- ^ Gatz, Donald F.; Smith, Lutero (junio de 1995). "El error estándar de una concentración media ponderada: I. Bootstrapping frente a otros métodos". Ambiente Atmosférico . 29 (11): 1185-1193. Código Bib : 1995AtmEn..29.1185G. doi :10.1016/1352-2310(94)00210-C.- enlace pdf
- ^ Endlich, RM; Eymon, BP; Ferek, RJ; Valdés, AD; Maxwell, C. (1 de diciembre de 1988). "Análisis estadístico de las mediciones de la química de la precipitación en el este de los Estados Unidos. Parte I: patrones y correlaciones estacionales y regionales". Revista de Meteorología y Climatología Aplicadas . 27 (12): 1322-1333. Código Bib : 1988JApMe..27.1322E. doi : 10.1175/1520-0450(1988)027<1322:SAOPCM>2.0.CO;2 .
- ^ "Biblioteca científica GNU - Manual de referencia: muestras ponderadas". Gnu.org . Consultado el 22 de diciembre de 2017 .
- ^ "Error estándar ponderado y su impacto en las pruebas de significancia (WinCross vs. Quantum y SPSS), Dr. Albert Madansky" (PDF) . Grupo analítico.com . Consultado el 22 de diciembre de 2017 .
- ^ ab Price, George R. (abril de 1972). "Ampliación de las matemáticas de selección de covarianza" (PDF) . Anales de genética humana . 35 (4): 485–490. doi :10.1111/j.1469-1809.1957.tb01874.x. PMID 5073694. S2CID 37828617.
- ^ 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
- ^ James, Federico (2006). Métodos estadísticos en física experimental (2ª ed.). Singapur: World Scientific. pag. 324.ISBN 981-270-527-9.
- ^ GH Hardy, JE Littlewood y G. Pólya. Desigualdades (2ª ed.), Cambridge University Press, ISBN 978-0-521-35880-4 , 1988.
- ^ Jane Grossman, Michael Grossman, Robert Katz. Los primeros sistemas de cálculo diferencial e integral ponderado, ISBN 0-9771170-1-4 , 1980.
Lectura adicional
- Bevington, Philip R (1969). Reducción de Datos y Análisis de Errores para las Ciencias Físicas . Nueva York, Nueva York: McGraw-Hill. OCLC 300283069.
- Strutz, T. (2010). Ajuste e incertidumbre de datos (una introducción práctica a los mínimos cuadrados ponderados y más) . Vieweg+Teubner. ISBN 978-3-8348-1022-9.
Enlaces externos
- David Terr. "Media ponderada". MundoMatemático .
- Herramienta para calcular el promedio ponderado