1

Como podría generar histogramas con distinto ancho de intervalo (clase), considerando que tengo los datos del intervalo y su respectiva frecuencia.

    inf = [0,1,2,3,4,5,10,15,20,25,30] #limite inferior

    sup = [1,2,3,4,5,10,15,20,25,30,50] #limite superior

    frec = [18,165,112,87,43,117,51,32,10,8,8] # frecuencia
2
  • Has intentado algo?
    – Alfabravo
    Commented el 19 abr. 2020 a las 4:56
  • Mayúsculas == Gritar. Edita el título por favor, gracias :D
    – fredyfx
    Commented el 19 abr. 2020 a las 22:26

1 respuesta 1

0

plt.bar acepta unos parámetros útiles en este caso:

  • x: las posiciones en el eje x
  • height: las alturas de cada barra
  • width: la anchura de cada barra (o una anchura para cada una por separado, o una anchura igual para todas)
  • align: center (por defecto, centrado en las posiciones de x) o edge (las barras empiezan en las posiciones de x)

Para las anchuras, hay que calcular las diferencias entre sup e inf. Lo más fácil va con numpy, que lo hace en una sola operación.

import matplotlib.pyplot as plt
import numpy as np

inf = [0, 1, 2, 3, 4, 5, 10, 15, 20, 25, 30]  # limite inferior
sup = [1, 2, 3, 4, 5, 10, 15, 20, 25, 30, 50]  # limite superior
frec = [18, 165, 112, 87, 43, 117, 51, 32, 10, 8, 8]  # frecuencia

inf = np.array(inf)
sup = np.array(sup)
frec = np.array(frec)
plt.bar (inf, frec, width=sup-inf, align='edge', ec='white')

plt.xticks(np.concatenate((inf, sup)))
plt.show()

ejemplo de la histograma

PD: Para obtener una histograma con la densidad (cada barra con una superficie relativa a su participación al total, y todas sumando a uno), se puede hacer calculando frec/(sup-inf)/frec.sum().

También se puede aproximar un kde tomando una muestra artificial:

suma = frec.sum()
plt.bar(inf, frec/(sup-inf)/suma, width=sup-inf, align='edge', ec='white')

muestra = np.concatenate([np.linspace(i, s, f) for i, s, f in zip(inf, sup, frec)])

from scipy.stats import gaussian_kde
kde = gaussian_kde(muestra)
x = np.linspace(inf[0], sup[-1], 500)
plt.plot(x, kde(x), color='crimson')
plt.xticks(np.concatenate((inf, sup)))
plt.show()

densidad y kde

1

Tu Respuesta

By clicking “Publica tu respuesta”, you agree to our terms of service and acknowledge you have read our privacy policy.

¿No es la respuesta que buscas? Examina otras preguntas con la etiqueta o formula tu propia pregunta.