Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
97.4 519 1867 97 517 1864 105.4 510 1853 102.7 509 1866 98.1 501 1852 104.5 507 1823 87.4 569 1862 89.9 580 1859 109.8 578 1855 111.7 565 1871 98.6 547 1822 96.9 555 1865 95.1 562 1922 97 561 1960 112.7 555 1978 102.9 544 1975 97.4 537 1953 111.4 543 1931 87.4 594 1922 96.8 611 1904 114.1 613 1927 110.3 611 1943 103.9 594 1937 101.6 595 1960 94.6 591 1955 95.9 589 1958 104.7 584 1953 102.8 573 1983 98.1 567 1985 113.9 569 1980 80.9 621 1975 95.7 629 1961 113.2 628 1996 105.9 612 2008 108.8 595 2016 102.3 597 2025 99 593 2118 100.7 590 2111 115.5 580 2137 100.7 574 2110 109.9 573 2131 114.6 573 2111 85.4 620 2132 100.5 626 2140 114.8 620 2143 116.5 588 2139 112.9 566 2101 102 557 2089 106 561 2127 105.3 549 2084 118.8 532 1980 106.1 526 2001 109.3 511 1944 117.2 499 1783 92.5 555 1832 104.2 565 1827 112.5 542 1806 122.4 527 1818 113.3 510 1798 100 514 1824 110.7 517 1834 112.8 508 1851 109.8 493 1813 117.3 490 1788 109.1 469 1737 115.9 478 1699 96 528 1724 99.8 534 1686 116.8 518 1686 115.7 506 1682 99.4 502 1698 94.3 516 1742
Names of X columns:
Ind Werkl Man
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