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

Sun, 18 Dec 2011 09:46:25 -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/18/t1324219778eczl6jd6pl0c6tp.htm/, Retrieved Sun, 05 May 2024 20:10:17 +0000

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=156931, Retrieved Sun, 05 May 2024 20:10:17 +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

181

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

-     [Kendall tau Correlation Matrix] [WS10] [2011-12-09 08:45:50] [09e53a95f5780167f20e6b4304200573]
-    D    [Kendall tau Correlation Matrix] [Deel 3] [2011-12-18 14:46:25] [fd45a0bdad9cb518cfd91b75ef7bafe6] [Current]

Feedback Forum

Post a new message

Dataseries X:

Download CSV

Histogram

Boxplots

449290	240	36	278	151036	26258
312831	111	28	84	92556	14881
275326	91	26	129	57803	16957
266451	68	27	89	68701	18717
244190	57	29	60	41152	8464
233497	61	31	60	92596	26641
229437	77	26	67	35728	9724
221399	129	30	124	54002	20323
218900	217	33	100	62897	22557
215480	104	36	166	57021	14612
214975	117	25	302	58092	4362
204812	108	19	135	29245	6372
200652	99	25	255	47007	13448
200609	48	24	84	42009	14349
191923	58	27	58	36022	10881
191463	90	27	94	83700	20517
189007	87	21	114	46456	13275
186704	60	26	58	38844	8030
184024	40	20	90	40078	16318
183873	95	24	104	46609	11252
182328	61	24	94	73116	11219
178613	67	36	86	74990	19365
178228	67	21	94	49598	14426
176034	68	25	127	40400	2570
174876	54	33	49	29010	9863
174181	149	25	90	58534	15657
173779	80	31	113	66616	21224
172071	64	30	113	52143	20488
168465	52	20	74	33629	10265
166226	86	38	59	57476	7768
164516	42	20	111	31076	10462
161553	74	24	59	63369	5989
156916	85	18	109	20465	6200
155934	77	28	55	44663	15329
154806	38	23	70	27525	10698
153278	56	21	48	57786	6837
152455	65	25	88	41211	9188
151621	41	28	76	34711	9556
150416	54	30	54	31172	8247
148950	56	25	68	37477	7028
148531	32	26	56	23284	5753
147766	92	28	119	28549	6852
146552	55	26	171	41554	10615
146020	71	31	61	29351	9289
144677	77	19	81	59147	6185
144193	60	34	89	92901	12964
143910	46	26	71	48140	8895
143546	80	21	41	36205	5299
141868	37	29	97	44810	9435
140935	49	17	50	28735	5389
140303	89	28	69	75043	9970
139901	85	22	40	37403	13958
139780	25	21	160	30165	7233
139075	69	23	87	76542	13178
136333	75	29	72	66856	14884
136220	48	21	72	40715	7820
134379	53	23	89	33287	10165
134153	52	33	92	78950	24369
131263	52	18	37	100674	7642
130414	52	29	65	33277	13823
129419	45	32	267	53349	13073
128907	65	25	152	29653	9422
126905	54	27	111	34241	11580
126602	53	34	62	44093	15604
125760	49	20	61	26757	9831
125728	42	29	127	33994	12840
124626	57	15	81	53216	6616
124401	41	23	55	31032	7987
122294	82	34	142	46154	17327
122259	91	26	154	25629	8874
121593	72	24	158	49830	8058
121391	41	30	87	28113	8683
120414	95	17	87	74608	7192
120111	47	21	54	39644	13400
119070	50	11	74	23110	8374
115944	44	23	84	32665	8208
115572	37	28	59	65897	12112
115343	56	19	84	61281	10170
113870	61	21	51	30874	12396
113245	37	23	56	38610	6873
111940	59	23	32	35139	14146
111924	49	29	60	41194	4260
109807	56	21	65	32683	10566
108685	26	25	66	22527	5953
106742	65	21	57	37941	5322
104972	61	12	80	50008	5413
102378	36	24	48	64622	12310
102194	33	24	46	25820	7573
102141	71	19	81	31141	7230
101560	65	21	89	13310	3394
100798	43	25	49	28470	8191
98629	49	34	164	33797	5162
97181	32	25	132	28263	10226
95704	42	22	46	60132	11270
95276	62	28	60	52338	10837
94982	37	29	36	35378	9791
94341	53	23	66	19249	7488
94332	29	23	35	84205	9184
94127	18	22	50	23494	8181
93107	49	14	137	23333	8022
90385	54	24	49	41622	8620
89920	31	10	32	13018	1350
88874	39	17	46	22698	5141
86249	94	21	83	21055	7945
85610	31	25	58	98177	10623
84971	71	31	105	42249	10138
84315	36	26	44	61849	1661
82209	95	22	71	22883	6900
81293	31	35	49	24460	7162
80420	64	26	49	42005	14697
72904	37	12	27	15292	4574
72559	22	26	72	27084	7509
71921	37	32	79	10956	3495
71894	57	21	71	34545	1900
64149	52	15	72	13497	70
63952	23	17	48	8019	3080
60373	39	12	29	32689	443
60138	23	21	76	21152	4911
58280	20	20	44	31747	6562
56252	26	23	49	14399	4204
52197	52	15	40	10726	3149
51201	27	4	11	2325	593
50913	42	22	49	17215	3456
49176	34	13	30	40911	9570
48188	28	22	37	22346	2102
43410	19	1	63	2781	4
40248	16	4	34	5444	775
31970	15	21	40	4157	1283
27918	24	18	49	53249	3878
27676	22	16	59	5842	522
19764	12	2	19	5752	2416
19349	12	1	15	3895	786
11796	9	2	22	0	0
10674	9	0	7	0	0
8812	13	4	18	2179	548
7215	18	0	1	1423	603
7131	4	0	0	0	0
6836	3	1	5	0	0
5118	3	0	1	0	0
3616	5	0	5	0	0
0	1	0	0	0	0
0	0	0	0	0	0
0	0	0	0	0	0
0	0	0	0	0	0

Summary of computational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	1 seconds
R Server	'Gwilym Jenkins' @ jenkins.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 & 1 seconds \tabularnewline
R Server & 'Gwilym Jenkins' @ jenkins.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=156931&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]1 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Gwilym Jenkins' @ jenkins.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=156931&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=156931&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	1 seconds
R Server	'Gwilym Jenkins' @ jenkins.wessa.net

Correlations for all pairs of data series (method=pearson)
	Y	X1	X2	X3	X4	X5
Y	1	0.772	0.665	0.614	0.704	0.737
X1	0.772	1	0.548	0.603	0.613	0.621
X2	0.665	0.548	1	0.526	0.618	0.69
X3	0.614	0.603	0.526	1	0.439	0.458
X4	0.704	0.613	0.618	0.439	1	0.722
X5	0.737	0.621	0.69	0.458	0.722	1

\begin{tabular}{lllllllll}
\hline
Correlations for all pairs of data series (method=pearson) \tabularnewline
  & Y & X1 & X2 & X3 & X4 & X5 \tabularnewline
Y & 1 & 0.772 & 0.665 & 0.614 & 0.704 & 0.737 \tabularnewline
X1 & 0.772 & 1 & 0.548 & 0.603 & 0.613 & 0.621 \tabularnewline
X2 & 0.665 & 0.548 & 1 & 0.526 & 0.618 & 0.69 \tabularnewline
X3 & 0.614 & 0.603 & 0.526 & 1 & 0.439 & 0.458 \tabularnewline
X4 & 0.704 & 0.613 & 0.618 & 0.439 & 1 & 0.722 \tabularnewline
X5 & 0.737 & 0.621 & 0.69 & 0.458 & 0.722 & 1 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=156931&T=1

[TABLE]
[ROW][C]Correlations for all pairs of data series (method=pearson)[/C][/ROW]
[ROW][C] [/C][C]Y[/C][C]X1[/C][C]X2[/C][C]X3[/C][C]X4[/C][C]X5[/C][/ROW]
[ROW][C]Y[/C][C]1[/C][C]0.772[/C][C]0.665[/C][C]0.614[/C][C]0.704[/C][C]0.737[/C][/ROW]
[ROW][C]X1[/C][C]0.772[/C][C]1[/C][C]0.548[/C][C]0.603[/C][C]0.613[/C][C]0.621[/C][/ROW]
[ROW][C]X2[/C][C]0.665[/C][C]0.548[/C][C]1[/C][C]0.526[/C][C]0.618[/C][C]0.69[/C][/ROW]
[ROW][C]X3[/C][C]0.614[/C][C]0.603[/C][C]0.526[/C][C]1[/C][C]0.439[/C][C]0.458[/C][/ROW]
[ROW][C]X4[/C][C]0.704[/C][C]0.613[/C][C]0.618[/C][C]0.439[/C][C]1[/C][C]0.722[/C][/ROW]
[ROW][C]X5[/C][C]0.737[/C][C]0.621[/C][C]0.69[/C][C]0.458[/C][C]0.722[/C][C]1[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=156931&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=156931&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)
	Y	X1	X2	X3	X4	X5
Y	1	0.772	0.665	0.614	0.704	0.737
X1	0.772	1	0.548	0.603	0.613	0.621
X2	0.665	0.548	1	0.526	0.618	0.69
X3	0.614	0.603	0.526	1	0.439	0.458
X4	0.704	0.613	0.618	0.439	1	0.722
X5	0.737	0.621	0.69	0.458	0.722	1

Correlations for all pairs of data series with p-values
pair	Pearson r	Spearman rho	Kendall tau
Y;X1	0.7716	0.7354	0.5684
p-value	(0)	(0)	(0)
Y;X2	0.6652	0.584	0.4282
p-value	(0)	(0)	(0)
Y;X3	0.6142	0.645	0.4794
p-value	(0)	(0)	(0)
Y;X4	0.7036	0.6521	0.4906
p-value	(0)	(0)	(0)
Y;X5	0.7372	0.6973	0.5272
p-value	(0)	(0)	(0)
X1;X2	0.5482	0.4828	0.3525
p-value	(0)	(0)	(0)
X1;X3	0.603	0.64	0.4778
p-value	(0)	(0)	(0)
X1;X4	0.6129	0.6075	0.459
p-value	(0)	(0)	(0)
X1;X5	0.6211	0.5786	0.4248
p-value	(0)	(0)	(0)
X2;X3	0.5264	0.5083	0.3768
p-value	(0)	(0)	(0)
X2;X4	0.6177	0.5678	0.4237
p-value	(0)	(0)	(0)
X2;X5	0.69	0.6584	0.5019
p-value	(0)	(0)	(0)
X3;X4	0.4387	0.459	0.3269
p-value	(0)	(0)	(0)
X3;X5	0.4575	0.5421	0.3961
p-value	(0)	(0)	(0)
X4;X5	0.7216	0.7084	0.5493
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
Y;X1 & 0.7716 & 0.7354 & 0.5684 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
Y;X2 & 0.6652 & 0.584 & 0.4282 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
Y;X3 & 0.6142 & 0.645 & 0.4794 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
Y;X4 & 0.7036 & 0.6521 & 0.4906 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
Y;X5 & 0.7372 & 0.6973 & 0.5272 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
X1;X2 & 0.5482 & 0.4828 & 0.3525 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
X1;X3 & 0.603 & 0.64 & 0.4778 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
X1;X4 & 0.6129 & 0.6075 & 0.459 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
X1;X5 & 0.6211 & 0.5786 & 0.4248 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
X2;X3 & 0.5264 & 0.5083 & 0.3768 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
X2;X4 & 0.6177 & 0.5678 & 0.4237 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
X2;X5 & 0.69 & 0.6584 & 0.5019 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
X3;X4 & 0.4387 & 0.459 & 0.3269 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
X3;X5 & 0.4575 & 0.5421 & 0.3961 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
X4;X5 & 0.7216 & 0.7084 & 0.5493 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=156931&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]Y;X1[/C][C]0.7716[/C][C]0.7354[/C][C]0.5684[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]Y;X2[/C][C]0.6652[/C][C]0.584[/C][C]0.4282[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]Y;X3[/C][C]0.6142[/C][C]0.645[/C][C]0.4794[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]Y;X4[/C][C]0.7036[/C][C]0.6521[/C][C]0.4906[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]Y;X5[/C][C]0.7372[/C][C]0.6973[/C][C]0.5272[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]X1;X2[/C][C]0.5482[/C][C]0.4828[/C][C]0.3525[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]X1;X3[/C][C]0.603[/C][C]0.64[/C][C]0.4778[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]X1;X4[/C][C]0.6129[/C][C]0.6075[/C][C]0.459[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]X1;X5[/C][C]0.6211[/C][C]0.5786[/C][C]0.4248[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]X2;X3[/C][C]0.5264[/C][C]0.5083[/C][C]0.3768[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]X2;X4[/C][C]0.6177[/C][C]0.5678[/C][C]0.4237[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]X2;X5[/C][C]0.69[/C][C]0.6584[/C][C]0.5019[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]X3;X4[/C][C]0.4387[/C][C]0.459[/C][C]0.3269[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]X3;X5[/C][C]0.4575[/C][C]0.5421[/C][C]0.3961[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]X4;X5[/C][C]0.7216[/C][C]0.7084[/C][C]0.5493[/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=156931&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=156931&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
Y;X1	0.7716	0.7354	0.5684
p-value	(0)	(0)	(0)
Y;X2	0.6652	0.584	0.4282
p-value	(0)	(0)	(0)
Y;X3	0.6142	0.645	0.4794
p-value	(0)	(0)	(0)
Y;X4	0.7036	0.6521	0.4906
p-value	(0)	(0)	(0)
Y;X5	0.7372	0.6973	0.5272
p-value	(0)	(0)	(0)
X1;X2	0.5482	0.4828	0.3525
p-value	(0)	(0)	(0)
X1;X3	0.603	0.64	0.4778
p-value	(0)	(0)	(0)
X1;X4	0.6129	0.6075	0.459
p-value	(0)	(0)	(0)
X1;X5	0.6211	0.5786	0.4248
p-value	(0)	(0)	(0)
X2;X3	0.5264	0.5083	0.3768
p-value	(0)	(0)	(0)
X2;X4	0.6177	0.5678	0.4237
p-value	(0)	(0)	(0)
X2;X5	0.69	0.6584	0.5019
p-value	(0)	(0)	(0)
X3;X4	0.4387	0.459	0.3269
p-value	(0)	(0)	(0)
X3;X5	0.4575	0.5421	0.3961
p-value	(0)	(0)	(0)
X4;X5	0.7216	0.7084	0.5493
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