R version 2.12.0 (2010-10-15) Copyright (C) 2010 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(1418 + ,210907 + ,79 + ,94 + ,24188 + ,2 + ,2172 + ,179321 + ,108 + ,103 + ,32287 + ,4 + ,1583 + ,149061 + ,43 + ,93 + ,27101 + ,0 + ,1764 + ,237213 + ,78 + ,123 + ,19716 + ,0 + ,1495 + ,173326 + ,86 + ,148 + ,17753 + ,-4 + ,1373 + ,133131 + ,44 + ,90 + ,9028 + ,4 + ,2187 + ,258873 + ,104 + ,124 + ,18653 + ,4 + ,4041 + ,324799 + ,158 + ,168 + ,29498 + ,0 + ,1706 + ,230964 + ,102 + ,115 + ,27563 + ,-1 + ,2152 + ,236785 + ,77 + ,71 + ,18293 + ,0 + ,2242 + ,344297 + ,80 + ,108 + ,16116 + ,1 + ,2515 + ,174724 + ,123 + ,120 + ,26569 + ,0 + ,2147 + ,174415 + ,73 + ,114 + ,24785 + ,3 + ,1638 + ,223632 + ,105 + ,120 + ,23825 + ,-1 + ,2452 + ,294424 + ,107 + ,124 + ,34461 + ,4 + ,2662 + ,325107 + ,84 + ,126 + ,24919 + ,3 + ,865 + ,106408 + ,33 + ,37 + ,12558 + ,1 + ,1793 + ,96560 + ,42 + ,38 + ,7784 + ,0 + ,2527 + ,265769 + ,96 + ,120 + ,28522 + ,-2 + ,2747 + ,269651 + ,106 + ,93 + ,22265 + ,-3 + ,1324 + ,149112 + ,56 + ,95 + ,14459 + ,-4 + ,1383 + ,152871 + ,59 + ,90 + ,22240 + ,2 + ,4308 + ,362301 + ,76 + ,110 + ,11912 + ,2 + ,1831 + ,183167 + ,91 + ,138 + ,18220 + ,-4 + ,3373 + ,277965 + ,115 + ,133 + ,19199 + ,3 + ,2352 + ,218946 + ,76 + ,96 + ,25239 + ,2 + ,2144 + ,244052 + ,101 + ,164 + ,29801 + ,2 + ,4691 + ,341570 + ,94 + ,78 + ,18450 + ,0 + ,2694 + ,233328 + ,92 + ,102 + ,34861 + ,5 + ,1769 + ,206161 + ,75 + ,99 + ,16688 + ,-2 + ,3148 + ,311473 + ,128 + ,129 + ,24683 + ,0 + ,1954 + ,207176 + ,56 + ,114 + ,21436 + ,-2 + ,1226 + ,196553 + ,41 + ,99 + ,30546 + ,-3 + ,1496 + ,143246 + ,67 + ,104 + ,15977 + ,2 + ,1943 + ,182192 + ,77 + ,138 + ,14251 + ,2 + ,1762 + ,194979 + ,66 + ,151 + ,16851 + ,2 + ,1403 + ,167488 + ,69 + ,72 + ,21113 + ,0 + ,1425 + ,143756 + ,105 + ,120 + ,17401 + ,4 + ,1857 + ,275541 + ,116 + ,115 + ,23958 + ,4 + ,1420 + ,152299 + ,62 + ,98 + ,14587 + ,2 + ,1644 + ,193339 + ,100 + ,71 + ,20537 + ,2 + ,1054 + ,130585 + ,67 + ,107 + ,30495 + ,-4 + ,937 + ,112611 + ,46 + ,73 + ,7117 + ,3 + ,2547 + ,148446 + ,135 + ,129 + ,33473 + ,3 + ,1626 + ,182079 + ,124 + ,118 + ,21115 + ,2 + ,1964 + ,243060 + ,58 + ,104 + ,32902 + ,-1 + ,1381 + ,162765 + ,68 + ,107 + ,25131 + ,-3 + ,1290 + ,85574 + ,37 + ,36 + ,6943 + ,0 + ,1982 + ,225060 + ,93 + ,139 + ,31808 + ,1 + ,1590 + ,133328 + ,56 + ,56 + ,17014 + ,-3 + ,1281 + ,100750 + ,83 + ,93 + ,6440 + ,3 + ,1272 + ,101523 + ,59 + ,87 + ,18647 + ,0 + ,1944 + ,243511 + ,133 + ,110 + ,20556 + ,0 + ,1605 + ,152474 + ,106 + ,83 + ,22392 + ,0 + ,1386 + ,132487 + ,71 + ,98 + ,8388 + ,3 + ,2395 + ,317394 + ,116 + ,82 + ,22120 + ,-3 + ,2699 + ,244749 + ,98 + ,115 + ,20923 + ,0 + ,1606 + ,184510 + ,64 + ,140 + ,20237 + ,-4 + ,1204 + ,128423 + ,32 + ,120 + ,3769 + ,2 + ,1138 + ,97839 + ,25 + ,66 + ,12252 + ,-1 + ,1111 + ,172494 + ,46 + ,139 + ,21721 + ,3 + ,2186 + ,229242 + ,63 + ,119 + ,17939 + ,2 + ,3604 + ,351619 + ,95 + ,141 + ,23436 + ,5 + ,3261 + ,324598 + ,113 + ,133 + ,34538 + ,2 + ,1641 + ,195838 + ,111 + ,98 + ,25515 + ,-2 + ,2312 + ,254488 + ,120 + ,117 + ,29402 + ,0 + ,2201 + ,199476 + ,87 + ,105 + ,28732 + ,3 + ,961 + ,92499 + ,25 + ,55 + ,5250 + ,-2 + ,1900 + ,224330 + ,131 + ,132 + ,28608 + ,0 + ,1645 + ,181633 + ,47 + ,73 + ,14817 + ,6 + ,2429 + ,271856 + ,109 + ,86 + ,16714 + ,-3 + ,872 + ,95227 + ,37 + ,48 + ,1669 + ,3 + ,1018 + ,98146 + ,15 + ,48 + ,7768 + ,0 + ,1403 + ,118612 + ,54 + ,43 + ,7936 + ,-2 + ,616 + ,65475 + ,16 + ,46 + ,7294 + ,1 + ,1232 + ,108446 + ,22 + ,65 + ,13275 + ,0 + ,1255 + ,121848 + ,37 + ,52 + ,5401 + ,2 + ,995 + ,76302 + ,29 + ,68 + ,8702 + ,2 + ,2048 + ,98104 + ,55 + ,47 + ,8030 + ,-3 + ,301 + ,30989 + ,5 + ,41 + ,1278 + ,-2 + ,628 + ,31774 + ,0 + ,47 + ,1574 + ,1 + ,1597 + ,150580 + ,27 + ,71 + ,9653 + ,-4 + ,717 + ,54157 + ,37 + ,30 + ,7067 + ,0 + ,652 + ,59382 + ,29 + ,24 + ,1514 + ,1 + ,733 + ,84105 + ,17 + ,63 + ,5432 + ,0) + ,dim=c(6 + ,85) + ,dimnames=list(c('pageviews' + ,'time_in_rfc' + ,'blogged_comp.' + ,'feedb.mess.long' + ,'tot.revisions' + ,'Result_test') + ,1:85)) > y <- array(NA,dim=c(6,85),dimnames=list(c('pageviews','time_in_rfc','blogged_comp.','feedb.mess.long','tot.revisions','Result_test'),1:85)) > 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/www/rcomp/tmp/1pyui1324664656.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/www/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/rcomp/createtable") > > n <- length(y[,1]) > n [1] 6 > 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/www/rcomp/tmp/2e3bs1324664656.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) Pearson's product-moment correlation data: y[1, ] and y[2, ] t = 15.1361, df = 83, p-value < 2.2e-16 alternative hypothesis: true correlation is not equal to 0 95 percent confidence interval: 0.7874352 0.9047007 sample estimates: cor 0.8567742 > 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 14 warnings (use warnings() to see them) > a<-table.end(a) > table.save(a,file="/var/www/rcomp/tmp/3ik0c1324664656.tab") > > try(system("convert tmp/1pyui1324664656.ps tmp/1pyui1324664656.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.30 0.05 1.35