Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
161 452 25.60 7.40 149 462 23.70 7.10 139 455 22.00 6.80 135 461 21.30 6.90 130 461 20.70 7.20 127 463 20.40 7.40 122 462 20.30 7.30 117 456 20.40 6.90 112 455 19.80 6.90 113 456 19.50 6.80 149 472 23.10 7.10 157 472 23.50 7.20 157 471 23.50 7.10 147 465 22.90 7.00 137 459 21.90 6.90 132 465 21.50 7.10 125 468 20.50 7.30 123 467 20.20 7.50 117 463 19.40 7.50 114 460 19.20 7.50 111 462 18.80 7.30 112 461 18.80 7.00 144 476 22.60 6.70 150 476 23.30 6.50 149 471 23.00 6.50 134 453 21.40 6.50 123 443 19.90 6.60 116 442 18.80 6.80 117 444 18.60 6.90 111 438 18.40 6.90 105 427 18.60 6.80 102 424 19.90 6.80 95 416 19.20 6.50 93 406 18.40 6.10 124 431 21.10 6.10 130 434 20.50 5.90 124 418 19.10 5.70 115 412 18.10 5.90 106 404 17.00 5.90 105 409 17.10 6.10 105 412 17.40 6.30 101 406 16.80 6.20 95 398 15.30 5.90 93 397 14.30 5.70 84 385 13.40 5.40 87 390 15.30 5.60 116 413 22.10 6.20 120 413 23.70 6.30 117 401 22.20 6.00 109 397 19.50 5.60 105 397 16.60 5.50 107 409 17.30 5.90 109 419 19.80 6.50 109 424 21.20 6.80 108 428 21.50 6.80 107 430 20.60 6.50 99 424 19.10 6.20 103 433 19.60 6.20 131 456 23.50 6.50 137 459 24.00 6.70
Names of X columns:
wklhd<25 wklhd>25 wkz<25 wkz>25
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