Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
10519.20 1154.80 10414.90 1206.70 12476.80 1199.00 12384.60 1265.00 12266.70 1247.10 12919.90 1116.50 11497.30 1153.90 12142.00 1077.40 13919.40 1132.50 12656.80 1058.80 12034.10 1195.10 13199.70 1263.40 10881.30 1023.10 11301.20 1141.00 13643.90 1116.30 12517.00 1135.60 13981.10 1210.50 14275.70 1230.00 13425.00 1136.50 13565.70 1068.70 16216.30 1372.50 12970.00 1049.90 14079.90 1302.20 14235.00 1305.90 12213.40 1173.50 12581.00 1277.40 14130.40 1238.60 14210.80 1508.60 14378.50 1423.40 13142.80 1375.10 13714.70 1344.10 13621.90 1287.50 15379.80 1446.90 13306.30 1451.00 14391.20 1604.40 14909.90 1501.50 14025.40 1522.80 12951.20 1328.00 14344.30 1420.50 16093.40 1648.00 15413.60 1631.10 14705.70 1396.60 15972.80 1663.40 16241.40 1283.00 16626.40 1582.40 17136.20 1785.20 15622.90 1853.60 18003.90 1994.10 16136.10 2042.80 14423.70 1586.10 16789.40 1942.40 16782.20 1763.60 14133.80 1819.90 12607.00 1836.00 12004.50 1447.50 12175.40 1509.50 13268.00 1661.20 12299.30 1456.20 11800.60 1310.90 13873.30 1542.10 12315.00 1537.70
Names of X columns:
invoerEU invoerAM
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