Home » date » 2008 » Dec » 12 »

Trivariate scatterplot

*The author of this computation has been verified*
R Software Module: rwasp_cloud.wasp (opens new window with default values)
Title produced by software: Trivariate Scatterplots
Date of computation: Fri, 12 Dec 2008 06:05:11 -0700
 
Cite this page as follows:
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL http://www.freestatistics.org/blog/date/2008/Dec/12/t12290872622calny2fiedkc1v.htm/, Retrieved Fri, 12 Dec 2008 14:07:52 +0100
 
BibTeX entries for LaTeX users:
@Manual{KEY,
    author = {{YOUR NAME}},
    publisher = {Office for Research Development and Education},
    title = {Statistical Computations at FreeStatistics.org, URL http://www.freestatistics.org/blog/date/2008/Dec/12/t12290872622calny2fiedkc1v.htm/},
    year = {2008},
}
@Manual{R,
    title = {R: A Language and Environment for Statistical Computing},
    author = {{R Development Core Team}},
    organization = {R Foundation for Statistical Computing},
    address = {Vienna, Austria},
    year = {2008},
    note = {{ISBN} 3-900051-07-0},
    url = {http://www.R-project.org},
}
 
 
Original text written by user:
 
IsPrivate?
No (this computation is public)
 
User-defined keywords:
 
Dataseries X:
» Textbox « » Textfile « » CSV «
21863 20403 18792 17931 16475 16205 25134 31896 26537 22801 20200 19666 19809 18799 17884 17512 16327 16880 26537 31867 29427 25800 22041 21759 21333 20462 19594 18564 17640 18614 32562 35640 31865 28117 25508 25006 24452 22643 21474 20500 19505 21769 36062 38633 34629 30184 27271 26841 26482 25538 23789 22386 21087 22891 36192 38922 34669 30197 27001 25891 24879 23662 22741 21615 20305 21877 35369 37941 33480 29757 26323 25359 22207 21763 19944 19662 18624 19902 31726 32860 28894 22949 19758 18420 18245 16761 15341 14271 13418 15218 26485 27457 21402 17879 15607 15626 15303 14296 13686 12948 11609 14602 23629 24680
 
Dataseries Y:
» Textbox « » Textfile « » CSV «
60804 57907 54355 52536 49081 48877 64599 75314 71209 65210 59829 57656 57428 55315 52790 51050 48519 48354 65333 73990 72755 67424 59214 57427 56681 55437 53600 51641 49478 50124 71313 76208 74387 69520 64735 63413 62553 60109 57764 55667 53103 55301 76795 80928 79213 72759 67802 66940 66396 67539 67776 68014 68251 68488 68725 68962 69200 69437 68212 65444 63181 61198 59010 56388 53723 55340 75352 79817 78289 71892 66448 64167 61250 59580 56417 54662 53349 55385 73546 77683 74995 67282 60742 57283 57314 54704 51578 49962 46252 47234 64708 68753 62970 57474 52494 51831 51663 49637 46679 45557 41630 44417 60070 63157
 
Dataseries Z:
» Textbox « » Textfile « » CSV «
30811 29877 28303 27605 26074 26112 32350 35804 36574 34486 32158 30965 30505 29629 28169 26972 25752 25027 31530 34705 35223 33471 29239 27954 27727 27314 26576 25775 24669 24480 30834 33218 33783 32546 30661 30070 29722 29075 28136 27315 26125 26057 32601 34214 35232 33565 31931 31779 31626 31230 29574 28312 27186 27397 33387 34996 36251 34284 32349 30991 29916 29067 27978 26719 25544 25703 31703 33733 35121 32714 31111 29977 30375 29323 28193 27222 26904 27952 33512 36215 36856 35341 32624 30885 31108 30267 28645 28474 25805 24756 30437 33177 33069 31342 28912 28373 28599 27884 25727 25393 23147 23164 29286 31008
 
Output produced by software:


Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time3 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135
 
Charts produced by software:
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Dec/12/t12290872622calny2fiedkc1v/1236t1229087108.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Dec/12/t12290872622calny2fiedkc1v/1236t1229087108.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Dec/12/t12290872622calny2fiedkc1v/2afm11229087108.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Dec/12/t12290872622calny2fiedkc1v/2afm11229087108.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Dec/12/t12290872622calny2fiedkc1v/3edrh1229087108.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Dec/12/t12290872622calny2fiedkc1v/3edrh1229087108.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Dec/12/t12290872622calny2fiedkc1v/4bfkn1229087108.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Dec/12/t12290872622calny2fiedkc1v/4bfkn1229087108.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Dec/12/t12290872622calny2fiedkc1v/57f9l1229087108.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Dec/12/t12290872622calny2fiedkc1v/57f9l1229087108.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Dec/12/t12290872622calny2fiedkc1v/6fsf11229087108.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Dec/12/t12290872622calny2fiedkc1v/6fsf11229087108.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Dec/12/t12290872622calny2fiedkc1v/7tft91229087108.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Dec/12/t12290872622calny2fiedkc1v/7tft91229087108.ps (open in new window)


 
Parameters (Session):
par1 = 50 ; par2 = 50 ; par3 = Y ; par4 = Y ; par5 = Werkloosheid Vlaams Gewest ; par6 = Totale werkloosheid ; par7 = Werkloosheid Waals Gewest ;
 
Parameters (R input):
par1 = 50 ; par2 = 50 ; par3 = Y ; par4 = Y ; par5 = Werkloosheid Vlaams Gewest ; par6 = Totale werkloosheid ; par7 = Werkloosheid Waals Gewest ;
 
R code (references can be found in the software module):
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()
 





Copyright

Creative Commons License

This work is licensed under a Creative Commons Attribution-Noncommercial-Share Alike 3.0 License.

Software written by Ed van Stee & Patrick Wessa


Disclaimer

Information provided on this web site is provided "AS IS" without warranty of any kind, either express or implied, including, without limitation, warranties of merchantability, fitness for a particular purpose, and noninfringement. We use reasonable efforts to include accurate and timely information and periodically update the information, and software without notice. However, we make no warranties or representations as to the accuracy or completeness of such information (or software), and we assume no liability or responsibility for errors or omissions in the content of this web site, or any software bugs in online applications. Your use of this web site is AT YOUR OWN RISK. Under no circumstances and under no legal theory shall we be liable to you or any other person for any direct, indirect, special, incidental, exemplary, or consequential damages arising from your access to, or use of, this web site.


Privacy Policy

We may request personal information to be submitted to our servers in order to be able to:

  • personalize online software applications according to your needs
  • enforce strict security rules with respect to the data that you upload (e.g. statistical data)
  • manage user sessions of online applications
  • alert you about important changes or upgrades in resources or applications

We NEVER allow other companies to directly offer registered users information about their products and services. Banner references and hyperlinks of third parties NEVER contain any personal data of the visitor.

We do NOT sell, nor transmit by any means, personal information, nor statistical data series uploaded by you to third parties.

We carefully protect your data from loss, misuse, alteration, and destruction. However, at any time, and under any circumstance you are solely responsible for managing your passwords, and keeping them secret.

We store a unique ANONYMOUS USER ID in the form of a small 'Cookie' on your computer. This allows us to track your progress when using this website which is necessary to create state-dependent features. The cookie is used for NO OTHER PURPOSE. At any time you may opt to disallow cookies from this website - this will not affect other features of this website.

We examine cookies that are used by third-parties (banner and online ads) very closely: abuse from third-parties automatically results in termination of the advertising contract without refund. We have very good reason to believe that the cookies that are produced by third parties (banner ads) do NOT cause any privacy or security risk.

FreeStatistics.org is safe. There is no need to download any software to use the applications and services contained in this website. Hence, your system's security is not compromised by their use, and your personal data - other than data you submit in the account application form, and the user-agent information that is transmitted by your browser - is never transmitted to our servers.

As a general rule, we do not log on-line behavior of individuals (other than normal logging of webserver 'hits'). However, in cases of abuse, hacking, unauthorized access, Denial of Service attacks, illegal copying, hotlinking, non-compliance with international webstandards (such as robots.txt), or any other harmful behavior, our system engineers are empowered to log, track, identify, publish, and ban misbehaving individuals - even if this leads to ban entire blocks of IP addresses, or disclosing user's identity.


FreeStatistics.org is powered by