MATLAB: Can anyone debug the following code

digital image processingimage analysisstatisticsStatistics and Machine Learning Toolbox

The code is given below. The problem I am facing is while bringing two plots together using hold on function the result is not the exact result what I am getting in the seperate window. variable K is attached. Thanks

figure,histogram(K);hold on;
pd= fitdist(K','Kernel');
x_values = -30:0.1:40;
y = pdf(pd,x_values);
plot(x_values,y,'LineWidth',2)

Best Answer

becuase two plots have different y-axis scale! you have to modify y-axis scale!

 %code
figure,hist(K);hold on;
pd= fitdist(K','Kernel');
x_values = -30:0.1:40;
y = pdf(pd,x_values);
plot(x_values,y*(max(hist(K))/max(y)),'r','LineWidth',2)

Related Solutions

MATLAB: Get parameters for normal distributions in Kernel fit of probability distribution

Hello Darin,

I understand that you would like to get information on the normal distributions that make up the Kernal distribution. I assume that the Kernal distribution is calculated using the fitdist function.

If you check out the "Kernal Smoothing Function" example on the Kernal distribution documentation page, you can see that they actually plot the normal distributions that make up that particular Kernal distribution.

The mean of each normal distribution will be one of the data points, and the sigma of each normal distribution will be the 'BandWidth' of the Kernal distribution. Each normal distribution has the same standard deviation.

For example, if you wanted to modify that example, you could allow the default 'BandWidth' to be chosen, and add together the calculated normal distributions to show that they are the same as the plotted Kernal distribution:

% Plot the Kernal distribution
SixMPG = [13;15;23;29;32;34];
figure;
pdSix = fitdist(SixMPG,'Kernel');
x = 0:.1:45;
ySix = pdf(pdSix,x);
plot(x,ySix,'k-','LineWidth',2);
% Plot each individual normal distribution and scale its appearance on the plot
hold on;
yNormSum = zeros(size(ySix));
for i=1:6
    pd = makedist('Normal','mu',SixMPG(i),'sigma',pdSix.BandWidth);
    y = pdf(pd,x);
    y = y/6;
    plot(x,y,'b:');
    yNormSum = yNormSum+y;
end
plot(x,yNormSum,'g--','LineWidth',1.5)
hold off

MATLAB: Creating a Bell Curve of Data

Just looking at your data, histogram() function is useful.

histogram(tmax01,30);%however many bins you want

You could overlay:

pd = fitdist(tmax01,'Normal');
x_values = min(tmax01):.1:max(tmax01);%not sure what step to use
y = pdf(pd,x_values);
hold on;
plot(x_values,y,'LineWidth',2);
histogram(tmax01,30);

Best Answer

Related Solutions

MATLAB: Get parameters for normal distributions in Kernel fit of probability distribution

MATLAB: Creating a Bell Curve of Data

Related Question