Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
153.24 146.54 133.39 98.86 184.48 120.13 133.74 100.83 191.81 131.67 129.67 117.15 168.19 131.97 126.70 106.96 163.81 145.92 126.84 101.25 190.57 177.02 128.45 115.80 163.81 149.56 129.84 90.85 129.62 171.58 128.83 100.62 173.90 173.95 129.22 118.61 198.76 190.39 132.14 114.66 135.52 183.46 137.40 108.00 179.81 165.44 141.78 105.61 137.05 186.32 138.74 98.34 142.57 223.29 137.63 99.69 187.52 198.99 139.61 108.84 220.48 191.05 136.82 106.86 208.76 178.42 134.24 101.98 210.19 187.85 128.64 118.40 232.57 183.51 126.43 84.10 173.81 252.94 127.25 99.48 218.86 213.51 126.72 117.67 226.76 185.53 124.18 110.08 196.67 215.48 121.73 113.10 237.43 214.39 122.34 106.34 173.14 229.21 124.74 102.91 207.62 183.55 122.81 104.68 234.67 206.71 123.40 120.06 204.10 186.23 125.68 104.68 230.76 217.46 130.78 114.24 210.19 214.69 129.41 119.13 194.76 202.06 129.49 88.77 172.10 225.57 130.51 104.47 221.90 220.70 129.01 119.33 225.24 246.32 127.33 121.10 228.00 273.51 129.41 117.36 198.76 220.66 132.06 106.03 199.05 295.88 129.35 110.19 235.43 215.35 129.47 109.46 270.76 230.83 130.92 123.49 234.10 220.00 133.39 110.29 237.24 232.06 132.48 113.62 239.43 237.68 130.86 121.83 239.24 294.39 132.75 96.15 197.33 295.35 131.81 108.32 217.43 267.68 133.56 116.94 242.19 274.12 136.38 127.23 207.52 246.84 139.22 117.78 232.76 249.34 137.23 103.95 222.10 303.25 136.16 115.07 202.48 236.14 135.03 117.26 228.10 233.38 140.33 114.14 319.52 260.96 140.85 121.93 236.95 281.18 138.58 113.41 252.00 281.54 137.08 120.48 262.29 288.95 137.75 99.79 172.10 332.68 131.10 103.74 243.90 345.96 125.54 121.41 235.62 414.96 116.36 120.27 216.95 285.35 111.10 103.33 236.29 288.03 117.81 98.02
Names of X columns:
ENC IVC WK IP
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
1 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation