Python - вычисление линий тренда с ошибками

Итак, у меня есть некоторые данные, хранящиеся в виде двух списков, и я построил их с помощью

plot(datasetx, datasety)

Затем я установил линию тренда

trend = polyfit(datasetx, datasety)
trendx = []
trendy = []

for a in range(datasetx[0], (datasetx[-1]+1)):
    trendx.append(a)
    trendy.append(trend[0]*a**2 + trend[1]*a + trend[2])

plot(trendx, trendy)

Но у меня есть третий список данных , что является ошибкой в ​​исходном наборе данных. Я нормально строю диаграммы ошибок, но я не знаю, как найти ошибку в коэффициентах полиномиальной линии тренда.

Так скажем, моя линия тренда получилось 5x ^ 2 + 3x + 4 = y, должна быть какая-то ошибка для значений 5, 3 и 4.

Есть ли инструмент, использующий NumPy, который рассчитает это для меня?