Articulo de referencia

Varianza combinada

En estadística , la varianza combinada (también conocida como varianza combinada , varianza compuesta o varianza general , y escrita) σ 2 {\displaystyle \sigma ^{2}} Es un métod...

En estadística , la varianza combinada (también conocida como varianza combinada , varianza compuesta o varianza general , y escrita)σ2{\displaystyle \sigma ^{2}}Es un método para estimar la varianza de varias poblaciones diferentes cuando la media de cada población puede ser distinta, pero se puede suponer que la varianza de cada población es la misma. La estimación numérica resultante de este método también se denomina varianza combinada.

Bajo el supuesto de varianzas poblacionales iguales, la varianza de la muestra combinada proporciona una estimación de la varianza con mayor precisión que las varianzas de las muestras individuales. Esta mayor precisión puede aumentar la potencia estadística cuando se utiliza en pruebas estadísticas que comparan poblaciones, como la prueba t .

La raíz cuadrada de un estimador de varianza agrupada se conoce como desviación estándar agrupada (también conocida como desviación estándar combinada , desviación estándar compuesta o desviación estándar general ).

Motivación

En estadística , a menudo se recopilan datos de una variable dependiente , y , para un rango de valores de la variable independiente , x . Por ejemplo, se puede estudiar el consumo de combustible en función de la velocidad del motor, manteniendo constante la carga. Si para obtener una varianza pequeña en y se requieren numerosas pruebas repetidas para cada valor de x , el costo de las pruebas puede resultar prohibitivo. Se pueden obtener estimaciones razonables de la varianza utilizando el principio de varianza combinada tras repetir cada prueba para un valor de x determinado solo unas pocas veces.

Definición y cálculo

La varianza combinada es una estimación de la varianza común fija.σ2{\displaystyle \sigma ^{2}}subyacentes a diversas poblaciones que tienen diferentes medios.

Se nos proporciona un conjunto de varianzas muestrales.si2{\displaystyle s_{i}^{2}}donde las poblaciones están indexadasi=1,,metro{\displaystyle i=1,\ldots ,m},

si2{\displaystyle s_{i}^{2}}=1nortei1j=1nortei(yi,jy¯i)2.{\displaystyle {\frac {1}{n_{i}-1}}\sum _{j=1}^{n_{i}}\left(y_{i,j}-{\overline {y}}_{i}\right)^{2}.}

Suponiendo tamaños de muestra uniformes ,nortei=norte{\displaystyle n_{i}=n}, entonces la varianza combinadaspag2{\displaystyle s_{p}^{2}}se puede calcular mediante la media aritmética :

spag2=i=1metrosi2metro=s12+s22++smetro2metro.{\displaystyle s_{p}^{2}={\frac {\sum _{i=1}^{m}s_{i}^{2}}{m}}={\frac {s_{1}^{2}+s_{2}^{2}+\cdots +s_{m}^{2}}{m}}.}

Si los tamaños de muestra no son uniformes, entonces la varianza combinadaspag2{\displaystyle s_{p}^{2}}se puede calcular mediante el promedio ponderado , utilizando como ponderacioneswi=nortei1{\displaystyle w_{i}=n_{i}-1}los respectivos grados de libertad (véase también: corrección de Bessel ):

spag2=i=1metro(nortei1)si2i=1metro(nortei1)=(norte11)s12+(norte21)s22++(nortemetro1)smetro2norte1+norte2++nortemetrometro.{\displaystyle s_{p}^{2}={\frac {\sum _{i=1}^{m}(n_{i}-1)s_{i}^{2}}{\sum _{i=1}^{m}(n_{i}-1)}}={\frac {(n_{1}-1)s_{1}^{2}+(n_{2}-1)s_{2}^{2}+\cdots +(n_{m}-1)s_{m}^{2}}{n_{1}+n_{2}+\cdots +n_{m}-m}}.}

La distribución despag2(inorteimetro)σ2{\displaystyle {\frac {s_{p}^{2}(\sum _{i}n_{i}-m)}{\sigma ^{2}}}}esχ2(inorteimetro){\displaystyle \chi ^{2}(\sum _{i}n_{i}-m)}.

Demostración. Cuando hay una sola media, la distribución de(y1y¯,,ynortey¯){\displaystyle (y_{1}-{\bar {y}},\dots ,y_{n}-{\bar {y}})}es una gaussiana enΔnorte1{\displaystyle \Delta _{n-1}}, el(norte1){\displaystyle (n-1)}Simplex -dimensional, con desviación estándarσ{\displaystyle \sigma }. Donde hay múltiples medios, la distribución de(y1,1y¯1,,y1,norte1y¯1,,ymetro,1y¯metro,,ymetro,nortemetroy¯metro){\displaystyle (y_{1,1}-{\bar {y}}_{1},\dots ,y_{1,n_{1}}-{\bar {y}}_{1},\dots ,y_{m,1}-{\bar {y}}_{m},\dots ,y_{m,n_{m}}-{\bar {y}}_{m})}es una gaussiana enΔnorte11××Δnortemetro1{\displaystyle \Delta _{n_{1}-1}\times \dots \times \Delta _{n_{m}-1}}.

Variantes

La estimación insesgada de mínimos cuadrados deσ2{\displaystyle \sigma ^{2}}(como se presenta arriba), y la estimación de máxima verosimilitud sesgada que se muestra a continuación:

spag2=i=1norte(nortei1)si2i=1nortenortei,{\displaystyle s_{p}^{2}={\frac {\sum _{i=1}^{N}(n_{i}-1)s_{i}^{2}}{\sum _{i=1}^{N}n_{i}}},}

se utilizan en diferentes contextos. El primero puede dar una opinión imparcial.spag2{\displaystyle s_{p}^{2}}estimarσ2{\displaystyle \sigma ^{2}}cuando los dos grupos comparten una varianza poblacional igual. Este último puede dar una respuesta más eficiente.spag2{\displaystyle s_{p}^{2}}estimarσ2{\displaystyle \sigma ^{2}}, aunque sujeto a sesgo. Nótese que las cantidadessi2{\displaystyle s_{i}^{2}}En el lado derecho de ambas ecuaciones se encuentran las estimaciones insesgadas.

Ejemplo

Considere el siguiente conjunto de datos para y obtenidos en varios niveles de la variable independiente x . 

En la siguiente tabla se presentan el número de ensayos, la media, la varianza y la desviación estándar.

Estas estadísticas representan la varianza y la desviación estándar para cada subconjunto de datos en los distintos niveles de x . Si podemos suponer que los mismos fenómenos generan un error aleatorio en cada nivel de x , los datos anteriores se pueden "agrupar" para expresar una única estimación de la varianza y la desviación estándar. En cierto sentido, esto sugiere encontrar una varianza media o una desviación estándar entre los cinco resultados anteriores. Esta varianza media se calcula ponderando los valores individuales con el tamaño del subconjunto para cada nivel de x . Por lo tanto, la varianza agrupada se define por

spag2=(norte11)s12+(norte21)s22++(nortek1)sk2(norte11)+(norte21)++(nortek1){\displaystyle s_{p}^{2}={\frac {(n_{1}-1)s_{1}^{2}+(n_{2}-1)s_{2}^{2}+\cdots +(n_{k}-1)s_{k}^{2}}{(n_{1}-1)+(n_{2}-1)+\cdots +(n_{k}-1)}}}

donde n 1 , n 2 , . . ., n k son los tamaños de los subconjuntos de datos en cada nivel de la variable x , y s 1 2 , s 2 2 , . . ., s k 2 son sus respectivas varianzas.

Por lo tanto, la varianza combinada de los datos mostrados anteriormente es:

spag2=2.764{\displaystyle s_{p}^{2}=2.764\,}

Efecto sobre la precisión

La varianza combinada es una estimación que se obtiene cuando existe correlación entre los conjuntos de datos combinados o cuando el promedio de los conjuntos de datos no es idéntico. La varianza combinada es menos precisa cuanto mayor sea la correlación (distinta de cero) o mayor la diferencia entre los promedios de los conjuntos de datos.

La variación de datos para conjuntos de datos no superpuestos es:

σincógnita2=i[(norteincógnitai1)σincógnitai2+norteincógnitaiμincógnitai2][inorteincógnitai]μincógnita2inorteincógnitai1{\displaystyle \sigma _{X}^{2}={\frac {\sum _{i}\left[(N_{X_{i}}-1)\sigma _{X_{i}}^{2}+N_{X_{i}}\mu _{X_{i}}^{2}\right]-\left[\sum _{i}N_{X_{i}}\right]\mu _{X}^{2}}{\sum _{i}N_{X_{i}}-1}}}

donde la media se define como:

μincógnita=inorteincógnitaiμincógnitaiinorteincógnitai{\displaystyle \mu _{X}={\frac {\sum _{i}N_{X_{i}}\mu _{X_{i}}}{\sum _{i}N_{X_{i}}}}}

Dada una máxima verosimilitud sesgada definida como:

spag2=i=1k(nortei1)si2i=1knortei,{\displaystyle s_{p}^{2}={\frac {\sum _{i=1}^{k}(n_{i}-1)s_{i}^{2}}{\sum _{i=1}^{k}n_{i}}},}

Entonces, el error en la estimación de máxima verosimilitud sesgada es:

Error=spag2σincógnita2=i(norteincógnitai1)si2inorteincógnitai1inorteincógnitai1(i[(norteincógnitai1)σincógnitai2+norteincógnitaiμincógnitai2][inorteincógnitai]μincógnita2){\displaystyle {\begin{aligned}{\text{Error}}&=s_{p}^{2}-\sigma _{X}^{2}\\[6pt]&={\frac {\sum _{i}(N_{X_{i}}-1)s_{i}^{2}}{\sum _{i}N_{X_{i}}}}-{\frac {1}{\sum _{i}N_{X_{i}}-1}}\left(\sum _{i}\left[(N_{X_{i}}-1)\sigma _{X_{i}}^{2}+N_{X_{i}}\mu _{X_{i}}^{2}\right]-\left[\sum _{i}N_{X_{i}}\right]\mu _{X}^{2}\right)\end{aligned}}}

Suponiendo que N es grande de tal manera que:

inorteincógnitaiinorteincógnitai1{\displaystyle \sum _{i}N_{X_{i}}\approx \sum _{i}N_{X_{i}}-1}

Entonces, el error en la estimación se reduce a:

mi=(i[norteincógnitaiμincógnitai2][inorteincógnitai]μincógnita2)inorteincógnitai=μincógnita2i[norteincógnitaiμincógnitai2]inorteincógnitai{\displaystyle {\begin{aligned}E&=-{\frac {\left(\sum _{i}\left[N_{X_{i}}\mu _{X_{i}}^{2}\right]-\left[\sum _{i}N_{X_{i}}\right]\mu _{X}^{2}\right)}{\sum _{i}N_{X_{i}}}}\\[3pt]&=\mu _{X}^{2}-{\frac {\sum _{i}\left[N_{X_{i}}\mu _{X_{i}}^{2}\right]}{\sum _{i}N_{X_{i}}}}\end{aligned}}}

O alternativamente:

mi=[inorteincógnitaiμincógnitaiinorteincógnitai]2i[norteincógnitaiμincógnitai2]inorteincógnitai=[inorteincógnitaiμincógnitai]2inorteincógnitaii[norteincógnitaiμincógnitai2][inorteincógnitai]2{\displaystyle {\begin{aligned}E&=\left[{\frac {\sum _{i}N_{X_{i}}\mu _{X_{i}}}{\sum _{i}N_{X_{i}}}}\right]^{2}-{\frac {\sum _{i}\left[N_{X_{i}}\mu _{X_{i}}^{2}\right]}{\sum _{i}N_{X_{i}}}}\\[3pt]&={\frac {\left[\sum _{i}N_{X_{i}}\mu _{X_{i}}\right]^{2}-\sum _{i}N_{X_{i}}\sum _{i}\left[N_{X_{i}}\mu _{X_{i}}^{2}\right]}{\left[\sum _{i}N_{X_{i}}\right]^{2}}}\end{aligned}}}

Agregación de datos de desviación estándar

En lugar de estimar la desviación estándar combinada, la siguiente es la forma de agregar con precisión la desviación estándar cuando se dispone de más información estadística.

Estadísticas basadas en la población

Las poblaciones de conjuntos que pueden superponerse se pueden calcular simplemente de la siguiente manera:

norteincógnitaY=norteincógnita+norteYnorteincógnitaY{\displaystyle {\begin{aligned}&&N_{X\cup Y}&=N_{X}+N_{Y}-N_{X\cap Y}\\\end{aligned}}}

Las poblaciones de conjuntos que no se superponen se pueden calcular simplemente de la siguiente manera:

incógnitaY=norteincógnitaY=0norteincógnitaY=norteincógnita+norteY{\displaystyle {\begin{aligned}X\cap Y=\varnothing &\Rightarrow &N_{X\cap Y}&=0\\&\Rightarrow &N_{X\cup Y}&=N_{X}+N_{Y}\end{aligned}}}

Las desviaciones estándar de subpoblaciones no superpuestas ( XY = ∅ ) se pueden agregar de la siguiente manera si se conocen el tamaño (real o relativo entre sí) y las medias de cada una:

μincógnitaY=norteincógnitaμincógnita+norteYμYnorteincógnita+norteYσincógnitaY=norteincógnitaσincógnita2+norteYσY2norteincógnita+norteY+norteincógnitanorteY(norteincógnita+norteY)2(μincógnitaμY)2{\displaystyle {\begin{aligned}\mu _{X\cup Y}&={\frac {N_{X}\mu _{X}+N_{Y}\mu _{Y}}{N_{X}+N_{Y}}}\\[3pt]\sigma _{X\cup Y}&={\sqrt {{\frac {N_{X}\sigma _{X}^{2}+N_{Y}\sigma _{Y}^{2}}{N_{X}+N_{Y}}}+{\frac {N_{X}N_{Y}}{(N_{X}+N_{Y})^{2}}}(\mu _{X}-\mu _{Y})^{2}}}\end{aligned}}}

Por ejemplo, supongamos que se sabe que el hombre estadounidense promedio tiene una altura media de 70  pulgadas con una desviación estándar de tres pulgadas y que la mujer estadounidense promedio tiene una altura media de 65  pulgadas con una desviación estándar de dos pulgadas. Supongamos también que el número de hombres, N , es igual al número de mujeres. Entonces, la media y la desviación estándar de las alturas de los adultos estadounidenses podrían calcularse como:

μ=norte70+norte65norte+norte=70+652=67,5σ=32+222+(7065)222=12,753.57{\displaystyle {\begin{aligned}\mu &={\frac {N\cdot 70+N\cdot 65}{N+N}}={\frac {70+65}{2}}=67.5\\[3pt]\sigma &={\sqrt {{\frac {3^{2}+2^{2}}{2}}+{\frac {(70-65)^{2}}{2^{2}}}}}={\sqrt {12.75}}\approx 3.57\end{aligned}}}

Para el caso más general de M poblaciones no superpuestas, X 1 a X M , y la población agregadaincógnita=iincógnitai{\textstyle X\,=\,\bigcup _{i}X_{i}},

μincógnita=inorteincógnitaiμincógnitaiinorteincógnitaiσincógnita=inorteincógnitaiσincógnitai2inorteincógnitai+i<jnorteincógnitainorteincógnitaj(μincógnitaiμincógnitaj)2(inorteincógnitai)2{\displaystyle {\begin{aligned}\mu _{X}&={\frac {\sum _{i}N_{X_{i}}\mu _{X_{i}}}{\sum _{i}N_{X_{i}}}}\\[3pt]\sigma _{X}&={\sqrt {{\frac {\sum _{i}N_{X_{i}}\sigma _{X_{i}}^{2}}{\sum _{i}N_{X_{i}}}}+{\frac {\sum _{i<j}N_{X_{i}}N_{X_{j}}(\mu _{X_{i}}-\mu _{X_{j}})^{2}}{{\big (}\sum _{i}N_{X_{i}}{\big )}^{2}}}}}\end{aligned}}},

dónde

incógnitaiincógnitaj=, i<j.{\displaystyle X_{i}\cap X_{j}=\varnothing ,\quad \forall \ i<j.}

Si se conocen el tamaño (real o relativo entre sí), la media y la desviación estándar de dos poblaciones superpuestas, así como su intersección, entonces la desviación estándar de la población total aún se puede calcular de la siguiente manera:

μincógnitaY=1norteincógnitaY(norteincógnitaμincógnita+norteYμYnorteincógnitaYμincógnitaY)σincógnitaY=1norteincógnitaY(norteincógnita[σincógnita2+μincógnita2]+norteY[σY2+μY2]norteincógnitaY[σincógnitaY2+μincógnitaY2])μincógnitaY2{\displaystyle {\begin{aligned}\mu _{X\cup Y}&={\frac {1}{N_{X\cup Y}}}\left(N_{X}\mu _{X}+N_{Y}\mu _{Y}-N_{X\cap Y}\mu _{X\cap Y}\right)\\[3pt]\sigma _{X\cup Y}&={\sqrt {{\frac {1}{N_{X\cup Y}}}\left(N_{X}[\sigma _{X}^{2}+\mu _{X}^{2}]+N_{Y}[\sigma _{Y}^{2}+\mu _{Y}^{2}]-N_{X\cap Y}[\sigma _{X\cap Y}^{2}+\mu _{X\cap Y}^{2}]\right)-\mu _{X\cup Y}^{2}}}\end{aligned}}}

Si se suman dos o más conjuntos de datos punto por punto, la desviación estándar del resultado se puede calcular si se conoce la desviación estándar de cada conjunto de datos y la covarianza entre cada par de conjuntos de datos:

σincógnita=iσincógnitai2+2i,jcobertura(incógnitai,incógnitaj){\displaystyle \sigma _{X}={\sqrt {\sum _{i}{\sigma _{X_{i}}^{2}}+2\sum _{i,j}\operatorname {cov} (X_{i},X_{j})}}}

En el caso especial en que no existe correlación entre ningún par de conjuntos de datos, la relación se reduce a la raíz cuadrada de la suma de los cuadrados:

cobertura(incógnitai,incógnitaj)=0,i<jσincógnita=iσincógnitai2.{\displaystyle {\begin{aligned}&\operatorname {cov} (X_{i},X_{j})=0,\quad \forall i<j\\\Rightarrow &\;\sigma _{X}={\sqrt {\sum _{i}{\sigma _{X_{i}}^{2}}}}.\end{aligned}}}

Estadísticas basadas en muestras

Las desviaciones estándar de submuestras no superpuestas ( XY = ∅ ) se pueden agregar de la siguiente manera si se conocen el tamaño real y las medias de cada una:

μincógnitaY=1norteincógnitaY(norteincógnitaμincógnita+norteYμY)σincógnitaY=1norteincógnitaY1([norteincógnita1]σincógnita2+norteincógnitaμincógnita2+[norteY1]σY2+norteYμY2[norteincógnita+norteY]μincógnitaY2){\displaystyle {\begin{aligned}\mu _{X\cup Y}&={\frac {1}{N_{X\cup Y}}}\left(N_{X}\mu _{X}+N_{Y}\mu _{Y}\right)\\[3pt]\sigma _{X\cup Y}&={\sqrt {{\frac {1}{N_{X\cup Y}-1}}\left([N_{X}-1]\sigma _{X}^{2}+N_{X}\mu _{X}^{2}+[N_{Y}-1]\sigma _{Y}^{2}+N_{Y}\mu _{Y}^{2}-[N_{X}+N_{Y}]\mu _{X\cup Y}^{2}\right)}}\end{aligned}}}

Para el caso más general de M conjuntos de datos no superpuestos, X 1 a X M , y el conjunto de datos agregadoincógnita=iincógnitai{\textstyle X\,=\,\bigcup _{i}X_{i}},

μincógnita=1inorteincógnitai(inorteincógnitaiμincógnitai)σincógnita=1inorteincógnitai1(i[(norteincógnitai1)σincógnitai2+norteincógnitaiμincógnitai2][inorteincógnitai]μincógnita2){\displaystyle {\begin{aligned}\mu _{X}&={\frac {1}{\sum _{i}{N_{X_{i}}}}}\left(\sum _{i}{N_{X_{i}}\mu _{X_{i}}}\right)\\[3pt]\sigma _{X}&={\sqrt {{\frac {1}{\sum _{i}{N_{X_{i}}-1}}}\left(\sum _{i}{\left[(N_{X_{i}}-1)\sigma _{X_{i}}^{2}+N_{X_{i}}\mu _{X_{i}}^{2}\right]}-\left[\sum _{i}{N_{X_{i}}}\right]\mu _{X}^{2}\right)}}\end{aligned}}}

dónde

incógnitaiincógnitaj=,i<j.{\displaystyle X_{i}\cap X_{j}=\varnothing ,\quad \forall i<j.}

Si se conocen el tamaño, la media y la desviación estándar de dos muestras superpuestas, así como su intersección, entonces aún se puede calcular la desviación estándar de la muestra agregada. En general,

μincógnitaY=1norteincógnitaY(norteincógnitaμincógnita+norteYμYnorteincógnitaYμincógnitaY)σincógnitaY=[norteincógnita1]σincógnita2+norteincógnitaμincógnita2+[norteY1]σY2+norteYμY2[norteincógnitaY1]σincógnitaY2norteincógnitaYμincógnitaY2[norteincógnita+norteYnorteincógnitaY]μincógnitaY2norteincógnitaY1{\displaystyle {\begin{aligned}\mu _{X\cup Y}&={\frac {1}{N_{X\cup Y}}}\left(N_{X}\mu _{X}+N_{Y}\mu _{Y}-N_{X\cap Y}\mu _{X\cap Y}\right)\\[3pt]\sigma _{X\cup Y}&={\sqrt {\frac {[N_{X}-1]\sigma _{X}^{2}+N_{X}\mu _{X}^{2}+[N_{Y}-1]\sigma _{Y}^{2}+N_{Y}\mu _{Y}^{2}-[N_{X\cap Y}-1]\sigma _{X\cap Y}^{2}-N_{X\cap Y}\mu _{X\cap Y}^{2}-[N_{X}+N_{Y}-N_{X\cap Y}]\mu _{X\cup Y}^{2}}{N_{X\cup Y}-1}}}\end{aligned}}}

Véase también

Referencias

  • Killeen PR (mayo de 2005). "Una alternativa a las pruebas de significación de hipótesis nula" . Psychol Sci . 16 (5): 345– 53. doi : 10.1111/j.0956-7976.2005.01538.x . PMC 1473027. PMID 15869691 .  
  • Libro de Oro de la IUPAC: desviación estándar combinada
  • – haciendo referencia también a la d de Cohen (en la página 6)