Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
152 467 98.8 148 460 100.5 148 448 110.4 157 443 96.4 157 436 101.9 162 431 106.2 174 484 81 158 510 94.7 166 513 101 161 503 109.4 147 471 102.3 160 471 90.7 157 476 96.2 147 475 96.1 144 470 106 151 461 103.1 154 455 102 143 456 104.7 144 517 86 137 525 92.1 127 523 106.9 122 519 112.6 124 509 101.7 121 512 92 122 519 97.4 112 517 97 114 510 105.4 109 509 102.7 117 501 98.1 123 507 104.5 126 569 87.4 127 580 89.9 117 578 109.8 109 565 111.7 105 547 98.6 107 555 96.9 106 562 95.1 111 561 97 106 555 112.7 79 544 102.9 112 537 97.4 103 543 111.4 106 594 87.4 123 611 96.8 118 613 114.1 125 611 110.3 122 594 103.9 126 595 101.6 126 591 94.6 132 589 95.9 134 584 104.7 135 573 102.8 135 567 98.1 133 569 113.9 133 621 80.9 127 629 95.7 129 628 113.2 123 612 105.9 121 595 108.8 122 597 102.3 122 593 99 120 590 100.7 118 580 115.5 117 574 100.7 121 573 109.9 116 573 114.6 118 620 85.4 115 626 100.5 115 620 114.8 114 588 116.5 114 566 112.9 111 557 102
Names of X columns:
Man Werklh Prod
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