Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
613 0 4.9 0 611 25.7 4.6 25.7 594 24.7 4.3 24.7 595 24.2 4.3 24.2 591 23.6 4.6 23.6 589 24.4 5.1 24.4 584 22.5 4.8 22.5 573 19.4 4.5 19.4 567 18.1 4.9 18.1 569 18.1 5.1 18.1 621 20.7 5.1 20.7 629 19.1 5.2 19.1 628 18.3 4.5 18.3 612 16.9 4.6 16.9 595 17.9 4.9 17.9 597 20.2 4.6 20.2 593 21.2 4.4 21.2 590 23.8 3.7 23.8 580 24 4 24 574 26.6 4.2 26.6 573 25.3 3.9 25.3 573 27.6 3.6 27.6 620 24.7 3.6 24.7 626 26.6 3.2 26.6 620 24.4 3.2 24.4 588 24.6 3.5 24.6 566 26 3.6 26 557 24.8 3.7 24.8 561 24 3.8 24 549 22.7 3.8 22.7 532 23 3.8 23 526 24.1 3.3 24.1 511 24 3.3 24 499 22.7 3.4 22.7 555 22.6 3.1 22.6 565 23.1 3.5 23.1 542 24.4 4.2 24.4 527 23 4.9 23 510 22 5.1 22 514 21.3 5.5 21.3 517 21.5 5.6 21.5 508 21.3 6.4 21.3 493 23.2 6.2 23.2 490 21.8 7.2 21.8 469 23.3 7.8 23.3 478 21 7.9 21 528 22.4 7.4 22.4 534 20.4 7.5 20.4 518 19.9 6.7 19.9 506 21.3 5.1 21.3 502 18.9 4.6 18.9 516 15.6 4.3 15.6 528 12.5 3.9 12.5 533 7.8 2.6 7.8 536 5.5 2.6 5.5 537 4 1.6 4 524 3.3 0.9 3.3 536 3.7 0.3 3.7 587 3.1 1.2 3.1 597 5 0.8 5 581 6.3 1 6.3
Names of X columns:
WLH ICONS HICP INDPRO
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