Probability Density Function using ksdensity is not normalized

15 次查看（过去 30 天）
Ali
Ali 2014-6-28
回答： VladTheInstaller 2017-1-15
I have a vector "columnA" of N data points. I want to find the PDF. I use:
xi = min(columnA):1e-9:max(columnA);
f = ksdensity(columnA,xi);
plot(xi,f)
But when I use trapz to integrate f:
trapz(f)/length(xi)
the value is too far from 1. Even when increasing the range of xi, I still do not get reasonable value.

请先登录，再回答此问题。

回答（3 个）

VladTheInstaller
VladTheInstaller 2017-1-15
Actually, the output from ksdensity is normalized, but you will have to use numerical integration along the appropriate space. In your case,
trapz(xi,f)
should be close to 1.
Image Analyst
Image Analyst 2014-8-21
Why not use hist() or histc() to get the histogram? The histogram is essentially the probability density function.
Youssef Khmou
Youssef Khmou 2014-8-21
The ksdensity produces a Probability density function, no need to divide by the length of the x vector :
x=randn(200,1);
y=[min(x):0.1:max(x)];
p=ksdensity(x,y);
sum(p)
% plot(y,p)

类别

在 Help Center 和 File Exchange 中查找有关 Exploration and Visualization 的更多信息

标签

产品

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by