library(MASS)
data(iris)
## univariate example
ir <- iris[,1]
ir.gr <- iris[,5]
kda.fhat <- kda(x=ir, x.group=ir.gr, xmin=3, xmax=9)
plot(kda.fhat, xlab="Sepal length")
## bivariate example
ir <- iris[,1:2]
ir.gr <- iris[,5]
kda.fhat <- kda(x=ir, x.group=ir.gr)
plot(kda.fhat)
## trivariate example
ir <- iris[,1:3]
ir.gr <- iris[,5]
H <- Hkda(x=ir, x.group=ir.gr, bw="plugin", pilot="dscalar")
kda.fhat <- kda(x=ir, x.group=ir.gr, Hs=H)
plot(kda.fhat) ## colour=species, transparency=density heights
Run the code above in your browser using DataLab