ksdensity default bandwidth value

Hi,
Does any one knows how matlab computes the default bandwidth in the function ksdensity?
Thank you very much,

回答(3 个)

FK
FK 2022-1-22
The 'ksdensity' function uses "Silverman's rule" as given in the post of Ghada. The difference between the calculations of Alexandra and the calculations done in Matlab is caused by a different estimation of the the standard deviation. Instead of estimating it by it is done by applying a more robust estimator with (see this Wikipedia article for further information).
With m=[1,3,2,5] you can derive with Silverman's rule that the default bandwith is given by .

2 个评论

thank you, actually I was study how to obtain the default bandwith too
Hi FK, I can reproduce the results if the MAD is > 0. If it is zero, MATLAB seems to use something else. Any idea what that could be?

请先登录,再进行评论。

Hello Alexandra,
In MATLAB, the default bandwidth in 'ksdensity' function is calculated using the following formula:
BandWidth = sig * (4/(3*N))^(1/5);
where 'sig' is an estimated value of the standard deviation of the samples.
I hope this helps,
Ghada

1 个评论

I have used this formula to test the calculations of the bandwidth and I don't think it is correct.
Take a simple example: m=[1,3,2,5] Then from [a,b,c]=ksdensity(m,m,'function','cdf') c=1.1901 and sig*(4/(3*4))^(1/5)= 1.3709 So they don't yield the same answers.

请先登录,再进行评论。

you can try this: [f, xi, bw] = ksdensity(X); bw = 'gives you the default bandwidth'

标签

评论:

2024-9-19

Community Treasure Hunt

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

Start Hunting!

Translated by