Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- EmbedContour <- function(e, bw=1, grid=c(201,201),
- sub.ratio=.5, size.ratio=.5, cq=.666, add=T, ...) {
- krnn <- KernSmooth::bkde2D(e,bandwidth=bw*sub.ratio,gridsize=grid)
- krnw <- KernSmooth::bkde2D(e, bandwidth=bw, gridsize=grid,
- range.x=list(
- c(min(krnn$x1),max(krnn$x1)),
- c(min(krnn$x2),max(krnn$x2))))
- con <- list(x=krnw$x1, y=krnw$x2, z=krnn$fhat-sub.ratio*krnw$fhat)
- contour(con, add=add, drawlabels=F, levels=c(quantile(con$z, cq)), ...)
- }
Add Comment
Please, Sign In to add comment