scipyで最小二乗法

python用の科学計算ライプラリscipyを使った最小二乗法用のテンプレを貼って置きます。
SciPy Cookbook — SciPy Cookbook documentation
を参考にしました。

    from scipy.optimize import *
    fitfunc = lambda p, x: p[0] * x + p[1]  #フィットさせる方程式
    errfunc = lambda p, x, y: pow((fitfunc(p, x) - y), 2) #方程式とデーター間の距離
    p0 = [0, 0]  #パラメーターの初期値
    p1, success = leastsq(errfunc, p0, args = (x, y))  #最適化するようなパラメーターを受けとる
    plot(x, y, x, fitfunc(p1, x))