x <- array(list(140824 ,96 ,275243 ,151 ,71 ,32033 ,110459 ,71 ,180141 ,133 ,66 ,20654 ,105079 ,70 ,228925 ,170 ,75 ,16346 ,112098 ,134 ,218443 ,148 ,104 ,35926 ,43929 ,67 ,171533 ,88 ,52 ,10621 ,76173 ,8 ,70849 ,129 ,28 ,10024 ,187326 ,160 ,532492 ,128 ,125 ,43068 ,22807 ,1 ,33186 ,67 ,19 ,1271 ,144408 ,83 ,217320 ,132 ,59 ,34416 ,66485 ,82 ,213274 ,120 ,44 ,20318 ,79089 ,92 ,309323 ,169 ,111 ,24409 ,81625 ,117 ,242739 ,218 ,122 ,20648 ,68788 ,56 ,194882 ,122 ,76 ,12347 ,103297 ,139 ,364569 ,191 ,81 ,21857 ,69446 ,80 ,255231 ,162 ,86 ,11034 ,114948 ,176 ,391748 ,227 ,181 ,33433 ,167949 ,114 ,334118 ,156 ,75 ,35902 ,125081 ,105 ,374892 ,144 ,163 ,22355 ,125818 ,103 ,176082 ,111 ,56 ,31219 ,136588 ,135 ,266736 ,199 ,87 ,21983 ,112431 ,123 ,278265 ,187 ,70 ,40085 ,103037 ,87 ,442703 ,144 ,104 ,18507 ,82317 ,66 ,180393 ,89 ,42 ,16278 ,118906 ,103 ,189897 ,208 ,55 ,24662 ,83515 ,149 ,242176 ,165 ,125 ,31452 ,104581 ,113 ,241883 ,146 ,129 ,32580 ,103129 ,99 ,267318 ,158 ,129 ,22883 ,83243 ,117 ,277501 ,154 ,88 ,27652 ,37110 ,57 ,155915 ,117 ,52 ,9845 ,113344 ,141 ,365373 ,158 ,72 ,20190 ,139165 ,130 ,291907 ,195 ,80 ,46201 ,86652 ,50 ,247188 ,186 ,81 ,10971 ,112302 ,140 ,342856 ,197 ,95 ,34811 ,69652 ,43 ,101014 ,141 ,39 ,3029 ,119442 ,141 ,409551 ,168 ,100 ,38941 ,69867 ,83 ,273950 ,159 ,56 ,4958 ,101629 ,112 ,425253 ,161 ,122 ,32344 ,70168 ,79 ,227636 ,139 ,84 ,19433 ,31081 ,33 ,115658 ,55 ,33 ,12558 ,103925 ,137 ,375992 ,166 ,201 ,36524 ,92622 ,125 ,334004 ,151 ,82 ,26041 ,79011 ,76 ,186648 ,148 ,71 ,16637 ,93487 ,78 ,206196 ,115 ,76 ,28395 ,64520 ,68 ,182286 ,181 ,64 ,16747 ,93473 ,50 ,153778 ,73 ,82 ,9105 ,114360 ,101 ,455401 ,147 ,153 ,11941 ,33032 ,20 ,78800 ,82 ,42 ,7935 ,96125 ,101 ,208277 ,201 ,77 ,19499 ,151911 ,149 ,358127 ,193 ,121 ,22938 ,89256 ,108 ,182715 ,164 ,60 ,25314 ,95676 ,95 ,224649 ,158 ,76 ,28527 ,5950 ,8 ,24188 ,12 ,24 ,2694 ,149695 ,88 ,380576 ,163 ,326 ,20867 ,32551 ,21 ,65029 ,67 ,17 ,3597 ,31701 ,30 ,101097 ,52 ,64 ,5296 ,100087 ,97 ,279128 ,134 ,60 ,32982 ,169707 ,130 ,328024 ,230 ,86 ,38975 ,150491 ,132 ,359138 ,145 ,198 ,42721 ,120192 ,161 ,369539 ,153 ,145 ,41455 ,95893 ,89 ,266240 ,155 ,86 ,23923 ,151715 ,160 ,389724 ,198 ,147 ,26719 ,176225 ,139 ,305351 ,101 ,119 ,53405 ,59900 ,104 ,281929 ,169 ,118 ,12526 ,104767 ,92 ,264889 ,163 ,89 ,26584 ,114799 ,52 ,229585 ,139 ,68 ,37062 ,72128 ,117 ,235052 ,116 ,68 ,25696 ,143592 ,93 ,209030 ,145 ,135 ,24634 ,89626 ,85 ,240070 ,167 ,151 ,27269 ,131072 ,143 ,243937 ,135 ,84 ,25270 ,126817 ,143 ,175816 ,102 ,70 ,24634 ,81351 ,99 ,239337 ,173 ,70 ,17828 ,22618 ,22 ,73566 ,88 ,32 ,3007 ,88977 ,78 ,242622 ,175 ,86 ,20065 ,92059 ,83 ,196387 ,133 ,55 ,24648 ,81897 ,131 ,209049 ,148 ,64 ,21588 ,108146 ,140 ,363250 ,169 ,88 ,25217 ,126372 ,148 ,362030 ,143 ,96 ,30927 ,249771 ,80 ,217036 ,154 ,108 ,18487 ,71154 ,133 ,212643 ,164 ,66 ,18050 ,71571 ,97 ,202011 ,146 ,67 ,17696 ,55918 ,62 ,153654 ,109 ,48 ,17326 ,160141 ,161 ,475859 ,175 ,128 ,39361 ,38692 ,30 ,145943 ,99 ,69 ,9648 ,102812 ,118 ,287359 ,122 ,106 ,26759 ,56622 ,49 ,80953 ,28 ,25 ,7905 ,15986 ,52 ,150216 ,101 ,55 ,4527 ,123534 ,76 ,179317 ,139 ,60 ,41517 ,108535 ,85 ,394546 ,143 ,213 ,21261 ,93879 ,146 ,345373 ,206 ,120 ,36099 ,144551 ,165 ,179811 ,171 ,106 ,39039 ,56750 ,84 ,283459 ,150 ,100 ,13841 ,127654 ,165 ,274449 ,154 ,81 ,23841 ,65594 ,48 ,190312 ,114 ,65 ,8589 ,59938 ,149 ,280628 ,140 ,77 ,15049 ,146975 ,75 ,329007 ,156 ,85 ,39038 ,166616 ,83 ,194106 ,156 ,40 ,36918 ,168553 ,110 ,259365 ,127 ,59 ,39932 ,183500 ,164 ,297464 ,141 ,83 ,43840 ,165986 ,154 ,399815 ,251 ,155 ,43146 ,184923 ,165 ,402133 ,126 ,113 ,50099 ,140358 ,121 ,291970 ,198 ,138 ,40312 ,149959 ,150 ,296670 ,155 ,67 ,32616 ,57224 ,73 ,189276 ,138 ,188 ,11338 ,43750 ,13 ,43287 ,71 ,14 ,7409 ,48029 ,89 ,185969 ,84 ,83 ,18213 ,104978 ,105 ,250254 ,167 ,148 ,45873 ,100046 ,129 ,268391 ,155 ,50 ,39844 ,101047 ,169 ,314131 ,162 ,92 ,28317 ,197426 ,28 ,156967 ,112 ,81 ,24797 ,160902 ,118 ,161884 ,168 ,97 ,7471 ,147172 ,79 ,334485 ,157 ,74 ,27259 ,109432 ,147 ,300526 ,164 ,89 ,23201 ,1168 ,12 ,23623 ,0 ,11 ,238 ,83248 ,146 ,195817 ,155 ,73 ,28830 ,25162 ,23 ,61857 ,32 ,25 ,3913 ,45724 ,83 ,163871 ,169 ,50 ,9935 ,110529 ,163 ,428191 ,140 ,120 ,27738 ,855 ,4 ,21054 ,0 ,16 ,338 ,101382 ,81 ,252805 ,111 ,52 ,13326 ,14116 ,18 ,31961 ,25 ,22 ,3988 ,89506 ,114 ,333052 ,152 ,119 ,24347 ,135356 ,76 ,243223 ,183 ,68 ,27111 ,116066 ,55 ,177264 ,181 ,93 ,3938 ,144244 ,44 ,152043 ,107 ,54 ,17416 ,8773 ,16 ,38214 ,27 ,34 ,1888 ,102153 ,81 ,224597 ,163 ,43 ,18700 ,117440 ,137 ,357602 ,198 ,82 ,36809 ,104128 ,50 ,198104 ,205 ,61 ,24959 ,134238 ,142 ,423741 ,187 ,85 ,37343 ,134047 ,157 ,338606 ,187 ,99 ,21849 ,279488 ,141 ,417175 ,210 ,129 ,49809 ,79756 ,71 ,190323 ,151 ,37 ,21654 ,66089 ,42 ,102424 ,131 ,43 ,8728 ,102070 ,94 ,302158 ,171 ,347 ,20920 ,146760 ,117 ,437141 ,172 ,182 ,27195 ,154771 ,63 ,146250 ,164 ,56 ,1037 ,165933 ,127 ,395382 ,172 ,134 ,42570 ,64593 ,55 ,162575 ,143 ,79 ,17672 ,92280 ,117 ,278077 ,151 ,49 ,34245 ,67150 ,110 ,282410 ,158 ,96 ,16786 ,128692 ,39 ,219493 ,125 ,115 ,20954 ,124089 ,95 ,384177 ,169 ,122 ,16378 ,125386 ,128 ,246963 ,145 ,92 ,31852 ,37238 ,41 ,173260 ,79 ,63 ,2805 ,140015 ,146 ,333967 ,190 ,103 ,38086 ,150047 ,147 ,168994 ,192 ,58 ,21166 ,154451 ,119 ,253330 ,132 ,87 ,34672 ,156349 ,185 ,305217 ,168 ,109 ,36171 ,0 ,0 ,1 ,0 ,0 ,0 ,6023 ,4 ,14688 ,0 ,10 ,2065 ,0 ,0 ,98 ,0 ,1 ,0 ,0 ,0 ,455 ,0 ,2 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,84601 ,75 ,260345 ,133 ,88 ,19354 ,68946 ,157 ,409163 ,204 ,162 ,22124 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,203 ,0 ,4 ,0 ,1644 ,7 ,7199 ,0 ,5 ,556 ,6179 ,12 ,46660 ,15 ,20 ,2089 ,3926 ,0 ,17547 ,4 ,5 ,2658 ,52789 ,37 ,116969 ,152 ,45 ,1813 ,0 ,0 ,969 ,0 ,2 ,0 ,100350 ,59 ,229447 ,125 ,70 ,17372) ,dim=c(6 ,164) ,dimnames=list(c('#characters' ,'bloggedcomputations' ,'totalseconds' ,'PeerReviews' ,'#logins' ,'#revisions') ,1:164)) y <- array(NA,dim=c(6,164),dimnames=list(c('#characters','bloggedcomputations','totalseconds','PeerReviews','#logins','#revisions'),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 = '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/1ivzt1324374445.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/2ww151324374445.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/39qyg1324374445.tab") try(system("convert tmp/1ivzt1324374445.ps tmp/1ivzt1324374445.png",intern=TRUE))