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