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 computationTue, 28 Dec 2010 20:25:55 +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/28/t1293567824ufd1tky1ip8gstr.htm/, Retrieved Sun, 05 May 2024 02:14:31 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=116543, Retrieved Sun, 05 May 2024 02:14:31 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact110
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Kendall tau Correlation Matrix] [Kendell tau corre...] [2009-11-12 14:06:21] [21324e9cdf3569788a3d630236984d87]
-   PD    [Kendall tau Correlation Matrix] [] [2010-12-28 20:25:55] [1d208f56d63f78e3037c4c685f0bba30] [Current]
Feedback Forum

Post a new message
Dataseries X:
112,3	112,9	88,7	105,1
117,3	130,5	94,6	114,9
111,1	137,9	98,7	106,4
102,2	115	84,2	104,5
104,3	116,8	87,7	121,6
122,9	140,9	103,3	141,4
107,6	120,7	88,2	99
121,3	134,2	93,4	126,7
131,5	147,3	106,3	134,1
89	112,4	73,1	81,3
104,4	107,1	78,6	88,6
128,9	128,4	101,6	132,7
135,9	137,7	101,4	132,9
133,3	135	98,5	134,4
121,3	151	99	103,7
120,5	137,4	89,5	119,7
120,4	132,4	83,5	115
137,9	161,3	97,4	132,9
126,1	139,8	87,8	108,5
133,2	146	90,4	113,9
151,1	166,5	101,6	142
105	143,3	80	97,7
119	121	81,7	92,2
140,4	152,6	96,4	128,8
156,6	154,4	110,2	134,9
137,1	154,6	101,1	128,2
122,7	158	89,3	114,8
125,8	142,6	90	117,9
139,3	153,4	95,4	119,1
134,9	163,4	100,3	120,7
149,2	167,3	99,5	129,1
132,3	154,8	93,9	117,6
149	165,7	100,6	129,2
117,2	144,7	84,7	100
119,6	120,9	81,6	87
152	152,8	109	128
149,4	160,2	99	127,7
127,3	128,3	81,1	93,4
114,1	150,5	81,8	84,1
102,1	117	66,5	71,7
107,7	116	66,4	83,2
104,4	133,3	86,3	89,1
102,1	116,4	73,6	79,6
96	104	71,5	62,8
109,3	126,6	87,2	95,1
90	92,9	65,3	63,6
83,9	83,6	69,7	61,4
112	112,8	95,5	98,2
114,3	113,2	86,3	95,3
103,6	118,5	81	81,5
91,7	125,5	88,7	85,5
80,8	91,3	71,9	71,1
87,2	105,4	78,6	78,1
109,2	121,3	96	103
102,7	106,9	81,1	86
95,1	109,4	77,5	86,2
117,5	132,6	97,3	105,7
85,1	96,8	78,6	57,2
92,1	100,3	79	73,7
113,5	119,2	93,4	120,5





Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time1 seconds
R Server'RServer@AstonUniversity' @ vre.aston.ac.uk
R Framework error message
Warning: there are blank lines in the 'Data X' field.
Please, use NA for missing data - blank lines are simply
 deleted and are NOT treated as missing values.

\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 & 'RServer@AstonUniversity' @ vre.aston.ac.uk \tabularnewline
R Framework error message & 
Warning: there are blank lines in the 'Data X' field.
Please, use NA for missing data - blank lines are simply
 deleted and are NOT treated as missing values.
\tabularnewline \hline \end{tabular} %Source: https://freestatistics.org/blog/index.php?pk=116543&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]'RServer@AstonUniversity' @ vre.aston.ac.uk[/C][/ROW]
[ROW][C]R Framework error message[/C][C]
Warning: there are blank lines in the 'Data X' field.
Please, use NA for missing data - blank lines are simply
 deleted and are NOT treated as missing values.
[/C][/ROW] [/TABLE] Source: https://freestatistics.org/blog/index.php?pk=116543&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=116543&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 time1 seconds
R Server'RServer@AstonUniversity' @ vre.aston.ac.uk
R Framework error message
Warning: there are blank lines in the 'Data X' field.
Please, use NA for missing data - blank lines are simply
 deleted and are NOT treated as missing values.







Kendall tau rank correlations for all pairs of data series
pairtaup-value
tau( ProMet , MaApWe )0.694939463994084.44089209850063e-15
tau( ProMet , ElElAp )0.6201823027771292.96762614482304e-12
tau( ProMet , Transp )0.6832580278445111.31006316905768e-14
tau( MaApWe , ElElAp )0.5647143804002562.01694883017467e-10
tau( MaApWe , Transp )0.581520226662575.27282661977324e-11
tau( ElElAp , Transp )0.7263474743516942.22044604925031e-16

\begin{tabular}{lllllllll}
\hline
Kendall tau rank correlations for all pairs of data series \tabularnewline
pair & tau & p-value \tabularnewline
tau( ProMet , MaApWe ) & 0.69493946399408 & 4.44089209850063e-15 \tabularnewline
tau( ProMet , ElElAp ) & 0.620182302777129 & 2.96762614482304e-12 \tabularnewline
tau( ProMet , Transp ) & 0.683258027844511 & 1.31006316905768e-14 \tabularnewline
tau( MaApWe , ElElAp ) & 0.564714380400256 & 2.01694883017467e-10 \tabularnewline
tau( MaApWe , Transp ) & 0.58152022666257 & 5.27282661977324e-11 \tabularnewline
tau( ElElAp , Transp ) & 0.726347474351694 & 2.22044604925031e-16 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=116543&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( ProMet , MaApWe )[/C][C]0.69493946399408[/C][C]4.44089209850063e-15[/C][/ROW]
[ROW][C]tau( ProMet , ElElAp )[/C][C]0.620182302777129[/C][C]2.96762614482304e-12[/C][/ROW]
[ROW][C]tau( ProMet , Transp )[/C][C]0.683258027844511[/C][C]1.31006316905768e-14[/C][/ROW]
[ROW][C]tau( MaApWe , ElElAp )[/C][C]0.564714380400256[/C][C]2.01694883017467e-10[/C][/ROW]
[ROW][C]tau( MaApWe , Transp )[/C][C]0.58152022666257[/C][C]5.27282661977324e-11[/C][/ROW]
[ROW][C]tau( ElElAp , Transp )[/C][C]0.726347474351694[/C][C]2.22044604925031e-16[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=116543&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=116543&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( ProMet , MaApWe )0.694939463994084.44089209850063e-15
tau( ProMet , ElElAp )0.6201823027771292.96762614482304e-12
tau( ProMet , Transp )0.6832580278445111.31006316905768e-14
tau( MaApWe , ElElAp )0.5647143804002562.01694883017467e-10
tau( MaApWe , Transp )0.581520226662575.27282661977324e-11
tau( ElElAp , Transp )0.7263474743516942.22044604925031e-16



Parameters (Session):
par1 = pearson ;
Parameters (R input):
par1 = pearson ;
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')