Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
97.4 519 291 97 517 289 105.4 510 292 102.7 509 283 98.1 501 283 104.5 507 274 87.4 569 281 89.9 580 271 109.8 578 277 111.7 565 271 98.6 547 272 96.9 555 282 95.1 562 281 97 561 284 112.7 555 282 102.9 544 279 97.4 537 278 111.4 543 280 87.4 594 287 96.8 611 287 114.1 613 283 110.3 611 314 103.9 594 315 101.6 595 306 94.6 591 309 95.9 589 307 104.7 584 302 102.8 573 297 98.1 567 303 113.9 569 304 80.9 621 307 95.7 629 310 113.2 628 300 105.9 612 295 108.8 595 300 102.3 597 300 99 593 297 100.7 590 291 115.5 580 288 100.7 574 290 109.9 573 293 114.6 573 282 85.4 620 287 100.5 626 275 114.8 620 282 116.5 588 282 112.9 566 277 102 557 275 106 561 284 105.3 549 280 118.8 532 270 106.1 526 265 109.3 511 255 117.2 499 246 92.5 555 254 104.2 565 258 112.5 542 263 122.4 527 243 113.3 510 230 100 514 236 110.7 517 249 112.8 508 236 109.8 493 231 117.3 490 229 109.1 469 227 115.9 478 227 96 528 233 99.8 534 232 116.8 518 227 115.7 506 217 99.4 502 217 94.3 516 220
Names of X columns:
Ind Werkl Vrouw
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