Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
93.9 89.8 93.4 101.5 110.4 105.9 108.4 113.9 86.1 69.4 101.2 100.5 98.0 106.6 90.1 96.9 125.9 112.0 100.0 123.9 79.8 83.4 113.6 112.9 104.0 109.9 99.0 106.3 128.9 111.1 102.9 130.0 87.0 87.5 117.6 103.4 110.8 112.6 102.5 112.4 135.6 105.1 127.7 137.0 91.0 90.5 122.4 123.3 124.3 120.0 118.1 119.0 142.7 123.6 129.6 151.6 110.4 99.3 129.1 134.1
Data Y:
103.8 100.8 110.6 104.0 112.6 107.3 98.9 109.8 104.9 102.2 123.9 124.9 112.7 121.9 100.6 104.3 120.4 107.5 102.9 125.6 107.5 108.8 128.4 121.1 119.5 128.7 108.7 105.5 119.8 111.3 110.6 120.1 97.5 107.7 127.3 117.2 119.8 116.2 111.0 112.4 130.6 109.1 118.8 123.9 101.6 112.8 128.0 129.6 125.8 119.5 115.7 113.6 129.7 112.0 116.8 127.0 112.1 113.3 120.5 127.7
Data Z:
101.2 93.1 84.2 85.8 91.8 92.4 80.3 79.7 62.5 57.1 100.8 100.7 86.2 83.2 71.7 77.5 89.8 80.3 78.7 93.8 57.6 60.6 91.0 85.3 77.4 77.3 68.3 69.9 81.7 75.1 69.9 84.0 54.3 60.0 89.9 77.0 85.3 77.6 69.2 75.5 85.7 72.2 79.9 85.3 52.2 61.2 82.4 85.4 78.2 70.2 70.2 69.3 77.5 66.1 69.0 79.2 56.2 64.5 77.4 88.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