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(140824
+ ,272545
+ ,96
+ ,32033
+ ,110459
+ ,179444
+ ,71
+ ,20654
+ ,105079
+ ,222373
+ ,70
+ ,16346
+ ,112098
+ ,218443
+ ,134
+ ,35926
+ ,43929
+ ,171533
+ ,67
+ ,10621
+ ,76173
+ ,70849
+ ,8
+ ,10024
+ ,187326
+ ,517243
+ ,149
+ ,43068
+ ,22807
+ ,33186
+ ,1
+ ,1271
+ ,144408
+ ,217320
+ ,83
+ ,34416
+ ,66485
+ ,213274
+ ,82
+ ,20318
+ ,79089
+ ,307153
+ ,92
+ ,24409
+ ,81625
+ ,239766
+ ,117
+ ,20648
+ ,68788
+ ,185384
+ ,53
+ ,12347
+ ,103297
+ ,364569
+ ,139
+ ,21857
+ ,69446
+ ,251622
+ ,80
+ ,11034
+ ,114948
+ ,384524
+ ,175
+ ,33433
+ ,167949
+ ,325587
+ ,114
+ ,35902
+ ,125081
+ ,343085
+ ,100
+ ,22355
+ ,125818
+ ,176082
+ ,103
+ ,31219
+ ,136588
+ ,266736
+ ,135
+ ,21983
+ ,112431
+ ,278265
+ ,123
+ ,40085
+ ,103037
+ ,442703
+ ,87
+ ,18507
+ ,82317
+ ,180393
+ ,66
+ ,16278
+ ,118906
+ ,189897
+ ,103
+ ,24662
+ ,83515
+ ,238434
+ ,144
+ ,31452
+ ,104581
+ ,238006
+ ,113
+ ,32580
+ ,103129
+ ,267268
+ ,99
+ ,22883
+ ,83243
+ ,270787
+ ,117
+ ,27652
+ ,37110
+ ,155915
+ ,57
+ ,9845
+ ,113344
+ ,359764
+ ,138
+ ,20190
+ ,139165
+ ,283441
+ ,123
+ ,46201
+ ,86652
+ ,216853
+ ,44
+ ,10971
+ ,112302
+ ,336277
+ ,140
+ ,34811
+ ,69652
+ ,101014
+ ,43
+ ,3029
+ ,119442
+ ,396684
+ ,138
+ ,38941
+ ,69867
+ ,273950
+ ,83
+ ,4958
+ ,101629
+ ,425253
+ ,112
+ ,32344
+ ,70168
+ ,227636
+ ,79
+ ,19433
+ ,31081
+ ,115658
+ ,33
+ ,12558
+ ,103925
+ ,369702
+ ,135
+ ,36524
+ ,92622
+ ,329639
+ ,123
+ ,26041
+ ,79011
+ ,186409
+ ,76
+ ,16637
+ ,93487
+ ,206196
+ ,78
+ ,28395
+ ,64520
+ ,182286
+ ,68
+ ,16747
+ ,93473
+ ,153778
+ ,50
+ ,9105
+ ,114360
+ ,455401
+ ,101
+ ,11941
+ ,33032
+ ,78800
+ ,20
+ ,7935
+ ,96125
+ ,208277
+ ,101
+ ,19499
+ ,151911
+ ,358127
+ ,149
+ ,22938
+ ,89256
+ ,176357
+ ,99
+ ,25314
+ ,95676
+ ,224649
+ ,95
+ ,28527
+ ,5950
+ ,24188
+ ,8
+ ,2694
+ ,149695
+ ,380576
+ ,88
+ ,20867
+ ,32551
+ ,65029
+ ,21
+ ,3597
+ ,31701
+ ,101097
+ ,30
+ ,5296
+ ,100087
+ ,279128
+ ,97
+ ,32982
+ ,169707
+ ,328024
+ ,130
+ ,38975
+ ,150491
+ ,345663
+ ,132
+ ,42721
+ ,120192
+ ,367112
+ ,161
+ ,41455
+ ,95893
+ ,261528
+ ,89
+ ,23923
+ ,151715
+ ,385286
+ ,159
+ ,26719
+ ,176225
+ ,304468
+ ,139
+ ,53405
+ ,59900
+ ,272297
+ ,102
+ ,12526
+ ,104767
+ ,264889
+ ,92
+ ,26584
+ ,114799
+ ,229585
+ ,52
+ ,37062
+ ,72128
+ ,225460
+ ,107
+ ,25696
+ ,143592
+ ,203663
+ ,93
+ ,24634
+ ,89626
+ ,239986
+ ,85
+ ,27269
+ ,131072
+ ,238482
+ ,137
+ ,25270
+ ,126817
+ ,175816
+ ,143
+ ,24634
+ ,81351
+ ,239337
+ ,99
+ ,17828
+ ,22618
+ ,73566
+ ,22
+ ,3007
+ ,88977
+ ,242622
+ ,78
+ ,20065
+ ,92059
+ ,194855
+ ,83
+ ,24648
+ ,81897
+ ,209049
+ ,131
+ ,21588
+ ,108146
+ ,363250
+ ,140
+ ,25217
+ ,126372
+ ,350641
+ ,142
+ ,30927
+ ,249771
+ ,217036
+ ,80
+ ,18487
+ ,71154
+ ,208804
+ ,133
+ ,18050
+ ,71571
+ ,195793
+ ,95
+ ,17696
+ ,55918
+ ,153613
+ ,62
+ ,17326
+ ,160141
+ ,475859
+ ,161
+ ,39361
+ ,38692
+ ,145943
+ ,30
+ ,9648
+ ,102812
+ ,287359
+ ,118
+ ,26759
+ ,56622
+ ,80953
+ ,49
+ ,7905
+ ,15986
+ ,150216
+ ,52
+ ,4527
+ ,123534
+ ,179317
+ ,76
+ ,41517
+ ,108535
+ ,384720
+ ,81
+ ,21261
+ ,93879
+ ,342522
+ ,146
+ ,36099
+ ,144551
+ ,179811
+ ,165
+ ,39039
+ ,56750
+ ,273486
+ ,83
+ ,13841
+ ,127654
+ ,262811
+ ,161
+ ,23841
+ ,65594
+ ,190312
+ ,48
+ ,8589
+ ,59938
+ ,276134
+ ,149
+ ,15049
+ ,146975
+ ,328875
+ ,75
+ ,39038
+ ,161100
+ ,189654
+ ,83
+ ,34974
+ ,168553
+ ,259191
+ ,110
+ ,39932
+ ,183500
+ ,297464
+ ,164
+ ,43840
+ ,165986
+ ,393266
+ ,152
+ ,43146
+ ,184923
+ ,399801
+ ,165
+ ,50099
+ ,140358
+ ,291970
+ ,121
+ ,40312
+ ,149959
+ ,296670
+ ,150
+ ,32616
+ ,57224
+ ,186856
+ ,73
+ ,11338
+ ,43750
+ ,43287
+ ,13
+ ,7409
+ ,48029
+ ,185468
+ ,89
+ ,18213
+ ,104978
+ ,250254
+ ,105
+ ,45873
+ ,100046
+ ,268391
+ ,129
+ ,39844
+ ,101047
+ ,314131
+ ,169
+ ,28317
+ ,197426
+ ,153059
+ ,28
+ ,24797
+ ,160902
+ ,158511
+ ,118
+ ,7471
+ ,147172
+ ,316505
+ ,76
+ ,27259
+ ,109432
+ ,300526
+ ,147
+ ,23201
+ ,1168
+ ,23623
+ ,12
+ ,238
+ ,83248
+ ,195817
+ ,146
+ ,28830
+ ,25162
+ ,61857
+ ,23
+ ,3913
+ ,45724
+ ,163871
+ ,83
+ ,9935
+ ,110529
+ ,428191
+ ,163
+ ,27738
+ ,855
+ ,21054
+ ,4
+ ,338
+ ,101382
+ ,252805
+ ,81
+ ,13326
+ ,14116
+ ,31961
+ ,18
+ ,3988
+ ,89506
+ ,329583
+ ,114
+ ,24347
+ ,135356
+ ,243164
+ ,76
+ ,27111
+ ,116066
+ ,175582
+ ,55
+ ,3938
+ ,144244
+ ,152043
+ ,44
+ ,17416
+ ,8773
+ ,38214
+ ,16
+ ,1888
+ ,102153
+ ,224597
+ ,81
+ ,18700
+ ,117440
+ ,357602
+ ,137
+ ,36809
+ ,104128
+ ,198104
+ ,50
+ ,24959
+ ,134238
+ ,417399
+ ,142
+ ,37343
+ ,134047
+ ,338606
+ ,157
+ ,21849
+ ,279488
+ ,410575
+ ,141
+ ,49809
+ ,79756
+ ,187992
+ ,71
+ ,21654
+ ,66089
+ ,102424
+ ,42
+ ,8728
+ ,102070
+ ,301282
+ ,94
+ ,20920
+ ,146760
+ ,425430
+ ,115
+ ,27195
+ ,154771
+ ,143860
+ ,63
+ ,1037
+ ,165933
+ ,395382
+ ,127
+ ,42570
+ ,64593
+ ,162422
+ ,55
+ ,17672
+ ,92280
+ ,278077
+ ,117
+ ,34245
+ ,67150
+ ,282410
+ ,110
+ ,16786
+ ,128692
+ ,217665
+ ,38
+ ,20954
+ ,124089
+ ,384177
+ ,95
+ ,16378
+ ,125386
+ ,246963
+ ,128
+ ,31852
+ ,37238
+ ,173260
+ ,41
+ ,2805
+ ,140015
+ ,325961
+ ,145
+ ,38086
+ ,150047
+ ,168994
+ ,147
+ ,21166
+ ,154451
+ ,253330
+ ,119
+ ,34672
+ ,156349
+ ,305217
+ ,185
+ ,36171
+ ,6023
+ ,14688
+ ,4
+ ,2065
+ ,84601
+ ,251752
+ ,72
+ ,19354
+ ,68946
+ ,409163
+ ,157
+ ,22124
+ ,1644
+ ,7199
+ ,7
+ ,556
+ ,6179
+ ,46660
+ ,12
+ ,2089
+ ,3926
+ ,17547
+ ,0
+ ,2658
+ ,52789
+ ,116969
+ ,37
+ ,1813
+ ,100350
+ ,206501
+ ,52
+ ,17372)
+ ,dim=c(4
+ ,156)
+ ,dimnames=list(c('TijdInRFC'
+ ,'Characters'
+ ,'Blogs'
+ ,'Revisions')
+ ,1:156))
>  y <- array(NA,dim=c(4,156),dimnames=list(c('TijdInRFC','Characters','Blogs','Revisions'),1:156))
>  for (i in 1:dim(x)[1])
+  {
+  	for (j in 1:dim(x)[2])
+  	{
+  		y[i,j] <- as.numeric(x[i,j])
+  	}
+  }
> par1 = 'pearson'
> 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/wessaorg/rcomp/tmp/1fkvh1324371354.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] 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/wessaorg/rcomp/tmp/2oaym1324371354.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)

	Pearson's product-moment correlation

data:  y[1, ] and y[2, ] 
t = 10.4084, df = 154, p-value < 2.2e-16
alternative hypothesis: true correlation is not equal to 0 
95 percent confidence interval:
 0.5400135 0.7264089 
sample estimates:
      cor 
0.6426225 

> 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
> a<-table.end(a)
> table.save(a,file="/var/wessaorg/rcomp/tmp/325z11324371354.tab") 
> 
> try(system("convert tmp/1fkvh1324371354.ps tmp/1fkvh1324371354.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  0.775   0.118   0.883