Sesgo para el estimador de densidad del núcleo (caso periódico)

8

Estimador de densidad de Kernel está dada por f ( x , h ) = 1

f^(x,h)=1nhi=1nK(xXih)
donde iid con alguna densidad desconocida,- ancho de banda,X1,...Xnfh

K - función del núcleo ( , , ). El sesgo se puede calcular utilizando la expansión de Taylor: K(x)dx=1K(x)xdx=0K(x)x2dx<=- K(y)(f(x)hy+1

1hK(xyh)f(y)dyf(x)=K(y)(f(xhy)f(x))dy
=-K(y)(F(X)hy+12F(X)(hy)2+o(h2))rey=12F(X)h2+o(h2)

Cómo lidiar con el núcleo periódico F ( 0 01K(X)reX=1 , 0 01K(X)XreX=0 0 , 0 01K(X)X2reX< )?

¿Cómo puedo usar la expansión taylor? ( -No puedo usar las propiedades del kernel)0 011hK(y-Xh)F(y)rey=-Xh1-XhK(y)F(X-yh)rey0 01K(y)F(X-yh)rey

¿Podría recomendar un buen libro sobre el suavizado del núcleo para datos circulares?

Katja
fuente

Respuestas:

5

Un rápido google saca a relucir esto, lo que indica que al trabajar con datos circulares necesitará una definición diferente de 'sesgo' para comenzar:

Sin embargo, cuando el uso de datos en el círculo, no podemos utilizar la distancia en el espacio euclidiano, por lo que todas las diferencias θ - θ i debería sustituirse por teniendo en cuenta el ángulo entre dos vectores:

reyoθ)=θ-θyo=min(El |θ-θyoEl |,2π-El |θ-θyoEl |).

- Charles C Taylor. Selección automática de ancho de banda para la estimación de densidad circular. Estadísticas computacionales y análisis de datos Volumen 52, número 7, 15 de marzo de 2008, páginas 3493-3500. doi: 10.1016 / j.csda.2007.11.003

Él hace referencia a estos libros:

S. Rao Jammalamadaka y A. SenGupta, Temas en estadísticas circulares , World Scientific, Singapur (2001).

KV Mardia y PE Jupp, Directional Statistics , John Wiley, Chichester (1999).

una parada
fuente