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, 18 Dec 2009 08:21:58 -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/Dec/18/t1261149781zc5yzzkd77g1x4n.htm/, Retrieved Sat, 27 Apr 2024 06:31:43 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=69392, Retrieved Sat, 27 Apr 2024 06:31:43 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact133
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Bivariate Data Series] [Bivariate dataset] [2008-01-05 23:51:08] [74be16979710d4c4e7c6647856088456]
F RMPD  [Univariate Explorative Data Analysis] [Colombia Coffee] [2008-01-07 14:21:11] [74be16979710d4c4e7c6647856088456]
F RMPD    [Univariate Data Series] [] [2009-10-14 08:30:28] [74be16979710d4c4e7c6647856088456]
- RMPD      [Partial Correlation] [Partial Correlation] [2009-12-16 14:03:17] [4d62210f0915d3a20cbf115865da7cd4]
-    D        [Partial Correlation] [Partial Correlation] [2009-12-16 14:13:22] [4d62210f0915d3a20cbf115865da7cd4]
- RMP           [Trivariate Scatterplots] [Trivariate Scatte...] [2009-12-17 11:10:11] [4d62210f0915d3a20cbf115865da7cd4]
- RMPD              [Kendall tau Correlation Matrix] [Kendall tau corre...] [2009-12-18 15:21:58] [91df150cd527c563f0151b3a845ecd72] [Current]
-   P                 [Kendall tau Correlation Matrix] [Pearson Correlati...] [2010-12-14 18:53:01] [8214fe6d084e5ad7598b249a26cc9f06]
Feedback Forum

Post a new message
Dataseries X:
14.3	15.89	113
14.2	16.93	110
15.9	20.28	107
15.3	22.52	103
15.5	23.51	98
15.1	22.59	98
15	23.51	137
12.1	24.76	148
15.8	26.08	147
16.9	25.29	139
15.1	23.38	130
13.7	25.29	128
14.8	28.42	127
14.7	31.85	123
16	30.1	118
15.4	25.45	114
15	24.95	108
15.5	26.84	111
15.1	27.52	151
11.7	27.94	159
16.3	25.23	158
16.7	26.53	148
15	27.21	138
14.9	28.53	137
14.6	30.35	136
15.3	31.21	133
17.9	32.86	126
16.4	33.2	120
15.4	35.73	114
17.9	34.53	116
15.9	36.54	153
13.9	40.1	162
17.8	40.56	161
17.9	46.14	149
17.4	42.85	139
16.7	38.22	135
16	40.18	130
16.6	42.19	127
19.1	47.56	122
17.8	47.26	117
17.2	44.03	112
18.6	49.83	113
16.3	53.35	149
15.1	58.9	157
19.2	59.64	157
17.7	56.99	147
19.1	53.2	137
18	53.24	132
17.5	57.85	125
17.8	55.69	123
21.1	55.64	117
17.2	62.52	114
19.4	64.4	111
19.8	64.65	112
17.6	67.71	144
16.2	67.21	150
19.5	59.37	149
19.9	53.26	134
20	52.42	123
17.3	55.03	116




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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=69392&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







Kendall tau rank correlations for all pairs of data series
pairtaup-value
tau( BuitlH , POil )0.4985868914949272.26351553056503e-08
tau( BuitlH , WK<25j )-0.02577319587628870.773860876575599
tau( POil , WK<25j )0.1263542122281660.156563196462220

\begin{tabular}{lllllllll}
\hline
Kendall tau rank correlations for all pairs of data series \tabularnewline
pair & tau & p-value \tabularnewline
tau( BuitlH , POil ) & 0.498586891494927 & 2.26351553056503e-08 \tabularnewline
tau( BuitlH , WK<25j ) & -0.0257731958762887 & 0.773860876575599 \tabularnewline
tau( POil , WK<25j ) & 0.126354212228166 & 0.156563196462220 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=69392&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( BuitlH , POil )[/C][C]0.498586891494927[/C][C]2.26351553056503e-08[/C][/ROW]
[ROW][C]tau( BuitlH , WK<25j )[/C][C]-0.0257731958762887[/C][C]0.773860876575599[/C][/ROW]
[ROW][C]tau( POil , WK<25j )[/C][C]0.126354212228166[/C][C]0.156563196462220[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=69392&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=69392&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( BuitlH , POil )0.4985868914949272.26351553056503e-08
tau( BuitlH , WK<25j )-0.02577319587628870.773860876575599
tau( POil , WK<25j )0.1263542122281660.156563196462220



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