Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
17823.2 462 17872 455 17420.4 461 16704.4 461 15991.2 463 16583.6 462 19123.5 456 17838.7 455 17209.4 456 18586.5 472 16258.1 472 15141.6 471 19202.1 465 17746.5 459 19090.1 465 18040.3 468 17515.5 467 17751.8 463 21072.4 460 17170 462 19439.5 461 19795.4 476 17574.9 476 16165.4 471 19464.6 453 19932.1 443 19961.2 442 17343.4 444 18924.2 438 18574.1 427 21350.6 424 18594.6 416 19823.1 406 20844.4 431 19640.2 434 17735.4 418 19813.6 412 22160 404 20664.3 409 17877.4 412 20906.5 406 21164.1 398 21374.4 397 22952.3 385 21343.5 390 23899.3 413 22392.9 413 18274.1 401 22786.7 397 22321.5 397 17842.2 409 16373.5 419 15993.8 424 16446.1 428 17729 430 16643 424 16196.7 433 18252.1 456 17570.4 459 15836.8 446
Names of X columns:
Uitvoer Werkloosh.
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