swh:1:snp:ffdd0a7d2c8ea15ad41d45b3b178f668bd942287
Tip revision: 80ce80871d172ce0bb60c804a03eeb3a37e1bbd4 authored by Derek Young on 30 May 2010, 00:00:00 UTC
version 0.4.4
version 0.4.4
Tip revision: 80ce808
wkde.R
wkde <- function(x, u=x, w=rep(1, length(x)), bw=bw.nrd0(as.vector(x)), sym=FALSE) {
if (sym) {
return((wkde(x, u, w, bw) + wkde(x, -u, w, bw))/2)
}
Km <- exp(outer(x/bw, u/bw, function(a,b) -(a-b)^2/2))
normw <- matrix(w/sum(w), nrow=1)
as.vector(normw %*% Km) / (bw * sqrt(2*pi))
}