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(33024
+ ,31086
+ ,19828
+ ,18932
+ ,32526
+ ,30839
+ ,19967
+ ,18927
+ ,31455
+ ,30051
+ ,19814
+ ,19124
+ ,31524
+ ,29976
+ ,20053
+ ,19066
+ ,31856
+ ,30463
+ ,20719
+ ,19971
+ ,32696
+ ,31422
+ ,21174
+ ,20165
+ ,32584
+ ,31588
+ ,20648
+ ,19705
+ ,33498
+ ,31900
+ ,20659
+ ,19718
+ ,34175
+ ,32878
+ ,20733
+ ,19938
+ ,34172
+ ,33010
+ ,21069
+ ,20039
+ ,34379
+ ,32954
+ ,20566
+ ,19721
+ ,34988
+ ,33076
+ ,20839
+ ,19777
+ ,36158
+ ,35057
+ ,21615
+ ,20505
+ ,37411
+ ,35906
+ ,22739
+ ,21763
+ ,38015
+ ,36100
+ ,23222
+ ,22404
+ ,37577
+ ,35824
+ ,23031
+ ,22038
+ ,36354
+ ,34579
+ ,23014
+ ,22038
+ ,36030
+ ,34484
+ ,22868
+ ,21874
+ ,35636
+ ,33920
+ ,22182
+ ,21269
+ ,35669
+ ,34059
+ ,22177
+ ,21127
+ ,34635
+ ,33812
+ ,21216
+ ,20609
+ ,35496
+ ,34594
+ ,21031
+ ,20565
+ ,36376
+ ,36083
+ ,20968
+ ,19791
+ ,37635
+ ,36563
+ ,21049
+ ,20672
+ ,38875
+ ,37416
+ ,21033
+ ,20938
+ ,38372
+ ,37953
+ ,21078
+ ,20675
+ ,38897
+ ,37517
+ ,20702
+ ,19992
+ ,38018
+ ,37467
+ ,20309
+ ,19801
+ ,37325
+ ,36963
+ ,20449
+ ,20050
+ ,36893
+ ,36019
+ ,20737
+ ,20427
+ ,36117
+ ,35232
+ ,20849
+ ,20815
+ ,37599
+ ,36857
+ ,21966
+ ,21666
+ ,39037
+ ,37978
+ ,23100
+ ,22720
+ ,40809
+ ,40160
+ ,23975
+ ,23650
+ ,42508
+ ,42165
+ ,24350
+ ,24244
+ ,44021
+ ,43069
+ ,24020
+ ,23669
+ ,44088
+ ,43021
+ ,24005
+ ,23881
+ ,44510
+ ,43376
+ ,23602
+ ,23857
+ ,45786
+ ,43978
+ ,24120
+ ,23999
+ ,47349
+ ,45911
+ ,24847
+ ,24780
+ ,48696
+ ,47107
+ ,25702
+ ,25426
+ ,50598
+ ,49168
+ ,26312
+ ,26229
+ ,50066
+ ,48390
+ ,25891
+ ,25973
+ ,49367
+ ,47678
+ ,25172
+ ,25375
+ ,48784
+ ,47822
+ ,25698
+ ,25966
+ ,47841
+ ,46695
+ ,25833
+ ,25391
+ ,48300
+ ,47185
+ ,25658
+ ,26046
+ ,47518
+ ,45684
+ ,25269
+ ,25572
+ ,46504
+ ,44884
+ ,24846
+ ,24900
+ ,45147
+ ,44256
+ ,24390
+ ,24744
+ ,44404
+ ,43637
+ ,23954
+ ,24526
+ ,43455
+ ,42368
+ ,23828
+ ,24274
+ ,42299
+ ,40892
+ ,23507
+ ,23774
+ ,42105
+ ,40616
+ ,23144
+ ,23414
+ ,40152
+ ,39026
+ ,22302
+ ,23002
+ ,39519
+ ,38921
+ ,23028
+ ,23137
+ ,39633
+ ,38512
+ ,22741
+ ,22947
+ ,39376
+ ,38884
+ ,23129
+ ,23733
+ ,38850
+ ,38406
+ ,22911
+ ,23234
+ ,39657
+ ,38804
+ ,22071
+ ,22969
+ ,34804
+ ,34871
+ ,16466
+ ,17708
+ ,34372
+ ,34660
+ ,16370
+ ,17377
+ ,32678
+ ,33104
+ ,15049
+ ,16273
+ ,28420
+ ,28952
+ ,13174
+ ,14342
+ ,25420
+ ,26488
+ ,12231
+ ,13522
+ ,27683
+ ,29418
+ ,13620
+ ,15210
+ ,29904
+ ,32315
+ ,14317
+ ,16493
+ ,30546
+ ,32885
+ ,14039
+ ,16701
+ ,29142
+ ,31565
+ ,13526
+ ,15662
+ ,27724
+ ,30782
+ ,12826
+ ,15526
+ ,27069
+ ,30442
+ ,12360
+ ,15413
+ ,26665
+ ,30851
+ ,12592
+ ,15805
+ ,26004
+ ,30432
+ ,12381
+ ,15802
+ ,25767
+ ,31260
+ ,12554
+ ,16753
+ ,24915
+ ,30737
+ ,12338
+ ,16906
+ ,23689
+ ,30129
+ ,11768
+ ,16891
+ ,20915
+ ,27672
+ ,10687
+ ,15703
+ ,19414
+ ,26469
+ ,9964
+ ,15429
+ ,17824
+ ,24895
+ ,9338
+ ,14762
+ ,16348
+ ,24427
+ ,8697
+ ,14426
+ ,15571
+ ,23252
+ ,8068
+ ,14250
+ ,13929
+ ,21815
+ ,7295
+ ,13267
+ ,12480
+ ,20837
+ ,6372
+ ,12397
+ ,10837
+ ,18537
+ ,5649
+ ,11586
+ ,9473
+ ,17237
+ ,4926
+ ,10888
+ ,8051
+ ,15476
+ ,4199
+ ,9841
+ ,5278
+ ,10709
+ ,2568
+ ,6443
+ ,3008
+ ,6776
+ ,1461
+ ,4019
+ ,2404
+ ,5810
+ ,1173
+ ,3449
+ ,2298
+ ,5765
+ ,1084
+ ,3179
+ ,2260
+ ,5775
+ ,978
+ ,3341
+ ,1938
+ ,5589
+ ,947
+ ,3325
+ ,1371
+ ,4687
+ ,679
+ ,2478
+ ,1009
+ ,3630
+ ,457
+ ,1982
+ ,686
+ ,2552
+ ,262
+ ,1405
+ ,493
+ ,1928
+ ,218
+ ,1059
+ ,285
+ ,1323
+ ,132
+ ,740
+ ,192
+ ,1005
+ ,70
+ ,533
+ ,129
+ ,678
+ ,44
+ ,366
+ ,60
+ ,397
+ ,24
+ ,224
+ ,54
+ ,286
+ ,20
+ ,147
+ ,26
+ ,166
+ ,4
+ ,75
+ ,11
+ ,80
+ ,4
+ ,54
+ ,3
+ ,53
+ ,1
+ ,23
+ ,0
+ ,32
+ ,0
+ ,16
+ ,2
+ ,11
+ ,0
+ ,6
+ ,1
+ ,6
+ ,0
+ ,7
+ ,0
+ ,4
+ ,0
+ ,2
+ ,0
+ ,2
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,1
+ ,0
+ ,0)
+ ,dim=c(4
+ ,111)
+ ,dimnames=list(c('MVG'
+ ,'VVG'
+ ,'MWG'
+ ,'VWG')
+ ,1:111))
>  y <- array(NA,dim=c(4,111),dimnames=list(c('MVG','VVG','MWG','VWG'),1:111))
>  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: Patrick Wessa
> #To cite this work: Patrick Wessa, (2010), Multivariate Correlation Matrix (v1.0.4) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/Patrick.Wessa/rwasp_pairs.wasp#output/
> #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/1ueg01292700671.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")
> 
> n <- length(y[,1])
> n
[1] 4
> 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)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/2xee61292700671.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 = 14.835, p-value < 2.2e-16
alternative hypothesis: true tau is not equal to 0 
sample estimates:
      tau 
0.9542626 

> 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)
+ }
+ }
Warning messages:
1: In cor.test.default(y[i, ], y[j, ], method = "spearman") :
  Cannot compute exact p-values with ties
2: In cor.test.default(y[i, ], y[j, ], method = "spearman") :
  Cannot compute exact p-values with ties
3: In cor.test.default(y[i, ], y[j, ], method = "spearman") :
  Cannot compute exact p-values with ties
4: In cor.test.default(y[i, ], y[j, ], method = "spearman") :
  Cannot compute exact p-values with ties
5: In cor.test.default(y[i, ], y[j, ], method = "spearman") :
  Cannot compute exact p-values with ties
6: In cor.test.default(y[i, ], y[j, ], method = "spearman") :
  Cannot compute exact p-values with ties
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/3b6cf1292700671.tab") 
> 
> try(system("convert tmp/1ueg01292700671.ps tmp/1ueg01292700671.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  0.602   0.235   3.423