MVAcontbank3 (R 2.9.1)

MVAcontbank3 gives a contour plot of the kernel density estimate of variables X4 , X5 and X6 of the Swiss bank notes ("bank2.dat").

Download File

Song Song

Applied Multivariate Statistical Analysis

Thu, April 19 2012 by Dedy Dwi Prastyo

MVAcontbank3

bank2.dat


-


rm(list=ls(all=TRUE))
install.packages(c("KernSmooth","misc3d"))
library(KernSmooth)
library(misc3d)
xx<-read.table("bank2.dat")
d <- kde3d(xx[, 4], xx[, 5], xx[, 6], n = 15)
contour3d(d$d, level = c(max(d$d[10,10,])*.02, max(d$d[10,10,])*.5, max(d$d[10,10,])*1.3), fill = c(FALSE, FALSE, TRUE), col.mesh = c("green","red","blue") , engine = "standard", screen=list(z=210,x=-40,y=-295), scale=TRUE)