# NOT RUN {
data_toy <- c(rnorm(100, -3, 1), rnorm(100, 3, 1))
grid <- seq(-7, 7, length.out = 50)
est_model <- PYdensity(y = data_toy, mcmc = list(niter = 200, nburn = 100),
output = list(grid = grid))
x <- 1.4
dBNPdens(est_model, x)
# }
Run the code above in your browser using DataLab