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

Fri, 16 Dec 2011 10:37:49 -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/16/t1324049914whyt3bia0td4oz5.htm/, Retrieved Sun, 05 May 2024 15:12:10 +0000

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

115

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

-     [Kendall tau Correlation Matrix] [] [2010-12-05 17:44:33] [b98453cac15ba1066b407e146608df68]
- RMPD  [Kendall tau Correlation Matrix] [] [2011-12-13 08:44:05] [8fcdd1f5b88bf5ac5d2a0b8a91219b89]
-   PD    [Kendall tau Correlation Matrix] [] [2011-12-16 14:27:21] [8fcdd1f5b88bf5ac5d2a0b8a91219b89]
-    D      [Kendall tau Correlation Matrix] [] [2011-12-16 14:43:50] [8fcdd1f5b88bf5ac5d2a0b8a91219b89]
-    D          [Kendall tau Correlation Matrix] [] [2011-12-16 15:37:49] [888ed98a09d01be7e0be9dfdea403736] [Current]
-    D            [Kendall tau Correlation Matrix] [] [2011-12-16 15:42:20] [8fcdd1f5b88bf5ac5d2a0b8a91219b89] 

Feedback Forum

Post a new message

Dataseries X:

Download CSV

Histogram

Boxplots

210907	94	144	145
120982	103	103	101
176508	93	98	98
179321	103	135	132
123185	51	61	60
52746	70	39	38
385534	91	150	144
33170	22	5	5
101645	38	28	28
149061	93	84	84
165446	60	80	79
237213	123	130	127
173326	148	82	78
133131	90	60	60
258873	124	131	131
180083	70	84	84
324799	168	140	133
230964	115	151	150
236785	71	91	91
135473	66	138	132
202925	134	150	136
215147	117	124	124
344297	108	119	118
153935	84	73	70
132943	156	110	107
174724	120	123	119
174415	114	90	89
225548	94	116	112
223632	120	113	108
124817	81	56	52
221698	110	115	112
210767	133	119	116
170266	122	129	123
260561	158	127	125
84853	109	27	27
294424	124	175	162
101011	39	35	32
215641	92	64	64
325107	126	96	92
7176	NA	0	0
167542	70	84	83
106408	37	41	41
96560	38	47	47
265769	120	126	120
269651	93	105	105
149112	95	80	79
175824	77	70	65
152871	90	73	70
111665	80	57	55
116408	31	40	39
362301	110	68	67
78800	66	21	21
183167	138	127	127
277965	133	154	152
150629	113	116	113
168809	100	102	99
24188	7	7	7
329267	140	148	141
65029	61	21	21
101097	41	35	35
218946	96	112	109
244052	164	137	133
341570	78	135	123
103597	49	26	26
233328	102	230	230
256462	124	181	166
206161	99	71	68
311473	129	147	147
235800	62	190	179
177939	73	64	61
207176	114	105	101
196553	99	107	108
174184	70	94	90
143246	104	116	114
187559	116	106	103
187681	91	143	142
119016	74	81	79
182192	138	89	88
73566	67	26	25
194979	151	84	83
167488	72	113	113
143756	120	120	118
275541	115	110	110
243199	105	134	129
182999	104	54	51
135649	108	96	93
152299	98	78	76
120221	69	51	49
346485	111	121	118
145790	99	38	38
193339	71	145	141
80953	27	59	58
122774	69	27	27
130585	107	91	91
112611	73	48	48
286468	107	68	63
241066	93	58	56
148446	129	150	144
204713	69	74	73
182079	118	181	168
140344	73	65	64
220516	119	97	97
243060	104	121	117
162765	107	99	100
182613	99	152	149
232138	90	188	187
265318	197	138	127
85574	36	40	37
310839	85	254	245
225060	139	87	87
232317	106	178	177
144966	50	51	49
43287	64	49	49
155754	31	73	73
164709	63	176	177
201940	92	94	94
235454	106	120	117
220801	63	66	60
99466	69	56	55
92661	41	39	39
133328	56	66	64
61361	25	27	26
125930	65	65	64
100750	93	58	58
224549	114	98	95
82316	38	25	25
102010	44	26	26
101523	87	77	76
243511	110	130	129
22938	0	11	11
41566	27	2	2
152474	83	101	101
61857	30	31	28
99923	80	36	36
132487	98	120	89
317394	82	195	193
21054	0	4	4
209641	60	89	84
22648	28	24	23
31414	9	39	39
46698	33	14	14
131698	59	78	78
91735	49	15	14
244749	115	106	101
184510	140	83	82
79863	49	24	24
128423	120	37	36
97839	66	77	75
38214	21	16	16
151101	124	56	55
272458	152	132	131
172494	139	144	131
108043	38	40	39
328107	144	153	144
250579	120	143	139
351067	160	220	211
158015	114	79	78
98866	39	50	50
85439	78	39	39
229242	119	95	90
351619	141	169	166
84207	101	12	12
120445	56	63	57
324598	133	134	133
131069	83	69	69
204271	116	119	119
165543	90	119	119
141722	36	75	65
116048	50	63	61
250047	61	55	49
299775	97	103	101
195838	98	197	196
173260	78	16	15
254488	117	140	136
104389	148	89	89
136084	41	40	40
199476	105	125	123
92499	55	21	21
224330	132	167	163
135781	44	32	29
74408	21	36	35
81240	50	13	13
14688	0	5	5
181633	73	96	96
271856	86	151	151
7199	0	6	6
46660	13	13	13
17547	4	3	3
133368	57	57	56
95227	48	23	23
152601	46	61	57
			
98146	48	21	14
79619	32	43	43
59194	68	20	20
139942	87	82	72
118612	43	90	87
72880	67	25	21
65475	46	60	56
99643	46	61	59
71965	56	85	82
77272	48	43	43
49289	44	25	25
135131	60	41	38
108446	65	26	25
89746	55	38	38
44296	38	12	12
77648	52	29	29
181528	60	49	47
134019	54	46	45
124064	86	41	40
92630	24	31	30
121848	52	41	41
52915	49	26	25
81872	61	23	23
58981	61	14	14
53515	81	16	16
60812	43	25	26
56375	40	21	21
65490	40	32	27
80949	56	9	9
76302	68	35	33
104011	79	42	42
98104	47	68	68
67989	57	32	32
30989	41	6	6
135458	29	68	67
73504	3	33	33
63123	60	84	77
61254	30	46	46
74914	79	30	30
31774	47	0	0
81437	40	36	36
87186	48	47	46
50090	36	20	18
65745	42	50	48
56653	49	30	29
158399	57	30	28
46455	12	34	34
73624	40	33	33
38395	43	34	34
91899	33	37	33
139526	77	83	80
52164	43	32	32
51567	45	30	30
70551	47	43	41
84856	43	41	41
102538	45	51	51
86678	50	19	18
85709	35	37	34
34662	7	33	31
150580	71	41	39
99611	67	54	54
19349	0	14	14
99373	62	25	24
86230	54	25	24
30837	4	8	8
31706	25	26	26
89806	40	20	19
62088	38	11	11
40151	19	14	14
27634	17	3	1
76990	67	40	39
37460	14	5	5
54157	30	38	37
49862	54	32	32
84337	35	41	38
64175	59	46	47
59382	24	47	47
119308	58	37	37
76702	42	51	51
103425	46	49	45
70344	61	21	21
43410	3	1	1
104838	52	44	42
62215	25	26	26
69304	40	21	21
53117	32	4	4
19764	4	10	10
86680	49	43	43
84105	63	34	34
77945	67	32	31
89113	32	20	19
91005	23	34	34
40248	7	6	6
64187	54	12	11
50857	37	24	24
56613	35	16	16
62792	51	72	72
72535	39	27	21

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	'AstonUniversity' @ aston.wessa.net
R Framework error message	Warning: there are blank lines in the 'Data X' field. Please, use NA for missing data - blank lines are simply deleted and are NOT treated as missing values.

\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 & 'AstonUniversity' @ aston.wessa.net \tabularnewline
R Framework error message & Warning: there are blank lines in the 'Data X' field.
Please, use NA for missing data - blank lines are simply
 deleted and are NOT treated as missing values. \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=156057&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]'AstonUniversity' @ aston.wessa.net[/C][/ROW]
[ROW][C]R Framework error message[/C][C]Warning: there are blank lines in the 'Data X' field.
Please, use NA for missing data - blank lines are simply
 deleted and are NOT treated as missing values.[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=156057&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=156057&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	'AstonUniversity' @ aston.wessa.net
R Framework error message	Warning: there are blank lines in the 'Data X' field. Please, use NA for missing data - blank lines are simply deleted and are NOT treated as missing values.

Correlations for all pairs of data series (method=pearson)
	a	b	c	d
a	1	0.767	0.817	0.814
b	0.767	1	0.731	0.728
c	0.817	0.731	1	0.998
d	0.814	0.728	0.998	1

\begin{tabular}{lllllllll}
\hline
Correlations for all pairs of data series (method=pearson) \tabularnewline
  & a & b & c & d \tabularnewline
a & 1 & 0.767 & 0.817 & 0.814 \tabularnewline
b & 0.767 & 1 & 0.731 & 0.728 \tabularnewline
c & 0.817 & 0.731 & 1 & 0.998 \tabularnewline
d & 0.814 & 0.728 & 0.998 & 1 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=156057&T=1

[TABLE]
[ROW][C]Correlations for all pairs of data series (method=pearson)[/C][/ROW]
[ROW][C] [/C][C]a[/C][C]b[/C][C]c[/C][C]d[/C][/ROW]
[ROW][C]a[/C][C]1[/C][C]0.767[/C][C]0.817[/C][C]0.814[/C][/ROW]
[ROW][C]b[/C][C]0.767[/C][C]1[/C][C]0.731[/C][C]0.728[/C][/ROW]
[ROW][C]c[/C][C]0.817[/C][C]0.731[/C][C]1[/C][C]0.998[/C][/ROW]
[ROW][C]d[/C][C]0.814[/C][C]0.728[/C][C]0.998[/C][C]1[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=156057&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=156057&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)
	a	b	c	d
a	1	0.767	0.817	0.814
b	0.767	1	0.731	0.728
c	0.817	0.731	1	0.998
d	0.814	0.728	0.998	1

Correlations for all pairs of data series with p-values
pair	Pearson r	Spearman rho	Kendall tau
a;b	0.7668	0.8034	0.602
p-value	(0)	(0)	(0)
a;c	0.8168	0.8591	0.6675
p-value	(0)	(0)	(0)
a;d	0.8145	0.8576	0.6652
p-value	(0)	(0)	(0)
b;c	0.7307	0.7677	0.5648
p-value	(0)	(0)	(0)
b;d	0.7284	0.7684	0.5648
p-value	(0)	(0)	(0)
c;d	0.9982	0.9986	0.9801
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.7668 & 0.8034 & 0.602 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
a;c & 0.8168 & 0.8591 & 0.6675 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
a;d & 0.8145 & 0.8576 & 0.6652 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
b;c & 0.7307 & 0.7677 & 0.5648 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
b;d & 0.7284 & 0.7684 & 0.5648 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
c;d & 0.9982 & 0.9986 & 0.9801 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=156057&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.7668[/C][C]0.8034[/C][C]0.602[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]a;c[/C][C]0.8168[/C][C]0.8591[/C][C]0.6675[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]a;d[/C][C]0.8145[/C][C]0.8576[/C][C]0.6652[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]b;c[/C][C]0.7307[/C][C]0.7677[/C][C]0.5648[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]b;d[/C][C]0.7284[/C][C]0.7684[/C][C]0.5648[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]c;d[/C][C]0.9982[/C][C]0.9986[/C][C]0.9801[/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=156057&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=156057&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.7668	0.8034	0.602
p-value	(0)	(0)	(0)
a;c	0.8168	0.8591	0.6675
p-value	(0)	(0)	(0)
a;d	0.8145	0.8576	0.6652
p-value	(0)	(0)	(0)
b;c	0.7307	0.7677	0.5648
p-value	(0)	(0)	(0)
b;d	0.7284	0.7684	0.5648
p-value	(0)	(0)	(0)
c;d	0.9982	0.9986	0.9801
p-value	(0)	(0)	(0)

Figure 1

PNG link

Postscript link

PDF link

Parameters (Session):

par1 = 2 ; par2 = Do not include Seasonal Dummies ; par3 = No Linear Trend ;

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