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

2	73	69
1	58	53
1	68	43
2	69	60
1	62	49
1	68	62
1	65	45
1	65	50
1	81	75
1	73	82
2	64	60
1	68	59
1	51	21
1	68	40
1	61	62
2	77	54
1	69	47
1	73	59
2	61	37
2	62	43
1	63	48
1	69	59
2	47	79
2	66	62
1	58	16
2	63	38
1	69	58
2	59	60
1	59	72
2	63	67
2	65	55
1	65	47
2	71	59
1	60	49
2	66	47
1	67	57
2	81	39
1	62	49
1	63	26
2	73	53
2	55	75
1	59	65
1	64	49
2	63	48
2	74	45
2	67	31
1	64	67
1	73	61
1	54	49
1	76	69
2	74	54
2	63	80
2	73	57
2	67	34
2	68	69
1	66	44
2	62	70
2	71	51
1	68	66
1	63	18
1	75	74
1	77	59
2	62	48
1	74	55
2	67	44
2	56	56
2	60	65
2	58	77
1	65	46
2	49	70
1	61	39
2	66	55
2	64	44
2	65	45
1	46	45
2	81	25
2	65	49
1	72	65
2	65	45
2	74	71
1	69	48
1	59	41
2	58	40
1	71	64
2	79	56
2	68	52
1	66	41
2	62	45
1	69	49
1	60	42
2	63	54
1	62	40
1	61	40
2	65	51
1	64	48
2	67	80
2	56	38
2	56	57
1	48	28
1	74	51
1	69	46
1	62	58
1	73	67
1	64	72
1	57	26
1	57	54
2	60	53
2	61	69
1	72	64
1	57	47
1	51	43
1	63	66
1	54	54
1	72	62
1	62	52
1	68	64
1	62	55
2	63	57
1	77	74
1	57	32
1	57	38
1	61	66
2	66	37
1	65	26
1	63	64
1	59	28
2	66	66
1	68	65
1	72	48
1	68	44
2	68	64
1	67	39
1	59	50
1	56	52
1	62	66
2	55	48
2	72	70
2	68	66
1	67	61
1	54	31
2	69	61
1	61	54
1	55	34
2	75	62
1	55	47
1	49	52
2	54	37
1	51	46
1	66	50
1	73	61
2	63	70
2	61	38
1	74	63
2	81	34
1	58	46
1	62	40
1	64	30
1	62	35
1	85	51
1	74	56
1	51	68
1	66	39
2	61	44
1	72	58

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