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

41	34	38	30
39	33	32	28
30	29	35	31
31	34	33	35
34	32	37	35
35	35	29	37
39	41	31	39
34	27	36	31
36	40	35	38
37	40	38	37
38	36	31	37
36	40	34	35
38	43	35	37
39	40	38	42
33	33	37	28
32	37	33	37
36	32	32	36
38	26	38	37
39	36	38	33
32	39	32	40
32	38	33	30
31	34	31	36
39	35	38	33
37	41	39	40
39	42	32	37
41	36	32	37
36	39	35	39
33	33	37	35
33	33	33	36
34	36	33	34
31	37	28	36
27	36	32	32
37	34	31	33
34	32	37	27
34	35	30	37
32	39	33	32
29	30	31	31
36	25	33	31
29	29	31	32
35	39	33	37
37	31	32	25
34	26	33	30
38	28	32	37
35	40	33	37
38	32	28	40
37	35	35	35
38	32	39	35
33	41	34	43
36	34	38	32
38	36	32	42
32	38	38	35
32	34	30	27
32	32	33	30
34	34	38	31
32	32	32	36
37	40	32	36
39	43	34	41
29	35	34	34
37	45	36	36
35	36	34	33
30	39	28	35
38	31	34	28
34	36	35	33
31	36	35	38
34	37	31	37
35	40	37	39
36	35	35	34
30	36	27	32
39	32	40	36
35	36	37	36
38	37	36	35
31	42	38	33
34	37	39	42
38	36	41	36
34	36	27	33
39	33	30	36
37	37	37	32
34	35	31	35
28	37	31	38
37	28	27	33
33	33	36	32
37	45	38	38
35	38	37	39
37	43	33	39
32	37	34	39
33	36	31	30
38	40	39	38
33	39	34	38
29	43	32	42
33	32	33	41
31	37	36	31
36	34	32	39
35	44	41	40
32	35	28	31
29	34	30	34
39	37	36	23
37	40	35	28
35	36	31	36
37	44	34	41
32	35	36	29
38	34	36	31
37	40	35	33
36	34	37	35
32	39	28	35
33	36	39	34
40	40	32	40
38	37	35	34
41	35	39	36
36	45	35	35
43	39	42	39

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