This function implements bivariant Gaussian kernel density estimation. It can be used to estimate bivariant probability density function (pdf), cumulative distribution function (cdf) and inversed cdf (icdf) from a set of random data. The code is programmed in a way to be suitable for beginners to learn kernel density estimation.
Cite As
Yi Cao (2020). Bivariant Kernel Density Estimation (V2.1) (https://www.mathworks.com/matlabcentral/fileexchange/19280-bivariant-kernel-density-estimation-v2-1), MATLAB Central File Exchange. Retrieved .
Comments and Ratings (18)
Updates
| 1.2.0.0 | a bug fixed |
|
| 1.1.0.0 | new functionality |
MATLAB Release Compatibility
Platform Compatibility
Windows macOS LinuxTags
Acknowledgements
Inspired by: Probability Density Function (PDF) Estimator (V3.2)
Discover Live Editor
Create scripts with code, output, and formatted text in a single executable document.
html/
Select a Web Site
Choose a web site to get translated content where available and see local events and offers. Based on your location, we recommend that you select: .
Select web siteYou can also select a web site from the following list:
Americas
- América Latina (Español)
- Canada (English)
- United States (English)
Europe
- Belgium (English)
- Denmark (English)
- Deutschland (Deutsch)
- España (Español)
- Finland (English)
- France (Français)
- Ireland (English)
- Italia (Italiano)
- Luxembourg (English)
- Netherlands (English)
- Norway (English)
- Österreich (Deutsch)
- Portugal (English)
- Sweden (English)
- Switzerland
- United Kingdom (English)
Takumi Sase (view profile)
Sorry, please ignore the comment below.
Takumi Sase (view profile)
How to optimaze the number of peaks in the code?
Larissa (view profile)
So the function works once you figure out how the data needs to be structured. My only problem is that, as written, the actual pdf values scale down as your data range gets larger. Should this really be happening?
Up to and including the line
z=(X(k+zeros(n,1),:)-x)./h(ones(n,1),:);
the values are the same regardless of data scale. However, the next line scales by the bandwidth AGAIN and messes that up. Is this accurate? If so, could you explain the reasoning for this?
Joseph (view profile)
examples do work, and so does the code. Really straightforward coding and simple to manipulate for use in your individual needs, even if you're new to KDE
Xinyao Tang (view profile)
It seems that your bandwidth formula doesn't correct.
Adwin555 (view profile)
I just delete % to run the examples, but I got "Error in gkde2"
Matthew (view profile)
documentation is false, misleading, and old.
Michael (view profile)
nice code. I like that you put your references for methods, and I need to double-check them. Wish you had a link to Wolfram or an online resource.
Internal nomenclature could be cleaned up. The command call in the example should be able to be copied, pasted to command-line and executed.
sebastiano (view profile)
Yong-ha (view profile)
yolgs Hoffmann (view profile)
ah ok: The functions does not get p.xy, the information in the file is misleading. But instead requires p.x and p.y separately. It´s in the check function around line 180 f.
yolgs Hoffmann (view profile)
Nice one. But to really use it I definitely need to put in my own grid locations. So I have the same problem as Frank. Please - is there a solution? That would be really good.
Frank (view profile)
Nice program with the default values. Put I'm not seeing how I can input my own grid locations. The optional inputs say "p.xy - Nx2 matrix for locations to make estimation", however, I don't see where this is done within the program. Also, I get errors on line 260 when I try to input my own bandwidth values.
Turkay YILDIZ (view profile)
Warrick (view profile)
Very easy to understand for a Matlab novice like myself. However, is the bandwidth formula correct? It doesn't seem to match the reference? Also, my total probability density is almost, but not exactly 1 (~0.96)
Tyler Jandreau (view profile)
V. Poor (view profile)
very handy, thank you