R version 2.6.0 (2007-10-03)
Copyright (C) 2007 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.065757
+ ,0.938300
+ ,100.600000
+ ,100.600000
+ ,1.084952
+ ,0.921700
+ ,100.600000
+ ,100.600000
+ ,1.099505
+ ,0.909500
+ ,100.600000
+ ,100.600000
+ ,1.121076
+ ,0.892000
+ ,100.600000
+ ,100.600000
+ ,1.143903
+ ,0.874200
+ ,100.600000
+ ,100.600000
+ ,1.172058
+ ,0.853200
+ ,101.300000
+ ,101.300000
+ ,1.161845
+ ,0.860700
+ ,99.300000
+ ,99.300000
+ ,1.110494
+ ,0.900500
+ ,99.200000
+ ,99.200000
+ ,1.097574
+ ,0.911100
+ ,99.200000
+ ,99.200000
+ ,1.103875
+ ,0.905900
+ ,99.600000
+ ,99.200000
+ ,1.125746
+ ,0.888300
+ ,103.500000
+ ,99.200000
+ ,1.120574
+ ,0.892400
+ ,104.900000
+ ,99.200000
+ ,1.132118
+ ,0.883300
+ ,105.300000
+ ,99.600000
+ ,1.149425
+ ,0.870000
+ ,107.700000
+ ,101.600000
+ ,1.141813
+ ,0.875800
+ ,109.200000
+ ,101.600000
+ ,1.128923
+ ,0.885800
+ ,108.000000
+ ,101.600000
+ ,1.090513
+ ,0.917000
+ ,108.400000
+ ,101.600000
+ ,1.046682
+ ,0.955400
+ ,108.400000
+ ,101.600000
+ ,1.007861
+ ,0.992200
+ ,106.800000
+ ,101.600000
+ ,1.022704
+ ,0.977800
+ ,109.100000
+ ,103.300000
+ ,1.019576
+ ,0.980800
+ ,114.100000
+ ,103.300000
+ ,1.019264
+ ,0.981100
+ ,118.900000
+ ,103.300000
+ ,0.998602
+ ,1.001400
+ ,112.400000
+ ,103.300000
+ ,0.982029
+ ,1.018300
+ ,112.400000
+ ,103.300000
+ ,0.941442
+ ,1.062200
+ ,109.700000
+ ,99.600000
+ ,0.928247
+ ,1.077300
+ ,111.100000
+ ,99.600000
+ ,0.925326
+ ,1.080700
+ ,108.600000
+ ,99.600000
+ ,0.921829
+ ,1.084800
+ ,107.200000
+ ,99.600000
+ ,0.863409
+ ,1.158200
+ ,105.200000
+ ,99.600000
+ ,0.857412
+ ,1.166300
+ ,105.200000
+ ,99.600000
+ ,0.879353
+ ,1.137200
+ ,104.900000
+ ,101.200000
+ ,0.897747
+ ,1.113900
+ ,106.100000
+ ,101.200000
+ ,0.891107
+ ,1.122200
+ ,106.400000
+ ,101.200000
+ ,0.855286
+ ,1.169200
+ ,103.400000
+ ,101.200000
+ ,0.854555
+ ,1.170200
+ ,105.400000
+ ,101.200000
+ ,0.813935
+ ,1.228600
+ ,106.200000
+ ,101.200000
+ ,0.792833
+ ,1.261300
+ ,104.200000
+ ,97.700000
+ ,0.790764
+ ,1.264600
+ ,102.800000
+ ,97.700000
+ ,0.815528
+ ,1.226200
+ ,102.700000
+ ,97.700000
+ ,0.834376
+ ,1.198500
+ ,103.700000
+ ,99.900000
+ ,0.832848
+ ,1.200700
+ ,103.100000
+ ,99.900000
+ ,0.823859
+ ,1.213800
+ ,103.300000
+ ,99.900000
+ ,0.815262
+ ,1.226600
+ ,105.100000
+ ,99.900000
+ ,0.821288
+ ,1.217600
+ ,108.200000
+ ,99.900000
+ ,0.818465
+ ,1.221800
+ ,105.700000
+ ,99.900000
+ ,0.800641
+ ,1.249000
+ ,105.900000
+ ,99.900000
+ ,0.769764
+ ,1.299100
+ ,108.500000
+ ,99.900000
+ ,0.745823
+ ,1.340800
+ ,107.200000
+ ,99.900000
+ ,0.762253
+ ,1.311900
+ ,105.900000
+ ,99.900000
+ ,0.768403
+ ,1.301400
+ ,107.300000
+ ,99.900000
+ ,0.757518
+ ,1.320100
+ ,107.900000
+ ,99.900000
+ ,0.772917
+ ,1.293800
+ ,105.800000
+ ,99.900000
+ ,0.787774
+ ,1.269400
+ ,105.300000
+ ,99.900000
+ ,0.822030
+ ,1.216500
+ ,107.100000
+ ,99.900000
+ ,0.830772
+ ,1.203700
+ ,107.300000
+ ,99.900000
+ ,0.813537
+ ,1.229200
+ ,106.500000
+ ,99.900000
+ ,0.815927
+ ,1.225600
+ ,107.300000
+ ,99.900000
+ ,0.832293
+ ,1.201500
+ ,107.200000
+ ,99.900000
+ ,0.848464
+ ,1.178600
+ ,106.700000
+ ,99.900000
+ ,0.843455
+ ,1.185600
+ ,107.200000
+ ,99.900000
+ ,0.826241
+ ,1.210300
+ ,111.700000
+ ,99.900000
+ ,0.837661
+ ,1.193800
+ ,111.700000
+ ,99.900000
+ ,0.831947
+ ,1.202000
+ ,110.400000
+ ,99.900000
+ ,0.814930
+ ,1.227100
+ ,113.500000
+ ,99.900000
+ ,0.783085
+ ,1.277000
+ ,113.600000
+ ,99.900000
+ ,0.790514
+ ,1.265000
+ ,114.900000
+ ,99.900000
+ ,0.788395
+ ,1.268400
+ ,113.500000
+ ,99.900000
+ ,0.780579
+ ,1.281100
+ ,113.400000
+ ,102.600000
+ ,0.785731
+ ,1.272700
+ ,113.200000
+ ,102.600000
+ ,0.792959
+ ,1.261100
+ ,113.200000
+ ,102.600000
+ ,0.776337
+ ,1.288100
+ ,113.700000
+ ,102.600000
+ ,0.756830
+ ,1.321300
+ ,115.100000
+ ,102.600000
+ ,0.769290
+ ,1.299900
+ ,115.400000
+ ,102.600000
+ ,0.764877
+ ,1.307400
+ ,116.400000
+ ,102.600000
+ ,0.755173
+ ,1.324200
+ ,117.300000
+ ,102.600000
+ ,0.739864
+ ,1.351600
+ ,118.500000
+ ,102.600000
+ ,0.740138
+ ,1.351100
+ ,120.700000
+ ,102.600000
+ ,0.745212
+ ,1.341900
+ ,120.900000
+ ,102.600000
+ ,0.729076
+ ,1.371600
+ ,123.100000
+ ,102.600000
+ ,0.734107
+ ,1.362200
+ ,120.500000
+ ,102.600000)
+ ,dim=c(4
+ ,80)
+ ,dimnames=list(c('w'
+ ,'x'
+ ,'y'
+ ,'z')
+ ,1:80))
>  y <- array(NA,dim=c(4,80),dimnames=list(c('w','x','y','z'),1:80))
>  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: Wessa P., (2007), Multivariate Correlation Matrix (v1.0.3) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_pairs.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description
> 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/1lbib1196265893.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 
> 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 = -13.1287, p-value < 2.2e-16
alternative hypothesis: true tau is not equal to 0 
sample estimates:
tau 
 -1 

> 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/2dw9k1196265893.tab") 
> 
> system("convert tmp/1lbib1196265893.ps tmp/1lbib1196265893.png")
> 
> 
> proc.time()
   user  system elapsed 
  0.746   0.217   0.822