Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
14.9 5.9 2253 418 18.6 6.2 2218 438 19.1 6.2 1855 389 18.8 6.4 2187 368 18.2 6.8 1852 386 18 6.9 1570 261 19 7 1851 294 20.7 7 1954 263 21.2 6.9 1828 293 20.7 6.7 2251 303 19.6 6.7 2277 326 18.6 6.5 2085 314 18.7 6.4 2282 332 23.8 6.5 2266 347 24.9 6.5 1878 290 24.8 6.5 2267 340 23.8 6.7 2069 371 22.3 6.8 1746 340 21.7 7.2 2299 376 20.7 7.6 2360 322 19.7 7.6 2214 364 18.4 7.2 2825 379 17.4 6.4 2355 343 17 6.1 2333 358 18 6.3 3016 433 23.8 7.1 2155 344 25.5 7.5 2172 357 25.6 7.4 2150 385 23.7 7.1 2533 392 22 6.8 2058 308 21.3 6.9 2160 294 20.7 7.2 2260 300 20.4 7.4 2498 302 20.3 7.3 2695 400 20.4 6.9 2799 392 19.8 6.9 2946 373 19.5 6.8 2930 379 23.1 7.1 2318 303 23.5 7.2 2540 324 23.5 7.1 2570 353 22.9 7 2669 392 21.9 6.9 2450 327 21.5 7.1 2842 376 20.5 7.3 3440 329 20.2 7.5 2678 359 19.4 7.5 2981 413 19.2 7.5 2260 338 18.8 7.3 2844 422 18.8 7 2546 390 22.6 6.7 2456 370 23.3 6.5 2295 367 23 6.5 2379 406 21.4 6.5 2479 418 19.9 6.6 2057 346 18.8 6.8 2280 350 18.6 6.9 2351 329 18.4 6.9 2276 318 18.6 6.8 2548 381 19.9 6.8 2311 337 19.2 6.5 2201 372 18.4 6.1 2725 422
Names of X columns:
<25 >25 wngb ntwngb
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