Geom_density_2d

set.seed(1001)
d <- data.frame(x=rnorm(1000),y=rnorm(1000))
getLevel <- function(x,y,prob=0.95) {
    kk <- MASS::kde2d(x,y)
    dx <- diff(kk$x[1:2])
    dy <- diff(kk$y[1:2])
    sz <- sort(kk$z)
    c1 <- cumsum(sz) * dx * dy
    approx(c1, sz, xout = 1 - prob)$y
}
L95 <- getLevel(d$x,d$y)
library(ggplot2); theme_set(theme_bw())
ggplot(d,aes(x,y)) +
   stat_density2d(geom="tile", aes(fill = ..density..),
                  contour = FALSE)+
   stat_density2d(colour="red",breaks=L95)

Source: stackoverflow.com

Similar codes

-1 vote227 views1 answer
Matlab invert image
Matlab
1 vote281 views1 answer
Matlab unix time to datetime
Matlab
0 vote223 views1 answer
BIDS json IntendedFor field examples
Matlab
0 vote241 views1 answer
Symbolic polynomial fraction to row vector
Matlab
7 votes290 views2 answers
To detect if a data frame has nan values
Matlab
1 vote259 views1 answer
Load value inside job matlab
Matlab
0 vote228 views1 answer
Octave get range length
Matlab
1 vote278 views1 answer
Matlab plot with motion
Matlab
1 vote265 views1 answer
Scilab create multi parameter function
Matlab
0 vote237 views2 answers
Octave plot legend
Matlab