Free Statistics

of Irreproducible Research!

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_pairs.wasp
Title produced by softwareKendall tau Correlation Matrix
Date of computationThu, 22 Dec 2011 09:58:18 -0500
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2011/Dec/22/t13245659859emw0jtmjgl5vtw.htm/, Retrieved Fri, 03 May 2024 12:55:34 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=159546, Retrieved Fri, 03 May 2024 12:55:34 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact83
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Kendall tau Correlation Matrix] [Paper Deel 3 Kend...] [2011-12-22 14:58:18] [312c935a345ba403c8b6fcfc5a5ec794] [Current]
Feedback Forum

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




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time2 seconds
R Server'Herman Ole Andreas Wold' @ wold.wessa.net

\begin{tabular}{lllllllll}
\hline
Summary of computational transaction \tabularnewline
Raw Input & view raw input (R code)  \tabularnewline
Raw Output & view raw output of R engine  \tabularnewline
Computing time & 2 seconds \tabularnewline
R Server & 'Herman Ole Andreas Wold' @ wold.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=159546&T=0

[TABLE]
[ROW][C]Summary of computational transaction[/C][/ROW]
[ROW][C]Raw Input[/C][C]view raw input (R code) [/C][/ROW]
[ROW][C]Raw Output[/C][C]view raw output of R engine [/C][/ROW]
[ROW][C]Computing time[/C][C]2 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Herman Ole Andreas Wold' @ wold.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=159546&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=159546&T=0

As an alternative you can also use a QR Code:  

The GUIDs for individual cells are displayed in the table below:

Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time2 seconds
R Server'Herman Ole Andreas Wold' @ wold.wessa.net







Correlations for all pairs of data series (method=kendall)
PVSCBCCRTSTH
PV10.2510.6270.5270.6030.584
SC0.25110.2070.2620.2190.223
BC0.6270.20710.5740.6850.747
CR0.5270.2620.57410.5540.557
TS0.6030.2190.6850.55410.734
TH0.5840.2230.7470.5570.7341

\begin{tabular}{lllllllll}
\hline
Correlations for all pairs of data series (method=kendall) \tabularnewline
  & PV & SC & BC & CR & TS & TH \tabularnewline
PV & 1 & 0.251 & 0.627 & 0.527 & 0.603 & 0.584 \tabularnewline
SC & 0.251 & 1 & 0.207 & 0.262 & 0.219 & 0.223 \tabularnewline
BC & 0.627 & 0.207 & 1 & 0.574 & 0.685 & 0.747 \tabularnewline
CR & 0.527 & 0.262 & 0.574 & 1 & 0.554 & 0.557 \tabularnewline
TS & 0.603 & 0.219 & 0.685 & 0.554 & 1 & 0.734 \tabularnewline
TH & 0.584 & 0.223 & 0.747 & 0.557 & 0.734 & 1 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=159546&T=1

[TABLE]
[ROW][C]Correlations for all pairs of data series (method=kendall)[/C][/ROW]
[ROW][C] [/C][C]PV[/C][C]SC[/C][C]BC[/C][C]CR[/C][C]TS[/C][C]TH[/C][/ROW]
[ROW][C]PV[/C][C]1[/C][C]0.251[/C][C]0.627[/C][C]0.527[/C][C]0.603[/C][C]0.584[/C][/ROW]
[ROW][C]SC[/C][C]0.251[/C][C]1[/C][C]0.207[/C][C]0.262[/C][C]0.219[/C][C]0.223[/C][/ROW]
[ROW][C]BC[/C][C]0.627[/C][C]0.207[/C][C]1[/C][C]0.574[/C][C]0.685[/C][C]0.747[/C][/ROW]
[ROW][C]CR[/C][C]0.527[/C][C]0.262[/C][C]0.574[/C][C]1[/C][C]0.554[/C][C]0.557[/C][/ROW]
[ROW][C]TS[/C][C]0.603[/C][C]0.219[/C][C]0.685[/C][C]0.554[/C][C]1[/C][C]0.734[/C][/ROW]
[ROW][C]TH[/C][C]0.584[/C][C]0.223[/C][C]0.747[/C][C]0.557[/C][C]0.734[/C][C]1[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=159546&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=159546&T=1

As an alternative you can also use a QR Code:  

The GUIDs for individual cells are displayed in the table below:

Correlations for all pairs of data series (method=kendall)
PVSCBCCRTSTH
PV10.2510.6270.5270.6030.584
SC0.25110.2070.2620.2190.223
BC0.6270.20710.5740.6850.747
CR0.5270.2620.57410.5540.557
TS0.6030.2190.6850.55410.734
TH0.5840.2230.7470.5570.7341







Correlations for all pairs of data series with p-values
pairPearson rSpearman rhoKendall tau
PV;SC0.29720.3410.2511
p-value(0)(0)(0)
PV;BC0.76060.82490.6272
p-value(0)(0)(0)
PV;CR0.66520.71850.5273
p-value(0)(0)(0)
PV;TS0.7430.7920.6032
p-value(0)(0)(0)
PV;TH0.70540.77910.5836
p-value(0)(0)(0)
SC;BC0.2380.27590.2074
p-value(0)(0)(0)
SC;CR0.38950.35110.2616
p-value(0)(0)(0)
SC;TS0.25370.29170.2192
p-value(0)(0)(0)
SC;TH0.26320.30140.2228
p-value(0)(0)(0)
BC;CR0.76230.77220.574
p-value(0)(0)(0)
BC;TS0.84010.86760.6852
p-value(0)(0)(0)
BC;TH0.8840.90640.7471
p-value(0)(0)(0)
CR;TS0.73210.76070.5536
p-value(0)(0)(0)
CR;TH0.72060.75470.557
p-value(0)(0)(0)
TS;TH0.88690.90320.7339
p-value(0)(0)(0)

\begin{tabular}{lllllllll}
\hline
Correlations for all pairs of data series with p-values \tabularnewline
pair & Pearson r & Spearman rho & Kendall tau \tabularnewline
PV;SC & 0.2972 & 0.341 & 0.2511 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
PV;BC & 0.7606 & 0.8249 & 0.6272 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
PV;CR & 0.6652 & 0.7185 & 0.5273 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
PV;TS & 0.743 & 0.792 & 0.6032 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
PV;TH & 0.7054 & 0.7791 & 0.5836 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
SC;BC & 0.238 & 0.2759 & 0.2074 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
SC;CR & 0.3895 & 0.3511 & 0.2616 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
SC;TS & 0.2537 & 0.2917 & 0.2192 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
SC;TH & 0.2632 & 0.3014 & 0.2228 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
BC;CR & 0.7623 & 0.7722 & 0.574 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
BC;TS & 0.8401 & 0.8676 & 0.6852 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
BC;TH & 0.884 & 0.9064 & 0.7471 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
CR;TS & 0.7321 & 0.7607 & 0.5536 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
CR;TH & 0.7206 & 0.7547 & 0.557 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
TS;TH & 0.8869 & 0.9032 & 0.7339 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=159546&T=2

[TABLE]
[ROW][C]Correlations for all pairs of data series with p-values[/C][/ROW]
[ROW][C]pair[/C][C]Pearson r[/C][C]Spearman rho[/C][C]Kendall tau[/C][/ROW]
[ROW][C]PV;SC[/C][C]0.2972[/C][C]0.341[/C][C]0.2511[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]PV;BC[/C][C]0.7606[/C][C]0.8249[/C][C]0.6272[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]PV;CR[/C][C]0.6652[/C][C]0.7185[/C][C]0.5273[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]PV;TS[/C][C]0.743[/C][C]0.792[/C][C]0.6032[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]PV;TH[/C][C]0.7054[/C][C]0.7791[/C][C]0.5836[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]SC;BC[/C][C]0.238[/C][C]0.2759[/C][C]0.2074[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]SC;CR[/C][C]0.3895[/C][C]0.3511[/C][C]0.2616[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]SC;TS[/C][C]0.2537[/C][C]0.2917[/C][C]0.2192[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]SC;TH[/C][C]0.2632[/C][C]0.3014[/C][C]0.2228[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]BC;CR[/C][C]0.7623[/C][C]0.7722[/C][C]0.574[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]BC;TS[/C][C]0.8401[/C][C]0.8676[/C][C]0.6852[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]BC;TH[/C][C]0.884[/C][C]0.9064[/C][C]0.7471[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]CR;TS[/C][C]0.7321[/C][C]0.7607[/C][C]0.5536[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]CR;TH[/C][C]0.7206[/C][C]0.7547[/C][C]0.557[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]TS;TH[/C][C]0.8869[/C][C]0.9032[/C][C]0.7339[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=159546&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=159546&T=2

As an alternative you can also use a QR Code:  

The GUIDs for individual cells are displayed in the table below:

Correlations for all pairs of data series with p-values
pairPearson rSpearman rhoKendall tau
PV;SC0.29720.3410.2511
p-value(0)(0)(0)
PV;BC0.76060.82490.6272
p-value(0)(0)(0)
PV;CR0.66520.71850.5273
p-value(0)(0)(0)
PV;TS0.7430.7920.6032
p-value(0)(0)(0)
PV;TH0.70540.77910.5836
p-value(0)(0)(0)
SC;BC0.2380.27590.2074
p-value(0)(0)(0)
SC;CR0.38950.35110.2616
p-value(0)(0)(0)
SC;TS0.25370.29170.2192
p-value(0)(0)(0)
SC;TH0.26320.30140.2228
p-value(0)(0)(0)
BC;CR0.76230.77220.574
p-value(0)(0)(0)
BC;TS0.84010.86760.6852
p-value(0)(0)(0)
BC;TH0.8840.90640.7471
p-value(0)(0)(0)
CR;TS0.73210.76070.5536
p-value(0)(0)(0)
CR;TH0.72060.75470.557
p-value(0)(0)(0)
TS;TH0.88690.90320.7339
p-value(0)(0)(0)



Parameters (Session):
par1 = kendall ;
Parameters (R input):
par1 = kendall ;
R code (references can be found in the software module):
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')