[lmfit] 9. 対数正規分布関数に基づくモデルによるカーブフィッティング

はじめに

lmfitは非線形最小二乗法を用いてカーブフィットするためのライブラリであり、Scipy.optimize.curve_fitの拡張版に位置する。ここでは、lmfitで対数正規分布をカーブフィッティングする方法について説明する。

rng = default_rng()とし、rng.lognormalで対数正規分布に従うデータを作成し、np.histogramでヒストグラム化することで対数正規分布を得る。図で示すと以下のようになる。

lmfit.models の LognormalModelをモデル関数として用いる。

model.guess(histo, x=bins)により上図の対数正規分布のフィッティングパラメータの初期値を推定する。パラメータ(params)は以下のようになる。

model.fit(histo, params, x=bins)により、カーブフィッティングを実行する。

print(result.fit_report())により、フィッティングの結果を見ることができる。

result.plot_fit()によりデータとフィッティングカーブが表示される。

result.plot()とすることで残差とともにフィッティング結果が表示される。