Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
101 106 101 104 82 88 99 114 108 105 118 116 107 105 104 111 105 110 117 103 105 119 107 107 127 123 124 128 112 109 135 104 102 132 122 125 136 108 102 143 94 101 142 120 116 153 118 114 145 117 115 138 113 119 148 106 108 152 108 110 169 122 120 169 115 113 161 110 111 174 120 121 179 104 99 191 96 104 190 121 117 182 111 108 175 120 122 181 114 122 197 107 111 194 108 111 197 127 131 216 105 108 221 119 118 218 121 119 230 106 104 227 97 105 204 119 118 197 122 124 199 121 123 208 106 114 191 114 119 202 112 116 211 127 129 224 109 112 224 118 123 231 123 124 244 115 117 235 105 110 250 116 118 266 131 135 288 121 127 283 104 117 295 127 137 312 126 130 334 124 132 348 132 142 383 117 122 407 123 126
Names of X columns:
Grond Uitvoer Invoer
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