1 votos

¿Cómo hacer una estimación de densidad de núcleo para datos almacenados en un sistema polar?

Quiero hacer una estimación de densidad kernel para un conjunto de datos de velocidad del viento. Los datos se almacenan en speed y angle así

speed direction
15 220
23 210
17 200

Las direcciones están sectorizadas en 10 grados.

Quiero saber cómo puedo hacer una estimación de la densidad del núcleo para esto.

Si cambio directamente los datos en X-Y panel, los puntos quedarían así: enter image description here

Y el resultado de KDE se vería así:

enter image description here

Obviamente, el peaks del resultado kde se deben al sectored angles . ¿Cómo puedo obtener un resultado kde suavizado, como éste?

enter image description here

ACTUALIZACIÓN: Después de aleatorizar los ángulos: enter image description here

0voto

MichaelCG8 Puntos 66

Pregunta antigua, pero para quien encuentre esto ahora, los KDE tienen una asociación ancho de banda que es algo que se puede ajustar. Aumentar el ancho de banda suaviza el gráfico. La forma de ajustar esto dependerá del paquete de trazado que esté utilizando, y muchos aplicarán alguna regla para elegir una mejor estimación si no la especifica usted mismo.

https://en.wikipedia.org/wiki/Kernel_density_estimation#Bandwidth_selection

Aquí tienes un ejemplo 1D del efecto de cambiar el ancho de banda, de esa página de Wikipedia. Comparison of kde of 1D distributions with various bandwidths

En el caso de la pregunta original, si se dispone de los datos originales aún no divididos en intervalos de 10 grados, se evitaría la creación de los picos. Esos serían los mejores datos para proporcionar a la función kde.

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X