Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
467 98.1 98.6 460 101.1 98 448 111.1 106.8 443 93.3 96.7 436 100 100.2 431 108 107.7 484 70.4 92 510 75.4 98.4 513 105.5 107.4 503 112.3 117.7 471 102.5 105.7 471 93.5 97.5 476 86.7 99.9 475 95.2 98.2 470 103.8 104.5 461 97 100.8 455 95.5 101.5 456 101 103.9 517 67.5 99.6 525 64 98.4 523 106.7 112.7 519 100.6 118.4 509 101.2 108.1 512 93.1 105.4 519 84.2 114.6 517 85.8 106.9 510 91.8 115.9 509 92.4 109.8 501 80.3 101.8 507 79.7 114.2 569 62.5 110.8 580 57.1 108.4 578 100.8 127.5 565 100.7 128.6 547 86.2 116.6 555 83.2 127.4
Names of X columns:
Werkl Duurzaam N-Duurzaam
Type of Correlation
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='kendall') 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') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Kendall tau rank correlations for all pairs of data series',3,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'pair',1,TRUE) a<-table.element(a,'tau',1,TRUE) a<-table.element(a,'p-value',1,TRUE) a<-table.row.end(a) n <- length(y[,1]) n cor.test(y[1,],y[2,],method='kendall') for (i in 1:(n-1)) { for (j in (i+1):n) { a<-table.row.start(a) dum <- paste('tau(',dimnames(t(x))[[2]][i]) dum <- paste(dum,',') dum <- paste(dum,dimnames(t(x))[[2]][j]) dum <- paste(dum,')') a<-table.element(a,dum,header=TRUE) r <- cor.test(y[i,],y[j,],method='kendall') a<-table.element(a,r$estimate) a<-table.element(a,r$p.value) a<-table.row.end(a) } } a<-table.end(a) table.save(a,file='mytable.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