【Python】線形フィットをする方法

import numpy as np
import matplotlib.pyplot as plt
a, b = np.polyfit(x, f, 1)
fx = np.array(range(xmax+1))
ff = a * fx + b
plt.scatter(x, f)
plt.plot(fx, ff)
plt.show()

【参考】
【NumPy】回帰分析で直線近似（線形フィッティング） | アルゴリズム雑記