En matemáticas y estadística, una media circular (o también media angular) es un promedio concebido para ser aplicado con ángulos y cantidades cíclicas similares, como las horas del día y la parte fraccionaria de los números reales.
Su uso se justifica teniendo en cuenta que la mayoría de las medias habituales pueden no ser adecuadas para cantidades en forma de ángulo. Por ejemplo, la media aritmética de 0° y 360° es 180°, lo cual es engañoso porque 360° equivale a 0° si se considera el módulo de un ciclo completo.[1] Otro ejemplo es la "hora promedio" entre las 11 p. m. y la 1 a. m., que puede ser medianoche o mediodía, dependiendo de si las dos horas son parte de una sola noche o de un solo día del calendario.
La media circular es uno de los ejemplos más simples de estadísticas direccionales y de estadísticas en espacios no euclídeos. Este cálculo produce un resultado diferente al de la media aritmética, siendo la diferencia mayor cuando los ángulos están ampliamente distribuidos. Por ejemplo, la media aritmética de los tres ángulos 0°, 0° y 90° es (0° + 0° + 90°) / 3 = 30°, pero la media vectorial (el ángulo del vector resultante de sumar tres vectores unitarios, uno vertical y los otros dos horizontales) es arctan(1/2) = 26,565°. Además, con la media aritmética la varianza circular solo se define en ±180°.
Definición
Dado que la media aritmética no siempre es apropiada para los ángulos, se puede utilizar el siguiente método para obtener tanto el valor medio como la medida de la varianza de los ángulos:
Conviértanse todos los ángulos a los puntos correspondientes de la circunferencia goniométrica, por ejemplo, a . Es decir, convertir las coordenadas polares a coordenadas cartesianas. Luego, calcular la media aritmética de estos puntos. El punto resultante estará dentro del disco unidad, pero generalmente no quedará sobre la círcunferencia unitaria. Por último, convertir ese punto nuevamente a coordenadas polares. El ángulo es una media razonable de los ángulos de entrada. El radio resultante será 1 si todos los ángulos son iguales. Si los ángulos están distribuidos uniformemente en la circunferencia, entonces el radio resultante será 0, y no existe una media circular (de hecho, es imposible definir una operación media continua en el círculo). En otras palabras, el radio mide la concentración de los ángulos.
Dados los ángulos , una fórmula común para obtener la media es una variante basada en la función trigonométrica inversa arcotangente de dos parámetros:
Usando aritmética compleja
Se puede formular una definición equivalente utilizando números complejos:
- .
Para hacer coincidir la fórmula anterior con la que utiliza la media aritmética de los puntos, las sumas tendrían que dividirse por . Sin embargo, la escala no influye en el resultado de las funciones y , por lo que se puede omitir.
Esto puede expresarse de manera más sucinta al darse cuenta de que los datos direccionales son, de hecho, vectores de longitud unitaria. En el caso de datos unidimensionales, estos puntos de datos se pueden representar convenientemente como números complejos de magnitud unitaria , donde es el ángulo medido. El vector resultante medio para la muestra es entonces:
El ángulo medio muestral es entonces el argumento del vector resultante medio:
La longitud del vector resultante medio muestral es:
y tendrá un valor entre 0 y 1. Por lo tanto, el vector resultante medio de la muestra se puede representar como:
También se utilizan cálculos similares para definir la varianza circular.
Propiedades
La media circular,
- La distancia es igual a la mitad de la distancia euclidiana al cuadrado entre los dos puntos de la circunferencia unitaria asociada con y .
Ejemplo
Una forma sencilla de calcular la media de una serie de ángulos (en el intervalo [0°, 360°)) es calcular la media de los cosenos y senos de cada ángulo, y obtener el ángulo calculando la tangente inversa. Considérense los siguientes tres ángulos como ejemplo: 10, 20 y 30 grados. Intuitivamente, calcular la media implicaría sumar estos tres ángulos y dividirlos por 3, lo que en este caso daría como resultado un ángulo medio correcto de 20 grados. Al girar este sistema en sentido antihorario 15 grados, los tres ángulos se convierten en 355 grados, 5 grados y 15 grados. La media aritmética ahora es 125 grados, que es la respuesta incorrecta, ya que debería ser 5 grados. La media vectorial se puede calcular de la siguiente manera, utilizando el seno medio y el coseno medio :
Implementación
En este código Python se usan las horas del día para encontrar su promedio circular:
import math
def circular_mean(hours):
# Convert hours to radians
# What is the 15?! (24*15=360)
radians= [math.radians(hour * 15) for hour in hours]
# Calculate the sum of sin and cos values
sin_sum= sum([math.sin(rad) for rad in radians])
cos_sum= sum([math.cos(rad) for rad in radians])
# Calculate the circular mean using arctan2
mean_rad= math.atan2(sin_sum, cos_sum)
# Convert the mean back to hours
mean_hour= (math.degrees(mean_rad) / 15) % 24
return mean_hour
# Example usage:
hours= [0, 12, 18]
mean_hour= circular_mean(hours)
print("First Circular mean:", round(mean_hour, 2))
hours= [0, 12]
mean_hour= circular_mean(hours)
print("Second Circular mean:", round(mean_hour, 2))
hours= [0, 0, 12, 12, 24]
mean_hour= circular_mean(hours)
print("Third Circular mean:", round(mean_hour, 2))
Generalizaciones
Una serie de N vectores unitarios independentes se extraen de una distribución de von Mises-Fisher.
La estimación de máxima verosimilitud de la dirección media es simplemente la media aritmética normalizada, un valor estadístico suficiente:[2]
Se puede definir una media esférica ponderada basándose en la interpolación lineal esférica.[3]
Véase también
Referencias
- ↑ Christopher M. Bishop: Pattern Recognition and Machine Learning (Information Science and Statistics), ISBN 0-387-31073-8
- ↑ Mardia, Kanti; Jupp, P. E. (1999). Directional Statistics. John Wiley & Sons Ltd. ISBN 978-0-471-95333-3.
- ↑ Buss, Samuel R.; Fillmore, Jay P. (2001). «Spherical averages and applications to spherical splines and interpolation». ACM Transactions on Graphics (Association for Computing Machinery (ACM)) 20 (2): 95-126. ISSN 0730-0301. doi:10.1145/502122.502124.
Lectura adicional
Enlaces externos