Minuit

Google Colaboratory
https://colab.research.google.com/drive/1jufurZXLxzkC1FzyOXpq6I7_-oDQZvCt

import numpy as np
from iminuit import Minuit

dx = np.arange(10)
dy = ((50.0 * dx) + 100.0) * (1.0 + (0.05 * np.random.randn(len(dx))))
dy_error = np.sqrt(dy)
model = lambda x, a, b : ((a * x) + b)
chisq = lambda a, b : sum(((dy - model(dx, a, b)) ** 2.00) / (dy_error ** 2.00))

fit = Minuit(chisq, a=1.00, b=1.00)
fit.migrad()
fit.hesse()
fit.minos()

fit.values
## <ValueView a=51.49166936856244 b=97.93756640441136>

fit.draw_mncontour('a', 'b', cl=[1.00, 2.00, 3.00])

fit.draw_mnmatrix(cl=[1.00, 2.00, 3.00])

https://hackmd.io/@tenoto/Skl_wArqD

Daiphys is a professional-service company for research and development of leading-edge technologies in physics and engineering.
Get started accelerating your business through our deep expertise in R&D with AI, quantum computing, and space development; please get in touch with Daiphys today!

Daiphys Technologies LLC - https://www.daiphys.com/