Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
100.70 97.90 96.50 96.60 96.60 95.50 91.80 89.30 87.00 85.90 88.00 87.90 89.20 90.90 91.60 90.20 89.10 87.50 86.30 86.00 84.40 86.10 91.00 92.70 88.00 84.30 82.20 80.80 79.40 80.20 82.20 82.20 81.20 82.10 88.10 88.50 92.10 98.60 100.90 100.60 101.10 102.10 103.60 102.80 108.30 104.00 106.10 106.30 109.00 111.00 113.70 112.70 110.30 114.50 119.30 121.80 125.40 129.70 129.40 134.50 141.20 141.40 152.20 167.70 173.30 168.70 172.60 169.80 172.00 179.40 174.60 172.50 172.60 176.30 178.90 179.60 179.90 180.30 180.90 177.70
Data Y:
101.30 97.60 96.40 97.00 96.40 94.70 89.30 85.90 83.30 81.50 85.00 84.80 87.50 89.00 90.00 89.60 87.40 84.80 81.90 81.10 79.10 80.50 88.50 90.90 84.90 80.00 76.50 75.40 73.50 74.30 77.70 77.90 76.70 77.20 86.00 86.90 92.00 101.70 104.50 101.70 100.60 100.30 102.50 101.00 108.60 103.40 106.40 106.60 108.90 110.50 114.00 112.80 109.60 116.00 124.60 129.00 131.50 138.60 138.10 146.30 157.60 158.40 176.30 199.90 210.40 202.60 207.10 202.00 203.40 216.30 207.30 203.50 204.40 203.70 205.70 208.00 209.30 208.70 206.50 204.50
Data Z:
99.80 103.50 103.10 105.60 105.70 106.60 107.00 105.20 105.70 105.00 105.10 105.90 105.30 104.90 103.20 103.40 104.40 104.50 105.90 110.60 112.40 111.80 111.00 111.00 109.10 107.80 107.20 108.40 107.50 106.40 106.20 104.90 106.20 107.60 107.00 104.50 105.10 104.70 103.70 104.90 105.90 106.10 106.10 106.80 106.40 107.80 107.60 107.60 108.40 109.50 109.20 109.10 110.00 109.00 109.00 111.90 109.30 112.10 112.10 112.50 113.60 112.90 114.00 116.10 116.50 117.10 117.10 117.10 116.50 116.50 116.30 116.50 119.20 118.60 117.50 117.10 117.60 118.30 118.60 116.70
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