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(108.35
+ ,98.68
+ ,100.70
+ ,104.38
+ ,97.72
+ ,15.38
+ ,31.27
+ ,109.87
+ ,99.21
+ ,99.62
+ ,103.97
+ ,98.01
+ ,15.03
+ ,35.83
+ ,111.30
+ ,99.36
+ ,99.83
+ ,103.32
+ ,97.78
+ ,15.21
+ ,37.12
+ ,115.50
+ ,100.72
+ ,100.74
+ ,105.01
+ ,98.04
+ ,15.20
+ ,36.77
+ ,116.22
+ ,102.27
+ ,100.84
+ ,104.88
+ ,98.54
+ ,14.60
+ ,35.17
+ ,116.63
+ ,102.62
+ ,100.85
+ ,104.46
+ ,98.39
+ ,13.79
+ ,37.25
+ ,116.84
+ ,102.97
+ ,99.71
+ ,104.71
+ ,98.58
+ ,14.54
+ ,33.77
+ ,116.63
+ ,102.88
+ ,100.80
+ ,106.09
+ ,98.91
+ ,14.31
+ ,30.59
+ ,117.03
+ ,102.90
+ ,100.06
+ ,106.54
+ ,98.68
+ ,13.93
+ ,33.59
+ ,117.00
+ ,103.01
+ ,100.57
+ ,104.36
+ ,98.59
+ ,14.82
+ ,37.24
+ ,117.14
+ ,103.02
+ ,99.79
+ ,105.31
+ ,99.13
+ ,14.46
+ ,34.81
+ ,116.64
+ ,103.73
+ ,99.90
+ ,105.07
+ ,98.70
+ ,14.85
+ ,34.94
+ ,117.24
+ ,104.18
+ ,100.12
+ ,105.39
+ ,99.00
+ ,14.95
+ ,34.47
+ ,117.52
+ ,103.73
+ ,100.40
+ ,105.65
+ ,98.80
+ ,14.43
+ ,30.48
+ ,117.83
+ ,103.78
+ ,100.51
+ ,108.25
+ ,98.80
+ ,14.84
+ ,30.94
+ ,119.79
+ ,103.61
+ ,100.70
+ ,107.71
+ ,99.29
+ ,14.39
+ ,30.60
+ ,120.86
+ ,103.84
+ ,100.62
+ ,108.58
+ ,99.69
+ ,15.70
+ ,28.42
+ ,120.75
+ ,103.86
+ ,99.70
+ ,108.27
+ ,100.01
+ ,15.34
+ ,25.89
+ ,120.63
+ ,104.14
+ ,99.48
+ ,107.62
+ ,99.85
+ ,13.98
+ ,26.32
+ ,120.89
+ ,104.05
+ ,99.36
+ ,108.80
+ ,99.66
+ ,14.75
+ ,27.18
+ ,120.23
+ ,104.01
+ ,99.39
+ ,109.26
+ ,101.18
+ ,14.81
+ ,25.85
+ ,121.19
+ ,104.49
+ ,99.45
+ ,108.58
+ ,101.47
+ ,14.67
+ ,26.32
+ ,120.79
+ ,104.83
+ ,99.28
+ ,107.05
+ ,101.28
+ ,15.03
+ ,23.07
+ ,120.09
+ ,104.78
+ ,99.40
+ ,109.20
+ ,101.80
+ ,14.34
+ ,20.19
+ ,120.86
+ ,104.95
+ ,99.10
+ ,109.52
+ ,102.48
+ ,12.54
+ ,18.65
+ ,121.10
+ ,105.28
+ ,99.48
+ ,111.12
+ ,102.32
+ ,11.37
+ ,17.74
+ ,121.47
+ ,105.28
+ ,99.74
+ ,108.74
+ ,102.30
+ ,12.58
+ ,17.26
+ ,122.01
+ ,105.91
+ ,100.42
+ ,110.53
+ ,102.84
+ ,13.06
+ ,16.01
+ ,123.94
+ ,106.81
+ ,100.80
+ ,110.44
+ ,102.36
+ ,12.50
+ ,17.94
+ ,125.78
+ ,106.39
+ ,100.66
+ ,111.02
+ ,102.16
+ ,11.11
+ ,15.53
+ ,125.31
+ ,107.02
+ ,101.03
+ ,111.13
+ ,102.57
+ ,12.39
+ ,14.49
+ ,125.79
+ ,106.92
+ ,101.22
+ ,110.90
+ ,102.49
+ ,12.34
+ ,15.35
+ ,126.12
+ ,107.01
+ ,101.23
+ ,111.32
+ ,104.11
+ ,11.54
+ ,14.67
+ ,125.57
+ ,106.79
+ ,100.10
+ ,109.37
+ ,104.78
+ ,10.22
+ ,12.95
+ ,125.44
+ ,107.41
+ ,99.98
+ ,110.18
+ ,104.13
+ ,8.50
+ ,8.81
+ ,126.12
+ ,107.13
+ ,99.91
+ ,110.74
+ ,104.22
+ ,9.06
+ ,9.33
+ ,126.01
+ ,107.54
+ ,99.84
+ ,111.70
+ ,104.73
+ ,9.28
+ ,9.31
+ ,126.50
+ ,108.48
+ ,99.68
+ ,111.33
+ ,104.99
+ ,7.24
+ ,9.03
+ ,126.13
+ ,108.50
+ ,99.74
+ ,110.86
+ ,104.70
+ ,7.58
+ ,10.96
+ ,126.66
+ ,108.27
+ ,99.71
+ ,109.48
+ ,104.69
+ ,7.81
+ ,14.26
+ ,126.33
+ ,109.42
+ ,99.35
+ ,108.77
+ ,104.85
+ ,8.54
+ ,14.20
+ ,126.61
+ ,110.09
+ ,99.21
+ ,109.81
+ ,104.24
+ ,9.27
+ ,13.70
+ ,126.36
+ ,109.98
+ ,99.21
+ ,109.15
+ ,104.74
+ ,10.11
+ ,17.46
+ ,126.83
+ ,109.99
+ ,99.16
+ ,109.63
+ ,104.20
+ ,9.21
+ ,18.73
+ ,125.90
+ ,109.54
+ ,99.20
+ ,111.32
+ ,105.62
+ ,10.71
+ ,20.37
+ ,126.29
+ ,108.85
+ ,99.08
+ ,109.75
+ ,106.08
+ ,10.85
+ ,18.72
+ ,126.37
+ ,106.76
+ ,98.16
+ ,110.37
+ ,105.46
+ ,11.77
+ ,21.60
+ ,125.11
+ ,107.56
+ ,98.00
+ ,108.30
+ ,105.42
+ ,11.81
+ ,22.75)
+ ,dim=c(7
+ ,48)
+ ,dimnames=list(c('Coffee'
+ ,'Tea'
+ ,'Sugar'
+ ,'Water'
+ ,'Soda'
+ ,'SaraLee'
+ ,'Starbucks')
+ ,1:48))
>  y <- array(NA,dim=c(7,48),dimnames=list(c('Coffee','Tea','Sugar','Water','Soda','SaraLee','Starbucks'),1:48))
>  for (i in 1:dim(x)[1])
+  {
+  	for (j in 1:dim(x)[2])
+  	{
+  		y[i,j] <- as.numeric(x[i,j])
+  	}
+  }
> par1 = 'kendall'
> main = 'Correlation Matrix'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Dr. Ian E. Holliday
> #To cite this work: Ian E. Holliday, 2009, 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: 
> #Technical 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=par1)
+ 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/1nofh1293008046.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)
There were 21 warnings (use warnings() to see them)
> 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")
> 
> n <- length(y[,1])
> n
[1] 7
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,paste('Correlations for all pairs of data series (method=',par1,')',sep=''),n+1,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,' ',header=TRUE)
> for (i in 1:n) {
+ a<-table.element(a,dimnames(t(x))[[2]][i],header=TRUE)
+ }
> a<-table.row.end(a)
> for (i in 1:n) {
+ a<-table.row.start(a)
+ a<-table.element(a,dimnames(t(x))[[2]][i],header=TRUE)
+ for (j in 1:n) {
+ r <- cor.test(y[i,],y[j,],method=par1)
+ a<-table.element(a,round(r$estimate,3))
+ }
+ a<-table.row.end(a)
+ }
There were 49 warnings (use warnings() to see them)
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/28oe51293008046.tab") 
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,'Correlations for all pairs of data series with p-values',4,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'pair',1,TRUE)
> a<-table.element(a,'Pearson r',1,TRUE)
> a<-table.element(a,'Spearman rho',1,TRUE)
> a<-table.element(a,'Kendall tau',1,TRUE)
> a<-table.row.end(a)
> cor.test(y[1,],y[2,],method=par1)

	Kendall's rank correlation tau

data:  y[1, ] and y[2, ] 
z = 8.3831, p-value < 2.2e-16
alternative hypothesis: true tau is not equal to 0 
sample estimates:
    tau 
0.83785 

Warning message:
In cor.test.default(y[1, ], y[2, ], method = par1) :
  Cannot compute exact p-value with ties
> for (i in 1:(n-1))
+ {
+ for (j in (i+1):n)
+ {
+ a<-table.row.start(a)
+ dum <- paste(dimnames(t(x))[[2]][i],';',dimnames(t(x))[[2]][j],sep='')
+ a<-table.element(a,dum,header=TRUE)
+ rp <- cor.test(y[i,],y[j,],method='pearson')
+ a<-table.element(a,round(rp$estimate,4))
+ rs <- cor.test(y[i,],y[j,],method='spearman')
+ a<-table.element(a,round(rs$estimate,4))
+ rk <- cor.test(y[i,],y[j,],method='kendall')
+ a<-table.element(a,round(rk$estimate,4))
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'p-value',header=T)
+ a<-table.element(a,paste('(',round(rp$p.value,4),')',sep=''))
+ a<-table.element(a,paste('(',round(rs$p.value,4),')',sep=''))
+ a<-table.element(a,paste('(',round(rk$p.value,4),')',sep=''))
+ a<-table.row.end(a)
+ }
+ }
There were 42 warnings (use warnings() to see them)
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/3bpcb1293008046.tab") 
> try(system("convert tmp/1nofh1293008046.ps tmp/1nofh1293008046.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  0.850   0.277   1.734