Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
16.4 17.8 22.3 22.8 18.3 22.4 23.9 21.3 23.0 21.4 21.2 20.9 17.9 20.7 22.2 19.8 17.7 19.6 20.8 19.8 18.6 21. 18.6 18.9 17.3 20.0 19.9 19.5 16.2 17.6 19.8 19.4 17.2 21.1 17.8 17.5 18.0 19.1 17.7 19.2 15.1 16.3 18.6 17.2 17.8 19.1 16.6 16.0 16.7 17.4 17.9 17.8 13.9 15.9 17.9 15.4 16.4 17.9 15.3 14.6 14.9 15.0 16.7 16.3 11.7 15.1 15.5 15.0 15.4 16.0 14.7 14.8
Data Y:
18.0 19.6 23.3 23.7 20.3 22.8 24.3 21.5 23.5 22.2 20.9 22.2 19.5 21.1 22.0 19.2 17.8 19.2 19.9 19.6 18.1 20.4 18.1 18.6 17.6 19.4 19.3 18.6 16.9 16.4 19.0 18.7 17.1 21.5 17.8 18.1 19.0 18.9 16.8 18.1 15.7 15.1 18.3 16.5 16.9 18.4 16.4 15.7 16.9 16.6 16.7 16.6 14.4 14.5 17.5 14.3 15.4 17.2 14.6 14.2 14.9 14.1 15.6 14.6 11.9 13.5 14.2 13.7 14.4 15.3 14.3 14.5
Data Z:
94.3 99.4 115.7 116.8 99.8 96.0 115.9 109.1 117.3 109.8 112.8 110.7 100.0 113.3 122.4 112.5 104.2 92.5 117.2 109.3 106.1 118.8 105.3 106.0 102.0 112.9 116.5 114.8 100.5 85.4 86.6 109.9 100.7 115.5 100.7 99.0 102.3 108.8 105.9 113.2 95.7 80.9 113.9 98.1 102.8 104.7 95.9 94.6 101.6 103.9 110.3 114.1 96.8 87.4 111.4 97.4 102.9 112.7 97.0 95.1 96.9 98.6 111.7 109.8 89.9 87.4 104.5 98.1 102.7 105.4 97.0 97.4
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
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