Re: Gaussian fit

*To*: mathgroup at smc.vnet.net*Subject*: [mg23208] Re: [mg23148] Gaussian fit*From*: BobHanlon at aol.com*Date*: Mon, 24 Apr 2000 01:12:15 -0400 (EDT)*Sender*: owner-wri-mathgroup at wolfram.com

Needs["Statistics`NormalDistribution`"] Needs["Graphics`Graphics`"] freq = {{-0.95, 0}, {-0.85, 0}, {-0.75, 0}, {-0.65, 1}, {-0.55, 2}, {-0.45, 9}, {-0.35, 26}, {-0.25, 36}, {-0.15, 23}, {-0.05, 15}, {0.05, 8}, {0.15, 7}, {0.25, 7}, {0.35, 1}, {0.45, 0}, {0.55, 2}, {0.65, 1}, {0.75, 0}, {0.85, 0}, {0.95, 0}}; step = freq[[2, 1]] - freq[[1, 1]]; data = Flatten[Table[#[[1]], {#[[2]]}] & /@ freq]; Estimates for the mean and standard deviation are {m, s} = {Mean[data], StandardDeviation[data]}; The maximum likelihood estimate (MLE) equations are eqns = ((D[ExpandAll[Times @@ (PDF[NormalDistribution[mu, sigma], #] & /@ data)], #] == 0) & /@ {mu, sigma}); The MLE for the distribution is dist = NormalDistribution[mu, sigma] /. FindRoot[eqns, {mu, m}, {sigma, s}]; bc = GeneralizedBarChart[Transpose[ {Table[-1. + step(k - 1/2), {k, Length[freq]}], Transpose[freq][[2]]/(step*Length[data]), Table[step, {Length[freq]}]}], Frame -> True, Axes -> False, PlotRange -> All, DisplayFunction -> Identity]; plt = Plot[PDF[dist, x], {x, -1, 1}, Frame -> True, Axes -> False, PlotStyle -> {{RGBColor[0, 0, 1], AbsoluteThickness[2]}}, DisplayFunction -> Identity]; Show[{bc, plt}, DisplayFunction -> $DisplayFunction, ImageSize -> {500, 310}]; The same approach would be used for your refined data. However, the order of your second set of frequency data is the reverse of the order for your first set Bob Hanlon In a message dated 4/20/2000 4:29:35 AM, slipstk at olemiss.edu writes: >I am trying to determine the best Gaussian fit to a series of data >points. > >Data={{-0.95,0},{-0.85,0},{-0.75,0},{-0.65,1},{-0.55,2},{-0.45,9},{-0.3 >5,26},{-0.25,36},{-0.15,23},{-0.05,15},{0.05,8},{0.15,7},{0.25,7},{0.35,1 >},{0.45,0},{0.55,2},{0.65,1},{0.75,0},{0.85,0},{0.95,0}} > >A more detailed histogram of this looks like > >{{0,-0.98},{0,-0.94},{0,-0.9},{0,-0.86},{0,-0.82},{0,-0.78},{0,-0.74},{0, >-0.7},{0,-0.66},{1,-0.62},{0,-0.58},{1,-0.54},{2,-0.5},{5,-0.46},{3,-0.42 >},{10,-0.38},{14,-0.34},{10,-0.3},{13,-0.26},{15,-0.22},{10,-0.18},{9,-0. >14},{7,-0.1},{3,-0.06},{9,-0.02},{3,0.02},{5,0.06},{1,0.1},{4,0.14},{2,0. >18},{3,0.22},{4,0.26},{0,0.3},{0,0.34},{1,0.38},{0,0.42},{0,0.46},{0,0.5} >,{1,0.54},{1,0.58},{0,0.62},{1,0.66},{0,0.7},{0,0.74},{0,0.78},{0,0.82},{ >0,0.86},{0,0.9},{0,0.94},{0,0.98}} > >I've tried using NonLinearFit, but I can't get the parameters to make >sense. There is also the possibility that this may be indicative of two > >processes superimposed upon each other, in which case, the fit would be > >two Gaussian distributions added together. I've also looked through the > >Archives, but there doesn't seem to be much on this that is directly >useful. Either that, or the solution is so simple that I'm just not >seeing it. >