R version 2.8.1 (2008-12-22)
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.

> x <- array(list(100.30
+ ,100.00
+ ,100.00
+ ,NA
+ ,101.90
+ ,102.83
+ ,97.07
+ ,NA
+ ,102.10
+ ,109.50
+ ,92.57
+ ,NA
+ ,103.20
+ ,115.91
+ ,87.69
+ ,NA
+ ,103.70
+ ,107.94
+ ,82.97
+ ,NA
+ ,106.20
+ ,110.86
+ ,78.80
+ ,100.00
+ ,107.70
+ ,118.89
+ ,75.16
+ ,101.42
+ ,109.90
+ ,123.38
+ ,72.07
+ ,97.84
+ ,111.70
+ ,113.33
+ ,68.06
+ ,100.00
+ ,114.90
+ ,116.38
+ ,65.13
+ ,95.64
+ ,116.00
+ ,122.04
+ ,63.94
+ ,98.94
+ ,118.30
+ ,125.47
+ ,63.29
+ ,98.39
+ ,120.40
+ ,115.62
+ ,62.91
+ ,99.98
+ ,126.00
+ ,117.91
+ ,62.96
+ ,96.14
+ ,128.10
+ ,122.40
+ ,65.24
+ ,99.33
+ ,130.10
+ ,125.05
+ ,69.96
+ ,97.82
+ ,130.80
+ ,114.18
+ ,76.14
+ ,102.85
+ ,133.60
+ ,114.74
+ ,84.06
+ ,103.93
+ ,134.20
+ ,120.63
+ ,90.56
+ ,107.88
+ ,135.50
+ ,123.68
+ ,94.69
+ ,110.82
+ ,136.20
+ ,112.84
+ ,98.54
+ ,115.53
+ ,139.10
+ ,115.64
+ ,99.57
+ ,117.09
+ ,139.00
+ ,122.32
+ ,101.52
+ ,123.40
+ ,139.60
+ ,124.59
+ ,106.07
+ ,126.57
+ ,138.70
+ ,116.33
+ ,109.00
+ ,130.34
+ ,140.90
+ ,117.45
+ ,104.83
+ ,81.79
+ ,141.30
+ ,125.64
+ ,103.96
+ ,83.26
+ ,141.80
+ ,128.38
+ ,103.04
+ ,83.65
+ ,142.00
+ ,119.87
+ ,99.19
+ ,83.88
+ ,144.50
+ ,121.22
+ ,97.89
+ ,84.73
+ ,144.60
+ ,128.98
+ ,93.82
+ ,86.27
+ ,145.50
+ ,131.35
+ ,89.48
+ ,87.71
+ ,146.80
+ ,121.35
+ ,87.74
+ ,88.52
+ ,149.50
+ ,123.72
+ ,85.79
+ ,89.32
+ ,149.90
+ ,131.06
+ ,85.47
+ ,90.61
+ ,150.10
+ ,134.55
+ ,82.00
+ ,90.95
+ ,150.90
+ ,125.93
+ ,83.51
+ ,90.90
+ ,152.80
+ ,128.90
+ ,82.38
+ ,91.02
+ ,153.10
+ ,136.19
+ ,79.88
+ ,91.66
+ ,154.00
+ ,140.34
+ ,76.03
+ ,91.23
+ ,154.90
+ ,130.48
+ ,70.44
+ ,90.81
+ ,156.90
+ ,134.68
+ ,68.71
+ ,90.90
+ ,158.40
+ ,141.05
+ ,65.08
+ ,91.27
+ ,159.70
+ ,145.44
+ ,62.09
+ ,90.45
+ ,160.20
+ ,136.21
+ ,59.82
+ ,89.60
+ ,163.20
+ ,139.85
+ ,58.35
+ ,88.49
+ ,163.70
+ ,147.13
+ ,59.38
+ ,88.31
+ ,164.40
+ ,151.44
+ ,58.57
+ ,87.83
+ ,163.70
+ ,143.62
+ ,58.68
+ ,87.62
+ ,165.50
+ ,148.55
+ ,57.48
+ ,87.60
+ ,165.60
+ ,153.54
+ ,55.69
+ ,88.08
+ ,166.80
+ ,159.79
+ ,54.39
+ ,88.08
+ ,167.50
+ ,152.55
+ ,53.52
+ ,87.41
+ ,170.60
+ ,155.84
+ ,52.01
+ ,87.53
+ ,170.90
+ ,160.38
+ ,49.51
+ ,88.31
+ ,172.00
+ ,164.22
+ ,49.46
+ ,89.04
+ ,171.80
+ ,156.40
+ ,48.21
+ ,89.34
+ ,173.90
+ ,160.05
+ ,47.83
+ ,90.22
+ ,174.00
+ ,165.60
+ ,48.86
+ ,91.02
+ ,173.80
+ ,171.15
+ ,49.51
+ ,91.69
+ ,173.90
+ ,161.90
+ ,49.95
+ ,91.62
+ ,176.00
+ ,167.21
+ ,49.57
+ ,91.32
+ ,176.60
+ ,171.34
+ ,51.52
+ ,91.82
+ ,178.20
+ ,176.83
+ ,48.48
+ ,91.75
+ ,179.20
+ ,166.27
+ ,50.43
+ ,91.57
+ ,181.30
+ ,172.30
+ ,48.26
+ ,92.19
+ ,181.80
+ ,176.71
+ ,48.70
+ ,92.81
+ ,182.90
+ ,182.99
+ ,47.72
+ ,50.99
+ ,183.80
+ ,172.07
+ ,45.55
+ ,50.92
+ ,186.30
+ ,178.17
+ ,45.82
+ ,50.37
+ ,187.40
+ ,182.20
+ ,44.41
+ ,50.32
+ ,189.20
+ ,188.49
+ ,45.44
+ ,49.77
+ ,189.70
+ ,176.88
+ ,45.17
+ ,49.45
+ ,191.90
+ ,182.13
+ ,45.39
+ ,49.15
+ ,192.60
+ ,185.32
+ ,46.10
+ ,49.10
+ ,193.70
+ ,192.86
+ ,50.27
+ ,48.85
+ ,194.20
+ ,180.27
+ ,50.54
+ ,48.76
+ ,197.60
+ ,184.92
+ ,53.04
+ ,48.58
+ ,199.30
+ ,187.82
+ ,53.04
+ ,48.90
+ ,201.40
+ ,194.94
+ ,52.82
+ ,48.97
+ ,203.00
+ ,184.36
+ ,52.98
+ ,48.78
+ ,206.30
+ ,188.80
+ ,51.52
+ ,48.64
+ ,207.10
+ ,193.42
+ ,50.70
+ ,48.64
+ ,209.80
+ ,199.76
+ ,49.51
+ ,48.16
+ ,211.10
+ ,188.78
+ ,51.08
+ ,48.12
+ ,215.30
+ ,191.49
+ ,53.69
+ ,48.05
+ ,217.40
+ ,194.87
+ ,58.30
+ ,48.25
+ ,215.50
+ ,198.28
+ ,64.05
+ ,47.77
+ ,210.90
+ ,183.24
+ ,72.45
+ ,46.72
+ ,212.60
+ ,204.87
+ ,80.91
+ ,NA)
+ ,dim=c(4
+ ,90)
+ ,dimnames=list(c('RP'
+ ,'FCE'
+ ,'UNEM'
+ ,'IS')
+ ,1:90))
>  y <- array(NA,dim=c(4,90),dimnames=list(c('RP','FCE','UNEM','IS'),1:90))
>  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/rcomp/tmp/1vyvs1257972987.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/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/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 = 12.2229, p-value < 2.2e-16
alternative hypothesis: true tau is not equal to 0 
sample estimates:
      tau 
0.8758742 

> 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/rcomp/tmp/2ywyd1257972987.tab") 
> 
> system("convert tmp/1vyvs1257972987.ps tmp/1vyvs1257972987.png")
> 
> 
> proc.time()
   user  system elapsed 
  0.480   0.320   0.815