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

Tue, 20 Dec 2011 17:08:57 -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/20/t1324418952y9bh5n0erhumk8h.htm/, Retrieved Mon, 06 May 2024 06:03:28 +0000

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=158310, Retrieved Mon, 06 May 2024 06:03:28 +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

103

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

-       [Kendall tau Correlation Matrix] [kendall] [2011-12-20 22:08:57] [935c692b8d0e827208dbfd6a4efb0528] [Current]

Feedback Forum

Post a new message

Dataseries X:

Download CSV

Histogram

Boxplots

158258	48	18	63	20465	23975
186930	53	20	56	33629	85634
7215	0	0	0	1423	1929
128162	51	27	63	25629	36294
226974	76	31	116	54002	72255
500344	125	36	138	151036	189748
171007	59	23	71	33287	61834
179835	80	30	107	31172	68167
154581	55	30	50	28113	38462
278960	67	26	79	57803	101219
121844	50	24	58	49830	43270
183086	77	30	91	52143	76183
98796	44	22	41	21055	31476
209322	79	25	91	47007	62157
157125	51	18	61	28735	46261
154565	54	22	74	59147	50063
134198	75	33	131	78950	64483
69128	2	15	45	13497	2341
150680	73	34	110	46154	48149
27997	13	18	41	53249	12743
69919	19	15	37	10726	18743
233044	93	30	84	83700	97057
195820	38	25	67	40400	17675
127994	48	34	69	33797	33106
145433	50	21	58	36205	53311
170864	48	21	60	30165	42754
199655	60	25	88	58534	59056
188633	81	31	75	44663	101621
354266	60	31	98	92556	118120
192399	52	20	67	40078	79572
165753	50	28	84	34711	42744
173721	60	20	58	31076	65931
126739	53	17	35	74608	38575
224762	76	25	74	58092	28795
219428	63	24	89	42009	94440
0	0	0	0	0	0
217267	54	27	75	36022	38229
99706	44	14	39	23333	31972
136733	36	35	101	53349	40071
249965	83	34	135	92596	132480
232951	105	22	76	49598	62797
143755	37	34	118	44093	40429
95734	25	23	76	84205	45545
191416	63	24	65	63369	57568
114820	55	26	97	60132	39019
157625	41	22	67	37403	53866
81293	23	35	63	24460	38345
210040	63	24	96	46456	50210
223771	54	31	112	66616	80947
160344	68	26	75	41554	43461
48188	12	22	39	22346	14812
145235	84	21	63	30874	37819
287839	66	27	93	68701	102738
235223	56	30	76	35728	54509
195583	67	33	117	29010	62956
145942	40	11	30	23110	55411
207309	53	26	65	38844	50611
93764	26	26	78	27084	26692
151985	67	23	87	35139	60056
190545	36	38	85	57476	25155
146414	50	30	111	33277	42840
130794	48	19	60	31141	39358
124234	46	19	53	61281	47241
112718	53	26	67	25820	49611
160817	27	26	90	23284	41833
99070	38	33	100	35378	48930
178653	68	36	135	74990	110600
138708	93	25	71	29653	52235
114408	59	24	75	64622	53986
31970	5	21	42	4157	4105
224494	53	19	42	29245	59331
123328	36	12	8	50008	47796
113504	72	30	86	52338	38302
105932	49	21	41	13310	14063
162203	81	34	118	92901	54414
100098	27	32	91	10956	9903
174768	94	28	102	34241	53987
156752	71	28	89	75043	88937
77269	18	21	46	21152	21928
84971	34	31	60	42249	29487
80522	54	26	69	42005	35334
276525	44	29	95	41152	57596
62974	26	23	17	14399	29750
120296	44	25	61	28263	41029
75555	35	22	55	17215	12416
157988	32	26	55	48140	51158
223247	55	33	124	62897	79935
115019	58	24	73	22883	26552
99602	44	24	73	41622	25807
151804	39	21	67	40715	50620
146005	49	28	66	65897	61467
163444	72	27	75	76542	65292
151517	39	25	83	37477	55516
133686	28	15	55	53216	42006
58128	24	13	27	40911	26273
234325	49	36	115	57021	90248
195576	96	24	76	73116	61476
19349	13	1	0	3895	9604
213189	32	24	83	46609	45108
151672	41	31	90	29351	47232
59117	24	4	4	2325	3439
71931	52	20	56	31747	30553
126653	57	23	63	32665	24751
113552	28	23	52	19249	34458
85338	36	12	24	15292	24649
27676	2	16	17	5842	2342
138522	80	29	105	33994	52739
122417	29	26	20	13018	6245
0	0	0	0	0	0
87592	46	25	51	98177	35381
107205	25	21	76	37941	19595
144664	51	23	59	31032	50848
136540	59	21	70	32683	39443
71894	36	21	38	34545	27023
3616	0	0	0	0	0
0	0	0	0	0	0
175055	38	23	81	27525	61022
144618	68	33	78	66856	63528
152826	28	28	67	28549	34835
113245	36	23	89	38610	37172
43410	7	1	3	2781	13
175762	70	29	87	41211	62548
93634	30	17	48	22698	31334
117426	59	31	66	41194	20839
60493	3	12	32	32689	5084
19764	10	2	4	5752	9927
164062	46	21	70	26757	53229
128144	34	26	94	22527	29877
154959	54	29	91	44810	37310
11796	1	2	1	0	0
10674	0	0	0	0	0
138547	35	18	39	100674	50067
6836	0	1	0	0	0
154135	48	21	45	57786	47708
5118	5	0	0	0	0
40248	8	4	7	5444	6012
0	0	0	0	0	0
120460	36	25	75	28470	27749
88837	21	26	52	61849	47555
7131	0	0	0	0	0
9056	0	4	1	2179	1336
68916	15	17	49	8019	11017
132697	50	21	69	39644	55184
100681	17	22	56	23494	43485

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	'AstonUniversity' @ aston.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 & 'AstonUniversity' @ aston.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=158310&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]'AstonUniversity' @ aston.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=158310&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=158310&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	'AstonUniversity' @ aston.wessa.net

Correlations for all pairs of data series (method=kendall)
	A	B	C	D	E	F
A	1	0.577	0.437	0.539	0.467	0.664
B	0.577	1	0.44	0.52	0.476	0.587
C	0.437	0.44	1	0.681	0.429	0.418
D	0.539	0.52	0.681	1	0.468	0.513
E	0.467	0.476	0.429	0.468	1	0.539
F	0.664	0.587	0.418	0.513	0.539	1

\begin{tabular}{lllllllll}
\hline
Correlations for all pairs of data series (method=kendall) \tabularnewline
  & A & B & C & D & E & F \tabularnewline
A & 1 & 0.577 & 0.437 & 0.539 & 0.467 & 0.664 \tabularnewline
B & 0.577 & 1 & 0.44 & 0.52 & 0.476 & 0.587 \tabularnewline
C & 0.437 & 0.44 & 1 & 0.681 & 0.429 & 0.418 \tabularnewline
D & 0.539 & 0.52 & 0.681 & 1 & 0.468 & 0.513 \tabularnewline
E & 0.467 & 0.476 & 0.429 & 0.468 & 1 & 0.539 \tabularnewline
F & 0.664 & 0.587 & 0.418 & 0.513 & 0.539 & 1 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=158310&T=1

[TABLE]
[ROW][C]Correlations for all pairs of data series (method=kendall)[/C][/ROW]
[ROW][C] [/C][C]A[/C][C]B[/C][C]C[/C][C]D[/C][C]E[/C][C]F[/C][/ROW]
[ROW][C]A[/C][C]1[/C][C]0.577[/C][C]0.437[/C][C]0.539[/C][C]0.467[/C][C]0.664[/C][/ROW]
[ROW][C]B[/C][C]0.577[/C][C]1[/C][C]0.44[/C][C]0.52[/C][C]0.476[/C][C]0.587[/C][/ROW]
[ROW][C]C[/C][C]0.437[/C][C]0.44[/C][C]1[/C][C]0.681[/C][C]0.429[/C][C]0.418[/C][/ROW]
[ROW][C]D[/C][C]0.539[/C][C]0.52[/C][C]0.681[/C][C]1[/C][C]0.468[/C][C]0.513[/C][/ROW]
[ROW][C]E[/C][C]0.467[/C][C]0.476[/C][C]0.429[/C][C]0.468[/C][C]1[/C][C]0.539[/C][/ROW]
[ROW][C]F[/C][C]0.664[/C][C]0.587[/C][C]0.418[/C][C]0.513[/C][C]0.539[/C][C]1[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=158310&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=158310&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)
	A	B	C	D	E	F
A	1	0.577	0.437	0.539	0.467	0.664
B	0.577	1	0.44	0.52	0.476	0.587
C	0.437	0.44	1	0.681	0.429	0.418
D	0.539	0.52	0.681	1	0.468	0.513
E	0.467	0.476	0.429	0.468	1	0.539
F	0.664	0.587	0.418	0.513	0.539	1

Correlations for all pairs of data series with p-values
pair	Pearson r	Spearman rho	Kendall tau
A;B	0.7763	0.7539	0.5772
p-value	(0)	(0)	(0)
A;C	0.6777	0.5866	0.4366
p-value	(0)	(0)	(0)
A;D	0.7419	0.7119	0.5391
p-value	(0)	(0)	(0)
A;E	0.691	0.6206	0.4675
p-value	(0)	(0)	(0)
A;F	0.8637	0.8286	0.6636
p-value	(0)	(0)	(0)
B;C	0.6753	0.59	0.4396
p-value	(0)	(0)	(0)
B;D	0.7305	0.6863	0.5201
p-value	(0)	(0)	(0)
B;E	0.6613	0.6363	0.4759
p-value	(0)	(0)	(0)
B;F	0.7645	0.7663	0.5874
p-value	(0)	(0)	(0)
C;D	0.8823	0.8318	0.6811
p-value	(0)	(0)	(0)
C;E	0.6193	0.5728	0.4286
p-value	(0)	(0)	(0)
C;F	0.6227	0.5592	0.4177
p-value	(0)	(0)	(0)
D;E	0.6461	0.6175	0.4676
p-value	(0)	(0)	(0)
D;F	0.7159	0.6849	0.513
p-value	(0)	(0)	(0)
E;F	0.7563	0.7017	0.5394
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
A;B & 0.7763 & 0.7539 & 0.5772 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
A;C & 0.6777 & 0.5866 & 0.4366 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
A;D & 0.7419 & 0.7119 & 0.5391 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
A;E & 0.691 & 0.6206 & 0.4675 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
A;F & 0.8637 & 0.8286 & 0.6636 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
B;C & 0.6753 & 0.59 & 0.4396 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
B;D & 0.7305 & 0.6863 & 0.5201 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
B;E & 0.6613 & 0.6363 & 0.4759 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
B;F & 0.7645 & 0.7663 & 0.5874 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
C;D & 0.8823 & 0.8318 & 0.6811 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
C;E & 0.6193 & 0.5728 & 0.4286 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
C;F & 0.6227 & 0.5592 & 0.4177 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
D;E & 0.6461 & 0.6175 & 0.4676 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
D;F & 0.7159 & 0.6849 & 0.513 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
E;F & 0.7563 & 0.7017 & 0.5394 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=158310&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]A;B[/C][C]0.7763[/C][C]0.7539[/C][C]0.5772[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]A;C[/C][C]0.6777[/C][C]0.5866[/C][C]0.4366[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]A;D[/C][C]0.7419[/C][C]0.7119[/C][C]0.5391[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]A;E[/C][C]0.691[/C][C]0.6206[/C][C]0.4675[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]A;F[/C][C]0.8637[/C][C]0.8286[/C][C]0.6636[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]B;C[/C][C]0.6753[/C][C]0.59[/C][C]0.4396[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]B;D[/C][C]0.7305[/C][C]0.6863[/C][C]0.5201[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]B;E[/C][C]0.6613[/C][C]0.6363[/C][C]0.4759[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]B;F[/C][C]0.7645[/C][C]0.7663[/C][C]0.5874[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]C;D[/C][C]0.8823[/C][C]0.8318[/C][C]0.6811[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]C;E[/C][C]0.6193[/C][C]0.5728[/C][C]0.4286[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]C;F[/C][C]0.6227[/C][C]0.5592[/C][C]0.4177[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]D;E[/C][C]0.6461[/C][C]0.6175[/C][C]0.4676[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]D;F[/C][C]0.7159[/C][C]0.6849[/C][C]0.513[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]E;F[/C][C]0.7563[/C][C]0.7017[/C][C]0.5394[/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=158310&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=158310&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
A;B	0.7763	0.7539	0.5772
p-value	(0)	(0)	(0)
A;C	0.6777	0.5866	0.4366
p-value	(0)	(0)	(0)
A;D	0.7419	0.7119	0.5391
p-value	(0)	(0)	(0)
A;E	0.691	0.6206	0.4675
p-value	(0)	(0)	(0)
A;F	0.8637	0.8286	0.6636
p-value	(0)	(0)	(0)
B;C	0.6753	0.59	0.4396
p-value	(0)	(0)	(0)
B;D	0.7305	0.6863	0.5201
p-value	(0)	(0)	(0)
B;E	0.6613	0.6363	0.4759
p-value	(0)	(0)	(0)
B;F	0.7645	0.7663	0.5874
p-value	(0)	(0)	(0)
C;D	0.8823	0.8318	0.6811
p-value	(0)	(0)	(0)
C;E	0.6193	0.5728	0.4286
p-value	(0)	(0)	(0)
C;F	0.6227	0.5592	0.4177
p-value	(0)	(0)	(0)
D;E	0.6461	0.6175	0.4676
p-value	(0)	(0)	(0)
D;F	0.7159	0.6849	0.513
p-value	(0)	(0)	(0)
E;F	0.7563	0.7017	0.5394
p-value	(0)	(0)	(0)

Figure 1

PNG link

Postscript link

PDF link

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

Free Statistics

Description of Statistical Computation

Tree of Dependent Computations

Dataset

Tables (Output of Computation)

Figures (Output of Computation)

Input Parameters & R Code