fitting a histogram with Gaussian and getting best fit data for mean and st dev from the fitted curve
6 次查看(过去 30 天)
显示 更早的评论
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: ![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/200331/image.jpeg)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/200331/image.jpeg)
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:
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/200332/image.jpeg)
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 个评论
Rik
2019-1-7
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));
回答(0 个)
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!