Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
95.1 93.3 121.8 49 97 97.3 127.6 40.5 112.7 127 129.9 73.7 102.9 111.7 128 52.6 97.4 96.4 123.5 42 111.4 133 124 51 87.4 72.2 127.4 43.8 96.8 95.8 127.6 44.1 114.1 124.1 128.4 57 110.3 127.6 131.4 56.3 103.9 110.7 135.1 69.7 101.6 104.6 134 72.3 94.6 112.7 144.5 55.1 95.9 115.3 147.3 31.8 104.7 139.4 150.9 54 102.8 119 148.7 66.4 98.1 97.4 141.4 40.8 113.9 154 138.9 58 80.9 81.5 139.8 59.7 95.7 88.8 145.6 49.7 113.2 127.7 147.9 94.3 105.9 105.1 148.5 59.4 108.8 114.9 151.1 94.6 102.3 106.4 157.5 94.8 99 104.5 167.5 71.7 100.7 121.6 172.3 71.4 115.5 141.4 173.5 68.6 100.7 99 187.5 52.9 109.9 126.7 205.5 85.9 114.6 134.1 195.1 73.9 85.4 81.3 204.5 52.4 100.5 88.6 204.5 51.6 114.8 132.7 201.7 70.4 116.5 132.9 207 64.1 112.9 134.4 206.6 68.4 102 103.7 210.6 74.9 106 119.7 211.1 70.7 105.3 115 215 43.2 118.8 132.9 223.9 62.1 106.1 108.5 238.2 34.1 109.3 113.9 238.9 33.9 117.2 142 229.6 43 92.5 97.7 232.2 45.6 104.2 92.2 222.1 39.7 112.5 128.8 221.6 36.2 122.4 134.9 227.3 81.1 113.3 128.2 221 55 100 114.8 213.6 102 110.7 117.9 243.4 31.9 112.8 119.1 253.8 37.5 109.8 120.7 265.3 43.6 117.3 129.1 268.2 34 109.1 117.6 268.5 32 115.9 129.2 266.9 52.6 96 100 268.4 36.8 99.8 87 250.8 59.3 116.8 128 231.2 38.8 115.7 127.7 192 52 99.4 93.4 171.4 32.2 94.3 84.1 160 132.2 91 71.7 148.1 25.4
Names of X columns:
IPtot IPtran IGpci IVtran
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