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(129404 + ,20 + ,18158 + ,5636 + ,22622 + ,30 + ,28 + ,130358 + ,38 + ,30461 + ,9079 + ,73570 + ,42 + ,39 + ,7215 + ,0 + ,1423 + ,603 + ,1929 + ,0 + ,0 + ,112861 + ,49 + ,25629 + ,8874 + ,36294 + ,54 + ,54 + ,219904 + ,76 + ,48758 + ,17988 + ,62378 + ,86 + ,80 + ,396382 + ,104 + ,129230 + ,21325 + ,167760 + ,157 + ,144 + ,117604 + ,37 + ,27376 + ,8325 + ,52443 + ,36 + ,36 + ,126737 + ,53 + ,26706 + ,7117 + ,57283 + ,48 + ,48 + ,99729 + ,42 + ,26505 + ,7996 + ,36614 + ,45 + ,42 + ,256310 + ,62 + ,49801 + ,14218 + ,93268 + ,77 + ,71 + ,113066 + ,50 + ,46580 + ,6321 + ,35439 + ,49 + ,49 + ,157228 + ,65 + ,48352 + ,19690 + ,72405 + ,77 + ,74 + ,69952 + ,28 + ,13899 + ,5659 + ,24044 + ,28 + ,27 + ,152673 + ,48 + ,39342 + ,11370 + ,55909 + ,84 + ,83 + ,130642 + ,42 + ,27465 + ,4778 + ,44689 + ,31 + ,31 + ,125769 + ,47 + ,55211 + ,5954 + ,49319 + ,28 + ,28 + ,123467 + ,71 + ,74098 + ,22924 + ,62075 + ,99 + ,98 + ,56232 + ,0 + ,13497 + ,70 + ,2341 + ,2 + ,2 + ,108330 + ,50 + ,38338 + ,14369 + ,40551 + ,41 + ,43 + ,22762 + ,12 + ,52505 + ,3706 + ,11621 + ,25 + ,24 + ,48554 + ,16 + ,10663 + ,3147 + ,18741 + ,16 + ,16 + ,182081 + ,77 + ,74484 + ,16801 + ,84202 + ,96 + ,95 + ,140857 + ,29 + ,28895 + ,2162 + ,15334 + ,23 + ,22 + ,93773 + ,38 + ,32827 + ,4721 + ,28024 + ,33 + ,33 + ,133398 + ,50 + ,36188 + ,5290 + ,53306 + ,46 + ,45 + ,113933 + ,33 + ,28173 + ,6446 + ,37918 + ,59 + ,59 + ,153851 + ,49 + ,54926 + ,14711 + ,54819 + ,72 + ,66 + ,140711 + ,59 + ,38900 + ,13311 + ,89058 + ,72 + ,70 + ,303804 + ,55 + ,88530 + ,13577 + ,103354 + ,62 + ,56 + ,161651 + ,40 + ,35482 + ,14634 + ,70239 + ,55 + ,55 + ,123344 + ,40 + ,26730 + ,6931 + ,33045 + ,27 + ,27 + ,157640 + ,51 + ,29806 + ,9992 + ,63852 + ,41 + ,37 + ,91279 + ,41 + ,41799 + ,6185 + ,30905 + ,51 + ,48 + ,189374 + ,73 + ,54289 + ,3445 + ,24242 + ,26 + ,26 + ,178768 + ,51 + ,36805 + ,12327 + ,78907 + ,65 + ,64 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,175403 + ,46 + ,33146 + ,9898 + ,36005 + ,28 + ,21 + ,92342 + ,44 + ,23333 + ,8022 + ,31972 + ,44 + ,44 + ,100023 + ,31 + ,47686 + ,10765 + ,35853 + ,36 + ,36 + ,178277 + ,71 + ,77783 + ,22717 + ,115301 + ,100 + ,89 + ,145062 + ,61 + ,36042 + ,10090 + ,47689 + ,104 + ,101 + ,110980 + ,28 + ,34541 + ,12385 + ,34223 + ,35 + ,31 + ,86039 + ,21 + ,75620 + ,8513 + ,43431 + ,69 + ,65 + ,125481 + ,42 + ,60610 + ,5508 + ,52220 + ,73 + ,71 + ,95535 + ,44 + ,55041 + ,9628 + ,33863 + ,106 + ,102 + ,126456 + ,40 + ,32087 + ,11872 + ,46879 + ,53 + ,53 + ,61554 + ,15 + ,16356 + ,4186 + ,23228 + ,43 + ,41 + ,164752 + ,46 + ,40161 + ,10877 + ,42827 + ,49 + ,46 + ,159121 + ,43 + ,55459 + ,17066 + ,65765 + ,38 + ,37 + ,129362 + ,47 + ,36679 + ,9175 + ,38167 + ,51 + ,51 + ,48188 + ,12 + ,22346 + ,2102 + ,14812 + ,14 + ,14 + ,95461 + ,46 + ,27377 + ,10807 + ,32615 + ,40 + ,40 + ,229864 + ,56 + ,50273 + ,13662 + ,82188 + ,79 + ,77 + ,191094 + ,47 + ,32104 + ,9224 + ,51763 + ,52 + ,51 + ,150640 + ,48 + ,27016 + ,9001 + ,59325 + ,44 + ,43 + ,111388 + ,35 + ,19715 + ,7204 + ,48976 + ,34 + ,33 + ,165098 + ,44 + ,33629 + ,6572 + ,43384 + ,47 + ,47 + ,63205 + ,25 + ,27084 + ,7509 + ,26692 + ,32 + ,31 + ,109102 + ,47 + ,32352 + ,12920 + ,53279 + ,31 + ,31 + ,137303 + ,28 + ,51845 + ,5438 + ,20652 + ,40 + ,40 + ,125304 + ,48 + ,26591 + ,11489 + ,38338 + ,42 + ,42 + ,85332 + ,32 + ,29677 + ,6661 + ,36735 + ,34 + ,35 + ,95808 + ,28 + ,54237 + ,7941 + ,42764 + ,40 + ,40 + ,83419 + ,31 + ,20284 + ,6173 + ,44331 + ,35 + ,30 + ,101723 + ,13 + ,22741 + ,5562 + ,41354 + ,11 + ,11 + ,94982 + ,38 + ,34178 + ,9492 + ,47879 + ,43 + ,41 + ,129700 + ,39 + ,69551 + ,17456 + ,103793 + ,53 + ,53 + ,113325 + ,68 + ,29653 + ,9422 + ,52235 + ,82 + ,82 + ,81518 + ,32 + ,38071 + ,10913 + ,49825 + ,41 + ,41 + ,31970 + ,5 + ,4157 + ,1283 + ,4105 + ,6 + ,6 + ,192268 + ,53 + ,28321 + ,6198 + ,58687 + ,82 + ,81 + ,91086 + ,33 + ,40195 + ,4501 + ,40745 + ,47 + ,47 + ,80820 + ,54 + ,48158 + ,9560 + ,33187 + ,108 + ,100 + ,83261 + ,36 + ,13310 + ,3394 + ,14063 + ,46 + ,46 + ,116290 + ,52 + ,78474 + ,9871 + ,37407 + ,38 + ,38 + ,56544 + ,0 + ,6386 + ,2419 + ,7190 + ,0 + ,0 + ,116173 + ,52 + ,31588 + ,10630 + ,49562 + ,45 + ,45 + ,111488 + ,45 + ,61254 + ,8536 + ,76324 + ,57 + ,56 + ,60138 + ,16 + ,21152 + ,4911 + ,21928 + ,20 + ,18 + ,73422 + ,33 + ,41272 + ,9775 + ,27860 + ,56 + ,54 + ,67751 + ,48 + ,34165 + ,11227 + ,28078 + ,38 + ,37 + ,213351 + ,33 + ,37054 + ,6916 + ,49577 + ,42 + ,40 + ,51185 + ,24 + ,12368 + ,3424 + ,28145 + ,37 + ,37 + ,97181 + ,37 + ,23168 + ,8637 + ,36241 + ,36 + ,36 + ,45100 + ,17 + ,16380 + ,3189 + ,10824 + ,34 + ,34 + ,115801 + ,32 + ,41242 + ,8178 + ,46892 + ,53 + ,49 + ,186310 + ,55 + ,48450 + ,16739 + ,61264 + ,85 + ,82 + ,71960 + ,39 + ,20790 + ,6094 + ,22933 + ,36 + ,36 + ,80105 + ,31 + ,34585 + ,7237 + ,20787 + ,33 + ,33 + ,103613 + ,26 + ,35672 + ,7355 + ,43978 + ,57 + ,55 + ,98707 + ,37 + ,52168 + ,9734 + ,51305 + ,50 + ,50 + ,136234 + ,66 + ,53933 + ,11225 + ,55593 + ,71 + ,71 + ,136781 + ,35 + ,34474 + ,6213 + ,51648 + ,32 + ,31 + ,105863 + ,24 + ,43753 + ,4875 + ,30552 + ,45 + ,42 + ,38775 + ,18 + ,36456 + ,8159 + ,23470 + ,33 + ,31 + ,179997 + ,37 + ,51183 + ,11893 + ,77530 + ,53 + ,51 + ,169406 + ,86 + ,52742 + ,10754 + ,57299 + ,64 + ,64 + ,19349 + ,13 + ,3895 + ,786 + ,9604 + ,14 + ,14 + ,153069 + ,21 + ,37076 + ,9706 + ,34684 + ,38 + ,37 + ,109510 + ,32 + ,24079 + ,7796 + ,41094 + ,39 + ,37 + ,43803 + ,8 + ,2325 + ,593 + ,3439 + ,8 + ,8 + ,47062 + ,38 + ,29354 + ,5600 + ,25171 + ,38 + ,38 + ,110845 + ,45 + ,30341 + ,7245 + ,23437 + ,24 + ,23 + ,92517 + ,24 + ,18992 + ,7360 + ,34086 + ,22 + ,22 + ,58660 + ,23 + ,15292 + ,4574 + ,24649 + ,18 + ,18 + ,27676 + ,2 + ,5842 + ,522 + ,2342 + ,3 + ,1 + ,98550 + ,52 + ,28918 + ,10905 + ,45571 + ,49 + ,48 + ,43646 + ,5 + ,3738 + ,999 + ,3255 + ,5 + ,5 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,67312 + ,43 + ,95352 + ,9016 + ,30002 + ,47 + ,46 + ,57359 + ,18 + ,37478 + ,5134 + ,19360 + ,33 + ,33 + ,104330 + ,44 + ,26839 + ,6608 + ,43320 + ,44 + ,41 + ,70369 + ,45 + ,26783 + ,8577 + ,35513 + ,56 + ,57 + ,65494 + ,29 + ,33392 + ,1543 + ,23536 + ,49 + ,49 + ,3616 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,143931 + ,32 + ,25446 + ,9803 + ,54438 + ,45 + ,45 + ,117946 + ,65 + ,59847 + ,12140 + ,56812 + ,78 + ,78 + ,131175 + ,26 + ,28162 + ,6678 + ,33838 + ,51 + ,46 + ,84336 + ,24 + ,33298 + ,6420 + ,32366 + ,25 + ,25 + ,43410 + ,7 + ,2781 + ,4 + ,13 + ,1 + ,1 + ,136250 + ,62 + ,37121 + ,7979 + ,55082 + ,62 + ,59 + ,79015 + ,30 + ,22698 + ,5141 + ,31334 + ,29 + ,29 + ,92937 + ,49 + ,27615 + ,1311 + ,16612 + ,26 + ,26 + ,57586 + ,3 + ,32689 + ,443 + ,5084 + ,4 + ,4 + ,19764 + ,10 + ,5752 + ,2416 + ,9927 + ,10 + ,10 + ,105757 + ,42 + ,23164 + ,8396 + ,47413 + ,43 + ,43 + ,97213 + ,18 + ,20304 + ,5462 + ,27389 + ,36 + ,36 + ,113402 + ,40 + ,34409 + ,7271 + ,30425 + ,43 + ,41 + ,11796 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,7627 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,121085 + ,29 + ,92538 + ,4423 + ,33510 + ,33 + ,32 + ,6836 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,139563 + ,46 + ,46037 + ,5331 + ,40389 + ,53 + ,53 + ,5118 + ,5 + ,0 + ,0 + ,0 + ,0 + ,0 + ,40248 + ,8 + ,5444 + ,775 + ,6012 + ,6 + ,6 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,95079 + ,21 + ,23924 + ,6676 + ,22205 + ,19 + ,18 + ,80763 + ,21 + ,52230 + ,1489 + ,17231 + ,26 + ,26 + ,7131 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,4194 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60378 + ,15 + ,8019 + ,3080 + ,11017 + ,16 + ,16 + ,109173 + ,47 + ,34542 + ,11409 + ,46741 + ,84 + ,84 + ,83484 + ,17 + ,21157 + ,6769 + ,39869 + ,28 + ,22) + ,dim=c(7 + ,144) + ,dimnames=list(c('timeRFC' + ,'bloggcomp' + ,'characters' + ,'revisions' + ,'seconds' + ,'inclhyperlinks' + ,'inclblogs') + ,1:144)) > y <- array(NA,dim=c(7,144),dimnames=list(c('timeRFC','bloggcomp','characters','revisions','seconds','inclhyperlinks','inclblogs'),1:144)) > 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/wessaorg/rcomp/tmp/1a9pt1323976965.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] 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) + } > a<-table.end(a) > table.save(a,file="/var/wessaorg/rcomp/tmp/2ha041323976965.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 = 10.5078, p-value < 2.2e-16 alternative hypothesis: true tau is not equal to 0 sample estimates: tau 0.5962786 > 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 21 warnings (use warnings() to see them) > a<-table.end(a) > table.save(a,file="/var/wessaorg/rcomp/tmp/3zz011323976965.tab") > > try(system("convert tmp/1a9pt1323976965.ps tmp/1a9pt1323976965.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.556 0.110 1.655