x <- array(list(0
,210907
,0
,2
,0
,179321
,NA
,4
,0
,149061
,0
,0
,0
,237213
,1
,0
,0
,173326
,NA
,-4
,0
,133131
,1
,4
,0
,258873
,NA
,4
,0
,324799
,1
,0
,0
,230964
,0
,-1
,0
,236785
,1
,0
,0
,344297
,1
,1
,0
,174724
,1
,0
,0
,174415
,1
,3
,0
,223632
,1
,-1
,0
,294424
,0
,4
,0
,325107
,1
,3
,0
,106408
,0
,1
,0
,96560
,0
,0
,0
,265769
,1
,-2
,0
,269651
,NA
,-3
,0
,149112
,0
,-4
,0
,152871
,0
,2
,0
,362301
,1
,2
,0
,183167
,0
,-4
,0
,277965
,NA
,3
,0
,218946
,1
,2
,0
,244052
,1
,2
,0
,341570
,1
,0
,0
,233328
,NA
,5
,0
,206161
,NA
,-2
,0
,311473
,NA
,0
,0
,207176
,NA
,-2
,0
,196553
,1
,-3
,0
,143246
,0
,2
,0
,182192
,NA
,2
,0
,194979
,NA
,2
,0
,167488
,NA
,0
,0
,143756
,0
,4
,0
,275541
,NA
,4
,0
,152299
,1
,2
,0
,193339
,1
,2
,0
,130585
,0
,-4
,0
,112611
,1
,3
,0
,148446
,1
,3
,0
,182079
,0
,2
,0
,243060
,1
,-1
,0
,162765
,1
,-3
,0
,85574
,1
,0
,0
,225060
,0
,1
,0
,133328
,1
,-3
,0
,100750
,1
,3
,0
,101523
,1
,0
,0
,243511
,1
,0
,0
,152474
,1
,0
,0
,132487
,1
,3
,0
,317394
,0
,-3
,0
,244749
,1
,0
,0
,184510
,NA
,-4
,0
,128423
,0
,2
,0
,97839
,0
,-1
,1
,172494
,NA
,3
,1
,229242
,1
,2
,1
,351619
,NA
,5
,1
,324598
,0
,2
,1
,195838
,0
,-2
,1
,254488
,0
,0
,1
,199476
,NA
,3
,1
,92499
,1
,-2
,1
,224330
,0
,0
,1
,181633
,1
,6
,1
,271856
,1
,-3
,1
,95227
,1
,3
,1
,98146
,0
,0
,1
,118612
,0
,-2
,1
,65475
,1
,1
,1
,108446
,0
,0
,1
,121848
,0
,2
,1
,76302
,1
,2
,1
,98104
,0
,-3
,1
,30989
,1
,-2
,1
,31774
,0
,1
,1
,150580
,1
,-4
,1
,54157
,NA
,0
,1
,59382
,0
,1
,1
,84105
,0
,0)
,dim=c(4
,85)
,dimnames=list(c('pop'
,'time_in_rfc'
,'gender'
,'total_tests')
,1:85))
 y <- array(NA,dim=c(4,85),dimnames=list(c('pop','time_in_rfc','gender','total_tests'),1:85))
 for (i in 1:dim(x)[1])
 {
 	for (j in 1:dim(x)[2])
 	{
 		y[i,j] <- as.numeric(x[i,j])
 	}
 }
par1 = 'pearson'
main = 'Correlation Matrix'
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', ...)
}
postscript(file="/var/wessaorg/rcomp/tmp/1irho1323610308.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
pairs(t(y),diag.panel=panel.hist, upper.panel=panel.smooth, lower.panel=panel.tau, main=main)
dev.off()

#Note: the /var/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
load(file="/var/wessaorg/rcomp/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="/var/wessaorg/rcomp/tmp/2jq841323610308.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="/var/wessaorg/rcomp/tmp/32jpt1323610308.tab") 

try(system("convert tmp/1irho1323610308.ps tmp/1irho1323610308.png",intern=TRUE))