Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
130 461 14.1 2 127 463 14.8 2.3 122 462 16.8 2.8 117 456 15.4 2.4 112 455 15.2 2.3 113 456 16.9 2.7 149 472 14.1 2.7 157 472 14.7 2.9 157 471 16.5 3 147 465 15.2 2.2 137 459 17.6 2.3 132 465 18 2.8 125 468 16.9 2.8 123 467 16.7 2.8 117 463 19.7 2.2 114 460 15.9 2.6 111 462 17.4 2.8 112 461 17.7 2.5 144 476 15.2 2.4 150 476 15.7 2.3 149 471 17.2 1.9 134 453 17.7 1.7 123 443 17.9 2 116 442 16.2 2.1 117 444 17.5 1.7 111 438 16.8 1.8 105 427 19.1 1.8 102 424 16.7 1.8 95 416 18.2 1.3 93 406 18.5 1.3 124 431 17.8 1.3 130 434 16.4 1.2 124 418 18 1.4 115 412 20.3 2.2 106 404 19.5 2.9 105 409 18 3.1 105 412 20.2 3.5 101 406 19 3.6 95 398 20.2 4.4 93 397 21.5 4.1 84 385 19.7 5.1 87 390 21.1 5.8 116 413 20.2 5.9 120 413 18.2 5.4 117 401 21.3 5.5 109 397 20.4 4.8 105 397 17.2 3.2 107 409 15.8 2.7 109 419 15.1 2.1 109 424 14.5 1.9 108 428 15.8 0.6 107 430 14.3 0.7 99 424 13.9 -0.2 103 433 15.5 -1 131 456 14.3 -1.7 137 459 13.6 -0.7 135 446 16.3 -1 124 441 16.8 -0.9 118 439 16 0 121 454 16.8 0.3 121 460 16 0.8
Names of X columns:
Werkloosheid-25 Werkloosheid+25 Invoer CPindex
Type of Correlation
kendall
pearson
spearman
kendall
Chart options
Title:
R Code
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', ...) } bitmap(file='test1.png') pairs(t(y),diag.panel=panel.hist, upper.panel=panel.smooth, lower.panel=panel.tau, main=main) dev.off() load(file='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='mytable.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='mytable1.tab')
Compute
Summary of computational transaction
Raw Input
view raw input (R code)
Raw Output
view raw output of R engine
Computing time
0 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation