Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
28.9 23.7 21.2 26.1 22.2 16.8 19.7 21.1 22.9 24.7 20.8 23 25.5 23.5 23.2 25.8 27 23.4 23.9 18.4 25.4 22.6 29.5 24.9 19.7 23.2 21.1 21.6 21.9 21.5 22.1 24.2 24.6 25.3 24.9 21.2 24.3 19 25.1 23.8 27.2 25.8 20.4 18.7 23.1 22.2 27.6 23.2 20.8
Data Y:
30.6 32.4 36.4 27.6 26 33.4 23.2 35.4 33.8 29.5 30.4 31.3 28.3 34.4 30 36.8 37.1 31.2 26.1 27.8 29.3 31.3 35.1 33.4 41.5 36.6 44 43.4 24.1 37.8 23 30.2 41.7 33.5 27 30.6 29.6 28 37 39 30.2 37.4 19.5 31.4 29.3 27.9 31.3 38.2 36.4
Data Z:
1.72 3.83 3.29 1.34 4.53 3.66 4.62 5.14 4.09 2.38 3.06 3.1 2 1.33 1.34 1.39 2.35 3.57 2.95 4.89 2.47 2.66 1.06 2.36 3.06 1.71 5.75 6.26 4.46 2.25 3.61 2.43 1.56 2.01 1.36 4.27 2.09 4.36 2.06 1.5 1.6 2.13 1.37 4.75 3.09 4.1 0.88 2.63 2
Sample Range:
(leave blank to include all observations)
From:
To:
gridsize on x-axis
(?)
gridsize on y-axis
(?)
plot contours
Y
Y
N
plot points
Y
Y
N
Name of dataset X
Name of dataset Y
Name of dataset Z
Chart options
R Code
par7 <- 'Variable Z' par6 <- 'Variable Y' par5 <- 'Variable X' par4 <- 'Y' par3 <- 'Y' par2 <- '50' par1 <- '50' x <- array(x,dim=c(length(x),1)) colnames(x) <- par5 y <- array(y,dim=c(length(y),1)) colnames(y) <- par6 z <- array(z,dim=c(length(z),1)) colnames(z) <- par7 d <- data.frame(cbind(z,y,x)) colnames(d) <- list(par7,par6,par5) par1 <- as.numeric(par1) par2 <- as.numeric(par2) if (par1>500) par1 <- 500 if (par2>500) par2 <- 500 if (par1<10) par1 <- 10 if (par2<10) par2 <- 10 library(GenKern) library(lattice) panel.hist <- function(x, ...) { usr <- par('usr'); on.exit(par(usr)) par(usr = c(usr[1:2], 0, 1.5) ) h <- hist(x, plot = FALSE) breaks <- h$breaks; nB <- length(breaks) y <- h$counts; y <- y/max(y) rect(breaks[-nB], 0, breaks[-1], y, col='black', ...) } bitmap(file='cloud1.png') cloud(z~x*y, screen = list(x=-45, y=45, z=35),xlab=par5,ylab=par6,zlab=par7) dev.off() bitmap(file='cloud2.png') cloud(z~x*y, screen = list(x=35, y=45, z=25),xlab=par5,ylab=par6,zlab=par7) dev.off() bitmap(file='cloud3.png') cloud(z~x*y, screen = list(x=35, y=-25, z=90),xlab=par5,ylab=par6,zlab=par7) dev.off() bitmap(file='pairs.png') pairs(d,diag.panel=panel.hist) dev.off() x <- as.vector(x) y <- as.vector(y) z <- as.vector(z) bitmap(file='bidensity1.png') op <- KernSur(x,y, xgridsize=par1, ygridsize=par2, correlation=cor(x,y), xbandwidth=dpik(x), ybandwidth=dpik(y)) image(op$xords, op$yords, op$zden, col=terrain.colors(100), axes=TRUE,main='Bivariate Kernel Density Plot (x,y)',xlab=par5,ylab=par6) if (par3=='Y') contour(op$xords, op$yords, op$zden, add=TRUE) if (par4=='Y') points(x,y) (r<-lm(y ~ x)) abline(r) box() dev.off() bitmap(file='bidensity2.png') op <- KernSur(y,z, xgridsize=par1, ygridsize=par2, correlation=cor(y,z), xbandwidth=dpik(y), ybandwidth=dpik(z)) op image(op$xords, op$yords, op$zden, col=terrain.colors(100), axes=TRUE,main='Bivariate Kernel Density Plot (y,z)',xlab=par6,ylab=par7) if (par3=='Y') contour(op$xords, op$yords, op$zden, add=TRUE) if (par4=='Y') points(y,z) (r<-lm(z ~ y)) abline(r) box() dev.off() bitmap(file='bidensity3.png') op <- KernSur(x,z, xgridsize=par1, ygridsize=par2, correlation=cor(x,z), xbandwidth=dpik(x), ybandwidth=dpik(z)) op image(op$xords, op$yords, op$zden, col=terrain.colors(100), axes=TRUE,main='Bivariate Kernel Density Plot (x,z)',xlab=par5,ylab=par7) if (par3=='Y') contour(op$xords, op$yords, op$zden, add=TRUE) if (par4=='Y') points(x,z) (r<-lm(z ~ x)) abline(r) box() dev.off()
Compute
Summary of computational transaction
Raw Input
view raw input (R code)
Raw Output
view raw output of R engine
Computing time
0 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation