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(41 + ,1785 + ,276257 + ,93 + ,492 + ,140824 + ,165 + ,34 + ,1227 + ,180480 + ,63 + ,436 + ,110459 + ,135 + ,44 + ,1944 + ,229345 + ,62 + ,697 + ,105079 + ,121 + ,38 + ,2683 + ,218443 + ,135 + ,1137 + ,112098 + ,148 + ,27 + ,1247 + ,171533 + ,48 + ,380 + ,43929 + ,73 + ,35 + ,631 + ,70849 + ,46 + ,179 + ,76173 + ,49 + ,33 + ,4845 + ,536497 + ,109 + ,2354 + ,187326 + ,185 + ,18 + ,381 + ,33186 + ,46 + ,111 + ,22807 + ,5 + ,34 + ,2066 + ,217320 + ,75 + ,740 + ,144408 + ,125 + ,33 + ,1759 + ,213274 + ,72 + ,595 + ,66485 + ,93 + ,46 + ,2157 + ,310843 + ,82 + ,809 + ,79089 + ,154 + ,57 + ,2211 + ,242788 + ,63 + ,693 + ,81625 + ,98 + ,37 + ,1967 + ,195022 + ,60 + ,738 + ,68788 + ,70 + ,55 + ,3005 + ,367785 + ,117 + ,1184 + ,103297 + ,148 + ,44 + ,2228 + ,261990 + ,50 + ,713 + ,69446 + ,100 + ,62 + ,5009 + ,392509 + ,132 + ,1729 + ,114948 + ,150 + ,40 + ,2353 + ,335528 + ,63 + ,844 + ,167949 + ,197 + ,39 + ,3278 + ,376673 + ,90 + ,1298 + ,125081 + ,114 + ,32 + ,1473 + ,181980 + ,44 + ,514 + ,125818 + ,169 + ,51 + ,2347 + ,266736 + ,64 + ,689 + ,136588 + ,200 + ,49 + ,2523 + ,278265 + ,103 + ,837 + ,112431 + ,148 + ,39 + ,2987 + ,461287 + ,101 + ,1330 + ,103037 + ,140 + ,25 + ,1523 + ,195786 + ,62 + ,491 + ,82317 + ,74 + ,56 + ,1761 + ,197058 + ,69 + ,622 + ,118906 + ,128 + ,45 + ,3660 + ,250284 + ,159 + ,1332 + ,83515 + ,140 + ,38 + ,2949 + ,250092 + ,74 + ,1046 + ,104581 + ,116 + ,45 + ,2914 + ,274121 + ,101 + ,1082 + ,103129 + ,147 + ,43 + ,2009 + ,278027 + ,73 + ,636 + ,83243 + ,132 + ,32 + ,1555 + ,165597 + ,53 + ,586 + ,37110 + ,70 + ,41 + ,3050 + ,371452 + ,140 + ,1170 + ,113344 + ,144 + ,50 + ,2438 + ,296386 + ,54 + ,976 + ,139165 + ,155 + ,50 + ,2092 + ,248320 + ,84 + ,721 + ,86652 + ,165 + ,51 + ,2404 + ,351980 + ,56 + ,863 + ,112302 + ,161 + ,37 + ,1018 + ,101014 + ,40 + ,343 + ,69652 + ,31 + ,44 + ,3462 + ,412722 + ,87 + ,1278 + ,119442 + ,199 + ,42 + ,2764 + ,273950 + ,87 + ,1186 + ,69867 + ,78 + ,44 + ,3710 + ,425531 + ,77 + ,1334 + ,101629 + ,121 + ,36 + ,1947 + ,231912 + ,67 + ,652 + ,70168 + ,112 + ,17 + ,947 + ,115658 + ,36 + ,284 + ,31081 + ,41 + ,43 + ,3609 + ,376008 + ,51 + ,1273 + ,103925 + ,158 + ,41 + ,3324 + ,335039 + ,44 + ,1518 + ,92622 + ,123 + ,41 + ,1842 + ,194127 + ,75 + ,715 + ,79011 + ,104 + ,38 + ,1869 + ,206947 + ,85 + ,671 + ,93487 + ,94 + ,49 + ,1813 + ,182286 + ,97 + ,486 + ,64520 + ,73 + ,45 + ,2496 + ,153778 + ,82 + ,1022 + ,93473 + ,52 + ,42 + ,5557 + ,457592 + ,860 + ,2084 + ,114360 + ,71 + ,26 + ,918 + ,78800 + ,57 + ,330 + ,33032 + ,21 + ,52 + ,2254 + ,208277 + ,80 + ,658 + ,96125 + ,155 + ,50 + ,4103 + ,359144 + ,120 + ,1385 + ,151911 + ,174 + ,45 + ,2241 + ,184648 + ,68 + ,930 + ,89256 + ,136 + ,40 + ,1943 + ,234078 + ,48 + ,620 + ,95676 + ,128 + ,4 + ,496 + ,24188 + ,20 + ,218 + ,5950 + ,7 + ,44 + ,2594 + ,380576 + ,81 + ,840 + ,149695 + ,165 + ,18 + ,744 + ,65029 + ,21 + ,255 + ,32551 + ,21 + ,14 + ,1161 + ,101097 + ,70 + ,454 + ,31701 + ,35 + ,38 + ,3121 + ,288327 + ,128 + ,1149 + ,100087 + ,137 + ,61 + ,2621 + ,334829 + ,86 + ,684 + ,169707 + ,174 + ,39 + ,3906 + ,359400 + ,222 + ,1190 + ,150491 + ,257 + ,42 + ,2722 + ,369577 + ,63 + ,1079 + ,120192 + ,207 + ,40 + ,2312 + ,269961 + ,86 + ,883 + ,95893 + ,103 + ,51 + ,4061 + ,389738 + ,70 + ,1331 + ,151715 + ,171 + ,28 + ,3195 + ,309474 + ,146 + ,1159 + ,176225 + ,279 + ,43 + ,3041 + ,282769 + ,72 + ,1217 + ,59900 + ,83 + ,42 + ,2644 + ,269324 + ,59 + ,946 + ,104767 + ,130 + ,37 + ,1496 + ,234773 + ,58 + ,579 + ,114799 + ,131 + ,30 + ,1835 + ,237561 + ,58 + ,474 + ,72128 + ,126 + ,39 + ,2090 + ,214916 + ,60 + ,630 + ,143592 + ,158 + ,44 + ,2772 + ,240376 + ,89 + ,843 + ,89626 + ,138 + ,36 + ,2440 + ,247447 + ,78 + ,893 + ,131072 + ,200 + ,28 + ,1718 + ,181550 + ,62 + ,633 + ,126817 + ,104 + ,47 + ,2566 + ,242152 + ,67 + ,873 + ,81351 + ,111 + ,23 + ,893 + ,73566 + ,39 + ,385 + ,22618 + ,26 + ,48 + ,2227 + ,246125 + ,60 + ,729 + ,88977 + ,115 + ,38 + ,1878 + ,199565 + ,94 + ,774 + ,92059 + ,127 + ,42 + ,2177 + ,222676 + ,67 + ,769 + ,81897 + ,140 + ,46 + ,2352 + ,363558 + ,96 + ,996 + ,108146 + ,121 + ,37 + ,2981 + ,365442 + ,51 + ,1194 + ,126372 + ,183 + ,41 + ,1926 + ,217036 + ,50 + ,575 + ,249771 + ,68 + ,42 + ,2400 + ,213466 + ,62 + ,725 + ,71154 + ,112 + ,41 + ,2034 + ,204477 + ,71 + ,706 + ,71571 + ,103 + ,36 + ,1800 + ,169080 + ,49 + ,665 + ,55918 + ,63 + ,45 + ,4168 + ,478396 + ,114 + ,1259 + ,160141 + ,166 + ,26 + ,1369 + ,145943 + ,45 + ,653 + ,38692 + ,38 + ,44 + ,2403 + ,288626 + ,57 + ,694 + ,102812 + ,163 + ,8 + ,870 + ,80953 + ,31 + ,437 + ,56622 + ,59 + ,27 + ,1968 + ,150221 + ,175 + ,822 + ,15986 + ,27 + ,38 + ,1569 + ,179317 + ,70 + ,458 + ,123534 + ,108 + ,38 + ,3962 + ,395368 + ,282 + ,1545 + ,108535 + ,88 + ,57 + ,2928 + ,349460 + ,91 + ,987 + ,93879 + ,92 + ,45 + ,3098 + ,180679 + ,72 + ,1051 + ,144551 + ,170 + ,40 + ,2557 + ,286578 + ,62 + ,838 + ,56750 + ,98 + ,42 + ,2329 + ,274477 + ,75 + ,703 + ,127654 + ,205 + ,31 + ,1858 + ,195623 + ,86 + ,613 + ,65594 + ,96 + ,36 + ,3146 + ,282361 + ,89 + ,1128 + ,59938 + ,107 + ,40 + ,2582 + ,329121 + ,138 + ,967 + ,146975 + ,150 + ,40 + ,1824 + ,198658 + ,68 + ,617 + ,165904 + ,138 + ,35 + ,1938 + ,264817 + ,72 + ,654 + ,169265 + ,177 + ,39 + ,2210 + ,300481 + ,61 + ,805 + ,183500 + ,213 + ,65 + ,4119 + ,403404 + ,130 + ,1355 + ,165986 + ,208 + ,33 + ,3844 + ,406613 + ,85 + ,1456 + ,184923 + ,307 + ,51 + ,2791 + ,294371 + ,83 + ,878 + ,140358 + ,125 + ,42 + ,2425 + ,313267 + ,86 + ,887 + ,149959 + ,208 + ,36 + ,2052 + ,189276 + ,116 + ,663 + ,57224 + ,73 + ,19 + ,602 + ,43287 + ,43 + ,214 + ,43750 + ,49 + ,25 + ,2195 + ,189520 + ,87 + ,733 + ,48029 + ,82 + ,44 + ,2413 + ,250254 + ,80 + ,830 + ,104978 + ,206 + ,40 + ,2632 + ,270832 + ,110 + ,1174 + ,100046 + ,112 + ,44 + ,2759 + ,314153 + ,59 + ,1068 + ,101047 + ,139 + ,30 + ,1268 + ,160308 + ,44 + ,413 + ,197426 + ,60 + ,45 + ,3233 + ,162843 + ,85 + ,946 + ,160902 + ,70 + ,42 + ,2025 + ,344925 + ,60 + ,657 + ,147172 + ,112 + ,45 + ,2310 + ,300526 + ,70 + ,690 + ,109432 + ,142 + ,1 + ,398 + ,23623 + ,9 + ,156 + ,1168 + ,11 + ,40 + ,2205 + ,195817 + ,54 + ,779 + ,83248 + ,130 + ,11 + ,530 + ,61857 + ,25 + ,192 + ,25162 + ,31 + ,45 + ,1611 + ,163931 + ,107 + ,461 + ,45724 + ,132 + ,38 + ,3153 + ,428191 + ,63 + ,1213 + ,110529 + ,219 + ,0 + ,387 + ,21054 + ,2 + ,146 + ,855 + ,4 + ,30 + ,2137 + ,252805 + ,67 + ,866 + ,101382 + ,102 + ,8 + ,492 + ,31961 + ,22 + ,200 + ,14116 + ,39 + ,41 + ,3674 + ,335888 + ,155 + ,1290 + ,89506 + ,125 + ,48 + ,2136 + ,246100 + ,79 + ,715 + ,135356 + ,121 + ,48 + ,1748 + ,180591 + ,112 + ,514 + ,116066 + ,42 + ,32 + ,1790 + ,163400 + ,50 + ,697 + ,144244 + ,111 + ,8 + ,568 + ,38214 + ,52 + ,276 + ,8773 + ,16 + ,43 + ,2191 + ,224597 + ,113 + ,752 + ,102153 + ,70 + ,52 + ,2792 + ,357602 + ,115 + ,1021 + ,117440 + ,162 + ,53 + ,1420 + ,202565 + ,78 + ,494 + ,104128 + ,173 + ,49 + ,3718 + ,424398 + ,135 + ,1626 + ,134238 + ,171 + ,48 + ,2554 + ,348017 + ,120 + ,884 + ,134047 + ,172 + ,56 + ,3460 + ,421610 + ,122 + ,1187 + ,279488 + ,254 + ,40 + ,1284 + ,192170 + ,49 + ,488 + ,79756 + ,90 + ,36 + ,1130 + ,102510 + ,56 + ,403 + ,66089 + ,50 + ,44 + ,2974 + ,302158 + ,155 + ,977 + ,102070 + ,113 + ,46 + ,4346 + ,444599 + ,155 + ,1525 + ,146760 + ,187 + ,43 + ,1785 + ,148707 + ,107 + ,551 + ,154771 + ,16 + ,46 + ,4350 + ,407736 + ,146 + ,1807 + ,165933 + ,175 + ,39 + ,1920 + ,164406 + ,76 + ,723 + ,64593 + ,90 + ,41 + ,1923 + ,278077 + ,83 + ,632 + ,92280 + ,140 + ,46 + ,2541 + ,282461 + ,192 + ,898 + ,67150 + ,145 + ,32 + ,1845 + ,219544 + ,69 + ,621 + ,128692 + ,141 + ,45 + ,4081 + ,384177 + ,131 + ,1606 + ,124089 + ,125 + ,39 + ,2339 + ,246963 + ,67 + ,811 + ,125386 + ,241 + ,21 + ,2035 + ,173260 + ,37 + ,716 + ,37238 + ,16 + ,49 + ,3108 + ,336715 + ,60 + ,1001 + ,140015 + ,175 + ,55 + ,1974 + ,176654 + ,127 + ,732 + ,150047 + ,132 + ,36 + ,2651 + ,253341 + ,52 + ,1024 + ,154451 + ,154 + ,48 + ,2702 + ,307133 + ,71 + ,831 + ,156349 + ,198 + ,0 + ,2 + ,1 + ,0 + ,0 + ,0 + ,0 + ,0 + ,207 + ,14688 + ,0 + ,85 + ,6023 + ,5 + ,0 + ,5 + ,98 + ,0 + ,0 + ,0 + ,0 + ,0 + ,8 + ,455 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,43 + ,2255 + ,260901 + ,67 + ,773 + ,84601 + ,125 + ,52 + ,3447 + ,409280 + ,123 + ,1128 + ,68946 + ,174 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,0 + ,4 + ,203 + ,0 + ,0 + ,0 + ,0 + ,0 + ,151 + ,7199 + ,0 + ,74 + ,1644 + ,6 + ,5 + ,474 + ,46660 + ,7 + ,259 + ,6179 + ,13 + ,1 + ,141 + ,17547 + ,3 + ,69 + ,3926 + ,3 + ,45 + ,1111 + ,118589 + ,102 + ,301 + ,52789 + ,35 + ,0 + ,29 + ,969 + ,0 + ,0 + ,0 + ,0 + ,34 + ,2011 + ,233108 + ,53 + ,668 + ,100350 + ,80) + ,dim=c(7 + ,164) + ,dimnames=list(c('reviews' + ,'pgviews' + ,'timerfc' + ,'prview' + ,'compview' + ,'char' + ,'hyper ') + ,1:164)) > y <- array(NA,dim=c(7,164),dimnames=list(c('reviews','pgviews','timerfc','prview','compview','char','hyper '),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 = '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/1wy451324654149.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) Warning messages: 1: In strwidth(labels, "user") : font width unknown for character 0xd 2: In text.default(x, y, txt, cex = cex, font = font) : font width unknown for character 0xd 3: In text.default(x, y, txt, cex = cex, font = font) : font metrics unknown for character 0xd > 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/2aeb81324654149.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 = 8.9405, p-value < 2.2e-16 alternative hypothesis: true tau is not equal to 0 sample estimates: tau 0.4774096 > 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/36jav1324654149.tab") > > try(system("convert tmp/1wy451324654149.ps tmp/1wy451324654149.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.676 0.119 1.785