Iterate Minimize Genfit [RATIONAL]

Iterate Minimize Genfit [RATIONAL] - Iterate Minimize Genfit [RATIONAL] - Messages

#1 Posted: 2/20/2017 5:33:28 PM
Jean Giraud

Jean Giraud

983 likes in 6866 posts.

Group: User

... another good example, observe Pearson correlation.

Jean

Genfit Rational INITIALIZE.sm (63 KiB) downloaded 68 time(s).
1 users liked this post
djlischer 3/19/2017 2:53:00 PM
#2 Posted: 3/19/2017 12:31:15 PM
djlischer

djlischer

0 likes in 2 posts.

Group: User

Jean,

Thank you for your file on Genfit curving fitting with a rational function. I'm new to SMath and I learned a lot from this file! I've made a few small changes in this new Rev2 version that you may or may not find useful.

The biggest change was that I used the methodology mentioned here to initialize all the fit terms rather than just the numerator terms. If there are 'k' terms in the fit equation, you simply enter 'k' data points that you want the initization to fit exactly. The program then solves linearly for the required fit parameters. This seems to work pretty well, but I found the results can be sensitive to the initially fitted data points.

Also, I changed how the rational equation is entered. Now you can just enter they desired polynomial orders for the numerator and the denominator and the program develops the equation itself. This makes it easier to experiment with different fits since you don't have to retype the equation each time. I also added a function for Pearson's correlation coefficient to make it easier to look at the correlation of various fits.

Regards,
Jeff

Genfit Rational INITIALIZE_Rev2.sm (91 KiB) downloaded 52 time(s).
#3 Posted: 3/20/2017 10:05:02 PM
Jean Giraud

Jean Giraud

983 likes in 6866 posts.

Group: User

Jeff,

Thanks for your attentive input. I'm familiar to setup the model num/den.
Your proposal is what we have developed in Mathcad "GenfitMatrix".
In this application, your initials is just great, further refined "Minimize".
Differently than the PT100 [Platinum Resistance Thermomether], I found no
relationship between the TC 'J' data set and an eventual XFR [Transfer Function].
Thus: found no model function. True: initializing is capital and there is no rule.
Fitting data set is an art [F.B. Hildebrand].

Oh ! NIST muffed Han, Kirby. These two have a model function [attached].
Are you familiar with Mathcad ? My Mathcad 11 still working.
Mathcad 11 "GenfitMatrix" ranks first in my tool box.
If you are interested: can mail to you in case Smath forum couldn't take.

Thanks, most appreciated.

Genfit Kirby_2 [Transmute].sm (33 KiB) downloaded 44 time(s).

Genfit Kirby_2 [Reduced].sm (37 KiB) downloaded 42 time(s).

Genfit Rational INITIALIZE_Rev3 [Jeff].sm (90 KiB) downloaded 50 time(s).

Genfit Hahn.sm (34 KiB) downloaded 44 time(s).



#4 Posted: 3/21/2017 11:37:43 PM
djlischer

djlischer

0 likes in 2 posts.

Group: User

Hi Jean,
Thanks for the files. I like what you did with the initialization in the Rational Rev 3 file. I'll take a look at the other files too. Yes, I do know Mathcad. I use Mathcad 15 at work, but I don't have it at home. As I was picking up SMath, I thought about what functions I was missing and Given/Find and Minerr/Minimize were two that came to mind. That's how I stumbled on your Minimize samples.
Regards, Jeff
  • New Posts New Posts
  • No New Posts No New Posts