How can I extend the normal distribution curve in histfit function
9 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
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;
0 commentaires
Réponse acceptée
the cyclist
le 3 Août 2014
Modifié(e) : the cyclist
le 3 Août 2014
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.
3 commentaires
the cyclist
le 5 Août 2014
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.
Plus de réponses (1)
the cyclist
le 2 Août 2014
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.
Voir également
Catégories
En savoir plus sur Descriptive Statistics and Visualization dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!