Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
4651.600942 4648.46719 4629.664677 4595.193403 4573.257137 4592.05965 4751.881013 4805.154801 4811.422306 4805.154801 4751.881013 4755.014766 4742.479757 4736.212252 4720.543491 4686.072217 4667.269704 4673.537208 4836.492323 4861.562341 4858.428589 4808.288553 4755.014766 4761.28227 4748.747261 4739.346005 4708.008482 4689.205969 4686.072217 4686.072217 4833.358571 4852.161084 4833.358571 4733.0785 4664.135951 4635.932181 4648.46719 4610.862164 4557.588376 4538.785863 4491.779579 4454.174553 4629.664677 4661.002199 4588.925898 4541.919615 4488.645827 4501.180836 4510.582093 4482.378323 4435.37204 4425.970783 4360.161986 4388.365756 4545.053367 4563.85588 4513.715845 4476.110818 4463.57581 4507.448341
Data Y:
52.034 58.146 79.035 63.969 34.831 42.602 23.589 21.974 29.651 39.547 27.957 10.511 63.044 50.898 63.547 55.028 40.073 56.080 28.091 29.211 33.415 35.336 31.502 1.007 75.840 69.091 88.380 61.147 59.750 50.533 33.467 28.944 32.053 42.095 27.840 -686 69.362 55.107 73.801 50.981 48.575 52.660 37.173 25.894 30.201 46.248 26.577 1.911 68.401 65.458 67.183 76.063 48.064 52.591 34.502 22.268 30.174 39.902 14.115 -2.715
Data Z:
562000 561000 555000 544000 537000 543000 594000 611000 613000 611000 594000 595000 591000 589000 584000 573000 567000 569000 621000 629000 628000 612000 595000 597000 593000 590000 580000 574000 573000 573000 620000 626000 620000 588000 566000 557000 561000 549000 532000 526000 511000 499000 555000 565000 542000 527000 510000 514000 517000 508000 493000 490000 469000 478000 528000 534000 518000 506000 502000 516000
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