how to find standart devuation of biological data

Question

olga 2013-8-21

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/85228-how-to-find-standart-devuation-of-biological-data

I have biological data with big noise(see plot down). On this plot you can see several peaks, in my interest are two bigest of them. I want to find standart deviation of each of this peaks. Previously I used http://www.mathworks.com/matlabcentral/fileexchange/23452-interactive-peak-fitter-version-9-2 this function (Peak Fitter ) , but since I have biological data(very variable) this function doesn't work always. The algoritm which comes to my mind is next: 1)smooth and find peaks 2)cut peaks from each other manualy 3)fit them with gaussian 4)find standart deviation But I don't like to cut manualy this 2 peaks fron each other, I want to do it automaticaly but I don't know how it's possiable... May be someone knows how I can cut them automaticaly or have any ideas of alternative algoritm to find standart deviation

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Jan 2013-8-21

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/85228-how-to-find-standart-devuation-of-biological-data#answer_94713

编辑：Jan 2013-8-21

You could fit a Gaussian Mixture Model (GMM) with two Gaussians to the data using the Expectation Maximziation (EM-) Algorithm and extract the standard deviations from the result. Use findpeaks to initialize the mean values and start with small standard deviations. Applying a bandpass filter prior to this will definitely improve the result.

2 个评论
显示无隐藏无

olga 2013-8-21

Thank you for your help! But I'm begginer of using MatLab...can you wright how can I extract the standard deviations from the result?

Jan 2013-8-21

在 MATLAB Online 中打开

Take a look at http://www.mathworks.de/de/help/stats/gmdistribution.fit.html

Here you will find a very good description on how to do GMM fitting in Matlab.

In short call on your signal X

obj = gmdistribution.fit(X,k)

In your case, you would choose k=2 as you want to fit two Gaussians. The fitting results (mean and standard deviation) will be accessible in obj.mu and obj.Sigma, respectively. Look into the help file on how to set parameter initialization values.

请先登录，再进行评论。