R version 2.8.0 (2008-10-20)
Copyright (C) 2008 The R Foundation for Statistical Computing
ISBN 3-900051-07-0

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> y <- c(18562.2,23184.8,22925.8,21490.2,23243.2,21688.7,21423.1,21211.2,17877.4,20664.3,22160.0,19813.6,17735.4,19640.2,20844.4,19823.1,18594.6,21350.6,18574.1,18924.2,17343.4,19961.2,19932.1,19464.6,16165.4,17574.9,19795.4,19439.5,17170.0,21072.4,17751.8,17515.5,18040.3,19090.1,17746.5,19202.1,15141.6,16258.1,18586.5,17209.4,17838.7,19123.5,16583.6,15991.2,16704.4,17420.4,17872.0,17823.2,13866.5,15912.8,17870.5,15420.3,16379.4,17903.9,15305.8,14583.3,14861.0,14968.9,16726.5,16283.6,11703.7,15101.8,15469.7,14956.9,15370.6,15998.1,14725.1,14768.9,13659.6,15070.3,16943.0,15761.3,12083.0,15023.6,15106.5,15498.6,15258.9,15859.4,14205.3,14291.1,12875.1,14755.3,15873.2,14803.0,12520.2,14568.3,15644.8,15454.9,14254.0,16754.8,14944.2,15044.5)
> x <- c(103.2,106.9,118.6,111.9,119.0,116.3,116.2,115.6,107.5,116.8,124.1,115.0,105.3,104.1,119.6,112.2,109.6,121.8,108.8,111.5,108.5,115.4,119.3,116.5,101.9,96.6,116.6,112.5,103.7,118.5,105.0,105.0,109.1,112.6,108.2,114.4,95.9,90.8,115.0,99.8,103.0,108.4,99.2,100.6,107.1,107.0,111.9,115.6,97.7,97.3,111.8,99.3,104.6,113.3,98.7,98.2,102.5,100.8,111.4,108.9,90.4,94.6,104.3,99.0,103.0,105.1,98.9,101.0,96.5,102.8,112.3,106.5,91.9,92.7,104.2,101.7,102.8,105.6,96.9,97.6,93.7,102.1,106.6,100.2,92.0,86.8,104.8,100.0,96.8,110.6,100.7,101.5)
> ylab = 'Totale Uitvoer van Belgie'
> xlab = 'Totale Productie (zonder bouwnijverheid)'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> library('Kendall')
> k <- Kendall(x,y)
> postscript(file="/var/www/html/freestat/rcomp/tmp/1yt7q1228679263.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> 
> plot(x,y,main='Scatterplot',xlab=xlab,ylab=ylab)
> grid()
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/24ep61228679263.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> 
> plot(rank(x),rank(y),main='Scatterplot of Ranks',xlab=xlab,ylab=ylab)
> grid()
> dev.off()
null device 
          1 
> 
> #Note: the /var/www/html/freestat/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/html/freestat/rcomp/createtable")
> 
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,'Kendall tau Rank Correlation',2,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Kendall tau',header=TRUE)
> a<-table.element(a,k$tau)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'2-sided p-value',header=TRUE)
> a<-table.element(a,k$sl)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Score',header=TRUE)
> a<-table.element(a,k$S)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Var(Score)',header=TRUE)
> a<-table.element(a,k$varS)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Denominator',header=TRUE)
> a<-table.element(a,k$D)
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/html/freestat/rcomp/tmp/3tgrc1228679263.tab") 
> 
> system("convert tmp/1yt7q1228679263.ps tmp/1yt7q1228679263.png")
> system("convert tmp/24ep61228679263.ps tmp/24ep61228679263.png")
> 
> 
> proc.time()
   user  system elapsed 
  0.862   0.523   0.994