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 computationFri, 13 Nov 2009 09:15:10 -0700
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2009/Nov/13/t1258128979y42lhcwrlpqrzgt.htm/, Retrieved Sun, 05 May 2024 13:03:23 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=56861, Retrieved Sun, 05 May 2024 13:03:23 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact132
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Kendall tau Correlation Matrix] [3/11/2009] [2009-11-02 21:25:00] [b98453cac15ba1066b407e146608df68]
-    D    [Kendall tau Correlation Matrix] [Kendall tau Corre...] [2009-11-13 16:15:10] [befe6dd6a614b6d3a2a74a47a0a4f514] [Current]
Feedback Forum

Post a new message
Dataseries X:
20.7	8.3	17.3	0.9
21.2	8.2	17.8	1.3
20.7	8	13.9	1.4
19.6	7.9	12.6	1.5
18.6	7.6	14.4	1.1
18.7	7.6	11.5	1.6
23.8	8.3	13.1	1.5
24.9	8.4	16.2	1.6
24.8	8.4	19.4	1.7
23.8	8.4	20.7	1.6
22.3	8.4	23.4	1.7
21.7	8.6	24.2	1.6
20.7	8.9	23.9	1.6
19.7	8.8	24.8	1.3
18.4	8.3	25.7	1.1
17.4	7.5	27.5	1.6
17	7.2	28.8	1.9
18	7.4	28.4	1.6
23.8	8.8	31.8	1.7
25.5	9.3	29.2	1.6
25.6	9.3	30	1.4
23.7	8.7	29.4	2.1
22	8.2	27.5	1.9
21.3	8.3	25.6	1.7
20.7	8.5	27	1.8
20.4	8.6	25.2	2
20.3	8.5	24.6	2.5
20.4	8.2	21.9	2.1
19.8	8.1	20.8	2.1
19.5	7.9	21.3	2.3
23.1	8.6	20.6	2.4
23.5	8.7	21.2	2.4
23.5	8.7	25.1	2.3
22.9	8.5	27.1	1.7
21.9	8.4	28.1	2
21.5	8.5	29	2.3
20.5	8.7	29.4	2
20.2	8.7	31	2
19.4	8.6	32.1	1.3
19.2	8.5	34.1	1.7
18.8	8.3	34.3	1.9
18.8	8	35.1	1.7
22.6	8.2	36.5	1.6
23.3	8.1	35.3	1.7
23	8.1	34.3	1.8
21.4	8	35.7	1.9
19.9	7.9	34.5	1.9
18.8	7.9	35.7	1.9
18.6	8	34.9	2
18.4	8	33.6	2.1
18.6	7.9	32.5	1.9
19.9	8	34.6	1.9
19.2	7.7	35.6	1.3
18.4	7.2	37.3	1.3
21.1	7.5	36	1.4
20.5	7.3	33.8	1.2
19.1	7	33.5	1.3
18.1	7	31.1	1.8
17	7	31.2	2.2
17.1	7.2	28.9	2.6
17.4	7.3	30.7	2.8
16.8	7.1	30.5	3.1
15.3	6.8	30.8	3.9
14.3	6.4	26.9	3.7
13.4	6.1	29.3	4.6
15.3	6.5	27.5	5.1
22.1	7.7	25.5	5.2
23.7	7.9	26	4.9
22.2	7.5	20.8	5.1
19.5	6.9	14.9	4.8
16.6	6.6	6.9	3.9
17.3	6.9	0	3.5




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time0 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135
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 & 0 seconds \tabularnewline
R Server & 'Gwilym Jenkins' @ 72.249.127.135 \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=56861&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]0 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Gwilym Jenkins' @ 72.249.127.135[/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=56861&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=56861&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 time0 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135
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'.







Kendall tau rank correlations for all pairs of data series
pairtaup-value
tau( WJ , TW )0.5502179300926052.38835617949462e-11
tau( WJ , CJ )-0.09572664123239460.237263375852732
tau( WJ , GI )-0.1564964520433230.058952182589252
tau( TW , CJ )-0.05035834469077680.539385875389011
tau( TW , GI )-0.1924157984698850.0219003234642248
tau( CJ , GI )-0.06805873187293250.409669813771907

\begin{tabular}{lllllllll}
\hline
Kendall tau rank correlations for all pairs of data series \tabularnewline
pair & tau & p-value \tabularnewline
tau( WJ , TW ) & 0.550217930092605 & 2.38835617949462e-11 \tabularnewline
tau( WJ , CJ ) & -0.0957266412323946 & 0.237263375852732 \tabularnewline
tau( WJ , GI
 ) & -0.156496452043323 & 0.058952182589252 \tabularnewline
tau( TW , CJ ) & -0.0503583446907768 & 0.539385875389011 \tabularnewline
tau( TW , GI
 ) & -0.192415798469885 & 0.0219003234642248 \tabularnewline
tau( CJ , GI
 ) & -0.0680587318729325 & 0.409669813771907 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=56861&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( WJ , TW )[/C][C]0.550217930092605[/C][C]2.38835617949462e-11[/C][/ROW]
[ROW][C]tau( WJ , CJ )[/C][C]-0.0957266412323946[/C][C]0.237263375852732[/C][/ROW]
[ROW][C]tau( WJ , GI
 )[/C][C]-0.156496452043323[/C][C]0.058952182589252[/C][/ROW]
[ROW][C]tau( TW , CJ )[/C][C]-0.0503583446907768[/C][C]0.539385875389011[/C][/ROW]
[ROW][C]tau( TW , GI
 )[/C][C]-0.192415798469885[/C][C]0.0219003234642248[/C][/ROW]
[ROW][C]tau( CJ , GI
 )[/C][C]-0.0680587318729325[/C][C]0.409669813771907[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=56861&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=56861&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
pairtaup-value
tau( WJ , TW )0.5502179300926052.38835617949462e-11
tau( WJ , CJ )-0.09572664123239460.237263375852732
tau( WJ , GI )-0.1564964520433230.058952182589252
tau( TW , CJ )-0.05035834469077680.539385875389011
tau( TW , GI )-0.1924157984698850.0219003234642248
tau( CJ , GI )-0.06805873187293250.409669813771907



Parameters (Session):
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')