R version 2.9.0 (2009-04-17)
Copyright (C) 2009 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.

> x <- array(list(1.01135
+ ,0.99108
+ ,0.96782
+ ,1.41683
+ ,1.01747
+ ,1.01106
+ ,0.93564
+ ,1.41879
+ ,1.03066
+ ,1.02303
+ ,0.88861
+ ,1.42140
+ ,1.04183
+ ,1.03884
+ ,0.87376
+ ,1.42270
+ ,1.05315
+ ,1.04008
+ ,0.86386
+ ,1.42596
+ ,1.06308
+ ,1.04561
+ ,0.90842
+ ,1.42792
+ ,1.07429
+ ,1.03687
+ ,0.87624
+ ,1.43053
+ ,1.08626
+ ,1.06073
+ ,0.83168
+ ,1.43770
+ ,1.10025
+ ,1.10986
+ ,0.78960
+ ,1.44097
+ ,1.12714
+ ,1.14308
+ ,0.79703
+ ,1.44162
+ ,1.14231
+ ,1.13704
+ ,0.80198
+ ,1.44684
+ ,1.14965
+ ,1.11881
+ ,0.76238
+ ,1.44814
+ ,1.16536
+ ,1.10102
+ ,0.80446
+ ,1.44945
+ ,1.18178
+ ,1.10427
+ ,0.85396
+ ,1.45140
+ ,1.19467
+ ,1.21031
+ ,0.83168
+ ,1.45401
+ ,1.21240
+ ,1.21034
+ ,0.82426
+ ,1.45532
+ ,1.22326
+ ,1.18008
+ ,0.86386
+ ,1.45727
+ ,1.23535
+ ,1.17951
+ ,0.90842
+ ,1.46119
+ ,1.25060
+ ,1.15236
+ ,0.97772
+ ,1.46836
+ ,1.26165
+ ,1.18474
+ ,0.99257
+ ,1.47097
+ ,1.27666
+ ,1.24201
+ ,0.99257
+ ,1.47358
+ ,1.29593
+ ,1.27613
+ ,1.00000
+ ,1.47619
+ ,1.30998
+ ,1.22767
+ ,0.97030
+ ,1.47750
+ ,1.32252
+ ,1.20806
+ ,0.94059
+ ,1.48076
+ ,1.33510
+ ,1.19442
+ ,0.95050
+ ,1.48271
+ ,1.34013
+ ,1.17686
+ ,0.93317
+ ,1.48598
+ ,1.35081
+ ,1.18862
+ ,0.94802
+ ,1.48793
+ ,1.36653
+ ,1.18155
+ ,1.00495
+ ,1.48989
+ ,1.37428
+ ,1.18306
+ ,1.01485
+ ,1.49054
+ ,1.38213
+ ,1.20374
+ ,0.99257
+ ,1.49250
+ ,1.39470
+ ,1.22945
+ ,1.04950
+ ,1.49380
+ ,1.40414
+ ,1.24972
+ ,1.07921
+ ,1.49446
+ ,1.41745
+ ,1.29155
+ ,1.14851
+ ,1.49576
+ ,1.43398
+ ,1.30633
+ ,1.14356
+ ,1.49772
+ ,1.39402
+ ,1.22020
+ ,1.09901
+ ,1.49967
+ ,1.40632
+ ,1.17153
+ ,1.08911
+ ,1.50033
+ ,1.41982
+ ,1.17140
+ ,1.09406
+ ,1.50163
+ ,1.43532
+ ,1.09528
+ ,1.06188
+ ,1.50554
+ ,1.44687
+ ,1.16604
+ ,1.09406
+ ,1.50685
+ ,1.46064
+ ,1.14549
+ ,1.05198
+ ,1.51142
+ ,1.47084
+ ,1.18018
+ ,1.02723
+ ,1.51272
+ ,1.48212
+ ,1.16872
+ ,1.02475
+ ,1.51990
+ ,1.49316
+ ,1.18306
+ ,1.08663
+ ,1.53359
+ ,1.50550
+ ,1.18789
+ ,1.12129
+ ,1.53947
+ ,1.51772
+ ,1.24369
+ ,1.19802
+ ,1.54534
+ ,1.31869
+ ,1.22060
+ ,1.20297
+ ,1.54664
+ ,1.33392
+ ,1.15531
+ ,1.13119
+ ,1.54925
+ ,1.34300
+ ,1.13510
+ ,1.13119
+ ,1.55382
+ ,1.35538
+ ,1.15243
+ ,1.09901
+ ,1.55838
+ ,1.21333
+ ,1.12756
+ ,1.03960
+ ,1.55969
+ ,1.18188
+ ,1.16054
+ ,0.95050
+ ,1.56556
+ ,1.10824
+ ,1.10329
+ ,1.00000
+ ,1.56686
+ ,1.04660
+ ,1.18980
+ ,1.00495
+ ,1.56882
+ ,1.05462
+ ,1.21024
+ ,0.99257
+ ,1.57012
+ ,1.06373
+ ,1.24295
+ ,0.96535
+ ,1.57273
+ ,0.95098
+ ,1.26390
+ ,0.97277
+ ,1.57469
+ ,0.95915
+ ,1.35493
+ ,0.99752
+ ,1.57991
+ ,0.97592
+ ,1.33673
+ ,0.93564
+ ,1.58447
+ ,0.99204
+ ,1.29075
+ ,0.89604
+ ,1.58643
+ ,1.00234
+ ,1.29929
+ ,0.89109
+ ,1.59100)
+ ,dim=c(4
+ ,60)
+ ,dimnames=list(c('leningbgez'
+ ,'saldozichtrek'
+ ,'rentevoet'
+ ,'buurloon')
+ ,1:60))
>  y <- array(NA,dim=c(4,60),dimnames=list(c('leningbgez','saldozichtrek','rentevoet','buurloon'),1:60))
>  for (i in 1:dim(x)[1])
+  {
+  	for (j in 1:dim(x)[2])
+  	{
+  		y[i,j] <- as.numeric(x[i,j])
+  	}
+  }
> main = 'Kendall tau Correlation Plot'
> #'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!)
> panel.tau <- function(x, y, digits=2, prefix='', cex.cor)
+ {
+ usr <- par('usr'); on.exit(par(usr))
+ par(usr = c(0, 1, 0, 1))
+ rr <- cor.test(x, y, method='kendall')
+ r <- round(rr$p.value,2)
+ txt <- format(c(r, 0.123456789), digits=digits)[1]
+ txt <- paste(prefix, txt, sep='')
+ if(missing(cex.cor)) cex <- 0.5/strwidth(txt)
+ text(0.5, 0.5, txt, cex = cex)
+ }
> 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='grey', ...)
+ }
> postscript(file="/var/www/html/rcomp/tmp/11a4y1258107499.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> pairs(t(y),diag.panel=panel.hist, upper.panel=panel.smooth, lower.panel=panel.tau, main=main)
> dev.off()
null device 
          1 
> 
> #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/html/rcomp/createtable")
> 
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,'Kendall tau rank correlations for all pairs of data series',3,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'pair',1,TRUE)
> a<-table.element(a,'tau',1,TRUE)
> a<-table.element(a,'p-value',1,TRUE)
> a<-table.row.end(a)
> n <- length(y[,1])
> n
[1] 4
> cor.test(y[1,],y[2,],method='kendall')

	Kendall's rank correlation tau

data:  y[1, ] and y[2, ] 
z = 1.9198, p-value = 0.05488
alternative hypothesis: true tau is not equal to 0 
sample estimates:
      tau 
0.1701046 

> for (i in 1:(n-1))
+ {
+ for (j in (i+1):n)
+ {
+ a<-table.row.start(a)
+ dum <- paste('tau(',dimnames(t(x))[[2]][i])
+ dum <- paste(dum,',')
+ dum <- paste(dum,dimnames(t(x))[[2]][j])
+ dum <- paste(dum,')')
+ a<-table.element(a,dum,header=TRUE)
+ r <- cor.test(y[i,],y[j,],method='kendall')
+ a<-table.element(a,r$estimate)
+ a<-table.element(a,r$p.value)
+ a<-table.row.end(a)
+ }
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/2clyk1258107499.tab") 
> 
> system("convert tmp/11a4y1258107499.ps tmp/11a4y1258107499.png")
> 
> 
> proc.time()
   user  system elapsed 
  0.483   0.224   0.591