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

4	30	79	146283	1418
NA	28	58	98364	869
1	38	60	86146	1530
NA	30	108	96933	2172
NA	22	49	79234	901
NA	26	0	42551	463
3	25	121	195663	3201
NA	18	1	6853	371
NA	11	20	21529	1192
3	26	43	95757	1583
4	25	69	85584	1439
4	38	78	143983	1764
NA	44	86	75851	1495
3	30	44	59238	1373
NA	40	104	93163	2187
NA	34	63	96037	1491
2	47	158	151511	4041
4	30	102	136368	1706
4	31	77	112642	2152
NA	23	82	94728	1036
NA	36	115	105499	1882
4	36	101	121527	1929
5	30	80	127766	2242
NA	25	50	98958	1220
NA	39	83	77900	1289
4	34	123	85646	2515
4	31	73	98579	2147
5	31	81	130767	2352
4	33	105	131741	1638
4	25	47	53907	1222
NA	33	105	178812	1812
2	35	94	146761	1677
3	42	44	82036	1579
NA	43	114	163253	1731
NA	30	38	27032	807
4	33	107	171975	2452
NA	13	30	65990	829
NA	32	71	86572	1940
2	36	84	159676	2662
4	0	0	1929	186
NA	28	59	85371	1499
2	14	33	58391	865
4	17	42	31580	1793
4	32	96	136815	2527
NA	30	106	120642	2747
4	35	56	69107	1324
3	20	57	50495	2702
4	28	59	108016	1383
5	28	39	46341	1179
NA	39	34	78348	2099
4	34	76	79336	4308
NA	26	20	56968	918
5	39	91	93176	1831
NA	39	115	161632	3373
NA	33	85	87850	1713
NA	28	76	127969	1438
NA	4	8	15049	496
2	39	79	155135	2253
NA	18	21	25109	744
NA	14	30	45824	1161
3	29	76	102996	2352
4	44	101	160604	2144
3	21	94	158051	4691
4	16	27	44547	1112
NA	28	92	162647	2694
4	35	123	174141	1973
NA	28	75	60622	1769
NA	38	128	179566	3148
NA	23	105	184301	2474
NA	36	55	75661	2084
NA	32	56	96144	1954
4	29	41	129847	1226
NA	25	72	117286	1389
5	27	67	71180	1496
NA	36	75	109377	2269
4	28	114	85298	1833
NA	23	118	73631	1268
NA	40	77	86767	1943
2	23	22	23824	893
NA	40	66	93487	1762
4	28	69	82981	1403
4	34	105	73815	1425
NA	33	116	94552	1857
4	28	88	132190	1840
NA	34	73	128754	1502
NA	30	99	66363	1441
4	33	62	67808	1420
NA	22	53	61724	1416
4	38	118	131722	2970
NA	26	30	68580	1317
4	35	100	106175	1644
NA	8	49	55792	870
2	24	24	25157	1654
2	29	67	76669	1054
4	20	46	57283	937
3	29	57	105805	3004
NA	45	75	129484	2008
4	37	135	72413	2547
NA	33	68	87831	1885
2	33	124	96971	1626
3	25	33	71299	1468
2	32	98	77494	2445
4	29	58	120336	1964
4	28	68	93913	1381
NA	28	81	136048	1369
4	31	131	181248	1659
3	52	110	146123	2888
4	21	37	32036	1290
NA	24	130	186646	2845
2	41	93	102255	1982
5	33	118	168237	1904
NA	32	39	64219	1391
NA	19	13	19630	602
NA	20	74	76825	1743
4	31	81	115338	1559
NA	31	109	109427	2014
NA	32	151	118168	2143
4	18	51	84845	2146
2	23	28	153197	874
4	17	40	29877	1590
4	20	56	63506	1590
2	12	27	22445	1210
NA	17	37	47695	2072
4	30	83	68370	1281
NA	31	54	146304	1401
NA	10	27	38233	834
3	13	28	42071	1105
5	22	59	50517	1272
4	42	133	103950	1944
2	1	12	5841	391
NA	9	0	2341	761
4	32	106	84396	1605
NA	11	23	24610	530
4	25	44	35753	1988
4	36	71	55515	1386
4	31	116	209056	2395
NA	0	4	6622	387
NA	24	62	115814	1742
2	13	12	11609	620
NA	8	18	13155	449
4	13	14	18274	800
4	19	60	72875	1684
NA	18	7	10112	1050
4	33	98	142775	2699
NA	40	64	68847	1606
NA	22	29	17659	1502
3	38	32	20112	1204
2	24	25	61023	1138
NA	8	16	13983	568
NA	35	48	65176	1459
4	43	100	132432	2158
3	43	46	112494	1111
4	14	45	45109	1421
5	41	129	170875	2833
NA	38	130	180759	1955
4	45	136	214921	2922
NA	31	59	100226	1002
NA	13	25	32043	1060
NA	28	32	54454	956
4	31	63	78876	2186
NA	40	95	170745	3604
NA	30	14	6940	1035
4	16	36	49025	1417
4	37	113	122037	3261
4	30	47	53782	1587
2	35	92	127748	1424
NA	32	70	86839	1701
NA	27	19	44830	1249
2	20	50	77395	946
3	18	41	89324	1926
NA	31	91	103300	3352
1	31	111	112283	1641
3	21	41	10901	2035
3	39	120	120691	2312
NA	41	135	58106	1369
NA	13	27	57140	1577
NA	32	87	122422	2201
3	18	25	25899	961
4	39	131	139296	1900
2	14	45	52678	1254
4	7	29	23853	1335
3	17	58	17306	1597
NA	0	4	7953	207
4	30	47	89455	1645
4	37	109	147866	2429
NA	0	7	4245	151
NA	5	12	21509	474
NA	1	0	7670	141
NA	16	37	66675	1639
2	32	37	14336	872
NA	24	46	53608	1318
4	17	15	30059	1018
NA	11	42	29668	1383
4	24	7	22097	1314
4	22	54	96841	1335
4	12	54	41907	1403
4	19	14	27080	910
3	13	16	35885	616
NA	17	33	41247	1407
NA	15	32	28313	771
NA	16	21	36845	766
NA	24	15	16548	473
2	15	38	36134	1376
5	17	22	55764	1232
NA	18	28	28910	1521
NA	20	10	13339	572
NA	16	31	25319	1059
4	16	32	66956	1544
4	18	32	47487	1230
NA	22	43	52785	1206
NA	8	27	44683	1205
2	17	37	35619	1255
NA	18	20	21920	613
3	16	32	45608	721
3	23	0	7721	1109
4	22	5	20634	740
NA	13	26	29788	1126
3	13	10	31931	728
NA	16	27	37754	689
NA	16	11	32505	592
2	20	29	40557	995
5	22	25	94238	1613
4	17	55	44197	2048
NA	18	23	43228	705
NA	17	5	4103	301
NA	12	43	44144	1803
NA	7	23	32868	799
4	17	34	27640	861
NA	14	36	14063	1186
5	23	35	28990	1451
4	17	0	4694	628
4	14	37	42648	1161
NA	15	28	64329	1463
NA	17	16	21928	742
4	21	26	25836	979
2	18	38	22779	675
2	18	23	40820	1241
NA	17	22	27530	676
4	17	30	32378	1049
NA	16	16	10824	620
NA	15	18	39613	1081
NA	21	28	60865	1688
NA	16	32	19787	736
2	14	21	20107	617
NA	15	23	36605	812
NA	17	29	40961	1051
4	15	50	48231	1656
2	15	12	39725	705
NA	10	21	21455	945
NA	6	18	23430	554
4	22	27	62991	1597
4	21	41	49363	982
NA	1	13	9604	222
4	18	12	24552	1212
4	17	21	31493	1143
4	4	8	3439	435
5	10	26	19555	532
NA	16	27	21228	882
NA	16	13	23177	608
NA	9	16	22094	459
NA	16	2	2342	578
NA	17	42	38798	826
NA	7	5	3255	509
4	15	37	24261	717
NA	14	17	18511	637
NA	14	38	40798	857
4	18	37	28893	830
2	12	29	21425	652
2	16	32	50276	707
NA	21	35	37643	954
NA	19	17	30377	1461
NA	16	20	27126	672
NA	1	7	13	778
NA	16	46	42097	1141
NA	10	24	24451	680
NA	19	40	14335	1090
NA	12	3	5084	616
NA	2	10	9927	285
NA	14	37	43527	1145
4	17	17	27184	733
NA	19	28	21610	888
NA	14	19	20484	849
NA	11	29	20156	1182
NA	4	8	6012	528
2	16	10	18475	642
NA	20	15	12645	947
NA	12	15	11017	819
NA	15	28	37623	757
NA	16	17	35873	894

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