Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
8 11.1 5560 33842 8.1 10.9 3922 39120 7.7 10 3759 33702 7.5 9.2 4138 25094 7.6 9.2 4634 51442 7.8 9.5 3996 45594 7.8 9.6 4308 52518 7.8 9.5 4143 48564 7.5 9.1 4429 41745 7.5 8.9 5219 49585 7.1 9 4929 32747 7.5 10.1 5755 33379 7.5 10.3 5592 35645 7.6 10.2 4163 37034 7.7 9.6 4962 35681 7.7 9.2 5208 20972 7.9 9.3 4755 58552 8.1 9.4 4491 54955 8.2 9.4 5732 65540 8.2 9.2 5731 51570 8.2 9 5040 51145 7.9 9 6102 46641 7.3 9 4904 35704 6.9 9.8 5369 33253 6.7 10 5578 35193 6.7 9.8 4619 41668 6.9 9.3 4731 34865 7 9 5011 21210 7.1 9 5299 56126 7.2 9.1 4146 49231 7.1 9.1 4625 59723 6.9 9.1 4736 48103 7 9.2 4219 47472 6.8 8.8 5116 50497 6.4 8.3 4205 40059 6.7 8.4 4121 34149 6.6 8.1 5103 36860 6.4 7.7 4300 46356 6.3 7.9 4578 36577 6.2 7.9 3809 23872 6.5 8 5526 57276 6.8 7.9 4247 56389 6.8 7.6 3830 57657 6.4 7.1 4394 62300 6.1 6.8 4826 48929 5.8 6.5 4409 51168 6.1 6.9 4569 39636 7.2 8.2 4106 33213 7.3 8.7 4794 38127 6.9 8.3 3914 43291 6.1 7.9 3793 30600 5.8 7.5 4405 21956 6.2 7.8 4022 48033 7.1 8.3 4100 46148 7.7 8.4 4788 50736 7.9 8.2 3163 48114 7.7 7.7 3585 38390 7.4 7.2 3903 44112 7.5 7.3 4178 36287 8 8.1 3863 30333 8.1 8.5 4187 30301
Names of X columns:
Wm Wv Bouwv Inspw
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