Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
1023.10 1141.00 1116.30 1135.60 1210.50 1230.00 1136.50 1068.70 1372.50 1049.90 1302.20 1305.90 1173.50 1277.40 1238.60 1508.60 1423.40 1375.10 1344.10 1287.50 1446.90 1451.00 1604.40 1501.50 1522.80 1328.00 1420.50 1648.00 1631.10 1396.60 1663.40 1283.00 1582.40 1785.20 1853.60 1994.10 2042.80 1586.10 1942.40 1763.60 1819.90 1836.00 1447.50 1509.50 1661.20 1456.20 1310.90 1542.10 1537.70
Data Y:
10881.30 11301.20 13643.90 12517.00 13981.10 14275.70 13425.00 13565.70 16216.30 12970.00 14079.90 14235.00 12213.40 12581.00 14130.40 14210.80 14378.50 13142.80 13714.70 13621.90 15379.80 13306.30 14391.20 14909.90 14025.40 12951.20 14344.30 16093.40 15413.60 14705.70 15972.80 16241.40 16626.40 17136.20 15622.90 18003.90 16136.10 14423.70 16789.40 16782.20 14133.80 12607.00 12004.50 12175.40 13268.00 12299.30 11800.60 13873.30 12315.00
Data Z:
1.20 1.23 1.23 1.20 1.18 1.19 1.21 1.19 1.20 1.23 1.28 1.27 1.27 1.28 1.27 1.26 1.29 1.32 1.30 1.31 1.32 1.35 1.35 1.34 1.37 1.36 1.39 1.42 1.47 1.46 1.47 1.47 1.55 1.58 1.56 1.56 1.58 1.50 1.44 1.33 1.27 1.34 1.32 1.28 1.31 1.32 1.37 1.40 1.41
Sample Range:
(leave blank to include all observations)
From:
To:
gridsize on x-axis
(?)
gridsize on y-axis
(?)
plot contours
Y
N
plot points
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