R version 2.12.0 (2010-10-15)
Copyright (C) 2010 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(97
+ ,197426
+ ,39
+ ,178377
+ ,490
+ ,146
+ ,187326
+ ,173
+ ,250931
+ ,2563
+ ,116
+ ,184923
+ ,165
+ ,226168
+ ,1538
+ ,113
+ ,183500
+ ,181
+ ,211381
+ ,898
+ ,75
+ ,176225
+ ,139
+ ,214738
+ ,1212
+ ,228
+ ,169707
+ ,166
+ ,210012
+ ,790
+ ,138
+ ,169265
+ ,116
+ ,163073
+ ,738
+ ,153
+ ,167949
+ ,114
+ ,164263
+ ,845
+ ,248
+ ,165986
+ ,155
+ ,189944
+ ,1369
+ ,161
+ ,165933
+ ,127
+ ,147581
+ ,1830
+ ,155
+ ,165904
+ ,107
+ ,127667
+ ,711
+ ,142
+ ,160902
+ ,126
+ ,106330
+ ,992
+ ,145
+ ,160141
+ ,161
+ ,175721
+ ,1272
+ ,159
+ ,156349
+ ,185
+ ,169216
+ ,852
+ ,153
+ ,154771
+ ,63
+ ,18284
+ ,575
+ ,130
+ ,154451
+ ,121
+ ,134969
+ ,1101
+ ,177
+ ,151911
+ ,150
+ ,191889
+ ,1410
+ ,181
+ ,151715
+ ,160
+ ,197765
+ ,1352
+ ,140
+ ,150491
+ ,132
+ ,194679
+ ,1208
+ ,196
+ ,150047
+ ,147
+ ,75767
+ ,739
+ ,140
+ ,149959
+ ,176
+ ,195894
+ ,926
+ ,175
+ ,149695
+ ,88
+ ,191179
+ ,865
+ ,155
+ ,147172
+ ,82
+ ,178303
+ ,677
+ ,147
+ ,146975
+ ,75
+ ,135599
+ ,971
+ ,177
+ ,146760
+ ,128
+ ,195791
+ ,1574
+ ,159
+ ,144551
+ ,165
+ ,81716
+ ,1051
+ ,132
+ ,144408
+ ,88
+ ,115466
+ ,763
+ ,94
+ ,144244
+ ,62
+ ,88229
+ ,724
+ ,140
+ ,143592
+ ,93
+ ,113963
+ ,652
+ ,130
+ ,140824
+ ,96
+ ,186099
+ ,504
+ ,176
+ ,140358
+ ,121
+ ,117495
+ ,893
+ ,153
+ ,140015
+ ,146
+ ,145758
+ ,1034
+ ,179
+ ,139165
+ ,143
+ ,184531
+ ,1111
+ ,197
+ ,136588
+ ,135
+ ,134163
+ ,692
+ ,163
+ ,135356
+ ,76
+ ,91502
+ ,740
+ ,170
+ ,134238
+ ,152
+ ,191469
+ ,1716
+ ,145
+ ,134047
+ ,163
+ ,231257
+ ,884
+ ,129
+ ,131072
+ ,154
+ ,114268
+ ,925
+ ,57
+ ,128692
+ ,48
+ ,100187
+ ,723
+ ,144
+ ,127654
+ ,168
+ ,105590
+ ,732
+ ,92
+ ,126817
+ ,143
+ ,94333
+ ,637
+ ,144
+ ,126372
+ ,156
+ ,165278
+ ,1266
+ ,95
+ ,125818
+ ,103
+ ,111669
+ ,527
+ ,126
+ ,125386
+ ,128
+ ,134218
+ ,811
+ ,97
+ ,125081
+ ,121
+ ,135213
+ ,1390
+ ,144
+ ,124089
+ ,96
+ ,130332
+ ,1613
+ ,137
+ ,123534
+ ,76
+ ,100922
+ ,459
+ ,155
+ ,120192
+ ,161
+ ,197680
+ ,1118
+ ,163
+ ,119442
+ ,141
+ ,189723
+ ,1293
+ ,227
+ ,118906
+ ,103
+ ,102509
+ ,636
+ ,187
+ ,117440
+ ,137
+ ,157384
+ ,1031
+ ,148
+ ,116066
+ ,55
+ ,24469
+ ,524
+ ,208
+ ,114948
+ ,176
+ ,165354
+ ,1775
+ ,136
+ ,114799
+ ,66
+ ,153242
+ ,669
+ ,134
+ ,114360
+ ,101
+ ,79367
+ ,2089
+ ,149
+ ,113344
+ ,155
+ ,230054
+ ,1230
+ ,160
+ ,112431
+ ,123
+ ,140303
+ ,847
+ ,187
+ ,112302
+ ,145
+ ,198299
+ ,906
+ ,146
+ ,112098
+ ,134
+ ,106194
+ ,1154
+ ,102
+ ,110529
+ ,164
+ ,232241
+ ,1251
+ ,143
+ ,110459
+ ,75
+ ,113854
+ ,510
+ ,115
+ ,109432
+ ,148
+ ,116938
+ ,698
+ ,151
+ ,108535
+ ,85
+ ,118845
+ ,1586
+ ,144
+ ,108146
+ ,140
+ ,100125
+ ,1001
+ ,118
+ ,105079
+ ,70
+ ,99776
+ ,710
+ ,98
+ ,104978
+ ,117
+ ,139292
+ ,906
+ ,142
+ ,104767
+ ,103
+ ,124527
+ ,1030
+ ,151
+ ,104581
+ ,116
+ ,116136
+ ,1092
+ ,171
+ ,104128
+ ,50
+ ,122975
+ ,511
+ ,156
+ ,103925
+ ,152
+ ,164808
+ ,1319
+ ,171
+ ,103297
+ ,139
+ ,101345
+ ,1186
+ ,142
+ ,103129
+ ,114
+ ,158376
+ ,1201
+ ,148
+ ,103037
+ ,110
+ ,150773
+ ,1443
+ ,96
+ ,102812
+ ,120
+ ,136323
+ ,703
+ ,151
+ ,102153
+ ,98
+ ,80716
+ ,862
+ ,173
+ ,102070
+ ,94
+ ,86480
+ ,1031
+ ,151
+ ,101629
+ ,112
+ ,188355
+ ,1348
+ ,77
+ ,101382
+ ,81
+ ,127097
+ ,866
+ ,135
+ ,101047
+ ,169
+ ,135848
+ ,1079
+ ,71
+ ,100350
+ ,62
+ ,75882
+ ,695
+ ,133
+ ,100087
+ ,102
+ ,129711
+ ,1229
+ ,139
+ ,100046
+ ,133
+ ,128602
+ ,1288
+ ,201
+ ,96125
+ ,107
+ ,106314
+ ,764
+ ,141
+ ,95893
+ ,90
+ ,81180
+ ,919
+ ,158
+ ,95676
+ ,99
+ ,160792
+ ,691
+ ,126
+ ,93879
+ ,152
+ ,170492
+ ,1099
+ ,119
+ ,93487
+ ,84
+ ,133252
+ ,766
+ ,65
+ ,93473
+ ,57
+ ,121850
+ ,1150
+ ,128
+ ,92622
+ ,126
+ ,134097
+ ,1566
+ ,147
+ ,92280
+ ,118
+ ,147341
+ ,668
+ ,90
+ ,92059
+ ,101
+ ,91313
+ ,910
+ ,169
+ ,89626
+ ,85
+ ,134904
+ ,894
+ ,150
+ ,89506
+ ,118
+ ,160501
+ ,1351
+ ,156
+ ,89256
+ ,129
+ ,104864
+ ,1187
+ ,179
+ ,88977
+ ,85
+ ,111563
+ ,784
+ ,149
+ ,86652
+ ,50
+ ,114198
+ ,758
+ ,94
+ ,84601
+ ,85
+ ,105406
+ ,816
+ ,154
+ ,83515
+ ,158
+ ,96785
+ ,1370
+ ,103
+ ,83248
+ ,146
+ ,106020
+ ,785
+ ,148
+ ,83243
+ ,150
+ ,153990
+ ,763
+ ,84
+ ,82317
+ ,77
+ ,111848
+ ,569
+ ,144
+ ,81897
+ ,131
+ ,89770
+ ,781
+ ,203
+ ,81625
+ ,132
+ ,94853
+ ,743
+ ,160
+ ,81351
+ ,107
+ ,102204
+ ,900
+ ,152
+ ,79756
+ ,80
+ ,122531
+ ,575
+ ,147
+ ,79089
+ ,114
+ ,169351
+ ,981
+ ,111
+ ,79011
+ ,97
+ ,80238
+ ,784
+ ,89
+ ,76173
+ ,8
+ ,47552
+ ,179
+ ,87
+ ,72128
+ ,163
+ ,145707
+ ,542
+ ,121
+ ,71571
+ ,102
+ ,75881
+ ,746
+ ,146
+ ,71154
+ ,137
+ ,80906
+ ,767
+ ,100
+ ,70168
+ ,79
+ ,104470
+ ,695
+ ,127
+ ,69867
+ ,83
+ ,100826
+ ,1186
+ ,153
+ ,69652
+ ,56
+ ,33750
+ ,456
+ ,87
+ ,69446
+ ,87
+ ,113713
+ ,724
+ ,129
+ ,68946
+ ,164
+ ,174586
+ ,1145
+ ,113
+ ,68788
+ ,57
+ ,72591
+ ,785
+ ,124
+ ,67150
+ ,110
+ ,114651
+ ,905
+ ,92
+ ,66485
+ ,104
+ ,110896
+ ,661
+ ,112
+ ,66089
+ ,65
+ ,61394
+ ,507
+ ,102
+ ,65594
+ ,48
+ ,92795
+ ,632
+ ,115
+ ,64593
+ ,60
+ ,72558
+ ,790
+ ,148
+ ,64520
+ ,68
+ ,54518
+ ,488
+ ,135
+ ,59938
+ ,149
+ ,82390
+ ,1128
+ ,97
+ ,59900
+ ,104
+ ,96252
+ ,1257
+ ,59
+ ,57224
+ ,86
+ ,80684
+ ,800
+ ,101
+ ,56750
+ ,89
+ ,115750
+ ,846
+ ,27
+ ,56622
+ ,49
+ ,55792
+ ,437
+ ,112
+ ,55918
+ ,74
+ ,83963
+ ,795
+ ,89
+ ,52789
+ ,37
+ ,15673
+ ,309
+ ,40
+ ,48029
+ ,120
+ ,88634
+ ,833
+ ,130
+ ,45724
+ ,87
+ ,74151
+ ,641
+ ,73
+ ,43929
+ ,83
+ ,100792
+ ,415
+ ,64
+ ,43750
+ ,13
+ ,19630
+ ,214
+ ,99
+ ,38692
+ ,30
+ ,68580
+ ,657
+ ,78
+ ,37238
+ ,41
+ ,10901
+ ,716
+ ,110
+ ,37110
+ ,67
+ ,64057
+ ,665)
+ ,dim=c(5
+ ,137)
+ ,dimnames=list(c('FbackMess'
+ ,'CompendiumCharacters'
+ ,'BloggedComputations'
+ ,'CompendiumSeconds'
+ ,'CourseViews
')
+ ,1:137))
>  y <- array(NA,dim=c(5,137),dimnames=list(c('FbackMess','CompendiumCharacters','BloggedComputations','CompendiumSeconds','CourseViews
'),1:137))
>  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'
> 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/rcomp/tmp/18txh1324666870.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)
Warning messages:
1: In strwidth(labels, "user") : font width unknown for character 0xd
2: In text.default(x, y, txt, cex = cex, font = font) :
  font width unknown for character 0xd
3: In text.default(x, y, txt, cex = cex, font = font) :
  font metrics unknown for character 0xd
> dev.off()
null device 
          1 
> 
> #Note: the /var/www/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/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/www/rcomp/tmp/2r1b21324666870.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 = 4.9096, p-value = 9.126e-07
alternative hypothesis: true tau is not equal to 0 
sample estimates:
      tau 
0.2842386 

> 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
> a<-table.end(a)
> table.save(a,file="/var/www/rcomp/tmp/3utvn1324666870.tab") 
> 
> try(system("convert tmp/18txh1324666870.ps tmp/18txh1324666870.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.120   0.090   1.198