Repository of Reproducible Computations

Free Statistics

of Irreproducible Research!

Author's title

Author

*The author of this computation has been verified*

R Software Module

rwasp_pairs.wasp

Title produced by software

Kendall tau Correlation Matrix

Date of computation

Wed, 21 Dec 2011 09:39:24 -0500

Cite this page as follows

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2011/Dec/21/t1324478399irvegwnylq7k6n5.htm/, Retrieved Tue, 07 May 2024 12:42:37 +0000

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=158757, Retrieved Tue, 07 May 2024 12:42:37 +0000

QR Codes:

Paste this QR Code to cite your computation.

Original text written by user:

IsPrivate?

No (this computation is public)

User-defined keywords

Estimated Impact

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)

-       [Kendall tau Correlation Matrix] [Paperdeel3] [2011-12-21 14:39:24] [20d72825ced8e4e20cafbe8bf70803f4] [Current]

Feedback Forum

Post a new message

Dataseries X:

Download CSV

Histogram

Boxplots

159261	1801	91	48	19	20465
189672	1717	59	53	20	33629
7215	192	18	0	0	1423
129098	2295	95	51	27	25629
230632	3450	136	76	31	54002
515038	6861	263	136	36	151036
180745	1795	56	62	23	33287
185559	1681	59	83	30	31172
154581	1897	44	55	30	28113
298001	2974	96	67	26	57803
121844	1946	75	50	24	49830
200907	2330	70	87	30	52143
101647	1839	100	46	22	21055
220269	3183	119	79	28	47007
170952	1486	61	56	18	28735
154647	1567	88	54	22	59147
142018	1756	57	81	33	78950
79030	1247	61	6	15	13497
167047	2779	87	74	34	46154
27997	726	24	13	18	53249
73019	1048	59	22	15	10726
241082	2805	100	99	30	83700
195820	1760	72	38	25	40400
142001	2266	54	59	34	33797
145433	1848	86	50	21	36205
183744	1665	32	50	21	30165
202357	2084	163	61	25	58534
201385	1448	94	90	31	44663
354924	2741	118	60	31	92556
192399	2112	44	52	20	40078
182286	1684	44	61	28	34711
181590	1616	45	60	22	31076
133801	2227	105	53	17	74608
233686	3088	123	76	25	58092
219428	2389	53	63	24	42009
0	1	1	0	0	0
223044	2099	63	54	28	36022
100129	1669	51	44	14	23333
145864	2137	49	42	35	53349
249965	2153	64	83	34	92596
242379	2390	71	105	22	49598
145794	1701	59	37	34	44093
103623	1049	33	25	23	84205
195891	2161	78	64	24	63369
117156	1276	50	55	26	60132
157787	1190	95	41	22	37403
81293	745	32	23	35	24460
243273	2374	103	77	24	46456
233155	2289	89	59	31	66616
160344	2639	59	68	26	41554
48188	658	28	12	22	22346
161922	1917	69	99	21	30874
307432	2557	74	78	27	68701
235223	2026	79	56	30	35728
195583	1911	59	67	33	29010
146061	1716	56	40	11	23110
208834	1852	67	53	26	38844
93764	981	24	26	26	27084
151985	1177	66	67	23	35139
195506	2849	97	36	38	57476
148922	1688	60	50	31	33277
142670	2162	81	51	20	31141
129561	1331	61	46	22	61281
122204	1307	38	57	26	25820
160930	1256	35	27	26	23284
99184	1294	41	38	33	35378
192811	2311	71	72	36	74990
138708	2897	65	93	25	29653
114408	1103	38	59	24	64622
31970	340	15	5	21	4157
225558	2791	112	53	19	29245
139220	1338	72	40	12	50008
113612	1441	68	72	30	52338
119537	1681	72	53	21	13310
162203	2650	67	81	34	92901
100098	1499	44	27	32	10956
174768	2302	60	94	28	34241
158459	2540	97	71	28	75043
80934	1000	30	20	21	21152
84971	1234	71	34	31	42249
80545	927	68	54	26	42005
287191	2176	64	49	29	41152
62974	957	28	26	23	14399
134091	1551	40	48	25	28263
75555	1014	46	35	22	17215
162154	1772	55	32	26	48140
227638	2630	229	55	33	62897
115367	1205	112	58	24	22883
115603	1392	63	44	24	41622
155537	1524	52	45	21	40715
153133	1829	41	49	28	65897
165618	2229	78	72	27	76542
151517	1233	57	39	25	37477
133686	1365	58	28	15	53216
61342	950	40	24	13	40911
245196	2319	117	52	36	57021
195576	1857	70	96	24	73116
19349	223	12	13	1	3895
225371	2390	105	38	24	46609
153213	1985	78	41	31	29351
59117	700	29	24	4	2325
91762	1062	24	54	21	31747
136769	1311	54	68	23	32665
114798	1157	61	28	23	19249
85338	823	40	36	12	15292
27676	596	22	2	16	5842
153535	1545	48	91	29	33994
122417	1130	37	29	26	13018
0	0	0	0	0	0
91529	1082	32	46	25	98177
107205	1135	67	25	21	37941
144664	1367	45	51	23	31032
146445	1506	63	60	21	32683
76656	870	60	36	21	34545
3616	78	5	0	0	0
0	0	0	0	0	0
183088	1130	44	40	23	27525
144677	1582	84	68	33	66856
159104	2034	98	28	30	28549
128944	970	39	41	23	38610
43410	778	19	7	1	2781
175774	1752	73	70	29	41211
95401	957	42	30	18	22698
134837	2098	55	69	33	41194
60493	731	40	3	12	32689
19764	285	12	10	2	5752
164062	1834	56	46	21	26757
132696	1148	33	34	28	22527
155367	1646	54	54	29	44810
11796	256	9	1	2	0
10674	98	9	0	0	0
142261	1404	57	39	18	100674
6836	41	3	0	1	0
162563	1824	63	48	21	57786
5118	42	3	5	0	0
40248	528	16	8	4	5444
0	0	0	0	0	0
122641	1073	47	38	25	28470
88837	1305	38	21	26	61849
7131	81	4	0	0	0
9056	261	14	0	4	2179
76611	934	24	15	17	8019
132697	1180	51	50	21	39644
100681	1148	20	17	22	23494

Summary of computational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	2 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=158757&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=158757&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=158757&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 Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	2 seconds
R Server	'Herman Ole Andreas Wold' @ wold.wessa.net

Correlations for all pairs of data series (method=pearson)
	x1	x2	x3	x4	x5
x1	1	0.891	0.768	0.788	0.679	0.69
x2	0.891	1	0.826	0.792	0.675	0.693
x3	0.768	0.826	1	0.641	0.557	0.617
x4	0.788	0.792	0.641	1	0.674	0.657
x5	0.679	0.675	0.557	0.674	1	0.617
	0.69	0.693	0.617	0.657	0.617	1

\begin{tabular}{lllllllll}
\hline
Correlations for all pairs of data series (method=pearson) \tabularnewline
  & x1 & x2 & x3 & x4 & x5 &  \tabularnewline
x1 & 1 & 0.891 & 0.768 & 0.788 & 0.679 & 0.69 \tabularnewline
x2 & 0.891 & 1 & 0.826 & 0.792 & 0.675 & 0.693 \tabularnewline
x3 & 0.768 & 0.826 & 1 & 0.641 & 0.557 & 0.617 \tabularnewline
x4 & 0.788 & 0.792 & 0.641 & 1 & 0.674 & 0.657 \tabularnewline
x5 & 0.679 & 0.675 & 0.557 & 0.674 & 1 & 0.617 \tabularnewline
 & 0.69 & 0.693 & 0.617 & 0.657 & 0.617 & 1 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=158757&T=1

[TABLE]
[ROW][C]Correlations for all pairs of data series (method=pearson)[/C][/ROW]
[ROW][C] [/C][C]x1[/C][C]x2[/C][C]x3[/C][C]x4[/C][C]x5[/C][C][/C][/ROW]
[ROW][C]x1[/C][C]1[/C][C]0.891[/C][C]0.768[/C][C]0.788[/C][C]0.679[/C][C]0.69[/C][/ROW]
[ROW][C]x2[/C][C]0.891[/C][C]1[/C][C]0.826[/C][C]0.792[/C][C]0.675[/C][C]0.693[/C][/ROW]
[ROW][C]x3[/C][C]0.768[/C][C]0.826[/C][C]1[/C][C]0.641[/C][C]0.557[/C][C]0.617[/C][/ROW]
[ROW][C]x4[/C][C]0.788[/C][C]0.792[/C][C]0.641[/C][C]1[/C][C]0.674[/C][C]0.657[/C][/ROW]
[ROW][C]x5[/C][C]0.679[/C][C]0.675[/C][C]0.557[/C][C]0.674[/C][C]1[/C][C]0.617[/C][/ROW]
[ROW][C][/C][C]0.69[/C][C]0.693[/C][C]0.617[/C][C]0.657[/C][C]0.617[/C][C]1[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=158757&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=158757&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=pearson)
	x1	x2	x3	x4	x5
x1	1	0.891	0.768	0.788	0.679	0.69
x2	0.891	1	0.826	0.792	0.675	0.693
x3	0.768	0.826	1	0.641	0.557	0.617
x4	0.788	0.792	0.641	1	0.674	0.657
x5	0.679	0.675	0.557	0.674	1	0.617
	0.69	0.693	0.617	0.657	0.617	1

Correlations for all pairs of data series with p-values
pair	Pearson r	Spearman rho	Kendall tau
x1;x2	0.8913	0.8662	0.7034
p-value	(0)	(0)	(0)
x1;x3	0.7678	0.7276	0.5602
p-value	(0)	(0)	(0)
x1;x4	0.7877	0.7646	0.5852
p-value	(0)	(0)	(0)
x1;x5	0.679	0.5873	0.4359
p-value	(0)	(0)	(0)
x1;	0.6903	0.6247	0.4708
p-value	(0)	(0)	(0)
x2;x3	0.8258	0.7969	0.6288
p-value	(0)	(0)	(0)
x2;x4	0.7922	0.7728	0.5978
p-value	(0)	(0)	(0)
x2;x5	0.6754	0.6257	0.4686
p-value	(0)	(0)	(0)
x2;	0.6929	0.6354	0.4742
p-value	(0)	(0)	(0)
x3;x4	0.6413	0.6652	0.4997
p-value	(0)	(0)	(0)
x3;x5	0.5573	0.4943	0.3643
p-value	(0)	(0)	(0)
x3;	0.6173	0.6042	0.4584
p-value	(0)	(0)	(0)
x4;x5	0.6744	0.5961	0.4437
p-value	(0)	(0)	(0)
x4;	0.6566	0.63	0.4723
p-value	(0)	(0)	(0)
x5;	0.6166	0.5727	0.4279
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
x1;x2 & 0.8913 & 0.8662 & 0.7034 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
x1;x3 & 0.7678 & 0.7276 & 0.5602 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
x1;x4 & 0.7877 & 0.7646 & 0.5852 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
x1;x5 & 0.679 & 0.5873 & 0.4359 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
x1; & 0.6903 & 0.6247 & 0.4708 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
x2;x3 & 0.8258 & 0.7969 & 0.6288 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
x2;x4 & 0.7922 & 0.7728 & 0.5978 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
x2;x5 & 0.6754 & 0.6257 & 0.4686 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
x2; & 0.6929 & 0.6354 & 0.4742 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
x3;x4 & 0.6413 & 0.6652 & 0.4997 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
x3;x5 & 0.5573 & 0.4943 & 0.3643 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
x3; & 0.6173 & 0.6042 & 0.4584 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
x4;x5 & 0.6744 & 0.5961 & 0.4437 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
x4; & 0.6566 & 0.63 & 0.4723 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
x5; & 0.6166 & 0.5727 & 0.4279 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=158757&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]x1;x2[/C][C]0.8913[/C][C]0.8662[/C][C]0.7034[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]x1;x3[/C][C]0.7678[/C][C]0.7276[/C][C]0.5602[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]x1;x4[/C][C]0.7877[/C][C]0.7646[/C][C]0.5852[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]x1;x5[/C][C]0.679[/C][C]0.5873[/C][C]0.4359[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]x1;[/C][C]0.6903[/C][C]0.6247[/C][C]0.4708[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]x2;x3[/C][C]0.8258[/C][C]0.7969[/C][C]0.6288[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]x2;x4[/C][C]0.7922[/C][C]0.7728[/C][C]0.5978[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]x2;x5[/C][C]0.6754[/C][C]0.6257[/C][C]0.4686[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]x2;[/C][C]0.6929[/C][C]0.6354[/C][C]0.4742[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]x3;x4[/C][C]0.6413[/C][C]0.6652[/C][C]0.4997[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]x3;x5[/C][C]0.5573[/C][C]0.4943[/C][C]0.3643[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]x3;[/C][C]0.6173[/C][C]0.6042[/C][C]0.4584[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]x4;x5[/C][C]0.6744[/C][C]0.5961[/C][C]0.4437[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]x4;[/C][C]0.6566[/C][C]0.63[/C][C]0.4723[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]x5;[/C][C]0.6166[/C][C]0.5727[/C][C]0.4279[/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=158757&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=158757&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
pair	Pearson r	Spearman rho	Kendall tau
x1;x2	0.8913	0.8662	0.7034
p-value	(0)	(0)	(0)
x1;x3	0.7678	0.7276	0.5602
p-value	(0)	(0)	(0)
x1;x4	0.7877	0.7646	0.5852
p-value	(0)	(0)	(0)
x1;x5	0.679	0.5873	0.4359
p-value	(0)	(0)	(0)
x1;	0.6903	0.6247	0.4708
p-value	(0)	(0)	(0)
x2;x3	0.8258	0.7969	0.6288
p-value	(0)	(0)	(0)
x2;x4	0.7922	0.7728	0.5978
p-value	(0)	(0)	(0)
x2;x5	0.6754	0.6257	0.4686
p-value	(0)	(0)	(0)
x2;	0.6929	0.6354	0.4742
p-value	(0)	(0)	(0)
x3;x4	0.6413	0.6652	0.4997
p-value	(0)	(0)	(0)
x3;x5	0.5573	0.4943	0.3643
p-value	(0)	(0)	(0)
x3;	0.6173	0.6042	0.4584
p-value	(0)	(0)	(0)
x4;x5	0.6744	0.5961	0.4437
p-value	(0)	(0)	(0)
x4;	0.6566	0.63	0.4723
p-value	(0)	(0)	(0)
x5;	0.6166	0.5727	0.4279
p-value	(0)	(0)	(0)

Figure 1

PNG link

Postscript link

PDF link

Parameters (Session):

par1 = pearson ;

Parameters (R input):

par1 = pearson ;

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

Free Statistics

Description of Statistical Computation

Tree of Dependent Computations

Dataset

Tables (Output of Computation)

Figures (Output of Computation)

Input Parameters & R Code