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

2	14	53
2	18	86
2	11	66
1	12	67
2	16	76
2	18	78
2	14	53
2	14	80
2	15	74
2	15	76
1	17	79
2	19	54
1	10	67
2	16	54
2	18	87
1	14	58
1	14	75
2	17	88
1	14	64
2	16	57
1	18	66
2	11	68
2	14	54
2	12	56
1	17	86
2	9	80
1	16	76
2	14	69
2	15	78
1	11	67
2	16	80
1	13	54
2	17	71
2	15	84
1	14	74
1	16	71
1	9	63
1	15	71
2	17	76
1	13	69
1	15	74
2	16	75
1	16	54
1	12	52
2	12	69
2	11	68
2	15	65
2	15	75
2	17	74
1	13	75
2	16	72
1	14	67
1	11	63
2	12	62
1	12	63
2	15	76
2	16	74
2	15	67
1	12	73
2	12	70
1	8	53
1	13	77
2	11	77
2	14	52
2	15	54
1	10	80
2	11	66
1	12	73
2	15	63
1	15	69
1	14	67
2	16	54
2	15	81
1	15	69
1	13	84
2	12	80
2	17	70
2	13	69
1	15	77
1	13	54
1	15	79
1	16	30
2	15	71
1	16	73
2	15	72
2	14	77
1	15	75
2	14	69
2	13	54
2	7	70
2	17	73
2	13	54
2	15	77
2	14	82
2	13	80
2	16	80
2	12	69
2	14	78
1	17	81
1	15	76
2	17	76
1	12	73
2	16	85
1	11	66
2	15	79
1	9	68
2	16	76
1	15	71
1	10	54
2	10	46
2	15	82
2	11	74
2	13	88
1	14	38
2	18	76
1	16	86
2	14	54
2	14	70
2	14	69
2	14	90
2	12	54
2	14	76
2	15	89
2	15	76
2	15	73
2	13	79
1	17	90
2	17	74
2	19	81
2	15	72
1	13	71
1	9	66
2	15	77
1	15	65
1	15	74
2	16	82
1	11	54
1	14	63
2	11	54
2	15	64
1	13	69
2	15	54
1	16	84
2	14	86
1	15	77
2	16	89
2	16	76
1	11	60
1	12	75
1	9	73
2	16	85
2	13	79
1	16	71
2	12	72
2	9	69
2	13	78
2	13	54
2	14	69
2	19	81
2	13	84
2	12	84
2	13	69

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