Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
14.5 14.8 392 127 14.3 14.7 394 123 15.3 16 392 118 14.4 15.4 396 114 13.7 15 392 108 14.2 15.5 396 111 13.5 15.1 419 151 11.9 11.7 421 159 14.6 16.3 420 158 15.6 16.7 418 148 14.1 15 410 138 14.9 14.9 418 137 14.2 14.6 426 136 14.6 15.3 428 133 17.2 17.9 430 126 15.4 16.4 424 120 14.3 15.4 423 114 17.5 17.9 427 116 14.5 15.9 441 153 14.4 13.9 449 162 16.6 17.8 452 161 16.7 17.9 462 149 16.6 17.4 455 139 16.9 16.7 461 135 15.7 16 461 130 16.4 16.6 463 127 18.4 19.1 462 122 16.9 17.8 456 117 16.5 17.2 455 112 18.3 18.6 456 113 15.1 16.3 472 149 15.7 15.1 472 157 18.1 19.2 471 157 16.8 17.7 465 147 18.9 19.1 459 137 19 18 465 132 18.1 17.5 468 125 17.8 17.8 467 123 21.5 21.1 463 117 17.1 17.2 460 114 18.7 19.4 462 111 19 19.8 461 112 16.4 17.6 476 144 16.9 16.2 476 150 18.6 19.5 471 149 19.3 19.9 453 134 19.4 20 443 123 17.6 17.3 442 116 18.6 18.9 444 117 18.1 18.6 438 111 20.4 21.4 427 105 18.1 18.6 424 102 19.6 19.8 416 95 19.9 20.8 406 93 19.2 19.6 431 124 17.8 17.7 434 130 19.2 19.8 418 124 22 22.2 412 115 21.1 20.7 404 106 19.5 17.9 409 105 22.2 20.9 412 105 20.9 21.2 406 101 22.2 21.4 398 95 23.5 23 397 93 21.5 21.3 385 84 24.3 23.9 390 87 22.8 22.4 413 116 20.3 18.3 413 120 23.7 22.8 401 117 23.3 22.3 397 109 19.6 17.8 397 105 18 16.4 409 107 17.3 16 419 109 16.8 16.4 424 109 18.2 17.7 428 108 16.5 16.6 430 107 16 16.2 424 99 18.4 18.3 433 103
Names of X columns:
invbe uitvbe werkl>25jr werkl<25jr
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