How can I extend the normal distribution curve in histfit function

Hello, I am trying to plot a normal distribution curve together with a kernel density and the histogram of the data using the code below. However, when I plot them I can see that the normal distribution curve is truncated up to some point as you can see in the image below. My question is How can I extend this curve a little bit more?

          figHandler=figure;
        h1=histfit(rho_fec(:,1),30,'normal');
        set(h1(1),'facecolor',[.8 .8 1]);
        set(h1(2), 'color' ,'b');
        hold on;
        h2=histfit(rho_fec(:,1),30,'kernel');
        set(h2(1),'facecolor',[.8 .8 1]);
        uistack(h1,'up');
        axis([0.3 0.6 0 110]);
        xlabel('$\hat{\rho}$','Interpreter','latex','FontSize',13)
        l=legend([h1(2) h2(2)],' Std. Normal density',' Kernel density',...
                'Location',[0.67,0.80,0.1,0.1]);
        set(l,'Interpreter','latex','FontSize',14)
        legend boxoff;
        hold off;

Answer 1

[the cyclist]

Hm. Definitely strange. I can't replicate it by just poking around randomly. It might depend on your data. You could try taking a look at the (x,y) data generated by histfit(), using

get(h1(2),'XData')
get(h1(2),'YData')

and see what they look like. Given your input argument, I would expect both XData and YData to be vectors of length 30. I wonder if there is a very tiny value or some weird NaN getting in there. If so, you could put a breakpoint inside the histfit() function and see what's going on.

[Alberto]

They are vectors of length 100, nothing weird xdata starts about 0.32, so that is exactly where the curve is starting...

[the cyclist]

Looking inside the histfit() function, I see the following code that defines the range of x:

% Find range for plotting
q = icdf(pd,[0.0013499 0.99865]); % three-sigma range for normal distribution
x = linspace(q(1),q(2));

It's semi-arbitrarily plotting a 3-sigma portion of the distribution.

If you want to do something different, you could modify and save your own version of that file.

[Alberto]

Effectively, increasing the upper and decreasing the lower limits the curve was plotted in more points. Thank you very much.

Answer 2

[the cyclist]

I don't think the issue is that the normal curve is truncated. I think the issue is that the histogram is plotted in front of the normal curve.

I expect you can use the uistack() command -- as you did earlier -- to move the normal curve to the front.

[Alberto]

Thank you for your answer, but apparently this is not the problem, sorry for not attaching a complete image of the plot. In the image below you can see how both curves are in front of the histogram.

Or maybe i am missing something like sending the histogram down, but I think I tried all I thought it would work with -uistack- and did not have success.

Thank you for your help...