Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
35810356.5 35492936.3 38937434.1 40059102.8 37708710.2 41570965.7 36333563 34181220.1 42593543.9 43119727.6 38497690.9 45473273.4 38399780.4 38882302.6 44051120.6 41677559.9 40699203.5 44150027.6 38225518.7 35447405.7 43075518.3 42302792 39743541.7 44670641.2 37123384 37668266.4 46117528.8 42273156.4 39404153.2 45799994.5 38602505.2 39454830.1 47427901.4 46497980.9 45057149.4 50615569.2 43033396.2 46013056.5 54222266.3 46417306.4 51046271.8 51201279.6 43475288.7 44968981.1 53939345.4 54549319.7 54072107.3 58434230.1 51158751 50039368 57872617.4 51642978.8 54534465.9 56094697.8 48189983.1 47492381 52987449.1 55719803.5 53922860.5 54931231.9
Data Y:
1202454.6 1201423.4 1505916 1513377.6 1977605.3 1873829.6 1424049.1 1322740 1584825.5 1680460.3 1648573.7 3095468.7 1307982.9 1367588.9 1572718.3 1611602.9 1641196.4 1845262.4 1464237.6 1402385.7 2077099.8 1691129.6 1729012.7 3347792.1 1365087.7 1545460 1844355.1 1775549.8 1721779.2 2128726.1 1664319.9 1769471.4 1904578.4 1872042.3 1802181 3222199.4 1491414.2 1658519.2 2079206.9 1748767.4 2084447.4 2067181.6 1718122.8 1782337.1 1958118.4 2028681.3 2076128.1 3383873 1870369 1654852.9 2074338.3 1888653.7 1991137.8 2168237.9 1867424.1 1842359.6 1927476.3 2065555.4 2455608.5 3336170.9
Data Z:
5.89 5.65 5.85 6.02 5.90 5.83 5.64 5.75 5.69 5.69 5.68 5.45 5.22 5.11 5.03 5.03 5.09 4.96 4.88 4.66 4.34 4.28 4.33 4.09 3.90 4.04 4.26 4.11 4.29 4.64 4.94 5.18 5.34 5.58 5.30 5.41 5.79 5.79 5.62 5.52 5.69 5.53 5.60 5.56 5.63 5.58 5.52 5.28 5.16 5.16 5.08 5.21 5.38 5.33 5.35 5.15 5.14 4.89 4.75 4.97
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