Wessa.net

Send output to:


Data X:
0.275 0.275 0.275 0.275 0.275 0.28 0.28 0.28 0.28 0.28 0.275 0.28 0.265 0.265 0.265 0.265 0.265 0.265 0.265 0.265 0.265 0.265 0.265 0.27 0.27 0.27 0.275 0.275 0.27 0.275 0.28 0.29 0.295 0.3 0.3 0.3 0.295 0.3 0.3 0.31 0.325 0.34 0.365 0.39 0.39 0.41 0.41 0.405 0.415 0.425 0.43 0.425 0.425 0.41 0.4 0.405 0.4 0.4 0.4 0.4

Data Y:
1.58 1.59 1.6 1.6 1.61 1.61 1.61 1.62 1.63 1.63 1.64 1.64 1.64 1.64 1.64 1.65 1.65 1.65 1.65 1.65 1.66 1.66 1.67 1.68 1.68 1.68 1.68 1.69 1.7 1.7 1.71 1.72 1.73 1.74 1.74 1.75 1.75 1.75 1.76 1.79 1.83 1.84 1.85 1.87 1.87 1.87 1.88 1.88 1.88 1.88 1.89 1.89 1.89 1.9 1.89 1.89 1.89 1.89 1.89 1.89

Data Z:
5.24 5.28 5.34 5.36 5.38 5.39 5.41 5.44 5.51 5.55 5.56 5.57 5.58 5.58 5.59 5.61 5.63 5.64 5.64 5.61 5.64 5.66 5.65 5.69 5.68 5.68 5.7 5.72 5.73 5.76 5.77 5.82 5.89 5.91 5.92 5.93 5.94 5.94 5.98 6.11 6.37 6.44 6.52 6.64 6.68 6.73 6.75 6.78 6.77 6.79 6.8 6.83 6.91 6.97 7 6.99 7.03 7.03 7.04 7.04


Sample Range: (leave blank to include all observations)
From:
To:

gridsize on x-axis (?)


gridsize on y-axis (?)


plot contours


plot points


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~xy, screen = list(x=-45, y=45, z=35),xlab=par5,ylab=par6,zlab=par7) dev.off() bitmap(file='cloud2.png') cloud(z~xy, 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()

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