Bivariate Guassian mixed model on 1D data

1 次查看（过去 30 天）

Trisha Kibaya 2018-8-24

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/416125-bivariate-guassian-mixed-model-on-1d-data

评论： Trisha Kibaya 2018-8-28

Hi, I am trying to plot a guassian curve on the histogram distribution of my 1D data. This is my code:

data = xlsread('csf.xlsx','ab'); 
     mu1 = mean('data');
     sigma1 = std('data'); 
     X = [mvnrnd(mu1,sigma1,1000)];
     numComponents=2;
paramEsts = fitgmdist(X,numComponents);  
MU=[paramEsts.mu(1);paramEsts.mu(2)];
SIGMA=cat(3,[paramEsts.Sigma(1)],[paramEsts.Sigma(2)]);
PPp=[paramEsts.PComponents(1),paramEsts.PComponents(2)];
gm = gmdistribution(MU,SIGMA,PPp); 
dimension = gm.NumVariables 
 rng(1); 
 min_val = min(data(:)-100); max_val = max (data(:)+1000);  
   % min_val = min('data'); max_val = max ('data');  
      histogram(data, 'Normalization', 'pdf', 'EdgeColor', 'none')
     x= min_val:0.01:max_val;
  xlim([min_val max_val])
  hold on; 
  xgrid = linspace (50,3000,500)'; 
  n1 = makedist('normal',paramEsts.mu(1),sqrt(paramEsts.Sigma(1)));
  n2 = makedist('normal',paramEsts.mu(2),sqrt(paramEsts.Sigma(2)));
  p = paramEsts.ComponentProportion;
    y = p(1)*pdf(n1,xgrid) + p(2)*pdf(n2,xgrid);
    hold on; plot(xgrid,y,'b'); hold off

And this is what I get:

However, I am trying to get something like this:

Please, how do I achieve this?

2 个评论
显示无隐藏无

the cyclist 2018-8-24

Can you upload either the Excel file or the resulting data variable in a MAT file (using that little paperclip icon)?

Trisha Kibaya 2018-8-27

data.mat

Sure. Thanks.

请先登录，再进行评论。

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

回答（1 个）

Image Analyst 2018-8-24

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/416125-bivariate-guassian-mixed-model-on-1d-data#answer_334034

fit_two_Gaussians.m

Try something like the attached demo that fits two Gaussians to data. Adapt as needed (like replace blue dots on plot with white bars or whatever).

7 个评论
显示 5更早的评论隐藏 5更早的评论

Image Analyst 2018-8-28

Make it so people can read it: http://www.mathworks.com/matlabcentral/answers/13205#answer_18099

Trisha Kibaya 2018-8-28

在 MATLAB Online 中打开

Thanks. I am still new at this, I had no idea that I posted my code in the wrong way. Sorry about that.

    xgrid = linspace (40,1000,112110'); 
    n1 = makedist('normal',paramEsts.mu(1),sqrt(paramEsts.Sigma(1)));
    n2 = makedist('normal',paramEsts.mu(2),sqrt(paramEsts.Sigma(2)));
    p = paramEsts.ComponentProportion;
    y1 = p(1)*pdf(n1,xgrid); hold on; plot (xgrid, y1,'b'); hold off
    y2 = p(2)*pdf(n2,xgrid); hold on; plot (xgrid, y2,'r'); hold off