Problem when plotting the figure

I have a vector called K [622x1] double I plot the PDF of K . I fitted my K vector to the Gaussian mixture distribution and ploted the PDF with a solid line :

>> GMModel = fitgmdist(k,2);
>> gm=gmdistribution(GMModel.mu,GMModel.Sigma);
>> pdfk=pdf(gm,k);
>> figure
>> plot(k,pdfk)

The plot has several lines and I dont know what is the problem and how to resolve it. Ca anyone help me with this? thank youuu

Answer 1

[Star Strider]

I do not have ‘k’, so I created my own, and was able to reproduce essentially the result you got.

The multiple lines disappear of you sort ‘k’ first:

k = sum([randn(622,1)+5; randn(622,1)+1],2);
k = sort(k);                                                % <— ADD ‘sort’ CALL
GMModel = fitgmdist(k,2);
gm=gmdistribution(GMModel.mu,GMModel.Sigma);
pdfk=pdf(gm,k);
figure
plot(k,pdfk)

.

[yusra Ch]

Thank you so much ! it works

[Star Strider]

As always, my pleasure!

[yusra Ch]

I have one more question, How to draw the PDF of the vector K ?

I have plotted the PDF fitted to the mixture Gaussian distribution and I want to compare the two plots.

[Star Strider]

I would do something like this (starting with my original code):

k = sum([randn(622,1)+5; 0.5*randn(622,1)+1],2);
k = sort(k);                                                % <— ADD ‘sort’ CALL
GMModel = fitgmdist(k,2);
gm=gmdistribution(GMModel.mu,GMModel.Sigma);
pdfk=pdf(gm,k);
[f,x] = ecdf(k);                                            % Empirical Cumulative Distribution Function
[fr,xr] = resample(f, x);                                   % Resample To Constant Sampling Interval
df = gradient(fr)./gradient(xr);                            % Calculate Derivative To Get PDF
dfs = smoothdata(df, 'gaussian', 150);                      % Smooth To Eliminate Noise

figure
plot(k,pdfk)
hold on
plot(xr, dfs)
hold off

This appears to produce a reasonable approximation. Experiment with the smoothdata function (introduced in R2017a) with your data to get different results. (I do not kinow what MATLAB version you are using. Other options, such as movmedian to do the smoothing were introduced in R2016a, and of course there are still other options.)

.