Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
0.527 0.472 0.000 0.052 0.313 0.364 0.363 -0.155 0.052 0.568 0.668 1.378 0.252 -0.402 -0.050 0.555 0.050 0.150 0.450 0.299 0.199 0.496 0.444 -0.393 -0.444 0.198 0.494 0.133 0.388 0.484 0.278 0.369 0.165 0.155 0.087 0.414 0.360 0.975 0.270 0.359 0.169 0.381 0.154 0.486 0.925 0.728 -0.014 0.046 -0.819 -1.674 -0.788 0.279 0.396 -0.141 -0.019 0.099 0.742 0.005 0.448
Data Y:
4.79 5.95 5.46 5.75 5.15 4.96 5.28 5.73 5.75 5.88 6.3 6.74 6.75 7.34 6.64 6.62 6.32 5.32 5.68 6.18 5.02 2.1 4 3 4.73 5.14 5.81 6.24 4.49 4.22 4.88 5.18 5.19 5.06 4.65 4.83 4.6 4.72 4.33 4.97 5.37 4.19 4.54 5.82 5.49 3.28 5.11 6.24 6.41 6.43 8.42 8.23 3.17 2.72 3 3.47 3.88 3.43 4.06
Data Z:
24710.92 23983.59 24434.12 23939.23 24290.02 24117.63 23724.64 22989.44 23716.86 25058.83 25059.00 23579.18 24209.03 24173.67 24706.39 24522.12 24766.15 25940.04 24985.78 24788.00 26544.56 28019.08 27285.71 29161.16 28357.73 27979.91 27543.95 27397.53 27623.59 27736.07 27803.79 27779.55 27524.13 27582.72 28638.95 28825.78 30132.61 29326.85 29075.62 28230.63 28118.36 28173.29 27396.91 24578.55 24504.77 27582.37 26920.31 25426.68 25390.80 25041.16 22769.42 22921.89 26267.63 27364.67 28382.59 29132.81 28214.51 28865.73 24405.35
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