species.py

df_Adelie = df[df['species'] == 'Adelie']
df_Gentoo = df[df['species'] == 'Gentoo']
df_Chinstrap = df[df['species'] == 'Chinstrap']

datasets = [df_Adelie,df_Gentoo,df_Chinstrap]
color = ['skyblue','red','orange']
zip_datasets_color = zip(datasets, color)
for d,c in zip_datasets_color:
    g = sns.lmplot(x = 'culmen_length_mm',
                   y = 'culmen_depth_mm', 
                   data = d, 
                   row='species',
                   scatter_kws={'color': c},
                   line_kws={'color': c})
    r = stats.pearsonr(d['culmen_length_mm'], d['culmen_depth_mm'])[0]
    ax = plt.gca()
    ax.text(.03, 1, 'r={:.2f}'.format(r),
            transform=ax.transAxes)