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

210907	56	30	145	0
120982	56	28	101	0
176508	54	38	98	0
179321	89	30	132	NA
123185	40	22	60	NA
52746	25	26	38	NA
385534	92	25	144	NA
33170	18	18	5	NA
101645	63	11	28	NA
149061	44	26	84	0
165446	33	25	79	NA
237213	84	38	127	0
173326	88	44	78	NA
133131	55	30	60	1
258873	60	40	131	NA
180083	66	34	84	NA
324799	154	47	133	0
230964	53	30	150	0
236785	119	31	91	1
135473	41	23	132	NA
202925	61	36	136	NA
215147	58	36	124	NA
344297	75	30	118	0
153935	33	25	70	NA
132943	40	39	107	NA
174724	92	34	119	0
174415	100	31	89	0
225548	112	31	112	NA
223632	73	33	108	0
124817	40	25	52	0
221698	45	33	112	NA
210767	60	35	116	NA
170266	62	42	123	NA
260561	75	43	125	NA
84853	31	30	27	NA
294424	77	33	162	0
101011	34	13	32	NA
215641	46	32	64	NA
325107	99	36	92	NA
7176	17	0	0	NA
167542	66	28	83	NA
106408	30	14	41	0
96560	76	17	47	0
265769	146	32	120	0
269651	67	30	105	NA
149112	56	35	79	1
175824	107	20	65	0
152871	58	28	70	0
111665	34	28	55	0
116408	61	39	39	NA
362301	119	34	67	0
78800	42	26	21	NA
183167	66	39	127	0
277965	89	39	152	NA
150629	44	33	113	NA
168809	66	28	99	NA
24188	24	4	7	NA
329267	259	39	141	0
65029	17	18	21	NA
101097	64	14	35	NA
218946	41	29	109	1
244052	68	44	133	0
341570	168	21	123	1
103597	43	16	26	0
233328	132	28	230	NA
256462	105	35	166	NA
206161	71	28	68	NA
311473	112	38	147	NA
235800	94	23	179	NA
177939	82	36	61	NA
207176	70	32	101	NA
196553	57	29	108	0
174184	53	25	90	NA
143246	103	27	114	0
187559	121	36	103	NA
187681	62	28	142	NA
119016	52	23	79	NA
182192	52	40	88	NA
73566	32	23	25	NA
194979	62	40	83	NA
167488	45	28	113	0
143756	46	34	118	0
		33	110	NA
		28	129	NA
		34	51	0
		30	93	NA
		33	76	0
		22	49	NA
		38	118	0
		26	38	NA
		35	141	0
		8	58	NA
		24	27	0
		29	91	0
		20	48	1
		29	63	0
		45	56	NA
		37	144	NA
		33	73	NA
		33	168	1
		25	64	1
		32	97	0
		29	117	1
		28	100	0
		28	149	NA
		31	187	NA
		52	127	NA
		21	37	0
		24	245	NA
		41	87	1
		33	177	NA
		32	49	NA
		19	49	NA
		20	73	NA
		31	177	0
		31	94	NA
		32	117	NA
		18	60	0
		23	55	0
		17	39	0
		20	64	1
		12	26	0
		17	64	NA
		30	58	1
		31	95	NA
		10	25	NA
		13	26	NA
		22	76	0
		42	129	0
		1	11	NA
		9	2	NA
		32	101	1
		11	28	NA
		25	36	0
		36	89	0
		31	193	NA
		0	4	NA
		24	84	NA
		13	23	0
		8	39	NA
		13	14	NA
		19	78	1
		18	14	NA
		33	101	1
		40	82	NA
		22	24	NA
		38	36	1
		24	75	1
		8	16	NA
		35	55	NA
		43	131	0
		43	131	0
		14	39	0
		41	144	0
		38	139	NA
		45	211	0
		31	78	NA
		13	50	NA
		28	39	NA
		31	90	0
		40	166	NA
		30	12	NA
		16	57	0
		37	133	0
		30	69	0
		35	119	NA
		32	119	NA
		27	65	NA
		20	61	1
		18	49	1
		31	101	0
		31	196	0
		21	15	0
		39	136	0
		41	89	NA
		13	40	NA
		32	123	NA
		18	21	0
		39	163	0
		14	29	0
		7	35	NA
		17	13	NA
		0	5	NA
		30	96	1
		37	151	NA
		0	6	NA
		5	13	NA
		1	3	NA
		16	56	NA
		32	23	NA
		24	57	NA
		17	14	1
		11	43	NA
		24	20	0
		22	72	NA
		12	87	1
		19	21	NA
		13	56	0
		17	59	NA
		15	82	NA
		16	43	NA
		24	25	NA
		15	38	NA
		17	25	0
		18	38	NA
		20	12	NA
		16	29	NA
		16	47	0
		18	45	0
		22	40	NA
		8	30	NA
		17	41	0
		18	25	NA
		16	23	NA
		23	14	0
		22	16	0
		13	26	NA
		13	21	0
		16	27	NA
		16	9	NA
		20	33	0
		22	42	0
		17	68	0
		18	32	NA
		17	6	NA
		12	67	NA
		7	33	NA
		17	77	0
		14	46	NA
		23	30	NA
		17	0	0
		14	36	0
		15	46	NA
		17	18	NA
		21	48	1
		18	29	0
		18	28	1
		17	34	NA
		17	33	0
		16	34	NA
		15	33	NA
		21	80	NA
		16	32	NA
		14	30	NA
		15	41	NA
		17	41	NA
		15	51	0
		15	18	0
		10	34	NA
		6	31	NA
		22	39	0
		21	54	0
		1	14	NA
		18	24	0
		17	24	0
		4	8	NA
		10	26	NA
		16	19	NA
		16	11	NA
		9	14	NA
		16	1	NA
		17	39	NA
		7	5	NA
		15	37	0
		14	32	NA
		14	38	NA
		18	47	0
		12	47	0
		16	37	0
		21	51	NA
		19	45	NA
		16	21	NA
		1	1	NA
		16	42	NA
		10	26	NA
		19	21	NA
		12	4	NA
		2	10	NA
		14	43	NA
		17	34	0
		19	31	NA
		14	19	NA
		11	34	NA
		4	6	NA
		16	11	0
		20	24	NA
		12	16	NA
		15	72	NA
		16	21	NA

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