Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
519 50556 1954 104 517 43901 1828 110 510 48572 2251 104 509 43899 2277 94 501 37532 2085 95 507 40357 2282 99 569 35489 2266 101 580 29027 1878 104 578 34485 2267 99 565 42598 2069 105 547 30306 1746 107 555 26451 2299 110 562 47460 2360 117 561 50104 2214 118 555 61465 2825 126 544 53726 2355 127 537 39477 2333 134 543 43895 3016 131 594 31481 2155 135 611 29896 2172 142 613 33842 2150 141 611 39120 2533 153 594 33702 2058 145 595 25094 2160 137 591 51442 2260 148 589 45594 2498 152 584 52518 2695 169 573 48564 2799 168 567 41745 2946 161 569 49585 2930 174 621 32747 2318 179 629 33379 2540 190 628 35645 2570 190 612 37034 2669 181 595 35681 2450 174 597 20972 2842 180 593 58552 3440 196 590 54955 2678 193 580 65540 2981 197 574 51570 2260 216 573 51145 2844 221 573 46641 2546 217 620 35704 2456 229 626 33253 2295 227 620 35193 2379 204 588 41668 2479 196 566 34865 2057 198 557 21210 2280 207 561 56126 2351 190 549 49231 2276 201 532 59723 2548 210 526 48103 2311 223 511 47472 2201 223 499 50497 2725 231 555 40059 2408 244 565 34149 2139 234 542 36860 1898 250 527 46356 2537 265 510 36577 2068 287 514 23872 2063 283 517 57276 2520 295 508 56389 2434 312 493 57657 2190 333 490 62300 2794 347 469 48929 2070 383 478 51168 2615 407 528 39636 2265 413 534 33213 2139 362 518 38127 2428 321 506 43291 2137 239 502 30600 1823 191 516 21956 2063 159
Names of X columns:
werkl inp bouwv prig
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