Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
252 467 98.8 247 460 100.5 218 448 110.4 226 443 96.4 212 436 101.9 219 431 106.2 227 484 81 221 510 94.7 222 513 101 221 503 109.4 202 471 102.3 211 471 90.7 217 476 96.2 202 475 96.1 199 470 106 195 461 103.1 271 455 102 249 456 104.7 259 517 86 254 525 92.1 222 523 106.9 216 519 112.6 211 509 101.7 203 512 92 200 519 97.4 192 517 97 192 510 105.4 188 509 102.7 179 501 98.1 192 507 104.5 205 569 87.4 203 580 89.9 189 578 109.8 168 565 111.7 169 547 98.6 170 555 96.9 177 562 95.1 177 561 97 166 555 112.7 153 544 102.9 159 537 97.4 163 543 111.4 167 594 87.4 195 611 96.8 179 613 114.1 245 611 110.3 243 594 103.9 240 595 101.6 240 591 94.6 239 589 95.9 236 584 104.7 235 573 102.8 226 567 98.1 223 569 113.9 225 621 80.9 218 629 95.7 215 628 113.2 206 612 105.9 201 595 108.8 200 597 102.3 196 593 99 193 590 100.7 193 580 115.5 184 574 100.7 188 573 109.9 187 573 114.6 185 620 85.4 185 626 100.5 186 620 114.8 179 588 116.5 175 566 112.9 173 557 102
Names of X columns:
Vrouw 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