Хорошо, вот, как я думаю, положение вещей сейчас такое:
Скажем, у нас 500 случайных дро из дистрибутива, который, по нашему мнению, может быть логнормальным:
X_lognorm = np.random.lognormal(mean=0.0, sigma=1.7, size=500)
Заговор
Импорт
import numpy as np
from scipy import stats
import plotly.graph_objs as go
from plotly.offline import init_notebook_mode, iplot
init_notebook_mode(connected=True)
Беги заговор
qq = stats.probplot(X_lognorm, dist='lognorm', sparams=(1))
x = np.array([qq[0][0][0],qq[0][0][-1]])
pts = go.Scatter(x=qq[0][0],
y=qq[0][1],
mode = 'markers',
showlegend=False
)
line = go.Scatter(x=x,
y=qq[1][1] + qq[1][0]*x,
showlegend=False,
mode='lines'
)
data = [pts, line]
layout = dict(xaxis = dict(zeroline = False,
linewidth = 1,
mirror = True),
yaxis = dict(zeroline = False,
linewidth = 1,
mirror = True),
)
fig = dict(data=data, layout=layout)
iplot(fig, show_link=False)