data(lansing)
maples <- lansing[as.character(lansing[,"species"]) == "maple", -3]
maplesBIC <- EMclust(maples)
maplesModel <- summary(maplesBIC, maples)
x <- grid1(100, range=c(0,1))
y <- x
xyDens <- do.call("dens", c(list(data=grid2(x, y)), maplesModel))
xyDens <- matrix(xyDens, ncol=100)
contour(xyDens)
points(maples, cex=.2, col="red")
image(xyDens)
points(maples, cex=.5)
Run the code above in your browser using DataLab