## univariate example
x <- rnorm.mixt(n=100, mus=1, sigmas=1, props=1)
fhat1 <- kdde(x=x, deriv.order=1)
plot(fhat1)
## bivariate example
data(unicef)
fhat2 <- kdde(x=unicef, deriv.order=2)
plot(fhat2, which.deriv.ind=2, display="persp", phi=20)
plot(fhat2, which.deriv.ind=2, display="filled.contour2", col.fun=topo.colors)
## d^2 f/(dx dy): purple=-ve, green=zero, beige=+ve
Run the code above in your browser using DataLab