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/wessaorg/rcomp/tmp/11opz1324668682.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() #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 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/2uifz1324668682.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) 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) } } a<-table.end(a) table.save(a,file="/var/wessaorg/rcomp/tmp/3xhnq1324668682.tab") try(system("convert tmp/11opz1324668682.ps tmp/11opz1324668682.png",intern=TRUE))