Curve Fitting in Python (2022)

[@sergioaltamirano5888]

There is a lot of content through the internet, particularly on KZbin, about 'data analysis' and 'machine learning' topics. But the content like this is not abundant. I have read many books too, and I kept struggling to find high quality content. Congratulations! I'm a new subscriber. Sending my best from Argentina.

[@andrealorenzetto4947]

With these videos my coding journey is just simpler and much more enjoyable, I am very grateful. Huge support from Italy

[@hakansert13]

Sometimes one may both know coding and physics but it becomes really exhausting to combine them together; you may not able to know where to start... You are directly attacking to that tricky intersection perfectly and that is exactly what is needed for lots of science students/academics! I just personally wanted to thank you for such great videos as a PhD student in physics who had worked with paper and pen for years but needed numerical computation eventually. Your videos are awesome and very informing! I hope you keep going.

[@skilz8098]

If you ever try your hand at 3D Graphics Programming and write the underlying core components of the 3D Rendering Engine where you implement a Physics Engine with basic kinematics, collision detection systems, A.I. Path Finding Algorithms, etc. then you get into doing Animations, Particle Effects, Rigid Bodies, Rag Dolls, Liquid or Water Animations, etc... by using advanced shaders for many various lighting, shading, shadow, and other special effects... That's one of the best ways to merge the two fields. Now, this kind of Physics programming is different from Data Science, Equation Solver type Physics system, but it's a great place to start and most of the topics and transitions would already be covered. A basic 3D Rendering Engine or 3D Game Engine uses mostly linear algebra and trigonometry, but there is still a bit of calculus here and there such as when basic integrals, performing interpolations of various kinds, using runge kutta as well as FFTs all come in handy. Especially FFTs when working with Audio Processing and 3D Sound Environments is another great benefit. You end up incorporating many mathematical and physics, (maybe even a little bit of chemistry, biology, and environmental systems depending on the context of your Engine and what you are going to build with it) concepts into your Engine that stem across many software engineering and program application development disciplines by using almost every kind of container and associated algorithm there are out there. Other key features would be file management and file parsers, you can end up writing your own scripting or run time interpreted language, if doing multiplayer over networking then multithreading and multiprocessing principles, compression and decompression algorithms, encryption and decryption algorithms come in handy. Then you have many other topics such as batch processing, Storage or Assets Memory Managers, Many Various Shaders with an accompanying Shader Manager, and so much more. This doesn't even include things such as GUIs, HUDs and Font Management... Now, once you have the Engine Built with all of its core components (and I don't mean using a 3D engine such as Unity or Unreal...) Then the fun begins when you start building your game that is seperate from then engine but requires it as a core underlying library to function properly. Now imagine building a game like Factorio or Minecraft... Both of those games are Turing Complete... or even a 2D with 3D layer slices such as in Oxygen Not Included or a semi-factory but defense tower game like Mindustry which are both Turing Complete... or something similar to those... Then I think you can build just about any kind of software application that you could imagine except maybe an actual robust compiler-debugger-linker suite, or an actual Operating System... Oh they could be done... but they're a challenge. I found it was actually a bit easier to learn how to build a Hardware Emulator in C++ with a focus on the 6502 for the NES... Now that was a fun project as well! So almost anything is possible to build but it does depend on the limitations of your hardware... This is just a summary of about a 20 year journey of independent studying purely self taught hobbyist and some of the many topics I've learned over the years! There's many other things that I didn't mention as this is starting to become quite long... It's been a great journey. Even watching Ben Eaters 8-Bit Breadboard CPU series was a delight! So much so, that I actually implemented it in Logisim.

[@mikhailbandurist8652]

That's the scientific channel that we need. Good sound (which is the most important), good picture, Chad-looking science guy (in movies smart people are often nerds with big glasses) and a good explanation - this is really nice! Thank you for video!

[@python4scientists]

Nice! Note that (in 19:40) you need to set absolute_sigma=True to get the absolute value errors on the parameters. A good example is using the known formulas to calculate the error on the parameters (A,B) of a strait line fit and compare with the results returning with absolute_sigma=False e absolute_sigma=True.

[@boseongcho62]

It's a good comment

[@crazyclown9914]

Two weeks ago I had a physics assignment involving this method, which was a absolute disaster. This is why I should watch your videos as soon as they show up. Anyway, Thanks for the amazing job!

[@moona5454]

Great video. Especially focusing and repeating multiple times that our goal is the parameters and not just having a good fit! Parameters = Context = Reason why we're doing the fit

[@otaviooliveira5157]

The kind of detailed explanation i needed! The examples nailed it! Thanks a lot :)

[@vadympasko100]

Absolutely love this intro to SciPy's 'curve_fit' function! Suggest you to make a video on spline fitting, which is included into SciPy's 'interpolate' module. With splines it is possible to find best fit to a really complex data, including those presented in this tutorial.

[@heera_ai]

Your initial guess was good enough 😂, great tutorial!

[@nathan6798]

Hi, first off EXCELLENT VIDEO THANK YOU SO MUCH Just for those who need this video, there are a couple of clarifications (sorry if you made them and I didnt notice), but its popt, pcov = scipy.optimize.curve_fit not just curve_fit () . Unless you made a shortcut beforehand Other than that , thank you again so much help !

[@rainycloud10]

This is exactly the tutorial I needed given I had to perform many counting experiments . I wish I had found this channel sooner before going insane over tutorials that made zero sense to me.

[@alejandraacosta248]

Wow! You explain very clearly! I was fighting with my code since I did not understand how to work with curve_fit and you fell out of the sky! Thank you so much!

[@j.abrahamhernandez3629]

dude your videos rock, i love the no-nonsense type of attitude when tackling things, cheers! (;

[@atmosagm]

This is excellent content. Clearly explained a problem I've been having trouble grasping for months. Thanks!

[@RutgerHaan]

Thanks a lot! This was a really thorough explanation which gave me everything I need to plot my data in a nice way :) Keep it up!

[@MrAman47]

Great video, just a little correction, as the physicist in me is happy to show himself every now and then: It's actually Lennard-Jones potential, sometimes also called Van der Waals potentials (dipole-dipole potential) big brother. I don't think there's a thing called Leonard-Weibeck potential reguarding atomic repulsion, but if there is and it's related to the LJ potential, let me know!

[@MrPSolver]

Oh I might have misspoke! Thank you for this correction 😂

[@swatibhargava4140]

Super helpful and so nicely explained! Thank you so much for making this video! The multiple examples helped me a lot and in fact covered what I was as looking for!

[@laurenceturpin1409]

Thank you for making this video I particularly liked that you gave more than one example.

[@loredo97]

Brilliant explanation of the covariance!

[@gloryths]

You give quality out there. Respect man.

[@yt-1161]

where is the csv file ? @13:47

[@siddharthpachpor5928]

superb video with respect to physics and maths! Best application based explanation. Period. Please I request you to post such concepts with such physics based applications visualisation. I just feel blessed. Thanks a lot. tons

[@chayanikarabha8556]

Thank you so much for making such videos. Specifically in physics. It is of so much help. Please keep doing such videos.

[@dylanskinner6815]

Your videos are sweet! I found you earlier today and have already learned a ton. Thanks!

[@boseongcho62]

Thank you now I understood how the errors behave!

[@abrahamdelacruz8389]

Dude! this video is so useful! Great explanation and great material!

[@nothingisreal6345]

Like 25 years ago i wrote code to fit multiple gauss function to some spectrum comming from Cherenkov radiation emitted from cristals. Those days everythibg had to be done in C/C++. I use the "Numerical recipies in C" book and lib and the Marquard Levenberg algorithm. Took 3 weeks to progam that (including some GUI and other stuff). Today with Scipy maybe one or two days... Some stuff actually get's better.

[@Bored_Trumpet]

Was expecting a good description of what's going on with the curve-fitting functions, albeit with boring data. But actually I found exactly what I was looking for. Trying to get the photopeaks and their resolutions in a similar spectroscopy. Also the walking through cell-by-cell helps for digestion. Used to begrudingly use python, but after a semester I much more prefer it over finicky excel sheets.

[@iurialmeida8979]

i found your videos about 10 minutes ago and mate, what an amazing content, you remind me of keith galli, he is awesome too

[@GabrieleIuzzolino]

Thank you! Your tutorials are very clear and useful!

[@DEChacker]

awesome video as always. Learned a lot - especially the cov-Matrix was super helpful One note: It is the Lennard-Jones potential, not the Leonard Weibeck potential

[@harveyjones5227]

Amazing! Thank you. This has been perfect help with plotting an exponential curve to my data.

[@MrSpaceboyy]

This is great! Glad to find your channel

[@ayushisuman4889]

Very helpful content. The error(as obtained from experiment) is sometimes not symmetric about the mean y position. Is it possible to encode that into curve_fit function?

[@JarredDavidson]

Excellent video. Very clear unpacking of the logic. Thank you! I've subbed.

[@skilz8098]

Have you done any videos on FFTs or any Runge Kutta methods?

[@ragibabsarramon7881]

Your Videos are really helping me. Thanks a lot.

[@antonXPS]

excellent explanation, very detailed

[@GabrieleIuzzolino]

I'd like to ask you a question about the standardisation: why didn't you just multiply the ordinates by 10^(-21), in order to have numbers in the order of mangitude of 1?

[@datastako156]

great videos! continue your greatwork bro! supporting from philippines here.

[@NaneRulz]

Great video! In the given case that you want to work multiple samples, and perform simultaneous fits. Is that possible? It is something I use often in XPS.

[@ender84]

Same question here. I have tried clustering without luck.

[@mis_llaneous]

I learnt a lot from thus vedio, thanks bro👍

[@chaitanyavarmamudunuri7270]

Thank you for the content, I would like to see a video on hysteresis curvefitting with multiple conditions

[@brucewernick6542]

The heatmap of the covariance is brilliant, I'm definitely going to experiment with that. You mentioned the importance of the diagonal, but what do the other values in the covariance matrix mean?

[@librealgerien]

That’s how one parameter affects another if at all.

[@chaitanyavarmamudunuri7270]

Thank you very much for such beautiful tutorials. Can you show fitting hysteresis models

[@danteng5651]

This is super clear!

[@Data_Handler-e6c]

How to come up with the equation for curve fit method? Suppose the data doesn't follow a specific distribution which can be determined by seeing the scatter plot.

[@andresderudder9950]

Great video ! I just finished a lab report last Saturday, this would been useful jajaja. I would like you to make a video about animations in python tool

[@rainycloud10]

Same. I wish I had found this tutorial a couple of weeks ago for my lab reports. But better late than never.

[@Heisenberg20023]

I had a question, what if the function had constants in its definition that were an array that take certain values at specific points of xdata?

[@Data_Handler-e6c]

How to come up with the equation for curve fit? Suppose the scatter plot is not showing a particular distribution.

[@raysu9853]

Great video, thank you. Hope to see you @ UBC some day :)

[@senolkurt7864]

Thanks for the great tutorial. Since in real life we have only the data, how can we find the best non-linear model function that fits to our data?

[@ncarl2335]

you gained a subscriber sir thanks

[@ahmedalshemi955]

Thanks a lot for this nice explanation! I wonder what you use to fix one or some of the fitted parameters using the Scipy curve fitting function?

[@andrewjolly319]

Please consider doing some stuff on Bayesian inference using Markov Chains, maybe using the PyMC library?

[@G.K.LAKSHMIChellu]

Hello sir How to put scientific notation in p0 values. Eg. 2.7778×10^-15

[@monenehmoneneh9561]

perfect please do more in physcis and more intuitive

[@preshitlimje6046]

Hello Mr. P Solver, First of all a great video on optimization. I have a question regarding the optimization and identification of values in a given curve. Since I am completely new to this, I am not much aware of reverse engineering. I have an experimental Time-Temperature curve of a heat transfer system and based on it I need to identify the convective heat transfer coefficient. Is it possible to make it? if yes could you please give me an example. Thanks in advance. With Regards, Preshit

[@frederic-louissauser6945]

pretty nice explanation, thanks! would it give some sense to use R or Scilab to achieve this, would it be more efficient?

[@ssj_hamood5364]

For the Gaussian fit again. If I add a +c factor at the end in the gauss_f definition, will the plot work if I expand popt with c_opt so that the measurement points are taken into account which do not look "parabolic"?

[@KunalSingh-my5nd]

My values in y axis are all less than 1. So taking the sq root considering Poison statistics of photons gives very high error. What should I do ? Should i take some kind of scaling factor?

[@odvutmanush3234]

could you please show how to add bound to parameter?

[@vinrai09]

Thanks for explaining this so well!! Super helpful!!

[@gerrievanstaden3416]

Is there a way to do R^2?

[@joaovictorbalieirodasilva1430]

Please, make a numba tutorial :) Your videos are awsome

[@carlosyordanocozdelacruz3370]

good video. How would I estimate the data based on those parameters?

[@yazito21]

Very helpful video...👍👍👍

[@dominikrodak4204]

Hello, Thanks a lot for yours tutorials! They help me a lot. I have a question. How to use curve_fit (or how to cope with the problem) if I have a system of ODEs. How to use it as a model for curve_Fit function?

[@rostamr4096]

This is so very helpful...thank you

[@sadaf6295]

Sir, can u make video on Heisenberg spin chain and correlation between spins.

[@kapfox]

You definitely making a good thing, killing a Matlab!

[@lookaway8496]

Say I have a histogram plot with bars. I would like to do a skew norm fit to it. How do I proceed with this? can anyone help? like what initial data should I have

[@03_bikramkesharipanda68]

Can we solve a multi-variable function through curve-fitting , for example - y = f(x1,x2,x3,x4)?

[@AJ-et3vf]

Great video! Thank you!

[@rinezaman492]

These videos are really great

[@bavneetkaur979]

Hello Mr. P solver, is there a way to solve for a complex function- meaning that the function that I have defrined returns complex(both real and imaginary parts) simultaneously, I tried doing it this with Scipy.optimize.curve_fit() but does not help If you have any idea pls share, thanks :)

[@lospaturno]

Can you tell something about fitting multiple gaussians in a signal?

[@philtoa334]

Very good ,very Useful very : )

[@sarasachandrikabhavanivajj6592]

I am getting the following error. I tried multiple methods like using numpy.asarray, np.array...... Is there any other solution? setting an array element with a sequence. The requested array has an inhomogeneous shape after 1 dimensions. The detected shape was (5,) + inhomogeneous part.

[@khairulfahim]

What is that algorith? I mean that is working in the background.

[@mohamadyusuf847]

how to get that sample data, if i need

[@ricardoferes9051]

Hey could you do a tutorial for a CFD simulation in python?

[@Raiden_Amani]

Thank you so much !! You've been a really big help.

[@colesmith2136]

Really great video!

[@minhangnguyen999]

OH GOD!!, thank you very much!!!!!

[@LOL-vt8jh]

Good job!

[@Science4Ever]

Great video as always! May I ask you to recommend some good books to dive further into this stuff? Maybe the ones you learned from, if any? It would be great! Keep doing these great videos, sir!

[@jay89boy]

ye im also wondering how did he reach all that knowledge

[@Science4Ever]

@@jay89boy Same here, he might and actually should do a short video about giving some advices on how to dig deeper into the topic, especially book recommendations. He is a great guy and I think he will read this and respond in one way or another.

[@amelieeee9692]

you left out a 2 in the gaussian formula

[@assassingamer9088]

but when i plug in 2 it gives wrong answer. i tried it on my dataset. without it it gives right answer. why?

[@romenaakter1042]

Simply WOW!

[@Path_k_pradeep]

Thanks 👏👏

[@bean217]

This isn't at all related to the video content, but I have the same exact scar on my forehead

[@obliquesealray2188]

great vid

[@STWNoman]

Perfect

[@ksrajavel]

Amazing

[@rishiraj2548]

👍💯

[@miszcz310]

Hmmm i was expecting something more clicking this video. When you have object returned by the fit method it has method called plot. It shows data with fitted model as well as plots the residuals.

[@MsJ0ni]

thank you..

[@TheGnarTube]

Your style annoys me but the content is good

[@opm6278]

this scipy shit hard to understand man.