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

210907	56	396	145	0
120982	56	297	101	0
176508	54	559	98	0
179321	89	967	132	NA
123185	40	270	60	NA
52746	25	143	38	NA
385534	92	1562	144	NA
33170	18	109	5	NA
101645	63	371	28	NA
149061	44	656	84	0
165446	33	511	79	NA
237213	84	655	127	0
173326	88	465	78	NA
133131	55	525	60	1
258873	60	885	131	NA
180083	66	497	84	NA
324799	154	1436	133	0
230964	53	612	150	0
236785	119	865	91	1
135473	41	385	132	NA
202925	61	567	136	NA
215147	58	639	124	NA
344297	75	963	118	0
153935	33	398	70	NA
132943	40	410	107	NA
174724	92	966	119	0
174415	100	801	89	0
225548	112	892	112	NA
223632	73	513	108	0
124817	40	469	52	0
221698	45	683	112	NA
210767	60	643	116	NA
170266	62	535	123	NA
260561	75	625	125	NA
84853	31	264	27	NA
294424	77	992	162	0
101011	34	238	32	NA
215641	46	818	64	NA
325107	99	937	92	NA
7176	17	70	0	NA
167542	66	507	83	NA
106408	30	260	41	0
96560	76	503	47	0
265769	146	927	120	0
269651	67	1269	105	NA
149112	56	537	79	1
175824	107	910	65	0
152871	58	532	70	0
111665	34	345	55	0
116408	61	918	39	NA
362301	119	1635	67	0
78800	42	330	21	NA
183167	66	557	127	0
277965	89	1178	152	NA
150629	44	740	113	NA
168809	66	452	99	NA
24188	24	218	7	NA
329267	259	764	141	0
65029	17	255	21	NA
101097	64	454	35	NA
218946	41	866	109	1
244052	68	574	133	0
341570	168	1276	123	1
103597	43	379	26	0
233328	132	825	230	NA
256462	105	798	166	NA
206161	71	663	68	NA
311473	112	1069	147	NA
235800	94	921	179	NA
177939	82	858	61	NA
207176	70	711	101	NA
196553	57	503	108	0
174184	53	382	90	NA
143246	103	464	114	0
187559	121	717	103	NA
187681	62	690	142	NA
119016	52	462	79	NA
182192	52	657	88	NA
73566	32	385	25	NA
194979	62	577	83	NA
167488	45	619	113	0
143756	46	479	118	0
			110	NA
			129	NA
			51	0
			93	NA
			76	0
			49	NA
			118	0
			38	NA
			141	0
			58	NA
			27	0
			91	0
			48	1
			63	0
			56	NA
			144	NA
			73	NA
			168	1
			64	1
			97	0
			117	1
			100	0
			149	NA
			187	NA
			127	NA
			37	0
			245	NA
			87	1
			177	NA
			49	NA
			49	NA
			73	NA
			177	0
			94	NA
			117	NA
			60	0
			55	0
			39	0
			64	1
			26	0
			64	NA
			58	1
			95	NA
			25	NA
			26	NA
			76	0
			129	0
			11	NA
			2	NA
			101	1
			28	NA
			36	0
			89	0
			193	NA
			4	NA
			84	NA
			23	0
			39	NA
			14	NA
			78	1
			14	NA
			101	1
			82	NA
			24	NA
			36	1
			75	1
			16	NA
			55	NA
			131	0
			131	0
			39	0
			144	0
			139	NA
			211	0
			78	NA
			50	NA
			39	NA
			90	0
			166	NA
			12	NA
			57	0
			133	0
			69	0
			119	NA
			119	NA
			65	NA
			61	1
			49	1
			101	0
			196	0
			15	0
			136	0
			89	NA
			40	NA
			123	NA
			21	0
			163	0
			29	0
			35	NA
			13	NA
			5	NA
			96	1
			151	NA
			6	NA
			13	NA
			3	NA
			56	NA
			23	NA
			57	NA
			14	1
			43	NA
			20	0
			72	NA
			87	1
			21	NA
			56	0
			59	NA
			82	NA
			43	NA
			25	NA
			38	NA
			25	0
			38	NA
			12	NA
			29	NA
			47	0
			45	0
			40	NA
			30	NA
			41	0
			25	NA
			23	NA
			14	0
			16	0
			26	NA
			21	0
			27	NA
			9	NA
			33	0
			42	0
			68	0
			32	NA
			6	NA
			67	NA
			33	NA
			77	0
			46	NA
			30	NA
			0	0
			36	0
			46	NA
			18	NA
			48	1
			29	0
			28	1
			34	NA
			33	0
			34	NA
			33	NA
			80	NA
			32	NA
			30	NA
			41	NA
			41	NA
			51	0
			18	0
			34	NA
			31	NA
			39	0
			54	0
			14	NA
			24	0
			24	0
			8	NA
			26	NA
			19	NA
			11	NA
			14	NA
			1	NA
			39	NA
			5	NA
			37	0
			32	NA
			38	NA
			47	0
			47	0
			37	0
			51	NA
			45	NA
			21	NA
			1	NA
			42	NA
			26	NA
			21	NA
			4	NA
			10	NA
			43	NA
			34	0
			31	NA
			19	NA
			34	NA
			6	NA
			11	0
			24	NA
			16	NA
			72	NA
			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')