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, 12 Dec 2010 10:18:09 +0000

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/2010/Dec/12/t1292149139l6608ddg6nmqh14.htm/, Retrieved Wed, 08 May 2024 00:10:46 +0000

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=108357, Retrieved Wed, 08 May 2024 00:10:46 +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

153

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

-       [Kendall tau Correlation Matrix] [Workshop 10(1)] [2010-12-12 10:18:09] [c9b1b69acb8f4b2b921fdfd5091a94b7] [Current]

Feedback Forum

Post a new message

Dataseries X:

Download CSV

Histogram

Boxplots

13	15	42	14	13
12	18	51	8	13
15	11	42	12	16
12	16	46	7	12
10	12	41	10	11
12	17	49	7	12
15	15	47	16	18
9	19	33	11	11
12	16	33	14	14
11	18	47	6	9
11	10	42	16	14
11	14	32	11	12
15	18	53	16	11
7	18	41	12	12
11	14	41	7	13
11	14	33	13	11
10	12	37	11	12
14	16	43	15	16
10	15	45	7	9
6	13	33	9	11
11	16	49	7	13
15	14	42	14	15
11	9	43	15	10
12	9	37	7	11
14	17	43	15	13
15	13	42	17	16
9	15	43	15	15
13	17	46	14	14
13	16	33	14	14
16	12	42	8	14
13	11	40	8	8
12	16	44	14	13
14	17	42	14	15
11	17	52	8	13
9	16	44	11	11
16	13	45	16	15
12	12	46	10	15
10	12	36	8	9
13	16	45	14	13
16	14	49	16	16
14	12	43	13	13
15	12	43	5	11
5	14	37	8	12
8	8	32	10	12
11	15	45	8	12
16	14	45	13	14
17	11	45	15	14
9	13	45	6	8
9	14	31	12	13
13	15	33	16	16
10	16	44	5	13
6	10	49	15	11
12	11	44	12	14
8	12	41	8	13
14	14	44	13	13
12	15	38	14	13
11	16	33	12	12
16	9	47	16	16
8	11	37	10	15
15	15	48	15	15
7	15	40	8	12
16	13	50	16	14
14	17	54	19	12
16	17	43	14	15
9	15	54	6	12
14	13	44	13	13
11	15	47	15	12
13	13	33	7	12
15	15	45	13	13
5	10	33	4	5
15	15	44	14	13
13	14	47	13	13
11	15	45	11	14
11	16	43	14	17
12	13	33	12	13
12	7	43	15	13
12	13	33	14	12
12	15	46	13	13
14	13	47	8	14
6	16	47	6	11
7	16	0	7	12
14	14	42	13	12
14	12	43	13	16
10	11	44	11	12
13	15	46	5	12
12	14	36	12	12
9	11	42	8	10
12	14	44	11	15
16	15	47	14	15
10	9	41	9	12
14	15	47	10	16
10	17	46	13	15
16	16	47	16	16
15	14	46	16	13
12	15	46	11	12
10	16	36	8	11
8	10	30	4	13
8	17	48	7	10
11	15	45	14	15
13	15	49	11	13
16	13	55	15	15
14	14	11	17	18
11	16	52	5	13
4	11	33	4	10
14	18	47	10	16
9	14	33	11	13
14	14	44	15	15
8	14	42	10	14
8	14	55	9	15
11	15	42	12	14
12	14	46	15	13
11	15	43	7	13
14	15	46	13	15
15	15	47	12	16
16	12	33	14	14
16	19	53	14	14
14	13	42	15	14
14	15	43	12	12
12	15	44	12	13
14	17	55	16	12
8	9	40	9	12
13	15	40	15	14
16	15	46	15	14
12	16	53	6	14
16	12	49	14	16
12	17	44	15	13
11	11	35	10	14
4	15	40	6	4
16	11	44	14	16
15	15	46	12	13
10	17	45	8	16
13	14	53	11	15
15	12	45	13	14
12	14	48	9	13
14	15	46	15	14
7	16	55	13	12
19	16	47	15	15
12	14	43	14	14
12	11	38	16	13
13	12	46	14	14
15	16	53	14	16
8	14	40	10	6
12	13	47	10	13
10	13	47	4	13
8	14	42	8	14
10	16	53	15	15
15	16	43	16	14
16	12	44	12	15
13	11	42	12	13
16	13	51	15	16
9	15	54	9	12
14	13	41	12	15
14	16	51	14	12
12	13	51	11	14

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	'RServer@AstonUniversity' @ vre.aston.ac.uk
R Framework error message	The field 'Names of X columns' contains a hard return which cannot be interpreted. Please, resubmit your request without hard returns in the 'Names of X columns'.

\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 & 'RServer@AstonUniversity' @ vre.aston.ac.uk \tabularnewline
R Framework error message & The field 'Names of X columns' contains a hard return which cannot be interpreted.
Please, resubmit your request without hard returns in the 'Names of X columns'. \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=108357&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]'RServer@AstonUniversity' @ vre.aston.ac.uk[/C][/ROW]
[ROW][C]R Framework error message[/C][C]The field 'Names of X columns' contains a hard return which cannot be interpreted.
Please, resubmit your request without hard returns in the 'Names of X columns'.[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=108357&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=108357&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	'RServer@AstonUniversity' @ vre.aston.ac.uk
R Framework error message	The field 'Names of X columns' contains a hard return which cannot be interpreted. Please, resubmit your request without hard returns in the 'Names of X columns'.

Correlations for all pairs of data series (method=kendall)
	Y1	X1	X2	X3	X4
Y1	1	0.008	0.229	0.444	0.437
X1	0.008	1	0.258	0.049	0.015
X2	0.229	0.258	1	0.104	0.166
X3	0.444	0.049	0.104	1	0.375
X4	0.437	0.015	0.166	0.375	1

\begin{tabular}{lllllllll}
\hline
Correlations for all pairs of data series (method=kendall) \tabularnewline
  & Y1 & X1 & X2 & X3 & X4
 \tabularnewline
Y1 & 1 & 0.008 & 0.229 & 0.444 & 0.437 \tabularnewline
X1 & 0.008 & 1 & 0.258 & 0.049 & 0.015 \tabularnewline
X2 & 0.229 & 0.258 & 1 & 0.104 & 0.166 \tabularnewline
X3 & 0.444 & 0.049 & 0.104 & 1 & 0.375 \tabularnewline
X4
 & 0.437 & 0.015 & 0.166 & 0.375 & 1 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=108357&T=1

[TABLE]
[ROW][C]Correlations for all pairs of data series (method=kendall)[/C][/ROW]
[ROW][C] [/C][C]Y1[/C][C]X1[/C][C]X2[/C][C]X3[/C][C]X4
[/C][/ROW]
[ROW][C]Y1[/C][C]1[/C][C]0.008[/C][C]0.229[/C][C]0.444[/C][C]0.437[/C][/ROW]
[ROW][C]X1[/C][C]0.008[/C][C]1[/C][C]0.258[/C][C]0.049[/C][C]0.015[/C][/ROW]
[ROW][C]X2[/C][C]0.229[/C][C]0.258[/C][C]1[/C][C]0.104[/C][C]0.166[/C][/ROW]
[ROW][C]X3[/C][C]0.444[/C][C]0.049[/C][C]0.104[/C][C]1[/C][C]0.375[/C][/ROW]
[ROW][C]X4
[/C][C]0.437[/C][C]0.015[/C][C]0.166[/C][C]0.375[/C][C]1[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=108357&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=108357&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)
	Y1	X1	X2	X3	X4
Y1	1	0.008	0.229	0.444	0.437
X1	0.008	1	0.258	0.049	0.015
X2	0.229	0.258	1	0.104	0.166
X3	0.444	0.049	0.104	1	0.375
X4	0.437	0.015	0.166	0.375	1

Correlations for all pairs of data series with p-values
pair	Pearson r	Spearman rho	Kendall tau
Y1;X1	0.0787	0.0177	0.008
p-value	(0.332)	(0.8276)	(0.8938)
Y1;X2	0.288	0.3073	0.2289
p-value	(3e-04)	(1e-04)	(1e-04)
Y1;X3	0.578	0.5841	0.444
p-value	(0)	(0)	(0)
Y1;X4	0.5686	0.5553	0.4365
p-value	(0)	(0)	(0)
X1;X2	0.2334	0.3467	0.2576
p-value	(0.0036)	(0)	(0)
X1;X3	0.0641	0.0578	0.0492
p-value	(0.4296)	(0.4765)	(0.41)
X1;X4	0.0855	0.0276	0.0153
p-value	(0.2915)	(0.7336)	(0.8026)
X2;X3	0.1194	0.1361	0.1037
p-value	(0.1404)	(0.0923)	(0.0746)
X2;X4	0.1498	0.2233	0.1664
p-value	(0.0637)	(0.0054)	(0.0053)
X3;X4	0.5042	0.4814	0.3748
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
Y1;X1 & 0.0787 & 0.0177 & 0.008 \tabularnewline
p-value & (0.332) & (0.8276) & (0.8938) \tabularnewline
Y1;X2 & 0.288 & 0.3073 & 0.2289 \tabularnewline
p-value & (3e-04) & (1e-04) & (1e-04) \tabularnewline
Y1;X3 & 0.578 & 0.5841 & 0.444 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
Y1;X4
 & 0.5686 & 0.5553 & 0.4365 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
X1;X2 & 0.2334 & 0.3467 & 0.2576 \tabularnewline
p-value & (0.0036) & (0) & (0) \tabularnewline
X1;X3 & 0.0641 & 0.0578 & 0.0492 \tabularnewline
p-value & (0.4296) & (0.4765) & (0.41) \tabularnewline
X1;X4
 & 0.0855 & 0.0276 & 0.0153 \tabularnewline
p-value & (0.2915) & (0.7336) & (0.8026) \tabularnewline
X2;X3 & 0.1194 & 0.1361 & 0.1037 \tabularnewline
p-value & (0.1404) & (0.0923) & (0.0746) \tabularnewline
X2;X4
 & 0.1498 & 0.2233 & 0.1664 \tabularnewline
p-value & (0.0637) & (0.0054) & (0.0053) \tabularnewline
X3;X4
 & 0.5042 & 0.4814 & 0.3748 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=108357&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]Y1;X1[/C][C]0.0787[/C][C]0.0177[/C][C]0.008[/C][/ROW]
[ROW][C]p-value[/C][C](0.332)[/C][C](0.8276)[/C][C](0.8938)[/C][/ROW]
[ROW][C]Y1;X2[/C][C]0.288[/C][C]0.3073[/C][C]0.2289[/C][/ROW]
[ROW][C]p-value[/C][C](3e-04)[/C][C](1e-04)[/C][C](1e-04)[/C][/ROW]
[ROW][C]Y1;X3[/C][C]0.578[/C][C]0.5841[/C][C]0.444[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]Y1;X4
[/C][C]0.5686[/C][C]0.5553[/C][C]0.4365[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]X1;X2[/C][C]0.2334[/C][C]0.3467[/C][C]0.2576[/C][/ROW]
[ROW][C]p-value[/C][C](0.0036)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]X1;X3[/C][C]0.0641[/C][C]0.0578[/C][C]0.0492[/C][/ROW]
[ROW][C]p-value[/C][C](0.4296)[/C][C](0.4765)[/C][C](0.41)[/C][/ROW]
[ROW][C]X1;X4
[/C][C]0.0855[/C][C]0.0276[/C][C]0.0153[/C][/ROW]
[ROW][C]p-value[/C][C](0.2915)[/C][C](0.7336)[/C][C](0.8026)[/C][/ROW]
[ROW][C]X2;X3[/C][C]0.1194[/C][C]0.1361[/C][C]0.1037[/C][/ROW]
[ROW][C]p-value[/C][C](0.1404)[/C][C](0.0923)[/C][C](0.0746)[/C][/ROW]
[ROW][C]X2;X4
[/C][C]0.1498[/C][C]0.2233[/C][C]0.1664[/C][/ROW]
[ROW][C]p-value[/C][C](0.0637)[/C][C](0.0054)[/C][C](0.0053)[/C][/ROW]
[ROW][C]X3;X4
[/C][C]0.5042[/C][C]0.4814[/C][C]0.3748[/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=108357&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=108357&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
Y1;X1	0.0787	0.0177	0.008
p-value	(0.332)	(0.8276)	(0.8938)
Y1;X2	0.288	0.3073	0.2289
p-value	(3e-04)	(1e-04)	(1e-04)
Y1;X3	0.578	0.5841	0.444
p-value	(0)	(0)	(0)
Y1;X4	0.5686	0.5553	0.4365
p-value	(0)	(0)	(0)
X1;X2	0.2334	0.3467	0.2576
p-value	(0.0036)	(0)	(0)
X1;X3	0.0641	0.0578	0.0492
p-value	(0.4296)	(0.4765)	(0.41)
X1;X4	0.0855	0.0276	0.0153
p-value	(0.2915)	(0.7336)	(0.8026)
X2;X3	0.1194	0.1361	0.1037
p-value	(0.1404)	(0.0923)	(0.0746)
X2;X4	0.1498	0.2233	0.1664
p-value	(0.0637)	(0.0054)	(0.0053)
X3;X4	0.5042	0.4814	0.3748
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