Repository of Reproducible Computations

Free Statistics

of Irreproducible Research!

Author's title

Author

*Unverified author*

R Software Module

rwasp_pairs.wasp

Title produced by software

Kendall tau Correlation Matrix

Date of computation

Wed, 28 Nov 2007 08:16:42 -0700

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/2007/Nov/28/t1196262400a5mx79wz1p7ngh3.htm/, Retrieved Thu, 02 May 2024 05:50:14 +0000

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=14501, Retrieved Thu, 02 May 2024 05:50:14 +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

Klaas Van Pelt

Estimated Impact

190

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

-       [Kendall tau Correlation Matrix] [Kendall Tau 1] [2007-11-28 15:16:42] [8fe24fc91a1c75911c9bd96fe8e98ab3] [Current]

Feedback Forum

Post a new message

Dataseries X:

Download CSV

Histogram

Boxplots

37702	124.9	99.9	92.1
30364	120.4	99.9	91.4
32609	141.2	99.9	90.7
30212	95.4	99.9	92.5
29965	111.4	99.9	94.4
28352	113.3	99.7	94.3
25814	78	99.6	87.3
22414	77.7	99.5	85.9
20506	110.1	100.6	89
28806	109.6	100.2	83.3
22228	111	100.1	78.6
13971	97	100.2	75.7
36845	90.9	99.1	79.6
35338	94.7	99.5	78.5
35022	113.5	99.5	82.6
34777	107.1	99.6	88.7
26887	103.7	99.5	88.5
23970	107.7	99.6	84.6
22780	75.2	99.8	83.4
17351	74.1	99.9	84.4
21382	115.7	100.5	94.1
24561	116.2	100.5	100.4
17409	95.4	100.5	93.1
11514	95.9	100.5	93.1
31514	89.9	99.5	82.1
27071	102.4	99.9	88.1
29462	130.2	100.4	87.7
26105	98.9	99.6	80.2
22397	102.1	99.5	73.8
23843	117.3	99.6	75.3
21705	87.8	98.4	77.3
18089	66.6	99.9	80.1
20764	106.3	100.3	81.3
25316	104.5	100.3	81.5
17704	113.1	101.3	83.2
15548	95.1	101	80.8
28029	96.9	99.7	81.3
29383	107.1	99.4	78.8
36438	135.4	99.9	82.8
32034	120.3	100.7	84.9
22679	105.5	99.8	93.2
24319	128.4	98.8	94.7
18004	78.2	99.6	94.8
17537	85.8	99.1	103.9
20366	128.7	100.3	107
22782	131.2	100.5	118.6
19169	128.9	100.8	112.2
13807	112.3	100.6	112.2
29743	115	99.1	93.8
25591	116	98.8	96.7
29096	140.2	99	108.7
26482	121.2	99.9	112.1
22405	104.1	99.5	107.2
27044	128.8	99.2	113.1
17970	76.1	99.6	120
18730	86.9	100.1	124.4
19684	127.6	99.8	139.5
19785	108	101.6	145.8
18479	113.5	101.7	135.6
10698	117.4	101.9	135.5
31956	102.4	100	141.6
29506	117.9	102	141.2
34506	132	102	141.6
27165	95	102.9	147.1
26736	120.8	102.7	146.5
23691	128.9	102.7	144.1
18157	88.4	102.7	148.5
17328	81.1	102.7	146.9
18205	120.7	102.6	135.7
20995	131	102.6	128.5
17382	123.6	102.5	128.7
9367	139	102.5	127.6
31124	109.6	102.1	122.5
26551	108	103.5	123.4
30651	130.3	103.2	129.4
25859	111.7	102.1	135.3
25100	105.9	103.8	138.5
25778	128.7	104	140.6
20418	89.9	103.9	144.2
18688	101	104.1	144.3

Summary of compuational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	1 seconds
R Server	'Gwilym Jenkins' @ 72.249.127.135

\begin{tabular}{lllllllll}
\hline
Summary of compuational 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' @ 72.249.127.135 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=14501&T=0

[TABLE]
[ROW][C]Summary of compuational 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' @ 72.249.127.135[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=14501&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=14501&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 compuational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	1 seconds
R Server	'Gwilym Jenkins' @ 72.249.127.135

Kendall tau rank correlations for all pairs of data series
pair	tau	p-value
tau( x , y )	0.170179225495645	0.0256660939808910
tau( x , z )	-0.188130696526411	0.0151990646263256
tau( x , a )	-0.119062721787171	0.118239231954348
tau( y , z )	0.157699483301410	0.0421359350004922
tau( y , a )	0.199460979121093	0.00896002434205045
tau( z , a )	0.442371072478323	1.16449809794972e-08

\begin{tabular}{lllllllll}
\hline
Kendall tau rank correlations for all pairs of data series \tabularnewline
pair & tau & p-value \tabularnewline
tau( x , y ) & 0.170179225495645 & 0.0256660939808910 \tabularnewline
tau( x , z ) & -0.188130696526411 & 0.0151990646263256 \tabularnewline
tau( x , a ) & -0.119062721787171 & 0.118239231954348 \tabularnewline
tau( y , z ) & 0.157699483301410 & 0.0421359350004922 \tabularnewline
tau( y , a ) & 0.199460979121093 & 0.00896002434205045 \tabularnewline
tau( z , a ) & 0.442371072478323 & 1.16449809794972e-08 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=14501&T=1

[TABLE]
[ROW][C]Kendall tau rank correlations for all pairs of data series[/C][/ROW]
[ROW][C]pair[/C][C]tau[/C][C]p-value[/C][/ROW]
[ROW][C]tau( x , y )[/C][C]0.170179225495645[/C][C]0.0256660939808910[/C][/ROW]
[ROW][C]tau( x , z )[/C][C]-0.188130696526411[/C][C]0.0151990646263256[/C][/ROW]
[ROW][C]tau( x , a )[/C][C]-0.119062721787171[/C][C]0.118239231954348[/C][/ROW]
[ROW][C]tau( y , z )[/C][C]0.157699483301410[/C][C]0.0421359350004922[/C][/ROW]
[ROW][C]tau( y , a )[/C][C]0.199460979121093[/C][C]0.00896002434205045[/C][/ROW]
[ROW][C]tau( z , a )[/C][C]0.442371072478323[/C][C]1.16449809794972e-08[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=14501&T=1

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

As an alternative you can also use a QR Code:

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

Kendall tau rank correlations for all pairs of data series
pair	tau	p-value
tau( x , y )	0.170179225495645	0.0256660939808910
tau( x , z )	-0.188130696526411	0.0151990646263256
tau( x , a )	-0.119062721787171	0.118239231954348
tau( y , z )	0.157699483301410	0.0421359350004922
tau( y , a )	0.199460979121093	0.00896002434205045
tau( z , a )	0.442371072478323	1.16449809794972e-08

Figure 1

PNG link

Postscript link

PDF link

Parameters (Session):

par1 = 60 ; par2 = 1 ; par3 = 0 ; par4 = 0 ; par5 = 12 ;

Parameters (R input):

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='kendall')
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')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Kendall tau rank correlations for all pairs of data series',3,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'pair',1,TRUE)
a<-table.element(a,'tau',1,TRUE)
a<-table.element(a,'p-value',1,TRUE)
a<-table.row.end(a)
n <- length(y[,1])
n
cor.test(y[1,],y[2,],method='kendall')
for (i in 1:(n-1))
{
for (j in (i+1):n)
{
a<-table.row.start(a)
dum <- paste('tau(',dimnames(t(x))[[2]][i])
dum <- paste(dum,',')
dum <- paste(dum,dimnames(t(x))[[2]][j])
dum <- paste(dum,')')
a<-table.element(a,dum,header=TRUE)
r <- cor.test(y[i,],y[j,],method='kendall')
a<-table.element(a,r$estimate)
a<-table.element(a,r$p.value)
a<-table.row.end(a)
}
}
a<-table.end(a)
table.save(a,file='mytable.tab')

Free Statistics

Description of Statistical Computation

Tree of Dependent Computations

Dataset

Tables (Output of Computation)

Figures (Output of Computation)

Input Parameters & R Code