• QR Code (size=2)
• QR Code (size=3)
• QR Code (size=4)
• QR Code (size=5)
• QR Code (size=6)

12	13	6
7	11	4
13	14	6
11	12	5
16	12	5
10	6	4
15	10	5
5	11	3
4	10	2
7	12	5
15	15	6
5	13	6
16	18	8
15	11	6
13	12	3
13	13	6
15	14	6
15	16	7
10	16	8
17	16	6
14	15	7
9	13	4
6	8	4
11	14	2
13	15	6
12	13	6
10	16	6
4	13	6
13	12	6
15	15	7
8	11	4
10	14	3
8	13	5
7	13	6
9	12	4
14	14	6
5	13	3
7	12	3
16	14	6
14	15	6
16	16	6
15	15	8
4	5	2
12	15	6
8	8	4
17	16	7
15	16	6
16	14	6
12	16	6
12	14	5
13	13	6
14	14	6
14	14	5
15	12	6
14	13	7
11	15	5
13	15	6
4	13	6
8	10	4
13	13	5
15	14	6
15	13	6
8	13	4
17	18	6
12	12	4
13	14	7
14	16	8
7	13	6
16	16	6
11	15	6
10	14	5
14	13	6
19	12	6
14	16	4
8	9	5
15	15	8
8	16	6
8	12	6
6	11	2
7	13	2
16	13	4
15	14	6
10	15	6
8	14	5
9	12	4
8	16	4
14	14	6
14	13	5
14	12	6
15	13	7
7	12	6
7	9	4
12	13	4
7	10	3
12	15	8
6	9	4
10	13	4
12	13	5
13	13	5
14	15	7
8	13	4
14	14	5
10	11	5
14	15	8
15	14	5
10	15	2
6	12	5
9	15	4
11	14	5
16	16	7
14	14	6
8	12	3
16	11	5
16	13	6
14	12	5
12	12	6
16	16	7
15	13	6
11	12	6
6	14	5
6	4	4
16	14	6
16	15	6
8	12	3
11	11	4
12	12	4
13	11	4
11	12	5
9	11	4
15	13	6
11	12	6
12	12	4
15	15	7
8	14	4
7	12	4
10	12	4
9	12	4
13	13	5
11	11	4
12	13	7
5	12	3
12	14	5
14	15	5
15	15	6
14	13	5
13	16	6
14	17	6
14	13	3
15	14	6
13	13	5
14	16	8
11	13	6
14	14	4
11	13	3
8	14	4
12	16	7

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')