Я попробовал ggplot и несколько других библиотек, но ни одна из них не смогла сделать то же самое.
from scipy import stats
from scipy.stats import norm, skew, kurtosis
def stat_kde_plot(input1,input2,input3):
f, (ax1,ax2,ax3) = plt.subplots(1,3,figsize=(15,5))
sns.kdeplot(df_protein[input1],ax = ax1,color='blue',shade=True,
label=("Skewness : %.2f"%(df_protein[input1].skew()),
"Kurtosis: %.2f"%(df_protein[input1].kurtosis())))