R version 2.13.0 (2011-04-13)
Copyright (C) 2011 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i486-pc-linux-gnu (32-bit)

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(264530
+ ,124252
+ ,25695
+ ,147
+ ,148
+ ,135248
+ ,98956
+ ,19967
+ ,126
+ ,124
+ ,207253
+ ,98073
+ ,14338
+ ,108
+ ,108
+ ,202898
+ ,106816
+ ,34117
+ ,145
+ ,142
+ ,145249
+ ,41449
+ ,9713
+ ,68
+ ,66
+ ,65295
+ ,76173
+ ,10024
+ ,49
+ ,47
+ ,442857
+ ,177551
+ ,39981
+ ,171
+ ,163
+ ,33186
+ ,22807
+ ,1271
+ ,5
+ ,5
+ ,183696
+ ,126938
+ ,30207
+ ,106
+ ,106
+ ,194006
+ ,61680
+ ,18035
+ ,88
+ ,87
+ ,287239
+ ,72117
+ ,21609
+ ,145
+ ,141
+ ,205260
+ ,79738
+ ,19836
+ ,93
+ ,88
+ ,141987
+ ,57793
+ ,9028
+ ,60
+ ,60
+ ,322679
+ ,91677
+ ,21750
+ ,145
+ ,145
+ ,199717
+ ,64631
+ ,10038
+ ,95
+ ,95
+ ,349227
+ ,106385
+ ,30276
+ ,144
+ ,137
+ ,278282
+ ,161961
+ ,34972
+ ,179
+ ,177
+ ,274382
+ ,112669
+ ,19954
+ ,102
+ ,102
+ ,157448
+ ,114029
+ ,28113
+ ,157
+ ,151
+ ,244170
+ ,124550
+ ,18830
+ ,170
+ ,156
+ ,256814
+ ,105416
+ ,37144
+ ,140
+ ,140
+ ,405942
+ ,72875
+ ,17916
+ ,133
+ ,130
+ ,161189
+ ,81964
+ ,16186
+ ,74
+ ,71
+ ,156389
+ ,104880
+ ,19195
+ ,120
+ ,116
+ ,200181
+ ,76302
+ ,29124
+ ,134
+ ,129
+ ,192645
+ ,96740
+ ,29813
+ ,108
+ ,107
+ ,249893
+ ,93071
+ ,20270
+ ,132
+ ,128
+ ,241171
+ ,78912
+ ,26105
+ ,125
+ ,119
+ ,143182
+ ,35224
+ ,9155
+ ,66
+ ,62
+ ,285266
+ ,90694
+ ,18113
+ ,130
+ ,124
+ ,243048
+ ,125369
+ ,40546
+ ,143
+ ,140
+ ,176062
+ ,80849
+ ,10096
+ ,150
+ ,144
+ ,305210
+ ,104434
+ ,32338
+ ,152
+ ,150
+ ,87995
+ ,65702
+ ,2871
+ ,28
+ ,28
+ ,343613
+ ,108179
+ ,36592
+ ,190
+ ,177
+ ,264159
+ ,63583
+ ,4914
+ ,73
+ ,73
+ ,394976
+ ,95066
+ ,30190
+ ,115
+ ,111
+ ,192718
+ ,62486
+ ,18153
+ ,101
+ ,98
+ ,114673
+ ,31081
+ ,12558
+ ,41
+ ,41
+ ,310170
+ ,94584
+ ,32894
+ ,147
+ ,139
+ ,292891
+ ,87408
+ ,24138
+ ,107
+ ,107
+ ,157518
+ ,68966
+ ,16628
+ ,103
+ ,102
+ ,180362
+ ,88766
+ ,26369
+ ,84
+ ,80
+ ,146175
+ ,57139
+ ,14171
+ ,68
+ ,66
+ ,140319
+ ,90586
+ ,8500
+ ,52
+ ,51
+ ,405293
+ ,109249
+ ,11940
+ ,70
+ ,69
+ ,78800
+ ,33032
+ ,7935
+ ,21
+ ,21
+ ,201970
+ ,96056
+ ,19456
+ ,155
+ ,155
+ ,312381
+ ,146648
+ ,21347
+ ,165
+ ,163
+ ,166270
+ ,80613
+ ,24095
+ ,124
+ ,121
+ ,199186
+ ,87026
+ ,26204
+ ,121
+ ,118
+ ,24188
+ ,5950
+ ,2694
+ ,7
+ ,7
+ ,348637
+ ,131106
+ ,20366
+ ,161
+ ,154
+ ,65029
+ ,32551
+ ,3597
+ ,21
+ ,21
+ ,101097
+ ,31701
+ ,5296
+ ,35
+ ,35
+ ,255082
+ ,91072
+ ,29463
+ ,125
+ ,122
+ ,287314
+ ,159803
+ ,35838
+ ,157
+ ,152
+ ,308944
+ ,143950
+ ,42590
+ ,256
+ ,255
+ ,280943
+ ,112368
+ ,38665
+ ,192
+ ,177
+ ,225816
+ ,82124
+ ,19442
+ ,86
+ ,83
+ ,348955
+ ,144068
+ ,25515
+ ,164
+ ,164
+ ,283283
+ ,162627
+ ,51318
+ ,213
+ ,202
+ ,199642
+ ,55062
+ ,11807
+ ,80
+ ,77
+ ,232791
+ ,95329
+ ,24130
+ ,122
+ ,118
+ ,212262
+ ,105612
+ ,34053
+ ,122
+ ,123
+ ,206826
+ ,62853
+ ,22641
+ ,113
+ ,109
+ ,180424
+ ,125976
+ ,18898
+ ,128
+ ,126
+ ,204450
+ ,79146
+ ,24539
+ ,117
+ ,114
+ ,200629
+ ,108461
+ ,21664
+ ,162
+ ,161
+ ,138731
+ ,99971
+ ,21577
+ ,87
+ ,85
+ ,219808
+ ,77826
+ ,16643
+ ,103
+ ,101
+ ,73566
+ ,22618
+ ,3007
+ ,26
+ ,25
+ ,222145
+ ,84892
+ ,18798
+ ,104
+ ,102
+ ,181728
+ ,92059
+ ,24648
+ ,127
+ ,126
+ ,150006
+ ,77993
+ ,20286
+ ,132
+ ,130
+ ,325723
+ ,104155
+ ,23999
+ ,112
+ ,112
+ ,268863
+ ,109840
+ ,26813
+ ,155
+ ,150
+ ,202410
+ ,238712
+ ,14718
+ ,57
+ ,54
+ ,173420
+ ,67486
+ ,16963
+ ,109
+ ,106
+ ,162366
+ ,68007
+ ,16673
+ ,92
+ ,90
+ ,136341
+ ,48194
+ ,14646
+ ,57
+ ,55
+ ,390163
+ ,134796
+ ,31772
+ ,145
+ ,139
+ ,145905
+ ,38692
+ ,9648
+ ,38
+ ,38
+ ,248834
+ ,93587
+ ,23096
+ ,152
+ ,148
+ ,80953
+ ,56622
+ ,7905
+ ,59
+ ,58
+ ,133301
+ ,15986
+ ,4527
+ ,27
+ ,27
+ ,138630
+ ,113402
+ ,37432
+ ,104
+ ,104
+ ,334082
+ ,97967
+ ,21082
+ ,80
+ ,75
+ ,277542
+ ,74844
+ ,30437
+ ,76
+ ,73
+ ,170849
+ ,136051
+ ,36288
+ ,163
+ ,157
+ ,236398
+ ,50548
+ ,12369
+ ,89
+ ,87
+ ,207178
+ ,112215
+ ,23774
+ ,199
+ ,186
+ ,157125
+ ,59591
+ ,8108
+ ,89
+ ,88
+ ,242395
+ ,59938
+ ,15049
+ ,107
+ ,107
+ ,273632
+ ,137639
+ ,36021
+ ,137
+ ,131
+ ,178489
+ ,143372
+ ,30391
+ ,123
+ ,123
+ ,210247
+ ,138599
+ ,30910
+ ,152
+ ,149
+ ,268066
+ ,174110
+ ,40656
+ ,202
+ ,201
+ ,351056
+ ,135062
+ ,35070
+ ,159
+ ,145
+ ,368833
+ ,175681
+ ,47250
+ ,282
+ ,273
+ ,247842
+ ,130307
+ ,36236
+ ,111
+ ,111
+ ,268118
+ ,139141
+ ,29601
+ ,197
+ ,195
+ ,174311
+ ,44244
+ ,10443
+ ,72
+ ,69
+ ,43287
+ ,43750
+ ,7409
+ ,49
+ ,49
+ ,182915
+ ,48029
+ ,18213
+ ,82
+ ,82
+ ,189021
+ ,95216
+ ,40856
+ ,192
+ ,193
+ ,237531
+ ,92288
+ ,36471
+ ,102
+ ,102
+ ,279589
+ ,94588
+ ,26077
+ ,127
+ ,124
+ ,106655
+ ,197426
+ ,24797
+ ,60
+ ,59
+ ,135798
+ ,151244
+ ,6816
+ ,61
+ ,61
+ ,292930
+ ,139206
+ ,25527
+ ,106
+ ,102
+ ,266805
+ ,106271
+ ,22139
+ ,139
+ ,138
+ ,23623
+ ,1168
+ ,238
+ ,11
+ ,11
+ ,174970
+ ,71764
+ ,24459
+ ,114
+ ,114
+ ,61857
+ ,25162
+ ,3913
+ ,31
+ ,28
+ ,151306
+ ,45635
+ ,9895
+ ,132
+ ,101
+ ,358662
+ ,101817
+ ,25902
+ ,210
+ ,208
+ ,21054
+ ,855
+ ,338
+ ,4
+ ,4
+ ,230091
+ ,100174
+ ,12937
+ ,98
+ ,93
+ ,31414
+ ,14116
+ ,3988
+ ,39
+ ,39
+ ,284519
+ ,85008
+ ,23370
+ ,119
+ ,114
+ ,209481
+ ,124254
+ ,24015
+ ,107
+ ,104
+ ,161691
+ ,105793
+ ,3870
+ ,41
+ ,40
+ ,137093
+ ,117129
+ ,14648
+ ,97
+ ,94
+ ,38214
+ ,8773
+ ,1888
+ ,16
+ ,16
+ ,166059
+ ,94747
+ ,16768
+ ,65
+ ,64
+ ,323179
+ ,107549
+ ,33400
+ ,156
+ ,154
+ ,186273
+ ,97392
+ ,23770
+ ,158
+ ,145
+ ,374269
+ ,126893
+ ,34762
+ ,160
+ ,150
+ ,280897
+ ,118850
+ ,18793
+ ,161
+ ,156
+ ,371680
+ ,234853
+ ,48186
+ ,238
+ ,229
+ ,179928
+ ,74783
+ ,20140
+ ,85
+ ,84
+ ,94381
+ ,66089
+ ,8728
+ ,50
+ ,49
+ ,273422
+ ,95684
+ ,19060
+ ,106
+ ,101
+ ,382564
+ ,139537
+ ,26880
+ ,184
+ ,179
+ ,118033
+ ,144253
+ ,415
+ ,15
+ ,15
+ ,370878
+ ,153824
+ ,38902
+ ,155
+ ,154
+ ,147989
+ ,63995
+ ,17375
+ ,90
+ ,90
+ ,236370
+ ,84891
+ ,31360
+ ,133
+ ,133
+ ,193456
+ ,61263
+ ,15051
+ ,136
+ ,133
+ ,189020
+ ,106221
+ ,16785
+ ,109
+ ,97
+ ,344751
+ ,113587
+ ,15886
+ ,118
+ ,116
+ ,224936
+ ,113864
+ ,28548
+ ,222
+ ,221
+ ,173260
+ ,37238
+ ,2805
+ ,16
+ ,15
+ ,291777
+ ,119906
+ ,34012
+ ,162
+ ,157
+ ,130908
+ ,135096
+ ,19215
+ ,108
+ ,105
+ ,209639
+ ,151611
+ ,34177
+ ,153
+ ,150
+ ,262412
+ ,144645
+ ,32990
+ ,181
+ ,177
+ ,1
+ ,0
+ ,0
+ ,0
+ ,0
+ ,14688
+ ,6023
+ ,2065
+ ,5
+ ,5
+ ,98
+ ,0
+ ,0
+ ,0
+ ,0
+ ,455
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,195822
+ ,77457
+ ,17428
+ ,113
+ ,111
+ ,347930
+ ,62464
+ ,19912
+ ,165
+ ,165
+ ,0
+ ,0
+ ,0
+ ,0
+ ,0
+ ,203
+ ,0
+ ,0
+ ,0
+ ,0
+ ,7199
+ ,1644
+ ,556
+ ,6
+ ,6
+ ,46660
+ ,6179
+ ,2089
+ ,13
+ ,13
+ ,17547
+ ,3926
+ ,2658
+ ,3
+ ,3
+ ,107465
+ ,42087
+ ,1801
+ ,33
+ ,33
+ ,969
+ ,0
+ ,0
+ ,0
+ ,0
+ ,182303
+ ,87656
+ ,16541
+ ,67
+ ,63)
+ ,dim=c(5
+ ,164)
+ ,dimnames=list(c('TimeRFCsec'
+ ,'CWCharacters'
+ ,'CWRevisions'
+ ,'CWHyperlinks'
+ ,'CWBlogs')
+ ,1:164))
>  y <- array(NA,dim=c(5,164),dimnames=list(c('TimeRFCsec','CWCharacters','CWRevisions','CWHyperlinks','CWBlogs'),1:164))
>  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 = 'tau 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/wessaorg/rcomp/tmp/12mj31323779407.ps",horizontal=F,onefile=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/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/wessaorg/rcomp/createtable")
> 
> n <- length(y[,1])
> n
[1] 5
> 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/wessaorg/rcomp/tmp/2xy171323779407.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 = 9.4569, p-value < 2.2e-16
alternative hypothesis: true tau is not equal to 0 
sample estimates:
      tau 
0.4981089 

> 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
7: In cor.test.default(y[i, ], y[j, ], method = "spearman") :
  Cannot compute exact p-values with ties
8: In cor.test.default(y[i, ], y[j, ], method = "spearman") :
  Cannot compute exact p-values with ties
9: In cor.test.default(y[i, ], y[j, ], method = "spearman") :
  Cannot compute exact p-values with ties
10: 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/wessaorg/rcomp/tmp/3ro881323779407.tab") 
> 
> try(system("convert tmp/12mj31323779407.ps tmp/12mj31323779407.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.074   0.109   1.171