Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
255198 276475 35907 104609 253353 272566 33335 101802 246057 264981 23988 94542 235372 263290 23099 93051 258556 296806 46390 124129 260993 303598 51588 130374 254663 286994 51579 123946 250643 276427 45390 114971 243422 266424 39215 105531 247105 267153 38433 104919 248541 268381 37676 104782 245039 262522 36055 101281 237080 255542 32986 94545 237085 253158 30953 93248 225554 243803 23558 84031 226839 250741 22487 87486 247934 280445 43528 115867 248333 285257 47913 120327 246969 270976 48621 117008 245098 261076 42169 108811 246263 255603 38444 104519 255765 260376 38692 106758 264319 263903 38124 109337 268347 264291 37886 109078 273046 263276 37310 108293 273963 262572 34689 106534 267430 256167 26450 99197 271993 264221 25565 103493 292710 293860 46562 130676 295881 300713 52653 137448 293299 287224 54807 134704
Names of X columns:
WLHmannen WLHvrouwen WLHjongeren_in_wachttijd WLHjonger_25J
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
1 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation