Hipótesis nula de equivalencia

11

Supongamos que son una muestra aleatoria simple de una distribución Normal . ( μ , σ 2 )X1,X2,...,Xn(μ,σ2)

Estoy interesado en hacer la siguiente prueba de hipótesis: para una constante dada .c > 0

H0:|μ|cH1:|μ|>c,
c>0

Estaba pensando en realizar dos pruebas unilaterales (TOST) de forma análoga a la situación habitual de las pruebas de bioequivalencia, donde nulo y es lugar, pero no sé si esto tiene sentido o es correcto.| μ | ct|μ|c

Mi idea es realizar las pruebas unilaterales y y rechaza la hipótesis nula global si uno de los valores es menor que un nivel de significancia .H 02 : μ - c

H01:μcH11:μ>c
p α
H02:μcH12:μ<c,
pα

¡Gracias por adelantado!

EDITAR:

He estado pensando un poco sobre esto, y creo que el enfoque que propuse no tiene un nivel de significancia .α

Suponga que se conoce el verdadero valor de es y .μ 0 σ 2μμ0σ2

La probabilidad de rechazar el valor nulo en la primera prueba es donde si el cdf estándar de la distribución Normal, y es un valor tal que .Φz1-αΦ(z1-α)=1-α

Pμ0(Rej.H01)=1Φ(z1α+cμ0σ/n),
Φz1αΦ(z1α)=1α

Si , . Entonces, si , . Alternativamente, si , .μ0=cPμ0(Rej.H01)=αμ0>cPμ0(Rej.H01)>αμ0<cPμ0(Rej.H01)<α

La probabilidad de rechazar el valor nulo en la segunda prueba es

Pμ0(Rej.H02)=Φ(z1αμ0+cσ/n).

Nuevamente, si tenemos . Del mismo modo, si , . Finalmente, si , .μ0=cPμ0(Rej.H02)=αμ0>cPμ0(Rej.H02)<αμ0<cPμ0(Rej.H02)>α

Dado que las regiones de rechazo de las dos pruebas son disjuntas, la probabilidad de rechazar es: H0

Pμ0(Rej.H0)=1Φ(z1α+cμ0σ/n)+Φ(z1αμ0+cσ/n)

Entonces, si , es un límite superior de la probabilidad de rechazar la hipótesis nula (global). Por lo tanto, el enfoque que propuse era demasiado liberal.μ[c,c]2α

Si no me equivoco, podemos alcanzar un nivel de significancia de haciendo las mismas dos pruebas y rechazando el valor nulo si el valor de una de ellas es menor que . Un argumento similar se mantiene cuando la varianza es desconocida y necesitamos aplicar la prueba .αpα/2t

Vic101
fuente
La edición está en el camino correcto :-).
whuber

Respuestas:

3

Muy interesante pregunta !!

Está utilizando la consecuencia lógica, es decir, la condición de vinculación. Esta condición de vinculación forma la base misma de la lógica clásica, garantiza la inferencia o deducción de un resultado de una premisa.

El razonamiento detrás de su propuesta es el siguiente:

Si implica , entonces los datos observados deberían generar más evidencia contra que .H0H0H0H0

En términos de sus hipótesis auxiliares y , tenemos , es decir, implica y también implica . Por lo tanto, de acuerdo con la condición de vinculación, debemos observar más evidencia contra que o . Luego, concluyó que si uno de los valores p calculados bajo o es suficientemente pequeño, el valor p calculado bajo será aún más pequeño.H01H02H0H01H02H0H01H0H02H0H01H02H01H02H0

Sin embargo, este razonamiento lógico no es válido para los valores p, es decir, los valores p no respetan la consecuencia lógica. Cada valor p se construye bajo una hipótesis nula específica, por lo tanto, los valores p para diferentes hipótesis nulas se calculan bajo diferentes métricas. Por esta razón, los valores p no pueden respetar el razonamiento lógico sobre el espacio del parámetro (o el espacio de las hipótesis nulas).

En Schervish (1996) y Patriota (2013) se presentan ejemplos de valores p que violan la condición de vinculación. El último artículo muestra ejemplos de una distribución normal bivariada y de un modelo de regresión (ver ejemplos 1.1 y 1.2 en las páginas 5 y 6, respectivamente). Eran Raviv proporciona un algoritmo en código R para el caso bivariado. El aprendizaje de estos ejemplos es: debe calcular el valor p directamente para la hipótesis nula de interés. Schervish (1996) proporciona una fórmula de valor p para su ejemplo cuando y , vea Fórmula (2) en la página 204. Si desea calcular un valor p, debe adecuar esa fórmula para Tu caso.n=1σ2=1

Patriota (2013) propone una nueva medida de evidencia para probar hipótesis nulas generales (hipótesis nulas compuestas o simples) que respeta la consecuencia lógica. Esta medida se llama valor s en el documento. El procedimiento es relativamente simple para su ejemplo:

  1. Encuentre un intervalo de confianza (1- ) para (uno asintótico): , donde es el promedio de la muestra, es la varianza de la muestra , es el cuantil de una distribución normal estándar es el tamaño de la muestra.αμI(μ,α)=[x¯zα/2s2n ; x¯+zα/2s2n]x¯s2zα/2α/2n

  2. Encuentre el valor para el cual la amplitud de es mínima y tiene al menos un elemento en común con (es decir, el borde de ). Este es el valor .αI(μ,α){c,c}[c,c]αs

  3. Por un lado, si , entonces la muestra observada está corroborando con la Hipótesis nula ; Si el valor es lo suficientemente pequeño, puede aceptar el valor nulo. Por otro lado, si , la muestra observada proporciona información contra la hipótesis nula ; Si el valor es lo suficientemente pequeño, puede rechazar el valor nulo. De lo contrario, no debe rechazar ni aceptar el nulo.x¯[c,c]H0:|μ|csx¯[c,c]H0s

Observe que, si y el valor respectivo es extremadamente pequeño, esto significa que la hipótesis alternativa está muy lejos del valor plausible máximo, . Si y el valor respectivo es extremadamente pequeño, esto significa que la hipótesis nula está muy lejos del valor plausible máximo, . Intente dibujar una imagen que represente el intervalo de confianza y la hipótesis nula de interés para comprender mejor las conclusiones. Para obtener más información, lea el documento original Patriota (2013).x¯[c,c]sx¯x¯[c,c]sx¯

Cómo encontrar umbrales objetivos para aceptar o rechazar el valor nulo utilizando este valor- sigue siendo un problema abierto. Este enfoque es bueno porque ahora podemos aceptar una hipótesis nula. Esto tiene sentido siempre que la muestra observada corrobora con el valor nulo y está muy lejos de la alternativa. En su ejemplo se puede observar para , , y . Es bastante simple ver que la densidad de datos está extremadamente concentrada en (diez veces el error estándar). Para tener una intersección no vacía con se requieren 99900 errores estándar. Por lo tanto, sería lo suficientemente justo como para aceptarc = 1000 ˉ x = 1 s 2 = 1 n = 10000 [ 0.9 , 1.1 ] [ - 1000 , 1000 ] H 0 : | μ | csc=1000x¯=1s2=1n=10000[0.9, 1.1][1000, 1000]H0:|μ|c en este caso.

Referencias

Patriota, AG (2013). Una medida clásica de evidencia para hipótesis nulas generales, Fuzzy Sets and Systems, 233, 74–88

Schervish, MJ (1996). Valores P: qué son y qué no son, The American Statistician, 50, 203–206.

Alexandre Patriota
fuente