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

119565	20	56	18158	30	28
127145	38	47	30461	42	39
7215	0	0	1423	0	0
112861	49	51	25629	54	54
197581	70	93	48758	86	80
377410	104	98	129230	157	144
117604	37	59	27376	36	36
120102	46	90	26706	48	48
96175	42	40	26505	45	42
253517	62	79	49801	77	71
108994	50	49	46580	49	49
156212	65	74	48352	77	74
68810	28	32	13899	28	27
149246	48	78	39342	84	83
125503	42	43	27465	31	31
125769	47	65	55211	28	28
123467	71	111	74098	99	98
56088	0	36	13497	2	2
108128	50	89	38338	41	43
22762	12	28	52505	25	24
48554	16	35	10663	16	16
159102	73	68	74484	96	95
139108	29	67	28895	23	22
92058	38	61	32827	33	33
126311	50	46	36188	46	45
113807	33	49	28173	59	59
138834	44	74	54926	72	66
135313	59	66	38900	72	70
279140	47	73	88530	62	56
158146	40	53	35482	55	55
107352	33	71	26730	27	27
149827	50	54	29806	41	37
90912	41	25	41799	51	48
171030	73	59	54289	26	26
162204	43	73	36805	65	64
0	0	0	0	0	0
163310	41	70	33146	28	21
92342	44	39	23333	44	44
98357	31	83	47686	36	36
178277	71	123	77783	100	89
134927	61	54	36042	104	101
104577	27	105	34541	35	31
81614	21	67	75620	69	65
119111	42	54	60610	73	71
83923	40	68	55041	106	102
109885	34	57	32087	53	53
52851	15	51	16356	43	41
153621	46	72	40161	49	46
152572	43	94	55459	38	37
120865	41	72	36679	51	51
48188	12	39	22346	14	14
90994	42	60	27377	40	40
215588	56	72	50273	79	77
176489	41	67	32104	52	51
147174	47	102	27016	44	43
104416	30	28	19715	34	33
156186	44	65	33629	47	47
60368	25	59	27084	32	31
95391	42	73	32352	31	31
126195	28	74	51845	40	40
96388	33	103	26591	42	42
77353	32	52	29677	34	35
89523	26	44	54237	40	40
79772	30	52	20284	35	30
96945	13	68	22741	11	11
88300	35	77	34178	43	41
126203	39	115	69551	53	53
110681	68	59	29653	82	82
81299	32	66	38071	41	41
31970	5	42	4157	6	6
185321	53	35	28321	82	81
87611	33	3	40195	47	47
73165	44	68	48158	108	100
82167	36	38	13310	46	46
113432	52	97	78474	38	38
49164	0	69	6386	0	0
105181	49	75	31588	45	45
105922	29	66	61254	57	56
60138	16	46	21152	20	18
73422	33	52	41272	56	54
67727	48	58	34165	38	37
193980	33	72	37054	42	40
51185	24	13	12368	37	37
84448	33	49	23168	36	36
41956	16	49	16380	34	34
110827	32	41	41242	53	49
167055	41	93	48259	83	80
63603	36	54	20790	36	36
64995	22	64	34585	33	33
93424	26	57	35672	57	55
97229	37	57	52168	50	50
112819	57	51	53933	71	71
111538	20	71	34474	32	31
102220	24	43	43753	45	42
38721	18	18	36456	33	31
168764	37	87	51183	53	51
161806	83	76	52742	64	64
19349	13	0	3895	14	14
128694	18	83	37076	38	37
101954	32	58	24079	39	37
43803	8	4	2325	8	8
47062	38	41	29354	38	38
104220	41	51	30341	24	23
85939	24	45	18992	22	22
58660	23	24	15292	18	18
27676	2	17	5842	3	1
93448	52	80	28918	49	48
43284	5	20	3738	5	5
0	0	0	0	0	0
64333	43	43	95352	47	46
57050	18	63	37478	33	33
96933	41	48	26839	44	41
70088	45	70	26783	56	57
65494	29	32	33392	49	49
3616	0	0	0	0	0
0	0	0	0	0	0
135104	32	69	25446	45	45
95554	41	56	60038	80	80
120307	17	62	28162	51	46
84336	24	77	33298	25	25
43410	7	3	2781	1	1
131452	62	65	37121	62	59
79015	30	37	22698	29	29
88043	49	49	27615	26	26
57586	3	32	32689	4	4
19764	10	4	5752	10	10
105757	42	55	23164	43	43
96410	18	81	20304	36	36
109491	40	90	34409	43	41
11796	1	1	0	0	0
7627	0	0	0	0	0
117413	29	38	92538	33	32
6836	0	0	0	0	0
134824	45	36	46037	53	53
5118	5	0	0	0	0
40248	8	7	5444	6	6
0	0	0	0	0	0
89017	21	67	23924	19	18
67140	16	45	52230	26	26
7131	0	0	0	0	0
4194	0	0	0	0	0
60378	15	45	8019	16	16
96971	40	60	34542	84	84
83484	17	48	21157	28	22

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