Note
Go to the end to download the full example code.
Model - savemodelresultΒΆ
[[Model]]
Model(gaussian)
[[Fit Statistics]]
# fitting method = leastsq
# function evals = 33
# data points = 101
# variables = 3
chi-square = 3.40883599
reduced chi-square = 0.03478404
Akaike info crit = -336.263713
Bayesian info crit = -328.418352
R-squared = 0.98533348
[[Variables]]
amplitude: 8.88021907 +/- 0.11359530 (1.28%) (init = 5)
center: 5.65866105 +/- 0.01030493 (0.18%) (init = 5)
sigma: 0.69765480 +/- 0.01030508 (1.48%) (init = 1)
fwhm: 1.64285148 +/- 0.02426660 (1.48%) == '2.3548200*sigma'
height: 5.07800563 +/- 0.06495769 (1.28%) == '0.3989423*amplitude/max(1e-15, sigma)'
[[Correlations]] (unreported correlations are < 0.100)
C(amplitude, sigma) = +0.5774
# <examples/doc_model_savemodelresult.py>
import numpy as np
from lmfit.model import save_modelresult
from lmfit.models import GaussianModel
data = np.loadtxt('model1d_gauss.dat')
x = data[:, 0]
y = data[:, 1]
gmodel = GaussianModel()
result = gmodel.fit(y, x=x, amplitude=5, center=5, sigma=1)
save_modelresult(result, 'gauss_modelresult.sav')
print(result.fit_report())
# <end examples/doc_model_savemodelresult.py>
Total running time of the script: (0 minutes 0.060 seconds)