Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
116.1 102.5 102.0 101.3 100.6 100.9 104.2 108.3 108.9 109.9 106.8 112.7 113.4 101.3 97.8 95.0 93.8 94.5 101.4 105.8 106.6 109.7 108.8 113.4 113.7 103.6 98.2 95.5 94.4 95.9 103.2 104.1 127.6 130.3 133.0 140.4 123.5 116.9 115.9 113.1 112.1 112.4 118.9 117.4 115.6 120.7 114.9 122.0 119.6 114.6 118.4 110.9 111.6 114.6 112.1 117.4 114.8 123.4 118.1 121.9 123.3
Data Y:
117.1 107.0 107.0 111.0 108.2 96.3 100.9 107.7 106.2 118.7 116.1 118.1 118.4 110.8 106.4 112.2 108.3 96.0 100.6 107.8 108.4 120.9 117.3 119.7 119.6 111.8 108.1 111.8 105.5 93.6 103.9 100.3 106.6 118.4 106.6 109.8 115.9 111.7 119.8 116.1 103.2 99.0 112.3 104.2 114.0 121.7 107.2 112.8 117.8 113.3 116.1 111.8 110.2 110.0 102.9 110.1 102.7 118.7 109.0 115.7 118.1
Data Z:
118.9 108.8 115.6 95.0 92.8 108.9 109.8 106.1 102.8 98.4 85.7 114.6 129.4 117.7 126.6 103.8 101.5 118.7 119.6 114.8 109.9 106.3 95.0 124.5 140.4 128.8 137.5 113.3 110.3 129.1 128.4 120.3 113.6 96.9 124.7 126.4 131.9 122.5 113.1 99.8 116.0 115.0 114.0 111.0 91.7 90.6 103.3 106.7 111.2 102.9 126.5 115.1 110.2 110.1 103.3 107.7 103.9 114.0 117.2 117.0 116.5
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