R version 2.15.2 (2012-10-26) -- "Trick or Treat" Copyright (C) 2012 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i686-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(10.81 + ,24563400 + ,-0.2643 + ,24.45 + ,2772.73 + ,0.0373 + ,115.7 + ,5.98 + ,9.12 + ,14163200 + ,-0.2643 + ,23.62 + ,2151.83 + ,0.0353 + ,109.2 + ,5.49 + ,11.03 + ,18184800 + ,-0.2643 + ,21.90 + ,1840.26 + ,0.0292 + ,116.9 + ,5.31 + ,12.74 + ,20810300 + ,-0.1918 + ,27.12 + ,2116.24 + ,0.0327 + ,109.9 + ,4.8 + ,9.98 + ,12843000 + ,-0.1918 + ,27.70 + ,2110.49 + ,0.0362 + ,116.1 + ,4.21 + ,11.62 + ,13866700 + ,-0.1918 + ,29.23 + ,2160.54 + ,0.0325 + ,118.9 + ,3.97 + ,9.40 + ,15119200 + ,-0.2246 + ,26.50 + ,2027.13 + ,0.0272 + ,116.3 + ,3.77 + ,9.27 + ,8301600 + ,-0.2246 + ,22.84 + ,1805.43 + ,0.0272 + ,114.0 + ,3.65 + ,7.76 + ,14039600 + ,-0.2246 + ,20.49 + ,1498.80 + ,0.0265 + ,97.0 + ,3.07 + ,8.78 + ,12139700 + ,0.3654 + ,23.28 + ,1690.20 + ,0.0213 + ,85.3 + ,2.49 + ,10.65 + ,9649000 + ,0.3654 + ,25.71 + ,1930.58 + ,0.019 + ,84.9 + ,2.09 + ,10.95 + ,8513600 + ,0.3654 + ,26.52 + ,1950.40 + ,0.0155 + ,94.6 + ,1.82 + ,12.36 + ,15278600 + ,0.0447 + ,25.51 + ,1934.03 + ,0.0114 + ,97.8 + ,1.73 + ,10.85 + ,15590900 + ,0.0447 + ,23.36 + ,1731.49 + ,0.0114 + ,95.0 + ,1.74 + ,11.84 + ,9691100 + ,0.0447 + ,24.15 + ,1845.35 + ,0.0148 + ,110.7 + ,1.73 + ,12.14 + ,10882700 + ,-0.0312 + ,20.92 + ,1688.23 + ,0.0164 + ,108.5 + ,1.75 + ,11.65 + ,10294800 + ,-0.0312 + ,20.38 + ,1615.73 + ,0.0118 + ,110.3 + ,1.75 + ,8.86 + ,16031900 + ,-0.0312 + ,21.90 + ,1463.21 + ,0.0107 + ,106.3 + ,1.75 + ,7.63 + ,13683600 + ,-0.0048 + ,19.21 + ,1328.26 + ,0.0146 + ,97.4 + ,1.73 + ,7.38 + ,8677200 + ,-0.0048 + ,19.65 + ,1314.85 + ,0.018 + ,94.5 + ,1.74 + ,7.25 + ,9874100 + ,-0.0048 + ,17.51 + ,1172.06 + ,0.0151 + ,93.7 + ,1.75 + ,8.03 + ,10725500 + ,0.0705 + ,21.41 + ,1329.75 + ,0.0203 + ,79.6 + ,1.75 + ,7.75 + ,8348400 + ,0.0705 + ,23.09 + ,1478.78 + ,0.022 + ,84.9 + ,1.34 + ,7.16 + ,8046200 + ,0.0705 + ,20.70 + ,1335.51 + ,0.0238 + ,80.7 + ,1.24 + ,7.18 + ,10862300 + ,-0.0134 + ,19.00 + ,1320.91 + ,0.026 + ,78.8 + ,1.24 + ,7.51 + ,8100300 + ,-0.0134 + ,19.04 + ,1337.52 + ,0.0298 + ,64.8 + ,1.26 + ,7.07 + ,7287500 + ,-0.0134 + ,19.45 + ,1341.17 + ,0.0302 + ,61.4 + ,1.25 + ,7.11 + ,14002500 + ,0.0812 + ,20.54 + ,1464.31 + ,0.0222 + ,81.0 + ,1.26 + ,8.98 + ,19037900 + ,0.0812 + ,19.77 + ,1595.91 + ,0.0206 + ,83.6 + ,1.26 + ,9.53 + ,10774600 + ,0.0812 + ,20.60 + ,1622.80 + ,0.0211 + ,83.5 + ,1.22 + ,10.54 + ,8960600 + ,0.1885 + ,21.21 + ,1735.02 + ,0.0211 + ,77.0 + ,1.01 + ,11.31 + ,7773300 + ,0.1885 + ,21.30 + ,1810.45 + ,0.0216 + ,81.7 + ,1.03 + ,10.36 + ,9579700 + ,0.1885 + ,22.33 + ,1786.94 + ,0.0232 + ,77.0 + ,1.01 + ,11.44 + ,11270700 + ,0.3628 + ,21.12 + ,1932.21 + ,0.0204 + ,81.7 + ,1.01 + ,10.45 + ,9492800 + ,0.3628 + ,20.77 + ,1960.26 + ,0.0177 + ,92.5 + ,1 + ,10.69 + ,9136800 + ,0.3628 + ,22.11 + ,2003.37 + ,0.0188 + ,91.7 + ,0.98 + ,11.28 + ,14487600 + ,0.2942 + ,22.34 + ,2066.15 + ,0.0193 + ,96.4 + ,1 + ,11.96 + ,10133200 + ,0.2942 + ,21.43 + ,2029.82 + ,0.0169 + ,88.5 + ,1.01 + ,13.52 + ,18659700 + ,0.2942 + ,20.14 + ,1994.22 + ,0.0174 + ,88.5 + ,1 + ,12.89 + ,15980700 + ,0.3036 + ,21.11 + ,1920.15 + ,0.0229 + ,93.0 + ,1 + ,14.03 + ,9732100 + ,0.3036 + ,21.19 + ,1986.74 + ,0.0305 + ,93.1 + ,1 + ,16.27 + ,14626300 + ,0.3036 + ,23.07 + ,2047.79 + ,0.0327 + ,102.8 + ,1.03 + ,16.17 + ,16904000 + ,0.3703 + ,23.01 + ,1887.36 + ,0.0299 + ,105.7 + ,1.26 + ,17.25 + ,13616700 + ,0.3703 + ,22.12 + ,1838.10 + ,0.0265 + ,98.7 + ,1.43 + ,19.38 + ,13772900 + ,0.3703 + ,22.40 + ,1896.84 + ,0.0254 + ,96.7 + ,1.61 + ,26.20 + ,28749200 + ,0.7398 + ,22.66 + ,1974.99 + ,0.0319 + ,92.9 + ,1.76 + ,33.53 + ,31408300 + ,0.7398 + ,24.21 + ,2096.81 + ,0.0352 + ,92.6 + ,1.93 + ,32.20 + ,26342800 + ,0.7398 + ,24.13 + ,2175.44 + ,0.0326 + ,102.7 + ,2.16 + ,38.45 + ,48909500 + ,0.6988 + ,23.73 + ,2062.41 + ,0.0297 + ,105.1 + ,2.28 + ,44.86 + ,41542400 + ,0.6988 + ,22.79 + ,2051.72 + ,0.0301 + ,104.4 + ,2.5 + ,41.67 + ,24857200 + ,0.6988 + ,21.89 + ,1999.23 + ,0.0315 + ,103.0 + ,2.63 + ,36.06 + ,34093700 + ,0.7478 + ,22.92 + ,1921.65 + ,0.0351 + ,97.5 + ,2.79 + ,39.76 + ,22555200 + ,0.7478 + ,23.44 + ,2068.22 + ,0.028 + ,103.1 + ,3 + ,36.81 + ,19067500 + ,0.7478 + ,22.57 + ,2056.96 + ,0.0253 + ,106.2 + ,3.04 + ,42.65 + ,19029100 + ,0.5651 + ,23.27 + ,2184.83 + ,0.0317 + ,103.6 + ,3.26 + ,46.89 + ,15223200 + ,0.5651 + ,24.95 + ,2152.09 + ,0.0364 + ,105.5 + ,3.5 + ,53.61 + ,21903700 + ,0.5651 + ,23.45 + ,2151.69 + ,0.0469 + ,87.5 + ,3.62 + ,57.59 + ,33306600 + ,0.6473 + ,23.42 + ,2120.30 + ,0.0435 + ,85.2 + ,3.78 + ,67.82 + ,23898100 + ,0.6473 + ,25.30 + ,2232.82 + ,0.0346 + ,98.3 + ,4 + ,71.89 + ,23279600 + ,0.6473 + ,23.90 + ,2205.32 + ,0.0342 + ,103.8 + ,4.16 + ,75.51 + ,40699800 + ,0.3441 + ,25.73 + ,2305.82 + ,0.0399 + ,106.8 + ,4.29 + ,68.49 + ,37646000 + ,0.3441 + ,24.64 + ,2281.39 + ,0.036 + ,102.7 + ,4.49 + ,62.72 + ,37277000 + ,0.3441 + ,24.95 + ,2339.79 + ,0.0336 + ,107.5 + ,4.59 + ,70.39 + ,39246800 + ,0.2415 + ,22.15 + ,2322.57 + ,0.0355 + ,109.8 + ,4.79 + ,59.77 + ,27418400 + ,0.2415 + ,20.85 + ,2178.88 + ,0.0417 + ,104.7 + ,4.94 + ,57.27 + ,30318700 + ,0.2415 + ,21.45 + ,2172.09 + ,0.0432 + ,105.7 + ,4.99 + ,67.96 + ,32808100 + ,0.3151 + ,22.15 + ,2091.47 + ,0.0415 + ,107.0 + ,5.24 + ,67.85 + ,28668200 + ,0.3151 + ,23.75 + ,2183.75 + ,0.0382 + ,100.2 + ,5.25 + ,76.98 + ,32370300 + ,0.3151 + ,25.27 + ,2258.43 + ,0.0206 + ,105.9 + ,5.25 + ,81.08 + ,24171100 + ,0.239 + ,26.53 + ,2366.71 + ,0.0131 + ,105.1 + ,5.25 + ,91.66 + ,25009100 + ,0.239 + ,27.22 + ,2431.77 + ,0.0197 + ,105.3 + ,5.25 + ,84.84 + ,32084300 + ,0.239 + ,27.69 + ,2415.29 + ,0.0254 + ,110.0 + ,5.24 + ,85.73 + ,50117500 + ,0.2127 + ,28.61 + ,2463.93 + ,0.0208 + ,110.2 + ,5.25 + ,84.61 + ,27522200 + ,0.2127 + ,26.21 + ,2416.15 + ,0.0242 + ,111.2 + ,5.26 + ,92.91 + ,26816800 + ,0.2127 + ,25.93 + ,2421.64 + ,0.0278 + ,108.2 + ,5.26 + ,99.80 + ,25136100 + ,0.273 + ,27.86 + ,2525.09 + ,0.0257 + ,106.3 + ,5.25 + ,121.19 + ,30295600 + ,0.273 + ,28.65 + ,2604.52 + ,0.0269 + ,108.5 + ,5.25 + ,122.04 + ,41526100 + ,0.273 + ,27.51 + ,2603.23 + ,0.0269 + ,105.3 + ,5.25 + ,131.76 + ,43845100 + ,0.3657 + ,27.06 + ,2546.27 + ,0.0236 + ,111.9 + ,5.26 + ,138.48 + ,39188900 + ,0.3657 + ,26.91 + ,2596.36 + ,0.0197 + ,105.6 + ,5.02 + ,153.47 + ,40496400 + ,0.3657 + ,27.60 + ,2701.50 + ,0.0276 + ,99.5 + ,4.94 + ,189.95 + ,37438400 + ,0.4643 + ,34.48 + ,2859.12 + ,0.0354 + ,95.2 + ,4.76 + ,182.22 + ,46553700 + ,0.4643 + ,31.58 + ,2660.96 + ,0.0431 + ,87.8 + ,4.49 + ,198.08 + ,31771400 + ,0.4643 + ,33.46 + ,2652.28 + ,0.0408 + ,90.6 + ,4.24 + ,135.36 + ,62108100 + ,0.5096 + ,30.64 + ,2389.86 + ,0.0428 + ,87.9 + ,3.94 + ,125.02 + ,46645400 + ,0.5096 + ,25.66 + ,2271.48 + ,0.0403 + ,76.4 + ,2.98 + ,143.50 + ,42313100 + ,0.5096 + ,26.78 + ,2279.10 + ,0.0398 + ,65.9 + ,2.61 + ,173.95 + ,38841700 + ,0.3592 + ,26.91 + ,2412.80 + ,0.0394 + ,62.3 + ,2.28 + ,188.75 + ,32650300 + ,0.3592 + ,26.82 + ,2522.66 + ,0.0418 + ,57.2 + ,1.98 + ,167.44 + ,34281100 + ,0.3592 + ,26.05 + ,2292.98 + ,0.0502 + ,50.4 + ,2 + ,158.95 + ,33096200 + ,0.7439 + ,24.36 + ,2325.55 + ,0.056 + ,51.9 + ,2.01 + ,169.53 + ,23273800 + ,0.7439 + ,25.94 + ,2367.52 + ,0.0537 + ,58.5 + ,2 + ,113.66 + ,43697600 + ,0.7439 + ,25.37 + ,2091.88 + ,0.0494 + ,61.4 + ,1.81 + ,107.59 + ,66902300 + ,0.139 + ,21.23 + ,1720.95 + ,0.0366 + ,38.8 + ,0.97 + ,92.67 + ,44957200 + ,0.139 + ,19.35 + ,1535.57 + ,0.0107 + ,44.9 + ,0.39 + ,85.35 + ,33800900 + ,0.139 + ,18.61 + ,1577.03 + ,0.0009 + ,38.6 + ,0.16 + ,90.13 + ,33487900 + ,0.1383 + ,16.37 + ,1476.42 + ,0.0003 + ,4.0 + ,0.15 + ,89.31 + ,27394900 + ,0.1383 + ,15.56 + ,1377.84 + ,0.0024 + ,25.3 + ,0.22 + ,105.12 + ,25963400 + ,0.1383 + ,17.70 + ,1528.59 + ,-0.0038 + ,26.9 + ,0.18 + ,125.83 + ,20952600 + ,0.2874 + ,19.52 + ,1717.30 + ,-0.0074 + ,40.8 + ,0.15 + ,135.81 + ,17702900 + ,0.2874 + ,20.26 + ,1774.33 + ,-0.0128 + ,54.8 + ,0.18 + ,142.43 + ,21282100 + ,0.2874 + ,23.05 + ,1835.04 + ,-0.0143 + ,49.3 + ,0.21 + ,163.39 + ,18449100 + ,0.0596 + ,22.81 + ,1978.50 + ,-0.021 + ,47.4 + ,0.16 + ,168.21 + ,14415700 + ,0.0596 + ,24.04 + ,2009.06 + ,-0.0148 + ,54.5 + ,0.16 + ,185.35 + ,17906300 + ,0.0596 + ,25.08 + ,2122.42 + ,-0.0129 + ,53.4 + ,0.15 + ,188.50 + ,22197500 + ,0.3201 + ,27.04 + ,2045.11 + ,-0.0018 + ,48.7 + ,0.12 + ,199.91 + ,15856500 + ,0.3201 + ,28.81 + ,2144.60 + ,0.0184 + ,50.6 + ,0.12 + ,210.73 + ,19068700 + ,0.3201 + ,29.86 + ,2269.15 + ,0.0272 + ,53.6 + ,0.12 + ,192.06 + ,30855100 + ,0.486 + ,27.61 + ,2147.35 + ,0.0263 + ,56.5 + ,0.11 + ,204.62 + ,21209000 + ,0.486 + ,28.22 + ,2238.26 + ,0.0214 + ,46.4 + ,0.13 + ,235.00 + ,19541600 + ,0.486 + ,28.83 + ,2397.96 + ,0.0231 + ,52.3 + ,0.16 + ,261.09 + ,21955000 + ,0.6129 + ,30.06 + ,2461.19 + ,0.0224 + ,57.7 + ,0.2 + ,256.88 + ,33725900 + ,0.6129 + ,25.51 + ,2257.04 + ,0.0202 + ,62.7 + ,0.2 + ,251.53 + ,28192800 + ,0.6129 + ,22.75 + ,2109.24 + ,0.0105 + ,54.3 + ,0.18 + ,257.25 + ,27377000 + ,0.6665 + ,25.52 + ,2254.70 + ,0.0124 + ,51.0 + ,0.18 + ,243.10 + ,16228100 + ,0.6665 + ,23.33 + ,2114.03 + ,0.0115 + ,53.2 + ,0.19 + ,283.75 + ,21278900 + ,0.6665 + ,24.34 + ,2368.62 + ,0.0114 + ,48.6 + ,0.19) + ,dim=c(8 + ,117) + ,dimnames=list(c('APPLE' + ,'VOLUME' + ,'REV.GROWTH' + ,'MICROSOFT' + ,'NASDAQ' + ,'INFLATION' + ,'CONS.CONF' + ,'FED.FUNDS.RATE') + ,1:117)) > y <- array(NA,dim=c(8,117),dimnames=list(c('APPLE','VOLUME','REV.GROWTH','MICROSOFT','NASDAQ','INFLATION','CONS.CONF','FED.FUNDS.RATE'),1:117)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > par20 = '' > par19 = '' > par18 = '' > par17 = '' > par16 = '' > par15 = '' > par14 = '' > par13 = '' > par12 = '' > par11 = '' > par10 = '' > par9 = '' > par8 = '' > par7 = '' > par6 = '' > par5 = '' > par4 = '' > par3 = '' > par2 = '' > par1 = 'pearson' > 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/fisher/rcomp/tmp/16ras1355268445.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/fisher/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/fisher/rcomp/createtable") > > n <- length(y[,1]) > n [1] 8 > 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/fisher/rcomp/tmp/2i6az1355268446.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 = 5.2726, df = 115, p-value = 6.381e-07 alternative hypothesis: true correlation is not equal to 0 95 percent confidence interval: 0.2823064 0.5765780 sample estimates: cor 0.4412274 > 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 25 warnings (use warnings() to see them) > a<-table.end(a) > table.save(a,file="/var/fisher/rcomp/tmp/31pq41355268446.tab") > > try(system("convert tmp/16ras1355268445.ps tmp/16ras1355268445.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.838 0.259 2.107