Free Statistics

of Irreproducible Research!

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_pairs.wasp
Title produced by softwareKendall tau Correlation Matrix
Date of computationWed, 15 Dec 2010 16:45:06 +0000
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2010/Dec/15/t1292431584m09ao0m5dflql66.htm/, Retrieved Fri, 03 May 2024 11:43:27 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=110558, Retrieved Fri, 03 May 2024 11:43:27 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact114
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Kendall tau Correlation Matrix] [Paper Kendall Cor...] [2010-12-15 16:45:06] [61e5ee05de011f44efa37f086a4e2271] [Current]
-    D    [Kendall tau Correlation Matrix] [paper Kendall tau...] [2010-12-25 20:18:46] [eeb33d252044f8583501f5ba0605ad6d]
Feedback Forum

Post a new message
Dataseries X:
130	461	14,1	2
127	463	14,8	2,3
122	462	16,8	2,8
117	456	15,4	2,4
112	455	15,2	2,3
113	456	16,9	2,7
149	472	14,1	2,7
157	472	14,7	2,9
157	471	16,5	3
147	465	15,2	2,2
137	459	17,6	2,3
132	465	18	2,8
125	468	16,9	2,8
123	467	16,7	2,8
117	463	19,7	2,2
114	460	15,9	2,6
111	462	17,4	2,8
112	461	17,7	2,5
144	476	15,2	2,4
150	476	15,7	2,3
149	471	17,2	1,9
134	453	17,7	1,7
123	443	17,9	2
116	442	16,2	2,1
117	444	17,5	1,7
111	438	16,8	1,8
105	427	19,1	1,8
102	424	16,7	1,8
95	416	18,2	1,3
93	406	18,5	1,3
124	431	17,8	1,3
130	434	16,4	1,2
124	418	18	1,4
115	412	20,3	2,2
106	404	19,5	2,9
105	409	18	3,1
105	412	20,2	3,5
101	406	19	3,6
95	398	20,2	4,4
93	397	21,5	4,1
84	385	19,7	5,1
87	390	21,1	5,8
116	413	20,2	5,9
120	413	18,2	5,4
117	401	21,3	5,5
109	397	20,4	4,8
105	397	17,2	3,2
107	409	15,8	2,7
109	419	15,1	2,1
109	424	14,5	1,9
108	428	15,8	0,6
107	430	14,3	0,7
99	424	13,9	-0,2
103	433	15,5	-1
131	456	14,3	-1,7
137	459	13,6	-0,7
135	446	16,3	-1
124	441	16,8	-0,9
118	439	16	0
121	454	16,8	0,3
121	460	16	0,8




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time2 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=110558&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 Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time2 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135







Correlations for all pairs of data series (method=kendall)
Werkloosheid-25Werkloosheid+25InvoerCPindex
Werkloosheid-2510.605-0.277-0.135
Werkloosheid+250.6051-0.396-0.107
Invoer-0.277-0.39610.363
CPindex-0.135-0.1070.3631

\begin{tabular}{lllllllll}
\hline
Correlations for all pairs of data series (method=kendall) \tabularnewline
  & Werkloosheid-25 & Werkloosheid+25 & Invoer & CPindex \tabularnewline
Werkloosheid-25 & 1 & 0.605 & -0.277 & -0.135 \tabularnewline
Werkloosheid+25 & 0.605 & 1 & -0.396 & -0.107 \tabularnewline
Invoer & -0.277 & -0.396 & 1 & 0.363 \tabularnewline
CPindex & -0.135 & -0.107 & 0.363 & 1 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=110558&T=1

[TABLE]
[ROW][C]Correlations for all pairs of data series (method=kendall)[/C][/ROW]
[ROW][C] [/C][C]Werkloosheid-25[/C][C]Werkloosheid+25[/C][C]Invoer[/C][C]CPindex[/C][/ROW]
[ROW][C]Werkloosheid-25[/C][C]1[/C][C]0.605[/C][C]-0.277[/C][C]-0.135[/C][/ROW]
[ROW][C]Werkloosheid+25[/C][C]0.605[/C][C]1[/C][C]-0.396[/C][C]-0.107[/C][/ROW]
[ROW][C]Invoer[/C][C]-0.277[/C][C]-0.396[/C][C]1[/C][C]0.363[/C][/ROW]
[ROW][C]CPindex[/C][C]-0.135[/C][C]-0.107[/C][C]0.363[/C][C]1[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=110558&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=110558&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)
Werkloosheid-25Werkloosheid+25InvoerCPindex
Werkloosheid-2510.605-0.277-0.135
Werkloosheid+250.6051-0.396-0.107
Invoer-0.277-0.39610.363
CPindex-0.135-0.1070.3631







Correlations for all pairs of data series with p-values
pairPearson rSpearman rhoKendall tau
Werkloosheid-25;Werkloosheid+250.78020.7870.6048
p-value(0)(0)(0)
Werkloosheid-25;Invoer-0.461-0.4058-0.2768
p-value(2e-04)(0.0012)(0.0019)
Werkloosheid-25;CPindex-0.2218-0.1905-0.1352
p-value(0.0859)(0.1413)(0.1299)
Werkloosheid+25;Invoer-0.6084-0.5517-0.3958
p-value(0)(0)(0)
Werkloosheid+25;CPindex-0.3793-0.2087-0.1066
p-value(0.0026)(0.1064)(0.2315)
Invoer;CPindex0.62180.49940.3633
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
Werkloosheid-25;Werkloosheid+25 & 0.7802 & 0.787 & 0.6048 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
Werkloosheid-25;Invoer & -0.461 & -0.4058 & -0.2768 \tabularnewline
p-value & (2e-04) & (0.0012) & (0.0019) \tabularnewline
Werkloosheid-25;CPindex & -0.2218 & -0.1905 & -0.1352 \tabularnewline
p-value & (0.0859) & (0.1413) & (0.1299) \tabularnewline
Werkloosheid+25;Invoer & -0.6084 & -0.5517 & -0.3958 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
Werkloosheid+25;CPindex & -0.3793 & -0.2087 & -0.1066 \tabularnewline
p-value & (0.0026) & (0.1064) & (0.2315) \tabularnewline
Invoer;CPindex & 0.6218 & 0.4994 & 0.3633 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=110558&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]Werkloosheid-25;Werkloosheid+25[/C][C]0.7802[/C][C]0.787[/C][C]0.6048[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]Werkloosheid-25;Invoer[/C][C]-0.461[/C][C]-0.4058[/C][C]-0.2768[/C][/ROW]
[ROW][C]p-value[/C][C](2e-04)[/C][C](0.0012)[/C][C](0.0019)[/C][/ROW]
[ROW][C]Werkloosheid-25;CPindex[/C][C]-0.2218[/C][C]-0.1905[/C][C]-0.1352[/C][/ROW]
[ROW][C]p-value[/C][C](0.0859)[/C][C](0.1413)[/C][C](0.1299)[/C][/ROW]
[ROW][C]Werkloosheid+25;Invoer[/C][C]-0.6084[/C][C]-0.5517[/C][C]-0.3958[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]Werkloosheid+25;CPindex[/C][C]-0.3793[/C][C]-0.2087[/C][C]-0.1066[/C][/ROW]
[ROW][C]p-value[/C][C](0.0026)[/C][C](0.1064)[/C][C](0.2315)[/C][/ROW]
[ROW][C]Invoer;CPindex[/C][C]0.6218[/C][C]0.4994[/C][C]0.3633[/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=110558&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=110558&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
pairPearson rSpearman rhoKendall tau
Werkloosheid-25;Werkloosheid+250.78020.7870.6048
p-value(0)(0)(0)
Werkloosheid-25;Invoer-0.461-0.4058-0.2768
p-value(2e-04)(0.0012)(0.0019)
Werkloosheid-25;CPindex-0.2218-0.1905-0.1352
p-value(0.0859)(0.1413)(0.1299)
Werkloosheid+25;Invoer-0.6084-0.5517-0.3958
p-value(0)(0)(0)
Werkloosheid+25;CPindex-0.3793-0.2087-0.1066
p-value(0.0026)(0.1064)(0.2315)
Invoer;CPindex0.62180.49940.3633
p-value(0)(0)(0)



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