Eu estou a tentar ajustar a distribuição log-normal a estes dados mas não sei como descobrir a média, o desvio-padrão e ajustar a curva da distribuição em Python (estou a usar o Jupyter no Anaconda). Já tinha tentado fazer isso com a distribuição normal, mas olhando melhor para o gráfico, parece que o que se ajusta mais aos dados é uma distribuição log-normal em vez de uma normal. Poderiam dizer-me como posso alterar o código e colocar uma distribuição log-normal em vez de normal, se faz favor, para o mesmo conjunto de dados?
#amostra 22A
import pandas as pd
import matplotlib.pyplot as plt
import numpy as np
from scipy.stats import norm
import matplotlib.ticker as tkr
import scipy, pylab
import locale
locale.setlocale(locale.LC_NUMERIC, "de_DE")
plt.rcParams['axes.formatter.use_locale'] = True
df=pd.DataFrame({'frequencia_relativa': {'0-1': 0.282, '1-2': 0.316, '2-3': 0.171, '3-4': 0.068, '4-5': 0.085, '5-6': 0.026, '6-7': 0.026, '7-8': 0.009, '8-9': 0.009, '9-10': 0.0, '10-11': 0.0, '11-12': 0.009}})
df=df.reindex(['0-1', '1-2', '2-3', '3-4', '4-5', '5-6', '6-7',
'7-8', '8-9', '9-10', '10-11', '11-12'])
plt.rcParams["figure.figsize"] = [15,8]
fig, ax =plt.subplots()
ax.bar(x=df.index, height=df.frequencia_relativa, alpha=0.5, width=0.9)
ax.set_xlabel('Tamanho lateral do triângulo ($\mu m$)', fontsize=22)
ax.set_ylabel('Frequência relativa', fontsize=22)
x_axis = np.arange(0, 10, 0.001)
plt.xticks(fontsize=18)
plt.yticks(fontsize=18)
#normal com média 2,232 e variância 1.88, conforme planilha
ax.plot(x_axis, norm.pdf(x_axis,2.232,1.888), linewidth=3)
#plt.show()
plt.savefig('output.png', dpi=500, bbox_inches='tight')
norm
porlognorm
no código. Veja documentação: docs.scipy.org/doc/scipy/reference/generated/…