Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
9 8 10 8 9 10 9 9 9 9 10 9 13 11 10 8 11 6 8 9 9 10 9 10 8 10 10 11 9 9 7 11 6 11 10 9 10 10 9 8 9 9 10 11 7 9 10 8 3 10 10 10 4 10 8 9 13 10 8 9 11 10 9 10 7 10 10 11 12 8 10 6 9 11 10 10 8 10 9 9 10 10 11 9 12 7 9 9 11 8 9 9 9 9 11 9 7 15 9 9 12 10 9 10 10 9 10 10 9 9 9 9 11 9 7 11 9 7 12 8 9 9 9 9 11 9 10 10 9 8 10 9 8 10 9 13 13 11 9 8 9 10 10 10 9 11 7 9 10 11 8 8 7 9 9 10 9 11 8 11 9 9 9 7 10 9 9
Data Y:
10 13 14 12 12 13 13 13 13 14 14 12 12 11 12 14 12 11 13 13 12 13 12 13 11 12 12 13 13 10 12 13 13 10 14 12 10 10 14 12 14 10 13 12 12 13 12 10 9 14 15 14 8 11 10 12 14 12 12 14 13 13 13 12 10 14 11 10 13 12 12 11 10 14 12 13 11 10 14 13 7 13 13 13 15 13 14 12 13 11 12 14 13 14 12 12 13 14 13 12 13 12 10 12 13 12 13 12 12 12 11 12 9 14 12 13 13 13 11 12 11 12 12 13 12 13 13 12 12 8 12 13 10 8 12 13 12 15 14 10 11 12 10 14 10 15 11 12 9 12 13 12 9 12 14 11 12 14 12 15 11 12 12 10 12 11 11
Data Z:
18 21 21 18 21 21 17 19 23 23 21 20 20 19 19 20 20 21 21 19 22 21 20 21 19 20 20 20 21 19 23 22 21 18 21 19 19 19 20 21 23 18 23 20 20 20 18 21 15 21 23 22 21 19 20 17 18 20 21 24 20 22 17 22 18 19 20 20 20 19 23 20 20 22 19 21 20 19 19 21 18 21 22 22 21 23 22 20 19 17 19 21 21 21 20 21 22 23 23 20 22 21 22 18 24 21 18 20 20 19 20 20 19 19 19 19 20 23 19 19 19 20 21 21 19 19 21 18 18 15 20 21 17 18 19 23 21 24 23 20 20 20 20 22 22 21 19 20 17 20 19 21 20 21 22 19 17 21 20 19 22 21 20 20 21 21 19
Sample Range:
(leave blank to include all observations)
From:
To:
gridsize on x-axis
(?)
gridsize on y-axis
(?)
plot contours
Y
N
plot points
Y
N
Name of dataset X
Name of dataset Y
Name of dataset Z
Chart options
R Code
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)) print(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)) print(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
1 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation