[SciPy-user] guassian_kde and kernel regression
Robert Kern
robert.kern@gmail....
Sun Nov 2 02:18:22 CST 2008
On Sat, Nov 1, 2008 at 18:51, Anne Archibald <aarchiba@physics.mcgill.ca> wrote:
> 2008/11/1 Frank Lagor <dfranci@seas.upenn.edu>:
>> This question is probably for Robert Kern, because I believe the he wrote
>> the gaussian_kde class in scipy.stats.kde, however I would very much
>> appreciate a response from anyone else who could help. My question is: Is
>> there currently any way to perform weighted kernel density estimation using
>> the gaussian_kde class? If not, what needs to be done, and how do I get
>> started?
>>
>> Just for clarity sake-- by weighted KDE I mean that I have more than just
>> the distribution of points for the density estimate. I also have an
>> associated probability with each point. In this case, I believe it becomes
>> a regression problem and I think is referred to as kernel regression. I
>> would very much like to use the class to perform both KDE and wKDE.
>
> The class does not support weights right now, but I don't think it
> would be very difficult to add them to most parts of the code,
> essentially just adding a "weights" optional argument. The automatic
> covariance selection would need some rethinking; you'd need to hunt
> down some research papers. (That method is only really appropriate for
> unimodal distributions anyway.) But it does seem valuable.
What Anne said.
--
Robert Kern
"I have come to believe that the whole world is an enigma, a harmless
enigma that is made terrible by our own mad attempt to interpret it as
though it had an underlying truth."
-- Umberto Eco
More information about the SciPy-user
mailing list