Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
1 6.60 6.3 2.0 8.3 4.5 42.0 3 1 3 2.547 4603.00 2.1 1.8 3.9 69.0 624.0 3 5 4 10.550 179.50 9.1 .7 9.8 27.0 180.0 4 4 4 0.023 0.30 15.8 3.9 19.7 19.0 35.0 1 1 1 160 169.00 5.2 1.0 6.2 30.4 392.0 4 5 4 3 25.60 10.9 3.6 14.5 28.0 63.0 1 2 1 52.160 440.00 8.3 1.4 9.7 50.0 230.0 1 1 1 0.425 6.40 11.0 1.5 12.5 7.0 112.0 5 4 4 465 423.00 3.2 .7 3.9 30.0 281.0 5 5 5 0.075 1.20 6.3 2.1 8.4 3.5 42.0 1 1 1 3 25.00 8.6 .0 8.6 50.0 28.0 2 2 2 0.785 3.50 6.6 4.1 10.7 6.0 42.0 2 2 2 0.2 5.00 9.5 1.2 10.7 10.4 120.0 2 2 2 27.660 115.00 3.3 .5 3.8 20.0 148.0 5 5 5 0.12 1.00 11.0 3.4 14.4 3.9 16.0 3 1 2 85 325.00 4.7 1.5 6.2 41.0 310.0 1 3 1 0.101 4.00 10.4 3.4 13.8 9.0 28.0 5 1 3 1.040 5.50 7.4 .8 8.2 7.6 68.0 5 3 4 521 655.00 2.1 .8 2.9 46.0 336.0 5 5 5 0.005 0.14 7.7 1.4 9.1 2.6 21.5 5 2 4 0.1 0.25 17.9 2.0 19.9 24.0 50.0 1 1 1 62 1320.00 6.1 1.9 8.0 100.0 267.0 1 1 1 0.023 0.40 11.9 1.3 13.2 3.2 19.0 4 1 3 0.048 0.33 10.8 2.0 12.8 2.0 30.0 4 1 3 1.700 6.30 13.8 5.6 19.4 5.0 12.0 2 1 1 3.500 10.80 14.3 3.1 17.4 6.5 120.0 2 1 1 0.48 15.50 15.2 1.8 17.0 12.0 140.0 2 2 2 10 115.00 10.0 .9 10.9 20.2 170.0 4 4 4 1.620 11.40 11.9 1.8 13.7 13.0 17.0 2 1 2 192 180.00 6.5 1.9 8.4 27.0 115.0 4 4 4 2.500 12.10 7.5 .9 8.4 18.0 31.0 5 5 5 0.28 1.90 10.6 2.6 13.2 4.7 21.0 3 1 3 4.235 50.40 7.4 2.4 9.8 9.8 52.0 1 1 1 6.800 179.00 8.4 1.2 9.6 29.0 164.0 2 3 2 0.75 12.30 5.7 .9 6.6 7.0 225.0 2 2 2 3.600 21.00 4.9 .5 5.4 6.0 225.0 3 2 3 55.500 175.00 3.2 .6 3.8 20.0 151.0 5 5 5 0.9 2.60 11.0 2.3 13.3 4.5 60.0 2 1 2 2 12.30 4.9 .5 5.4 7.5 200.0 3 1 3 0.104 2.50 13.2 2.6 15.8 2.3 46.0 3 2 2 4.190 58.00 9.7 .6 10.3 24.0 210.0 4 3 4 3.500 3.90 12.8 6.6 19.4 3.0 14.0 2 1 1
Names of X columns:
Wb Wbr SWS PS totS L Gt P S D
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