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

150596	535	18	20465	37
154801	396	20	33629	43
7215	72	0	1423	0
122139	617	26	25629	54
221399	1118	30	54002	86
441870	1755	34	151036	181
134379	498	23	33287	42
140428	355	30	31172	59
103255	413	30	28113	46
271630	891	26	57803	77
121593	629	24	49830	49
172071	611	30	52143	79
83707	564	19	21055	37
197412	964	25	47007	92
134398	362	17	28735	31
139224	442	19	59147	28
134153	391	33	78950	103
64149	305	15	13497	2
122294	721	34	46154	48
24889	206	15	53249	25
52197	310	15	10726	16
188915	686	27	83700	106
172874	590	25	40400	35
98575	558	34	33797	33
143546	569	21	36205	45
139780	513	21	30165	64
163784	602	25	58534	73
152479	276	28	44663	78
304108	791	26	92556	63
184024	815	20	40078	69
151621	427	28	34711	36
164516	496	20	31076	41
120289	655	17	74608	59
214701	857	25	58092	33
196865	736	24	42009	76
0	0	0	0	0
191678	884	27	36022	27
93107	483	14	23333	44
129352	495	32	53349	43
229143	749	31	92596	104
177063	627	21	49598	120
126602	597	34	44093	44
93742	348	23	84205	71
152153	711	24	63369	78
95704	322	22	60132	106
139793	280	22	37403	61
76348	205	35	24460	53
188980	648	21	46456	51
172100	580	31	66616	46
146552	875	26	41554	55
48188	205	22	22346	14
109185	363	21	30874	44
263652	757	27	68701	113
215609	647	26	35728	55
174876	584	33	29010	46
115124	457	11	23110	39
179712	438	26	38844	51
70369	235	26	27084	31
109215	312	21	35139	36
166096	877	38	57476	47
130414	454	29	33277	53
102057	668	19	31141	38
115310	346	19	61281	52
101181	377	24	25820	37
135228	365	26	23284	11
94982	391	29	35378	45
166919	476	34	74990	59
118169	747	25	29653	82
102361	246	24	64622	49
31970	101	21	4157	6
200413	901	19	29245	81
103381	334	12	50008	56
94940	404	28	52338	105
101560	442	21	13310	46
144176	627	34	92901	46
71921	345	32	10956	2
126905	538	27	34241	51
131184	741	26	75043	95
60138	253	21	21152	18
84971	395	31	42249	55
80420	211	26	42005	48
233569	670	26	41152	48
56252	244	23	14399	39
97181	438	25	28263	40
50800	255	22	17215	36
125941	434	26	48140	60
211032	613	33	62897	114
71960	233	22	22883	39
90379	360	24	41622	45
125650	486	21	40715	59
115572	535	28	65897	59
136266	585	22	76542	93
146715	402	22	37477	35
124626	466	15	53216	47
49176	291	13	40911	36
212926	691	36	57021	59
173884	515	24	73116	79
19349	67	1	3895	14
181141	712	24	46609	42
145502	770	31	29351	41
45448	247	4	2325	8
58280	240	20	31747	41
115944	360	23	32665	24
94341	249	23	19249	22
59090	138	12	15292	18
27676	194	16	5842	1
120586	285	28	33994	53
88011	227	10	13018	6
0	0	0	0	0
85610	306	25	98177	49
94530	355	21	37941	33
117769	397	21	31032	50
107653	369	21	32683	64
71894	287	21	34545	53
3616	14	0	0	0
0	0	0	0	0
154806	301	23	27525	48
136061	535	29	66856	90
141822	530	27	28549	46
106515	272	23	38610	29
43410	292	1	2781	1
146920	458	25	41211	64
88874	241	17	22698	29
111924	497	29	41194	27
60373	165	12	32689	4
19764	75	2	5752	10
121665	461	18	26757	47
108685	341	25	22527	44
124493	446	29	44810	51
11796	79	2	0	0
10674	33	0	0	0
131263	449	18	100674	38
6836	11	1	0	0
153278	606	21	57786	57
5118	6	0	0	0
40248	183	4	5444	6
0	0	0	0	0
100728	310	25	28470	22
84267	245	26	61849	34
7131	27	0	0	0
8812	97	4	2179	10
63952	247	17	8019	16
120111	273	21	39644	93
94127	386	22	23494	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')