¿Existe un estimador imparcial de la distancia de Hellinger entre dos distribuciones?


20

En un entorno donde se observa distribuido desde una distribución con densidad , me pregunto si hay un estimador imparcial (basado en 's) de la distancia de Hellinger a otra distribución con densidad , a saber, X1,,XnfXif0

H(f,f0)={1Xf(x)f0(x)dx}1/2.

44
Entonces f0 es conocido y fijo. Pero, ¿se conoce o proviene de una familia paramétrica o está haciendo esto en un marco no paramétrico con todo lo que sabe acerca de si proviene de su muestra? Creo que hace la diferencia al intentar una respuesta.
Michael R. Chernick

3
@MichaelChernick: suponga que todo lo que sabe sobre es la muestra . fX1,,Xn
Xi'an

2
No creo que se haya calculado (si existe). Si existe, entonces AIC tiene un hermano perdido.

44
Un ataque a este problema parece factible si asume que y son discretos. Esto lleva a un estimador obvio (calcule la distancia de Hellinger entre el FED y ). Bootstrapping (en teoría, ¡no a través de la simulación!) Nos dará una idea del posible sesgo, así como una forma de reducir (o incluso eliminar) el sesgo. Espero tener éxito con la distancia al cuadrado en lugar de la distancia en sí, porque es matemáticamente más manejable. El supuesto de una discreta no es un problema en las aplicaciones; el espacio de discreto es un subconjunto denso de todos modos. ff0f0ff
whuber

2
Me viene a la mente la prueba de Rosenblatt de que no existe un estimador imparcial de buena fe de . ¿Podemos superar eso y obtener un estimador no obsesionado de ? No lo sé. fH(f,f0)
Zen

Respuestas:


5

No existe un estimador imparcial de o de para de cualquier clase de distribución no paramétrica razonablemente amplia.HH2f

Podemos mostrar esto con el argumento maravillosamente simple de

Bickel y Lehmann (1969). Estimación imparcial en familias convexas . Los Anales de Estadística Matemática, 40 (5) 1523-1535. ( proyecto euclid )

algunas distribuciones , y , con las densidades correspondientes , y . Deje denotan , y dejar que ser alguna estimador de sobre la base de muestras iid .F0FGf0fgH(F)H(f,f0)H^(X)H(F)nXiF

Suponga que es imparcial para muestras de cualquier distribución de la forma Pero luego para que debe ser un polinomio enH^

Mα:=αF+(1α)G.
Q(α)=H(Mα)=x1xnH^(X)dMα(x1)dMα(xn)=x1xnH^(X)[αdF(x1)+(1α)dG(x1)][αdF(xn)+(1α)dG(xn)]=αnEXFn[H^(X)]++(1α)nEXGn[H^(X)],
Q(α)αde grado como máximo .n

Ahora, especializémonos en un caso razonable y demostremos que la correspondienteQ no es polinómica.

Supongamos que es una distribución que tiene densidad constante en [ - 1 , 1 ] : f 0 ( x ) = c para todos | x | 1 . (Su comportamiento fuera de ese rango no importa). Sea F una distribución admitida solo en [ - 1 , 0 ] , y G alguna distribución admitida solo en [ 0 , 1 ] .F0[1,1]f0(x)=c|x|1F[1,0]G[0,1]

Ahora dondee igualmente para. Tenga en cuenta que,para cualquier distribución,que tenga una densidad.

Q(α)=H(mα,f0)=1Rmα(x)f0(x)dx=110cαf(x)dx01c(1α)g(x)dx=1αBF1αBG,
BF:=Rf(x)f0(x)dxBGBF>0BG>0FG

1αBF1αBG no es un polinomio de ningún grado finito. Por lo tanto, ningún estimador puede ser imparcial para en todas las distribuciones con muchas muestras finitas.H^HMα

Del mismo modo, debido a que tampoco es un polinomio, no hay un estimador para que sea imparcial en todas las distribuciones con finitamente muchas muestras.1αBF1αBGH2Mα

Esto excluye prácticamente todas las clases de distribuciones no paramétricas razonables, excepto aquellas con densidades limitadas a continuación (una suposición que a veces hacen análisis no paramétricos). Probablemente podrías matar esas clases también con un argumento similar simplemente haciendo constantes las densidades o algo así.


13

No sé cómo construir (si existe) un estimador imparcial de la distancia Hellinger. Parece posible construir un estimador consistente. Tenemos una densidad fija conocida , y una muestra aleatoria de una densidad . Queremos estimar donde . Por el SLLN, sabemos que casi seguro, comof0X1,,Xnf>0

H(f,f0)=1Xf(x)f0(x)dx=1Xf0(x)f(x)f(x)dx
=1E[f0(X)f(X)],
Xf
11ni=1nf0(Xi)f(Xi)H(f,f0),
n. Por lo tanto, una forma razonable de estimar sería tomar un estimador de densidad (como un estimador de densidad de núcleo tradicional) de , y calcular H(f,f0)fn^f
H^=11ni=1nf0(Xi)fn^(Xi).

3
@ Zen: Buen punto! Considero esta respuesta como la respuesta porque me hizo darme cuenta de que suena muy parecido a una desviación estándar, para la cual no existe un estimador imparcial. En cuanto a la varianza de , no se preocupe: implica que este estimador tiene una varianza finita. HH^n2E[(f0(X)/f(X))2]=1
Xi'an

1
¡Gracias por la aclaración sobre la varianza del estimador, Xi'an!
Zen

2
Algunos trabajan en otros estimadores consistentes: (a) arxiv.org/abs/1707.03083 y trabajos relacionados basados ​​en estimadores de densidad -NN; (b) arxiv.org/abs/1402.2966 basado en la corrección de las estimaciones de densidad del núcleo; (c) ieeexplore.ieee.org/document/5605355 basado en una conexión a la clasificación. (Muchos de estos se basan en muestras de y , porque ese es el trabajo que conocía sobre la mano, pero creo que hay variantes para conocido .)kff0f0
Dougal
Al usar nuestro sitio, usted reconoce que ha leído y comprende nuestra Política de Cookies y Política de Privacidad.
Licensed under cc by-sa 3.0 with attribution required.