I used sns.distplot, but it is a deprecated function. Is there any good alternative that have 'fit'?
sns.distplot(y, kde=False, fit=stats.norm)
Right now I'm using this
s1, loc2, scale2 = lognorm.fit(y)
r2 = lognorm.rvs(s1, loc2, scale2, size=r_size)
sns.histplot(data=r2, kde='True', ax=axes[1], alpha = 0)
sns.histplot(data=y, ax=axes[1])
axes[1].set_title('Log Normal')