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, 12 Dec 2014 14:41: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/2014/Dec/12/t1418396085qlfv51nxdoyow7d.htm/, Retrieved Thu, 16 May 2024 20:03:40 +0000

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=266767, Retrieved Thu, 16 May 2024 20:03:40 +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] [] [2014-12-12 14:41:09] [37e054ac358b2aa7c2a1d0b751dfa890] [Current]

Feedback Forum

Post a new message

Dataseries X:

Download CSV

Histogram

Boxplots

23	24	13	27
22	18	23	30
26	17	24	24
41	21	22	16
23	17	18	27
33	19	17	18
31	19	22	24
35	24	27	24
28	17	15	18
31	20	19	22
23	24	23	25
25	18	11	16
30	22	19	18
30	17	14	24
19	22	20	24
32	26	21	29
50	23	23	22
27	20	19	21
36	20	17	23
31	18	24	24
26	25	25	23
32	25	16	19
35	19	21	24
30	26	28	20
38	24	25	24
41	24	22	30
27	20	16	17
28	21	22	22
24	18	14	24
21	14	10	20
39	14	18	23
33	26	18	19
28	16	17	22
47	23	19	24
26	21	23	20
25	11	18	24
34	21	17	26
30	16	19	24
30	21	21	24
25	24	21	24
19	14	10	21
28	18	18	22
39	29	22	29
20	12	9	23
30	19	24	22
31	21	21	25
19	20	15	23
25	19	17	24
52	28	29	30
33	21	26	24
22	21	18	24
32	16	20	20
17	13	11	16
31	22	19	27
20	30	14	13
29	18	24	29
37	22	14	27
21	20	17	24
23	19	15	24
30	19	22	23
21	21	18	22
24	19	23	26
40	23	27	26
20	17	13	21
33	19	21	23
20	15	13	20
26	20	20	28
22	24	25	29
32	17	13	16
13	10	10	25
28	23	23	28
32	22	18	24
27	21	15	24
32	18	25	24
23	19	20	12
28	23	17	22
23	19	23	22
29	24	21	24
26	20	16	26
15	17	17	24
14	20	14	26
19	22	21	22
19	15	13	23
26	23	17	29
33	18	13	16
35	20	25	18
28	20	16	22
25	19	12	23
41	24	22	30
28	13	17	24
25	19	25	21
26	16	16	23
41	22	24	14
28	21	17	25
26	16	15	17
24	19	16	24
32	14	14	23
25	19	20	22
22	20	15	16
29	19	18	22
36	23	24	30
40	20	18	25
27	13	15	21
35	17	23	22
18	22	21	23
36	20	20	24
27	19	24	22
31	21	23	21
16	16	15	26
26	24	21	24
20	23	21	27

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=266767&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=266767&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=266767&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=kendall)
	anderen	positief	negatief	organisatie
anderen	1	0.222	0.309	0.045
positief	0.222	1	0.341	0.201
negatief	0.309	0.341	1	0.176
organisatie	0.045	0.201	0.176	1

\begin{tabular}{lllllllll}
\hline
Correlations for all pairs of data series (method=kendall) \tabularnewline
  & anderen & positief & negatief & organisatie \tabularnewline
anderen & 1 & 0.222 & 0.309 & 0.045 \tabularnewline
positief & 0.222 & 1 & 0.341 & 0.201 \tabularnewline
negatief & 0.309 & 0.341 & 1 & 0.176 \tabularnewline
organisatie & 0.045 & 0.201 & 0.176 & 1 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=266767&T=1

[TABLE]
[ROW][C]Correlations for all pairs of data series (method=kendall)[/C][/ROW]
[ROW][C] [/C][C]anderen[/C][C]positief[/C][C]negatief[/C][C]organisatie[/C][/ROW]
[ROW][C]anderen[/C][C]1[/C][C]0.222[/C][C]0.309[/C][C]0.045[/C][/ROW]
[ROW][C]positief[/C][C]0.222[/C][C]1[/C][C]0.341[/C][C]0.201[/C][/ROW]
[ROW][C]negatief[/C][C]0.309[/C][C]0.341[/C][C]1[/C][C]0.176[/C][/ROW]
[ROW][C]organisatie[/C][C]0.045[/C][C]0.201[/C][C]0.176[/C][C]1[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=266767&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=266767&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)
	anderen	positief	negatief	organisatie
anderen	1	0.222	0.309	0.045
positief	0.222	1	0.341	0.201
negatief	0.309	0.341	1	0.176
organisatie	0.045	0.201	0.176	1

Correlations for all pairs of data series with p-values
pair	Pearson r	Spearman rho	Kendall tau
anderen;positief	0.3652	0.3018	0.2223
p-value	(1e-04)	(0.0013)	(0.0011)
anderen;negatief	0.4658	0.4225	0.309
p-value	(0)	(0)	(0)
anderen;organisatie	0.1066	0.0664	0.0454
p-value	(0.2656)	(0.4889)	(0.5072)
positief;negatief	0.5003	0.4544	0.3411
p-value	(0)	(0)	(0)
positief;organisatie	0.1984	0.2671	0.2011
p-value	(0.0369)	(0.0046)	(0.0039)
negatief;organisatie	0.2782	0.2383	0.1756
p-value	(0.0031)	(0.0118)	(0.0111)

\begin{tabular}{lllllllll}
\hline
Correlations for all pairs of data series with p-values \tabularnewline
pair & Pearson r & Spearman rho & Kendall tau \tabularnewline
anderen;positief & 0.3652 & 0.3018 & 0.2223 \tabularnewline
p-value & (1e-04) & (0.0013) & (0.0011) \tabularnewline
anderen;negatief & 0.4658 & 0.4225 & 0.309 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
anderen;organisatie & 0.1066 & 0.0664 & 0.0454 \tabularnewline
p-value & (0.2656) & (0.4889) & (0.5072) \tabularnewline
positief;negatief & 0.5003 & 0.4544 & 0.3411 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
positief;organisatie & 0.1984 & 0.2671 & 0.2011 \tabularnewline
p-value & (0.0369) & (0.0046) & (0.0039) \tabularnewline
negatief;organisatie & 0.2782 & 0.2383 & 0.1756 \tabularnewline
p-value & (0.0031) & (0.0118) & (0.0111) \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=266767&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]anderen;positief[/C][C]0.3652[/C][C]0.3018[/C][C]0.2223[/C][/ROW]
[ROW][C]p-value[/C][C](1e-04)[/C][C](0.0013)[/C][C](0.0011)[/C][/ROW]
[ROW][C]anderen;negatief[/C][C]0.4658[/C][C]0.4225[/C][C]0.309[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]anderen;organisatie[/C][C]0.1066[/C][C]0.0664[/C][C]0.0454[/C][/ROW]
[ROW][C]p-value[/C][C](0.2656)[/C][C](0.4889)[/C][C](0.5072)[/C][/ROW]
[ROW][C]positief;negatief[/C][C]0.5003[/C][C]0.4544[/C][C]0.3411[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]positief;organisatie[/C][C]0.1984[/C][C]0.2671[/C][C]0.2011[/C][/ROW]
[ROW][C]p-value[/C][C](0.0369)[/C][C](0.0046)[/C][C](0.0039)[/C][/ROW]
[ROW][C]negatief;organisatie[/C][C]0.2782[/C][C]0.2383[/C][C]0.1756[/C][/ROW]
[ROW][C]p-value[/C][C](0.0031)[/C][C](0.0118)[/C][C](0.0111)[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=266767&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=266767&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
anderen;positief	0.3652	0.3018	0.2223
p-value	(1e-04)	(0.0013)	(0.0011)
anderen;negatief	0.4658	0.4225	0.309
p-value	(0)	(0)	(0)
anderen;organisatie	0.1066	0.0664	0.0454
p-value	(0.2656)	(0.4889)	(0.5072)
positief;negatief	0.5003	0.4544	0.3411
p-value	(0)	(0)	(0)
positief;organisatie	0.1984	0.2671	0.2011
p-value	(0.0369)	(0.0046)	(0.0039)
negatief;organisatie	0.2782	0.2383	0.1756
p-value	(0.0031)	(0.0118)	(0.0111)

Meta Analysis of Correlation Tests
Number of significant by total number of Correlations
Type I error	Pearson r	Spearman rho	Kendall tau
0.01	0.67	0.67	0.67
0.02	0.67	0.83	0.83
0.03	0.67	0.83	0.83
0.04	0.83	0.83	0.83
0.05	0.83	0.83	0.83
0.06	0.83	0.83	0.83
0.07	0.83	0.83	0.83
0.08	0.83	0.83	0.83
0.09	0.83	0.83	0.83
0.1	0.83	0.83	0.83

\begin{tabular}{lllllllll}
\hline
Meta Analysis of Correlation Tests \tabularnewline
Number of significant by total number of Correlations \tabularnewline
Type I error & Pearson r & Spearman rho & Kendall tau \tabularnewline
0.01 & 0.67 & 0.67 & 0.67 \tabularnewline
0.02 & 0.67 & 0.83 & 0.83 \tabularnewline
0.03 & 0.67 & 0.83 & 0.83 \tabularnewline
0.04 & 0.83 & 0.83 & 0.83 \tabularnewline
0.05 & 0.83 & 0.83 & 0.83 \tabularnewline
0.06 & 0.83 & 0.83 & 0.83 \tabularnewline
0.07 & 0.83 & 0.83 & 0.83 \tabularnewline
0.08 & 0.83 & 0.83 & 0.83 \tabularnewline
0.09 & 0.83 & 0.83 & 0.83 \tabularnewline
0.1 & 0.83 & 0.83 & 0.83 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=266767&T=3

[TABLE]
[ROW][C]Meta Analysis of Correlation Tests[/C][/ROW]
[ROW][C]Number of significant by total number of Correlations[/C][/ROW]
[ROW][C]Type I error[/C][C]Pearson r[/C][C]Spearman rho[/C][C]Kendall tau[/C][/ROW]
[ROW][C]0.01[/C][C]0.67[/C][C]0.67[/C][C]0.67[/C][/ROW]
[ROW][C]0.02[/C][C]0.67[/C][C]0.83[/C][C]0.83[/C][/ROW]
[ROW][C]0.03[/C][C]0.67[/C][C]0.83[/C][C]0.83[/C][/ROW]
[ROW][C]0.04[/C][C]0.83[/C][C]0.83[/C][C]0.83[/C][/ROW]
[ROW][C]0.05[/C][C]0.83[/C][C]0.83[/C][C]0.83[/C][/ROW]
[ROW][C]0.06[/C][C]0.83[/C][C]0.83[/C][C]0.83[/C][/ROW]
[ROW][C]0.07[/C][C]0.83[/C][C]0.83[/C][C]0.83[/C][/ROW]
[ROW][C]0.08[/C][C]0.83[/C][C]0.83[/C][C]0.83[/C][/ROW]
[ROW][C]0.09[/C][C]0.83[/C][C]0.83[/C][C]0.83[/C][/ROW]
[ROW][C]0.1[/C][C]0.83[/C][C]0.83[/C][C]0.83[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=266767&T=3

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=266767&T=3

As an alternative you can also use a QR Code:

The GUIDs for individual cells are displayed in the table below:

Meta Analysis of Correlation Tests
Number of significant by total number of Correlations
Type I error	Pearson r	Spearman rho	Kendall tau
0.01	0.67	0.67	0.67
0.02	0.67	0.83	0.83
0.03	0.67	0.83	0.83
0.04	0.83	0.83	0.83
0.05	0.83	0.83	0.83
0.06	0.83	0.83	0.83
0.07	0.83	0.83	0.83
0.08	0.83	0.83	0.83
0.09	0.83	0.83	0.83
0.1	0.83	0.83	0.83

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')
ncorrs <- (n*n -n)/2
mycorrs <- array(0, dim=c(10,3))
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)
for (iii in 1:10) {
iiid100 <- iii / 100
if (rp$p.value < iiid100) mycorrs[iii, 1] = mycorrs[iii, 1] + 1
if (rs$p.value < iiid100) mycorrs[iii, 2] = mycorrs[iii, 2] + 1
if (rk$p.value < iiid100) mycorrs[iii, 3] = mycorrs[iii, 3] + 1
}
}
}
a<-table.end(a)
table.save(a,file='mytable1.tab')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Meta Analysis of Correlation Tests',4,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Number of significant by total number of Correlations',4,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Type I error',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)
for (iii in 1:10) {
iiid100 <- iii / 100
a<-table.row.start(a)
a<-table.element(a,round(iiid100,2),header=T)
a<-table.element(a,round(mycorrs[iii,1]/ncorrs,2))
a<-table.element(a,round(mycorrs[iii,2]/ncorrs,2))
a<-table.element(a,round(mycorrs[iii,3]/ncorrs,2))
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable2.tab')

Free Statistics

Description of Statistical Computation

Tree of Dependent Computations

Dataset

Tables (Output of Computation)

Figures (Output of Computation)

Input Parameters & R Code