Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
150.85 147.79 141.96 148.39 147.71 150.6 151.18 152.24 157.19 154.62 157.22 159.7 160.55 149.66 151.69 154.13 151.48 153.34 155.8 158.87 156.09 156.3 156.4 154.09 161.32 160.12 155.17 154.51 151.38 152.59 153.98 154.91 153.01 155.09 155.53 161.86 166.03 164.54 164.33 163.21 159.95 164.18 167.13 166.8 166.29 168.07 167.1 163.53 168.28 169.07 165.84 163.88 157.33 161 163.54 161.21 158.92 160.18 159.9 164.46
Data Y:
128.6 128.9 129.06 129.23 129.27 129.33 129.35 129.31 129.4 129.49 129.47 129.46 129.45 129.28 129.2 129.25 129.14 129.11 129.02 129.08 128.99 129.11 129.08 129.19 129.23 129.25 129.31 129.33 129.39 129.55 129.43 129.45 129.57 129.76 129.92 130.08 130.41 130.84 131.24 131.49 131.74 132.34 133.5 134.43 136.5 137.41 138.02 138.15 138.24 138.2 138.31 138.65 139.3 139.8 140.52 141.57 141.77 141.66 141.36 141.17
Data Z:
131.6 132.05 132.4 132.57 133.02 133.47 133.66 133.96 134.19 134.93 134.9 135.05 135.16 135.23 135.15 135.12 137.29 137.41 137.44 137.62 137.78 137.98 138.06 138.16 138.28 138.33 138.43 138.44 138.41 138.55 138.64 138.72 138.9 139.02 139.04 139.15 139.3 140.73 141.84 141.95 142.1 142.36 142.58 142.75 142.85 143.03 143.19 143.62 143.89 144.69 147.51 147.78 148.04 148.21 148.29 148.34 148.33 148.38 148.37 148.37
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