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

101,5	102,7	91,4	99,4
99,2	116,6	85,4	102,7
107,8	114,7	110,4	109,3
92,3	74,2	90,1	93,9
99,2	73,9	103,4	95,3
101,6	81,3	118,7	101,8
87,0	85,8	76,4	85,6
71,4	113,7	92,7	81,1
104,7	101,7	105,2	109,5
115,1	85,2	91,5	104,0
102,5	62,5	75,3	94,5
75,3	61,2	60,5	79,0
96,7	87,1	80,4	92,8
94,6	96,5	84,5	95,6
98,6	88,4	93,9	101,7
99,5	67,0	78,0	90,8
92,0	58,6	92,3	89,5
93,6	62,5	90,0	91,8
89,3	86,1	72,1	83,8
66,9	98,7	76,9	77,4
108,8	85,3	76,0	112,7
113,2	73,8	88,7	98,8
105,5	55,2	55,4	85,7
77,8	54,4	46,6	72,8
102,1	80,8	90,9	96,9
97,0	88,0	84,9	95,0
95,5	74,5	89,0	94,2
99,3	55,1	90,2	87,3
86,4	47,2	72,3	80,6
92,4	54,2	83,0	87,9
85,7	70,6	71,6	79,6
61,9	78,5	75,4	71,9
104,9	77,1	85,1	94,6
107,9	56,6	81,2	91,4
95,6	39,8	68,7	86,6
79,8	44,1	68,4	68,5
94,8	66,9	93,7	90,1
93,7	75,3	96,6	91,6
108,1	74,9	101,8	95,4
96,9	48,8	93,6	85,4
88,8	37,0	88,9	81,6
106,7	49,8	114,1	88,9
86,8	63,2	82,3	84,1
69,8	75,9	96,4	74,7
110,9	68,5	104,0	97,1
105,4	49,2	88,2	95,3
99,2	40,3	85,2	85,1
84,4	38,6	87,1	67,3
87,2	54,2	85,5	80,6
91,9	70,6	89,1	87,9
97,9	68,0	105,2	89,2
94,5	43,0	82,9	81,3
85,0	42,3	86,8	79,7
100,3	47,7	112,0	83,7
78,7	57,7	97,4	82,1
65,8	75,8	88,9	69,3
104,8	57,2	109,4	91,2
96,0	43,6	87,8	85,7
103,3	40,0	90,5	85,2
82,9	35,9	79,3	70,0
91,4	59,5	114,9	85,8
94,5	72,7	118,8	91,4
109,3	70,9	125,0	97,5
92,1	44,9	96,1	87,1
99,3	44,5	116,7	85,1
109,6	48,0	119,5	94,1
87,5	60,4	104,1	85,8
73,1	71,8	121,0	74,7
110,7	63,2	127,3	99,9
111,6	32,4	117,7	90,7
110,7	33,9	108,0	86,8
84,0	24,2	89,4	74,8
101,6	64,7	137,4	91,8
102,1	73,0	142,0	97,6
113,9	61,7	137,3	100,8
99,0	31,9	122,8	85,4
100,4	30,8	126,1	84,0
109,5	36,7	147,6	90,6
93,0	47,4	115,7	80,5
76,8	54,0	139,2	73,9
105,3	41,1	150,9	93,6

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='kendall')
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')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Kendall tau rank correlations for all pairs of data series',3,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'pair',1,TRUE)
a<-table.element(a,'tau',1,TRUE)
a<-table.element(a,'p-value',1,TRUE)
a<-table.row.end(a)
n <- length(y[,1])
n
cor.test(y[1,],y[2,],method='kendall')
for (i in 1:(n-1))
{
for (j in (i+1):n)
{
a<-table.row.start(a)
dum <- paste('tau(',dimnames(t(x))[[2]][i])
dum <- paste(dum,',')
dum <- paste(dum,dimnames(t(x))[[2]][j])
dum <- paste(dum,')')
a<-table.element(a,dum,header=TRUE)
r <- cor.test(y[i,],y[j,],method='kendall')
a<-table.element(a,r$estimate)
a<-table.element(a,r$p.value)
a<-table.row.end(a)
}
}
a<-table.end(a)
table.save(a,file='mytable.tab')