x <- array(list(1773 ,158258 ,90 ,48 ,6200 ,39 ,1704 ,186930 ,58 ,53 ,10265 ,46 ,192 ,7215 ,18 ,0 ,603 ,0 ,2295 ,129098 ,95 ,51 ,8874 ,54 ,3450 ,230632 ,136 ,76 ,20323 ,93 ,6813 ,508313 ,261 ,128 ,26258 ,198 ,1795 ,180745 ,56 ,62 ,10165 ,42 ,1681 ,185559 ,59 ,83 ,8247 ,59 ,1897 ,154581 ,44 ,55 ,8683 ,49 ,2917 ,290658 ,95 ,67 ,16957 ,83 ,1946 ,121844 ,75 ,50 ,8058 ,49 ,2148 ,184039 ,69 ,77 ,20488 ,83 ,1832 ,100324 ,98 ,46 ,7945 ,39 ,3138 ,215855 ,117 ,79 ,13448 ,93 ,1476 ,168265 ,58 ,56 ,5389 ,31 ,1567 ,154647 ,88 ,54 ,6185 ,29 ,1756 ,142018 ,57 ,81 ,24369 ,104 ,1247 ,79030 ,61 ,6 ,70 ,2 ,2779 ,167047 ,87 ,74 ,17327 ,46 ,726 ,27997 ,24 ,13 ,3878 ,27 ,1048 ,73019 ,59 ,22 ,3149 ,16 ,2805 ,241082 ,100 ,99 ,20517 ,108 ,1760 ,195820 ,72 ,38 ,2570 ,36 ,2261 ,141899 ,53 ,59 ,5162 ,33 ,1848 ,145433 ,86 ,50 ,5299 ,46 ,1665 ,183744 ,32 ,50 ,7233 ,65 ,2082 ,202232 ,161 ,61 ,15657 ,80 ,1440 ,199532 ,93 ,87 ,15329 ,81 ,2741 ,354924 ,118 ,60 ,14881 ,69 ,2112 ,192399 ,44 ,52 ,16318 ,69 ,1684 ,182286 ,44 ,61 ,9556 ,37 ,1616 ,181590 ,45 ,60 ,10462 ,45 ,2227 ,133801 ,105 ,53 ,7192 ,62 ,3088 ,233686 ,123 ,76 ,4362 ,33 ,2389 ,219428 ,53 ,63 ,14349 ,77 ,1 ,0 ,1 ,0 ,0 ,0 ,2099 ,223044 ,63 ,54 ,10881 ,34 ,1669 ,100129 ,51 ,44 ,8022 ,44 ,2095 ,136733 ,48 ,36 ,13073 ,43 ,2153 ,249965 ,64 ,83 ,26641 ,117 ,2390 ,242379 ,71 ,105 ,14426 ,125 ,1701 ,145794 ,59 ,37 ,15604 ,49 ,983 ,96404 ,32 ,25 ,9184 ,76 ,2161 ,195891 ,78 ,64 ,5989 ,81 ,1276 ,117156 ,50 ,55 ,11270 ,111 ,1189 ,157787 ,94 ,41 ,13958 ,61 ,745 ,81293 ,32 ,23 ,7162 ,56 ,2231 ,224049 ,100 ,67 ,13275 ,54 ,2242 ,223789 ,87 ,54 ,21224 ,47 ,2639 ,160344 ,59 ,68 ,10615 ,55 ,658 ,48188 ,28 ,12 ,2102 ,14 ,1917 ,161922 ,69 ,99 ,12396 ,44 ,2489 ,294283 ,73 ,74 ,18717 ,115 ,2026 ,235223 ,79 ,56 ,9724 ,57 ,1911 ,195583 ,59 ,67 ,9863 ,48 ,1716 ,146061 ,56 ,40 ,8374 ,40 ,1852 ,208834 ,67 ,53 ,8030 ,51 ,981 ,93764 ,24 ,26 ,7509 ,32 ,1177 ,151985 ,66 ,67 ,14146 ,36 ,2809 ,190545 ,95 ,36 ,7768 ,47 ,1688 ,148922 ,60 ,50 ,13823 ,51 ,2097 ,132856 ,80 ,48 ,7230 ,37 ,1309 ,126107 ,60 ,46 ,10170 ,52 ,1244 ,112718 ,37 ,53 ,7573 ,42 ,1256 ,160930 ,35 ,27 ,5753 ,11 ,1293 ,99184 ,40 ,38 ,9791 ,47 ,2303 ,192535 ,70 ,71 ,19365 ,59 ,2897 ,138708 ,65 ,93 ,9422 ,82 ,1103 ,114408 ,38 ,59 ,12310 ,49 ,340 ,31970 ,15 ,5 ,1283 ,6 ,2791 ,225558 ,112 ,53 ,6372 ,83 ,1333 ,137011 ,71 ,40 ,5413 ,56 ,1441 ,113612 ,68 ,72 ,10837 ,114 ,1623 ,108641 ,71 ,51 ,3394 ,46 ,2650 ,162203 ,67 ,81 ,12964 ,46 ,1499 ,100098 ,44 ,27 ,3495 ,2 ,2302 ,174768 ,60 ,94 ,11580 ,51 ,2540 ,158459 ,97 ,71 ,9970 ,96 ,1000 ,80934 ,30 ,20 ,4911 ,20 ,1234 ,84971 ,71 ,34 ,10138 ,57 ,927 ,80545 ,68 ,54 ,14697 ,49 ,2176 ,287191 ,64 ,49 ,8464 ,51 ,957 ,62974 ,28 ,26 ,4204 ,40 ,1551 ,134091 ,40 ,48 ,10226 ,40 ,1014 ,75555 ,46 ,35 ,3456 ,36 ,1771 ,162154 ,54 ,32 ,8895 ,64 ,2613 ,226638 ,227 ,55 ,22557 ,117 ,1203 ,115019 ,110 ,58 ,6900 ,40 ,1337 ,108749 ,62 ,44 ,8620 ,46 ,1524 ,155537 ,52 ,45 ,7820 ,61 ,1829 ,153133 ,41 ,49 ,12112 ,59 ,2229 ,165618 ,78 ,72 ,13178 ,94 ,1233 ,151517 ,57 ,39 ,7028 ,36 ,1365 ,133686 ,58 ,28 ,6616 ,51 ,950 ,61342 ,40 ,24 ,9570 ,39 ,2319 ,245196 ,117 ,52 ,14612 ,62 ,1857 ,195576 ,70 ,96 ,11219 ,79 ,223 ,19349 ,12 ,13 ,786 ,14 ,2390 ,225371 ,105 ,38 ,11252 ,45 ,1973 ,152796 ,76 ,41 ,9289 ,43 ,700 ,59117 ,29 ,24 ,593 ,8 ,1062 ,91762 ,24 ,54 ,6562 ,41 ,1311 ,136769 ,54 ,68 ,8208 ,25 ,1157 ,114798 ,61 ,28 ,7488 ,22 ,823 ,85338 ,40 ,36 ,4574 ,18 ,596 ,27676 ,22 ,2 ,522 ,3 ,1545 ,153535 ,48 ,91 ,12840 ,54 ,1130 ,122417 ,37 ,29 ,1350 ,6 ,0 ,0 ,0 ,0 ,0 ,0 ,1082 ,91529 ,32 ,46 ,10623 ,50 ,1135 ,107205 ,67 ,25 ,5322 ,33 ,1366 ,144664 ,44 ,51 ,7987 ,54 ,1452 ,136540 ,62 ,59 ,10566 ,63 ,870 ,76656 ,60 ,36 ,1900 ,56 ,78 ,3616 ,5 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,1127 ,183065 ,43 ,40 ,10698 ,49 ,1582 ,144677 ,84 ,68 ,14884 ,90 ,2034 ,159104 ,98 ,28 ,6852 ,51 ,919 ,113273 ,38 ,36 ,6873 ,29 ,778 ,43410 ,19 ,7 ,4 ,1 ,1752 ,175774 ,73 ,70 ,9188 ,68 ,957 ,95401 ,42 ,30 ,5141 ,29 ,1875 ,118893 ,54 ,59 ,4260 ,27 ,731 ,60493 ,40 ,3 ,443 ,4 ,285 ,19764 ,12 ,10 ,2416 ,10 ,1834 ,164062 ,56 ,46 ,9831 ,47 ,1147 ,132696 ,32 ,34 ,5953 ,44 ,1646 ,155367 ,54 ,54 ,9435 ,53 ,256 ,11796 ,9 ,1 ,0 ,0 ,98 ,10674 ,9 ,0 ,0 ,0 ,1404 ,142261 ,57 ,39 ,7642 ,40 ,41 ,6836 ,3 ,0 ,0 ,0 ,1824 ,162563 ,63 ,48 ,6837 ,57 ,42 ,5118 ,3 ,5 ,0 ,0 ,528 ,40248 ,16 ,8 ,775 ,6 ,0 ,0 ,0 ,0 ,0 ,0 ,1073 ,122641 ,47 ,38 ,8191 ,24 ,1305 ,88837 ,38 ,21 ,1661 ,34 ,81 ,7131 ,4 ,0 ,0 ,0 ,261 ,9056 ,14 ,0 ,548 ,10 ,934 ,76611 ,24 ,15 ,3080 ,16 ,1179 ,132697 ,50 ,50 ,13400 ,93 ,1147 ,100681 ,19 ,17 ,8181 ,28) ,dim=c(6 ,144) ,dimnames=list(c('Pages' ,'Seconds' ,'Logins' ,'Blogged' ,'Revisions' ,'Hyperlinks') ,1:144)) y <- array(NA,dim=c(6,144),dimnames=list(c('Pages','Seconds','Logins','Blogged','Revisions','Hyperlinks'),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' #'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/1xyaq1324473874.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/2lopp1324473874.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/3y6ll1324473874.tab") try(system("convert tmp/1xyaq1324473874.ps tmp/1xyaq1324473874.png",intern=TRUE))