Guest User

Untitled

a guest
Jan 19th, 2019
82
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.45 KB | None | 0 0
  1. EmbedContour <- function(e, bw=1, grid=c(201,201),
  2. sub.ratio=.5, size.ratio=.5, cq=.666, add=T, ...) {
  3. krnn <- KernSmooth::bkde2D(e,bandwidth=bw*sub.ratio,gridsize=grid)
  4. krnw <- KernSmooth::bkde2D(e, bandwidth=bw, gridsize=grid,
  5. range.x=list(
  6. c(min(krnn$x1),max(krnn$x1)),
  7. c(min(krnn$x2),max(krnn$x2))))
  8. con <- list(x=krnw$x1, y=krnw$x2, z=krnn$fhat-sub.ratio*krnw$fhat)
  9. contour(con, add=add, drawlabels=F, levels=c(quantile(con$z, cq)), ...)
  10. }
Add Comment
Please, Sign In to add comment