A pesar del nombre de la función (lmplot = linear model plot), también es posible utilizarla para entrenar un modelo de regresión polinomial, lo que se consigue añadiendo el parámetro order y estableciéndolo al orden del polinomio deseado. Si este parámetro toma el valor 2 o superior, la función sns.lmplot ejecuta la función numpy.polyfit para estimar el modelo:
sns.lmplot(x = "total_bill", y = "tip", data = tips, hue = "time", order = 3);
