Skip to main content

fitting - Using Manipulate to find best starting parameters to fit data



Often, curve fitting is very sensitive to starting values of parameters. It would be great, if I could to find such starting values of parameters "by hand" using Manipulate.


So, I plot in Manipulate my experimental data points and theoretical curve. I can change the shape of this curve by changing all parameters in manipulate in such way that they approximate my experimental data relatively well. Now I would like to run fitting using current parameters in Manipulate as starting values. Finally I would like to insert the parameters found by fitting back into Manipulate.


Here is example code for simple function. My data is more complex.


data = Table[{x, 8 x^3 - 7 x^2 - 10 x + 1 + RandomReal[{-5, 5}]}, {x, -2, 2, 0.1}];
Manipulate[
 Show[
  Plot[a x^3 + b x^2 + c x + d, {x, -2, 2}],
  ListPlot[data]
 ],
 {a, -10, 10},

 {b, -10, 10},
 {c, -10, 10},
 {d, -10, 10}
]

Answer



As I understand the question a curve fitting procedure that has the following properties is sought:



  1. Manually adjust the parameters to get an approximate fit.

  2. Use these parameters as the starting values for FindFit.

  3. Propagate the solution from FindFit back to the Manipulate parameters.


  4. Subsequently enable further editing of the Manipulate parameters and repeat the cycle.


The following code satisfies this criteria by wrapping Manipulate inside a DynamicModule and the use of a Button to indicate when FindFit should be run.


data = Table[{x, 8 x^3 - 7 x^2 - 10 x + 1 + RandomReal[{-5, 5}]}, {x, -2, 2, 0.1}];

DynamicModule[
 {
  sol
  },


 Manipulate[

  If[computeFlag == True,
   sol = FindFit[data, 
     aa x^3 + bb x^2 + cc x + 
      dd, {{aa, a}, {bb, b}, {cc, c}, {dd, d}}, x];
   {a, b, c, d} = {aa, bb, cc, dd} /. sol;
   computeFlag = False;
   ];


  Column[{
    Dynamic[
     Button["Compute",
      computeFlag = True
      ]
     ],
    Show[
     Plot[a x^3 + b x^2 + c x + d, {x, -2, 2},
      PlotStyle -> Black],
     ListPlot[data, PlotStyle -> Red],

     ImageSize -> 300,
     PlotRange -> {{-2.05, 2.05}, All}
     ]
    }],

  (*Manipulate variables*)
  {{computeFlag, False}, ControlType -> None},

  {{a, 1}, -10, 10, Appearance -> "Open"},
  {{b, 1}, -10, 10, Appearance -> "Open"},

  {{c, 1}, -10, 10, Appearance -> "Open"},
  {{d, 1}, -10, 10, Appearance -> "Open"}

  ] (* end of Manipulate *)
 ] (* end of DynamicModule *)

Below is a figure where the parameters have been manually adjusted.


Mathematica graphics


After clicking the Compute button FindFit propagates the solution back to the Manipulate parameters.


Mathematica graphics



The user is free to re-edit the Manipulate parameters and repeat the cycle.


Comments

Post a Comment

Popular posts from this blog

functions - Get leading series expansion term?

Given a function f[x] , I would like to have a function leadingSeries that returns just the leading term in the series around x=0 . For example: leadingSeries[(1/x + 2)/(4 + 1/x^2 + x)] x and leadingSeries[(1/x + 2 + (1 - 1/x^3)/4)/(4 + x)] -(1/(16 x^3)) Is there such a function in Mathematica? Or maybe one can implement it efficiently? EDIT I finally went with the following implementation, based on Carl Woll 's answer: lds[ex_,x_]:=( (ex/.x->(x+O[x]^2))/.SeriesData[U_,Z_,L_List,Mi_,Ma_,De_]:>SeriesData[U,Z,{L[[1]]},Mi,Mi+1,De]//Quiet//Normal) The advantage is, that this one also properly works with functions whose leading term is a constant: lds[Exp[x],x] 1 Answer Update 1 Updated to eliminate SeriesData and to not return additional terms Perhaps you could use: leadingSeries[expr_, x_] := Normal[expr /. x->(x+O[x]^2) /. a_List :> Take[a, 1]] Then for your examples: leadingSeries[(1/x + 2)/(4 + 1/x^2 + x), x] leadingSeries[Exp[x], x] leadingSeries[(1/x + 2 + (1 - 1/x...
Post a Comment
Read more

How to thread a list

I have data in format data = {{a1, a2}, {b1, b2}, {c1, c2}, {d1, d2}} Tableform: I want to thread it to : tdata = {{{a1, b1}, {a2, b2}}, {{a1, c1}, {a2, c2}}, {{a1, d1}, {a2, d2}}} Tableform: And I would like to do better then pseudofunction[n_] := Transpose[{data2[[1]], data2[[n]]}]; SetAttributes[pseudofunction, Listable]; Range[2, 4] // pseudofunction Here is my benchmark data, where data3 is normal sample of real data. data3 = Drop[ExcelWorkBook[[Column1 ;; Column4]], None, 1]; data2 = {a #, b #, c #, d #} & /@ Range[1, 10^5]; data = RandomReal[{0, 1}, {10^6, 4}]; Here is my benchmark code kptnw[list_] := Transpose[{Table[First@#, {Length@# - 1}], Rest@#}, {3, 1, 2}] &@list kptnw2[list_] := Transpose[{ConstantArray[First@#, Length@# - 1], Rest@#}, {3, 1, 2}] &@list OleksandrR[list_] := Flatten[Outer[List, List@First[list], Rest[list], 1], {{2}, {1, 4}}] paradox2[list_] := Partition[Riffle[list[[1]], #], 2] & /@ Drop[list, 1] RM[list_] := FoldList[Transpose[{First@li...
Post a Comment
Read more

front end - keyboard shortcut to invoke Insert new matrix

I frequently need to type in some matrices, and the menu command Insert > Table/Matrix > New... allows matrices with lines drawn between columns and rows, which is very helpful. I would like to make a keyboard shortcut for it, but cannot find the relevant frontend token command (4209405) for it. Since the FullForm[] and InputForm[] of matrices with lines drawn between rows and columns is the same as those without lines, it's hard to do this via 3rd party system-wide text expanders (e.g. autohotkey or atext on mac). How does one assign a keyboard shortcut for the menu item Insert > Table/Matrix > New... , preferably using only mathematica? Thanks! Answer In the MenuSetup.tr (for linux located in the $InstallationDirectory/SystemFiles/FrontEnd/TextResources/X/ directory), I changed the line MenuItem["&New...", "CreateGridBoxDialog"] to read MenuItem["&New...", "CreateGridBoxDialog", MenuKey["m", Modifiers-...
Post a Comment
Read more