fitting a histogram with Gaussian and getting best fit data for mean and st dev from the fitted curve

18 views (last 30 days)
I'm trying to obtain the mean (mu) and stand dev (sigma) for a Gaussian curve drawn to fit the histogram of a data set (see attached, "histogram sample.xlsx). I can generate a histogram with Guassian curve using, for instance,
>> pd = fitdist(x,'Normal')
pd =
NormalDistribution
Normal distribution
mu = 0.0266749 [0.0162315, 0.0371183]
sigma = 0.0945052 [0.0876766, 0.102496]
However, those values (mu, sigma) appear to be the values for the underlying complete data set, NOT for the Gaussian fit. The mean, stand deviation for the Gaussian fit are smaller, of course. Question: how to extract the mu and sigma for the fitted curve instead of the raw data.
For example, Graphpad Prism (a convenient statistics software) also gives mu = 0.0266749 and sigma = 0.0945052 for the same dataset:
Graphpad also fits a Gaussian and provies a best fit values for mu and sigma for the fitted curve, mu = 0.006117 and sigma = 0.04694:
I'd like to be able to extract these same data (best fit values) using MatLab but I don't know how to. Any help?
  7 Comments
Stephen
Stephen on 7 Jan 2019
Edited: Stephen on 7 Jan 2019
Rik, I can get the relevant toolbox if I don't already have it. But is there a MatLab function or script you can suggest rather than try to plug in something from Wikipedia?
Rik
Rik on 7 Jan 2019
If you have x-y data, you can fit it to the Gaussian function:
fun=@(x,scale,mu,sd) scale/sqrt(2*pi*sd^2)*exp(-(x-mu).^2/(2*sd^2));

Sign in to comment.

Answers (0)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!