Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
1 521 18308 185 4.041 7.2 1 367 1148 600 0.55 8.5 1 443 18068 372 3.665 5.7 1 365 7729 142 2.351 7.3 1 614 100484 432 29.76 7.5 1 385 16728 290 3.294 5 1 286 14630 346 3.287 6.7 1 397 4008 328 0.666 6.2 1 764 38927 354 12.938 7.3 1 427 22322 266 6.478 5 1 153 3711 320 1.108 2.8 1 231 3136 197 1.007 6.1 1 524 50508 266 11.431 7.1 1 328 28886 173 5.544 5.9 1 240 16996 190 2.777 4.6 1 286 13035 239 2.478 4.4 1 285 12973 190 3.685 7.4 1 569 16309 241 4.22 7.1 1 96 5227 189 1.228 7.5 1 498 19235 358 4.781 5.9 1 481 44487 315 6.016 9 1 468 44213 303 9.295 9.2 1 177 23619 228 4.375 5.1 1 198 9106 134 2.573 8.6 1 458 24917 189 5.117 6.6 1 108 3872 196 0.799 6.9 1 246 8945 183 1.578 2.7 1 291 2373 417 1.202 5.5 1 68 7128 233 1.109 7.2 1 311 23624 349 7.73 6.6 1 606 5242 284 1.515 6.9 1 512 92629 499 17.99 7.2 1 426 28795 231 6.629 5.8 1 47 4487 143 0.639 4.1 1 265 48799 249 10.847 6.4 1 370 14067 195 3.146 6.7 1 312 12693 288 2.842 6 1 222 62184 229 11.882 6.9 1 280 9153 287 1.003 8.5 1 759 14250 224 3.487 6.2 1 114 3680 161 0.696 3.4 1 419 18063 221 4.877 6.6 1 435 65112 237 16.987 6.6 1 186 11340 220 1.723 4.9 1 87 4553 185 0.563 6.4 1 188 28960 260 6.187 5.8 1 303 19201 261 4.867 6.3 1 102 7533 118 1.793 10.5 1 127 26343 268 4.892 5.4 1 251 1641 300 0.454 5.1 0 205 145360 237 10.379 6.8 0 453 9066420 240 82.422 5.6 0 320 1038933 185 16.491 3.8 0 405 2739420 201 60.876 8.2 0 89 61620 193 0.474 4.1 0 74 827530 254 7.523 2.8 0 101 534100 230 5.45 6.3 0 321 328755 197 10.605 11.4 0 315 1413895 248 40.397 19.4 0 229 2909136 258 60.607 5.8 0 302 3604246 206 58.133 6.9 0 216 917504 199 8.192 3.5
Names of X columns:
Pop Assaults BachDgr PoliceExp TotlPopul Unempl
Type of Correlation
kendall
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=par1) 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') n <- length(y[,1]) n a<-table.start() a<-table.row.start(a) a<-table.element(a,paste('Correlations for all pairs of data series (method=',par1,')',sep=''),n+1,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,' ',header=TRUE) for (i in 1:n) { a<-table.element(a,dimnames(t(x))[[2]][i],header=TRUE) } a<-table.row.end(a) for (i in 1:n) { a<-table.row.start(a) a<-table.element(a,dimnames(t(x))[[2]][i],header=TRUE) for (j in 1:n) { r <- cor.test(y[i,],y[j,],method=par1) a<-table.element(a,round(r$estimate,3)) } a<-table.row.end(a) } a<-table.end(a) table.save(a,file='mytable.tab') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Correlations for all pairs of data series with p-values',4,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'pair',1,TRUE) a<-table.element(a,'Pearson r',1,TRUE) a<-table.element(a,'Spearman rho',1,TRUE) a<-table.element(a,'Kendall tau',1,TRUE) a<-table.row.end(a) cor.test(y[1,],y[2,],method=par1) for (i in 1:(n-1)) { for (j in (i+1):n) { a<-table.row.start(a) dum <- paste(dimnames(t(x))[[2]][i],';',dimnames(t(x))[[2]][j],sep='') a<-table.element(a,dum,header=TRUE) rp <- cor.test(y[i,],y[j,],method='pearson') a<-table.element(a,round(rp$estimate,4)) rs <- cor.test(y[i,],y[j,],method='spearman') a<-table.element(a,round(rs$estimate,4)) rk <- cor.test(y[i,],y[j,],method='kendall') a<-table.element(a,round(rk$estimate,4)) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'p-value',header=T) a<-table.element(a,paste('(',round(rp$p.value,4),')',sep='')) a<-table.element(a,paste('(',round(rs$p.value,4),')',sep='')) a<-table.element(a,paste('(',round(rk$p.value,4),')',sep='')) a<-table.row.end(a) } } a<-table.end(a) table.save(a,file='mytable1.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