Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
528.00 4.256 65.1 478.00 5.008 72.5 469.00 4.181 70.6 490.00 4.666 81.9 493.00 4.476 72 508.00 4.912 70.9 517.00 4.340 71.1 514.00 3.830 65.1 510.00 4.285 78.1 527.00 5.537 92 542.00 3.813 77.8 565.00 4.610 63.3 555.00 4.207 56.2 499.00 5.099 79.2 511.00 4.112 69 526.00 4.199 66.1 532.00 5.110 77.5 549.00 4.218 69.3 561.00 4.736 70.2 557.00 4.624 70.2 566.00 4.145 78.2 588.00 5.299 85.4 620.00 5.011 82.4 626.00 4.730 61.2 620.00 4.619 52.2 573.00 5.578 85.3 573.00 5.369 79.9 574.00 4.902 72.2 580.00 6.102 85.7 590.00 5.024 75.5 593.00 5.731 69.2 597.00 5.732 77.6 595.00 4.491 85.3 612.00 4.755 77 628.00 5.208 89.9 629.00 4.962 60 621.00 4.163 54.3 569.00 5.592 84 567.00 5.754 69.9 573.00 4.929 75.1 584.00 5.219 81.7 589.00 4.429 69.9 591.00 4.143 68.3 595.00 4.308 77.3 594.00 3.996 77.4 611.00 4.634 85.3 613.00 4.138 91 611.00 3.759 60.6 594.00 3.922 57.6 543.00 5.560 93.8 537.00 4.004 78.7 544.00 3.937 80.3 555.00 5.250 89.8 561.00 3.908 77.5 562.00 4.814 71.7 555.00 4.407 83.2 547.00 3.243 86.2 565.00 3.740 100.7 578.00 3.949 100.8 580.00 3.711 57.1 569.00 3.796 62.5 507.00 4.145 79.7 501.00 3.499 80.3 509.00 4.164 92.4 510.00 3.902 91.8 517.00 3.186 85.8 519.00 3.353 84.2 512.00 3.475 93.1 509.00 3.032 101.2
Names of X columns:
W VPW PDC
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