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(100.00
+ ,100.00
+ ,100.00
+ ,100.00
+ ,100.00
+ ,98.04
+ ,100.78
+ ,102.58
+ ,97.56
+ ,95.10
+ ,101.24
+ ,84.79
+ ,95.12
+ ,92.16
+ ,102.35
+ ,74.23
+ ,92.68
+ ,90.20
+ ,102.41
+ ,82.99
+ ,92.68
+ ,89.22
+ ,102.67
+ ,93.30
+ ,97.56
+ ,94.12
+ ,102.54
+ ,98.45
+ ,107.32
+ ,100.00
+ ,102.74
+ ,91.24
+ ,112.20
+ ,100.00
+ ,103.07
+ ,65.21
+ ,112.20
+ ,98.04
+ ,103.39
+ ,57.22
+ ,112.20
+ ,97.06
+ ,103.72
+ ,73.45
+ ,114.63
+ ,97.06
+ ,103.91
+ ,71.65
+ ,117.07
+ ,97.06
+ ,104.37
+ ,58.76
+ ,117.07
+ ,95.10
+ ,104.89
+ ,58.25
+ ,114.63
+ ,93.14
+ ,105.61
+ ,69.85
+ ,114.63
+ ,92.16
+ ,106.00
+ ,71.39
+ ,114.63
+ ,91.18
+ ,106.13
+ ,71.39
+ ,112.20
+ ,91.18
+ ,106.20
+ ,68.04
+ ,121.95
+ ,97.06
+ ,106.26
+ ,65.98
+ ,131.71
+ ,102.94
+ ,106.33
+ ,53.35
+ ,134.15
+ ,103.92
+ ,106.98
+ ,59.79
+ ,136.59
+ ,103.92
+ ,107.44
+ ,55.67
+ ,136.59
+ ,102.94
+ ,107.50
+ ,57.47
+ ,141.46
+ ,103.92
+ ,107.57
+ ,61.86
+ ,146.34
+ ,105.88
+ ,107.63
+ ,73.20
+ ,148.78
+ ,105.88
+ ,108.15
+ ,71.39
+ ,148.78
+ ,104.90
+ ,108.35
+ ,75.52
+ ,146.34
+ ,103.92
+ ,109.07
+ ,75.00
+ ,146.34
+ ,103.92
+ ,109.39
+ ,69.33
+ ,148.78
+ ,105.88
+ ,109.78
+ ,61.34
+ ,158.54
+ ,111.76
+ ,109.98
+ ,66.49
+ ,173.17
+ ,119.61
+ ,110.18
+ ,82.22
+ ,180.49
+ ,121.57
+ ,110.31
+ ,72.68
+ ,180.49
+ ,121.57
+ ,110.57
+ ,70.10
+ ,182.93
+ ,120.59
+ ,110.63
+ ,65.21
+ ,185.37
+ ,121.57
+ ,110.70
+ ,69.59
+ ,190.24
+ ,122.55
+ ,110.76
+ ,62.37
+ ,190.24
+ ,122.55
+ ,111.15
+ ,64.43
+ ,187.80
+ ,121.57
+ ,111.48
+ ,59.54
+ ,185.37
+ ,120.59
+ ,112.13
+ ,62.11
+ ,182.93
+ ,119.61
+ ,112.13
+ ,65.98
+ ,178.05
+ ,118.63
+ ,112.20
+ ,71.13
+ ,185.37
+ ,123.53
+ ,112.20
+ ,69.85
+ ,195.12
+ ,129.41
+ ,112.46
+ ,62.89
+ ,195.12
+ ,131.37
+ ,112.85
+ ,63.40
+ ,192.68
+ ,129.41
+ ,113.31
+ ,54.64
+ ,190.24
+ ,126.47
+ ,113.37
+ ,51.29
+ ,187.80
+ ,125.49
+ ,113.37
+ ,47.94
+ ,190.24
+ ,124.51
+ ,113.44
+ ,48.45
+ ,187.80
+ ,123.53
+ ,113.89
+ ,46.91
+ ,182.93
+ ,121.57
+ ,114.16
+ ,44.85
+ ,178.05
+ ,118.63
+ ,114.74
+ ,44.07
+ ,173.17
+ ,117.65
+ ,114.81
+ ,35.57
+ ,170.73
+ ,116.67
+ ,114.22
+ ,32.73
+ ,178.05
+ ,122.55
+ ,114.55
+ ,30.67
+ ,190.24
+ ,129.41
+ ,114.74
+ ,32.99
+ ,192.68
+ ,131.37
+ ,115.26
+ ,30.67
+ ,192.68
+ ,130.39
+ ,114.87
+ ,31.44
+ ,190.24
+ ,127.45
+ ,114.87
+ ,37.89
+ ,190.24
+ ,126.47
+ ,114.94
+ ,37.63
+ ,192.68
+ ,127.45
+ ,115.00
+ ,50.52
+ ,190.24
+ ,126.47
+ ,115.98
+ ,48.45
+ ,185.37
+ ,123.53
+ ,116.44
+ ,52.32
+ ,180.49
+ ,121.57
+ ,117.03
+ ,52.58
+ ,175.61
+ ,118.63
+ ,117.09
+ ,48.97
+ ,168.29
+ ,116.67
+ ,117.16
+ ,46.39
+ ,173.17
+ ,120.59
+ ,117.22
+ ,49.48
+ ,182.93
+ ,127.45
+ ,117.22
+ ,49.48
+ ,185.37
+ ,127.45
+ ,117.29
+ ,50.77
+ ,180.49
+ ,123.53
+ ,117.35
+ ,63.40
+ ,178.05
+ ,119.61
+ ,117.55
+ ,60.82
+ ,175.61
+ ,118.63
+ ,117.68
+ ,65.21
+ ,178.05
+ ,117.65
+ ,117.68
+ ,59.54
+ ,175.61
+ ,115.69
+ ,118.26
+ ,51.03
+ ,173.17
+ ,113.73
+ ,118.66
+ ,37.63
+ ,170.73
+ ,111.76
+ ,119.05
+ ,32.47
+ ,168.29
+ ,109.80
+ ,119.18
+ ,40.72
+ ,165.85
+ ,109.80
+ ,119.44
+ ,44.85
+ ,175.61
+ ,115.69
+ ,119.77
+ ,48.71
+ ,185.37
+ ,122.55
+ ,119.83
+ ,47.68
+ ,187.80
+ ,123.53
+ ,119.90
+ ,41.75
+ ,185.37
+ ,121.57
+ ,119.96
+ ,33.51
+ ,182.93
+ ,118.63
+ ,120.09
+ ,36.60
+ ,182.93
+ ,117.65
+ ,120.29
+ ,29.64
+ ,185.37
+ ,117.65
+ ,120.35
+ ,10.82
+ ,185.37
+ ,116.67
+ ,121.27
+ ,19.07
+ ,185.37
+ ,115.69
+ ,121.72
+ ,26.29
+ ,182.93
+ ,112.75
+ ,122.37
+ ,38.92
+ ,178.05
+ ,110.78
+ ,122.50
+ ,47.94
+ ,175.61
+ ,109.80
+ ,122.57
+ ,40.98
+ ,180.49
+ ,113.73
+ ,122.64
+ ,26.55
+ ,195.12
+ ,119.61
+ ,122.83
+ ,11.34
+ ,200.00
+ ,119.61
+ ,122.90
+ ,21.13
+ ,195.12
+ ,114.71
+ ,122.96
+ ,22.16
+ ,187.80
+ ,109.80
+ ,122.96
+ ,14.95
+ ,187.80
+ ,107.84
+ ,123.03
+ ,15.21
+ ,190.24
+ ,106.86
+ ,123.03
+ ,24.48
+ ,190.24
+ ,105.88
+ ,123.55
+ ,25.26
+ ,187.80
+ ,102.94
+ ,123.94
+ ,31.70
+ ,182.93
+ ,100.00
+ ,125.18
+ ,30.15
+ ,178.05
+ ,98.04
+ ,125.57
+ ,21.65
+ ,173.17
+ ,97.06
+ ,125.90
+ ,19.07
+ ,173.17
+ ,100.98
+ ,126.22
+ ,16.75
+ ,175.61
+ ,104.90
+ ,126.48
+ ,23.45
+ ,165.85
+ ,101.96
+ ,126.68
+ ,30.67
+ ,160.98
+ ,99.02
+ ,127.14
+ ,33.51
+ ,156.10
+ ,95.10
+ ,127.14
+ ,39.43
+ ,156.10
+ ,92.16
+ ,127.14
+ ,50.00
+ ,158.54
+ ,87.25
+ ,127.14
+ ,46.13
+ ,153.66
+ ,82.35
+ ,127.53
+ ,50.26
+ ,143.90
+ ,79.41
+ ,127.85
+ ,58.25
+ ,134.15
+ ,81.37
+ ,127.98
+ ,52.58
+ ,126.83
+ ,79.41
+ ,127.98
+ ,55.67
+ ,119.51
+ ,78.43
+ ,127.98
+ ,70.88
+ ,131.71
+ ,85.29
+ ,127.98
+ ,71.91
+ ,141.46
+ ,90.20
+ ,128.51
+ ,74.23
+ ,139.02
+ ,88.24
+ ,128.96
+ ,86.60
+ ,136.59
+ ,87.25
+ ,129.16
+ ,76.55
+ ,134.15
+ ,83.33
+ ,129.29
+ ,79.90
+ ,131.71
+ ,79.41
+ ,129.48
+ ,64.18
+ ,131.71
+ ,73.53
+ ,129.55
+ ,56.70
+ ,131.71
+ ,69.61
+ ,130.14
+ ,57.99
+ ,134.15
+ ,67.65
+ ,130.46
+ ,53.87
+ ,141.46
+ ,69.61
+ ,131.31
+ ,71.91
+ ,139.02
+ ,68.63
+ ,131.90
+ ,80.93
+ ,131.71
+ ,65.69
+ ,132.35
+ ,75.52
+ ,136.59
+ ,68.63
+ ,132.75
+ ,68.30
+ ,141.46
+ ,71.57
+ ,132.75
+ ,68.81
+ ,151.22
+ ,75.49
+ ,133.07
+ ,58.25
+ ,165.85
+ ,82.35
+ ,133.14
+ ,60.57
+ ,163.41
+ ,82.35
+ ,133.27
+ ,54.90
+ ,163.41
+ ,86.27
+ ,133.40
+ ,56.19
+ ,156.10
+ ,89.22
+ ,133.59
+ ,74.74
+ ,153.66
+ ,88.24
+ ,133.79
+ ,67.78
+ ,153.66
+ ,84.31
+ ,134.18
+ ,68.81
+ ,156.10
+ ,77.45
+ ,134.70
+ ,46.65
+ ,153.66
+ ,75.49
+ ,134.96
+ ,34.28
+ ,146.34
+ ,76.47
+ ,135.09
+ ,22.68
+ ,153.66
+ ,90.20
+ ,135.55
+ ,32.99
+ ,153.66
+ ,92.16
+ ,135.68
+ ,32.47
+ ,160.98
+ ,90.20
+ ,136.01
+ ,32.47
+ ,182.93
+ ,85.29
+ ,136.07
+ ,33.25
+ ,190.24
+ ,82.35
+ ,136.33
+ ,28.35
+ ,192.68
+ ,84.31
+ ,136.40
+ ,35.31
+ ,190.24
+ ,88.24
+ ,136.79
+ ,31.19
+ ,185.37
+ ,89.22
+ ,136.86
+ ,44.85
+ ,182.93
+ ,85.29
+ ,136.92
+ ,45.36
+ ,185.37
+ ,80.39
+ ,136.99
+ ,38.14
+ ,182.93
+ ,77.45
+ ,137.51
+ ,26.80
+ ,178.05
+ ,77.45
+ ,137.90
+ ,41.75
+ ,185.37
+ ,89.22
+ ,138.10
+ ,38.40
+ ,182.93
+ ,92.16
+ ,138.29
+ ,46.13
+ ,185.37
+ ,92.16
+ ,138.42
+ ,46.39
+ ,192.68
+ ,89.22
+ ,138.55
+ ,40.72
+ ,192.68
+ ,88.24
+ ,138.88
+ ,47.94
+ ,197.56
+ ,91.18
+ ,139.20
+ ,44.85
+ ,200.00
+ ,97.06
+ ,139.40
+ ,40.98
+ ,195.12
+ ,96.08
+ ,139.60
+ ,32.47
+ ,182.93
+ ,91.18
+ ,139.79
+ ,29.12
+ ,165.85
+ ,81.37
+ ,140.12
+ ,49.48
+ ,158.54
+ ,78.43
+ ,140.25
+ ,67.27
+ ,160.98
+ ,83.33
+ ,140.83
+ ,58.25
+ ,185.37
+ ,101.96
+ ,141.16
+ ,62.11
+ ,195.12
+ ,108.82
+ ,141.49
+ ,58.25
+ ,197.56
+ ,106.86
+ ,141.68
+ ,52.32
+ ,187.80
+ ,98.04
+ ,141.88
+ ,73.71
+ ,182.93
+ ,90.20
+ ,142.14
+ ,65.72
+ ,185.37
+ ,90.20
+ ,142.27
+ ,58.51
+ ,190.24
+ ,93.14
+ ,142.60
+ ,58.25
+ ,190.24
+ ,94.12
+ ,142.79
+ ,66.24
+ ,190.24
+ ,93.14
+ ,143.05
+ ,79.12
+ ,182.93
+ ,89.22
+ ,143.77
+ ,71.13
+ ,182.93
+ ,87.25
+ ,144.10
+ ,64.69
+ ,173.17
+ ,88.24
+ ,144.16
+ ,73.97
+ ,182.93
+ ,99.02
+ ,144.68
+ ,80.93
+ ,182.93
+ ,100.98
+ ,144.81
+ ,80.15
+ ,185.37
+ ,100.00
+ ,144.94
+ ,81.44
+ ,187.80
+ ,94.12
+ ,145.14
+ ,63.66
+ ,187.80
+ ,90.20
+ ,145.40
+ ,66.24
+ ,192.68
+ ,91.18
+ ,145.53
+ ,74.48
+ ,197.56
+ ,92.16
+ ,145.73
+ ,67.78
+ ,200.00
+ ,92.16
+ ,146.12
+ ,61.34
+ ,200.00
+ ,90.20
+ ,146.84
+ ,43.56
+ ,200.00
+ ,88.24
+ ,147.10
+ ,50.52
+ ,192.68
+ ,88.24
+ ,147.36
+ ,56.44
+ ,178.05
+ ,88.24
+ ,147.62
+ ,48.20
+ ,168.29
+ ,96.08
+ ,147.75
+ ,41.24
+ ,160.98
+ ,98.04
+ ,148.08
+ ,42.01
+ ,163.41
+ ,96.08
+ ,148.27
+ ,31.44
+ ,168.29
+ ,91.18
+ ,148.60
+ ,31.19
+ ,170.73
+ ,88.24
+ ,148.79
+ ,38.40
+ ,173.17
+ ,88.24
+ ,148.99
+ ,42.27
+ ,175.61
+ ,89.22
+ ,149.05
+ ,42.78
+ ,173.17
+ ,89.22
+ ,149.25
+ ,45.62
+ ,168.29
+ ,89.22
+ ,149.38
+ ,46.91
+ ,170.73
+ ,90.20
+ ,149.45
+ ,45.88
+ ,165.85
+ ,86.27
+ ,149.58
+ ,32.99
+ ,156.10
+ ,81.37
+ ,149.77
+ ,33.25
+ ,163.41
+ ,82.35
+ ,149.97
+ ,35.31
+ ,160.98
+ ,79.41
+ ,150.03
+ ,28.87
+ ,156.10
+ ,75.49
+ ,150.16
+ ,38.92
+ ,153.66
+ ,77.45
+ ,150.55
+ ,57.73
+ ,151.22
+ ,77.45
+ ,150.68
+ ,75.77
+ ,158.54
+ ,78.43
+ ,151.14
+ ,79.64
+ ,165.85
+ ,77.45
+ ,151.27
+ ,89.18
+ ,165.85
+ ,74.51
+ ,151.99
+ ,93.81
+ ,156.10
+ ,69.61
+ ,153.36
+ ,113.14
+ ,148.78
+ ,66.67
+ ,153.95
+ ,106.96
+ ,141.46
+ ,63.73
+ ,154.53
+ ,134.28
+ ,148.78
+ ,67.65
+ ,154.66
+ ,149.48
+ ,175.61
+ ,80.39
+ ,154.92
+ ,152.32
+ ,178.05
+ ,85.29
+ ,155.38
+ ,138.92
+ ,168.29
+ ,81.37
+ ,155.84
+ ,140.72
+ ,148.78
+ ,77.45
+ ,155.97
+ ,121.65
+ ,141.46
+ ,73.53
+ ,156.56
+ ,80.93
+ ,151.22
+ ,76.47
+ ,156.69
+ ,67.78
+ ,173.17
+ ,81.37
+ ,156.88
+ ,59.79
+ ,187.80
+ ,82.35
+ ,157.01
+ ,49.74
+ ,192.68
+ ,80.39
+ ,157.27
+ ,15.98
+ ,187.80
+ ,75.49
+ ,157.47
+ ,15.46
+ ,180.49
+ ,70.59
+ ,157.99
+ ,-9.54
+ ,182.93
+ ,71.57
+ ,158.45
+ ,-28.35
+ ,195.12
+ ,79.41
+ ,158.64
+ ,-43.30
+ ,197.56
+ ,83.33
+ ,159.10
+ ,-20.10)
+ ,dim=c(4
+ ,224)
+ ,dimnames=list(c('Belg_Werkl_Men'
+ ,'Belg_Werkl_Women'
+ ,'Belg_Uurl'
+ ,'Belg_Infl')
+ ,1:224))
>  y <- array(NA,dim=c(4,224),dimnames=list(c('Belg_Werkl_Men','Belg_Werkl_Women','Belg_Uurl','Belg_Infl'),1:224))
>  for (i in 1:dim(x)[1])
+  {
+  	for (j in 1:dim(x)[2])
+  	{
+  		y[i,j] <- as.numeric(x[i,j])
+  	}
+  }
> par20 = ''
> par19 = ''
> par18 = ''
> par17 = ''
> par16 = ''
> par15 = ''
> par14 = ''
> par13 = ''
> par12 = ''
> par11 = ''
> par10 = ''
> par9 = ''
> par8 = ''
> par7 = ''
> par6 = ''
> par5 = ''
> par4 = ''
> par3 = ''
> par2 = ''
> par1 = ''
> ylab = ''
> xlab = ''
> main = ''
> #'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/1m14u1257934818.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 = 6.4847, p-value = 8.893e-11
alternative hypothesis: true tau is not equal to 0 
sample estimates:
      tau 
0.2984354 

> 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/2jieu1257934818.tab") 
> 
> system("convert tmp/1m14u1257934818.ps tmp/1m14u1257934818.png")
> 
> 
> proc.time()
   user  system elapsed 
  0.785   0.300   0.889