why my kernel density function never touches x axes?
显示 更早的评论
Hello:
I want plot kernel density function with the attached test data. I have used foll. code:
Xgrid1 = linspace(min(X),max(X),1000);
pdfEst_1 = ksdensity(X,Xgrid1,'function','pdf');
plot(Xgrid1,pdfEst_1,'Color','B','linew',1.8);
However, the generated density function never touches the x-axes & always moves up, without any closed bounds. Any possible ways that the generated PDFs touches x - axes?
4 个评论
the cyclist
2021-1-28
Can you upload a file with H_10?
Poulomi Ganguli
2021-1-28
the cyclist
2021-1-28
Is that what is in the file Test1? Just those 15 data points?
It will depend on the kernel you choose, but the estimated density will be zero (or approach zero) as you move away from where the function has support. You haven't defined any grid points that are far away from where X has support.
Poulomi Ganguli
2021-1-28
采纳的回答
the cyclist
2021-1-28
Now that we've clarified some things, let me capture my reply in an actual answer.
The issue is that you have defined your X range to cover only places where the estimated PDF has support, because you actually have some data there. If, instead, you extend the range far enough away from the data, then the estimated PDF will go to zero (or at least approach zero).
For example:
extendRange = (max(X)-min(X))/2;
Xgrid1 = linspace(min(X) - extendRange,max(X) + extendRange,1000);
pdfEst_1 = ksdensity(X,Xgrid1,'function','pdf');
figure
plot(Xgrid1,pdfEst_1,'Color','B','linew',1.8);
Exactly how far out you need to go will depend on the kernel you choose for ksdensity.
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Kernel Distribution 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
