Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
773.9 902.2 7.1 0.8833 795.9 891.9 7.7 0.87 836.3 874 8.2 0.8758 876.1 930.9 9.2 0.8858 851.7 944.2 10.2 0.917 692.4 935.9 10.5 0.9554 877.3 937.1 10.8 0.9922 536.8 885.1 10.8 0.9778 705.9 892.4 11.1 0.9808 951 987.3 11 0.9811 755.7 946.3 11.1 1.0014 695.5 799.6 10.2 1.0183 744.8 875.4 9.6 1.0622 672.1 846.2 9.3 1.0773 666.6 880.6 8.9 1.0807 760.8 885.7 8.8 1.0848 756 868.9 8.9 1.1582 604.4 882.5 9 1.1663 883.9 789.6 8.9 1.1372 527.9 773.3 8.8 1.1139 756.2 804.3 9 1.1222 812.9 817.8 8.8 1.1692 655.6 836.7 8.8 1.1702 707.6 721.8 9.7 1.2286 612.6 760.8 9.8 1.2613 659.2 841.4 10 1.2646 833.4 1045.6 10.1 1.2262 727.8 949.2 10.1 1.1985 797.2 850.1 10 1.2007 753 957.4 9 1.2138 762 851.8 8.8 1.2266 613.7 913.9 8.4 1.2176 759.2 888 8 1.2218 816.4 973.8 8.1 1.249 736.8 927.6 8 1.2991 680.1 833 7.8 1.3408 736.5 879.5 7.7 1.3119 637.2 797.3 7.3 1.3014 801.9 834.5 7.2 1.3201 772.3 735.1 7.2 1.2938 897.3 835 6.7 1.2694 792.1 892.8 6.7 1.2165 826.8 697.2 6.4 1.2037 666.8 821.1 6.4 1.2292 906.6 732.7 6.3 1.2256 871.4 797.6 6 1.2015 891 866.3 5.5 1.1786 739.2 826.3 4.7 1.1856 833.6 778.6 3.8 1.2103 715.6 779.2 3.6 1.1938 871.6 951 2.5 1.202 751.6 692.3 1.7 1.2271 1005.50 841.4 1.6 1.277 681.2 857.3 2.2 1.265 837.3 760.7 2.2 1.2684 674.7 841.2 2 1.2811 806.3 810.3 1.9 1.2727 860.2 1007.4 1.6 1.2611 689.8 931.3 1.9 1.2881 691.6 931.2 2.6 1.3213
Names of X columns:
EXP IMP HDB EURODOL
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