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, 09 Dec 2009 07:37:33 -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/09/t1260369493hums4ueu7p29g9u.htm/, Retrieved Mon, 29 Apr 2024 13:22:10 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=64971, Retrieved Mon, 29 Apr 2024 13:22:10 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact105
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Kendall tau Correlation Matrix] [] [2009-12-09 14:37:33] [9adf7044e3e2072a25a3bb76b79e4d2e] [Current]
Feedback Forum

Post a new message
Dataseries X:
95.1	98.2	121.8	113.4
97.0	98.7	127.6	101.3
112.7	113.3	129.9	97.8
102.9	104.6	128.0	95.0
97.4	99.3	123.5	93.8
111.4	111.8	124.0	94.5
87.4	97.3	127.4	101.4
96.8	97.7	127.6	105.8
114.1	115.6	128.4	106.6
110.3	111.9	131.4	109.7
103.9	107.0	135.1	108.8
101.6	107.1	134.0	113.4
94.6	100.6	144.5	113.7
95.9	99.2	147.3	103.6
104.7	108.4	150.9	98.2
102.8	103.0	148.7	95.5
98.1	99.8	141.4	94.4
113.9	115.0	138.9	95.9
80.9	90.8	139.8	103.2
95.7	95.9	145.6	104.1
113.2	114.4	147.9	127.6
105.9	108.2	148.5	130.3
108.8	112.6	151.1	133.0
102.3	109.1	157.5	140.4
99.0	105.0	167.5	123.5
100.7	105.0	172.3	116.9
115.5	118.5	173.5	115.9
100.7	103.7	187.5	113.1
109.9	112.5	205.5	112.1
114.6	116.6	195.1	112.4
85.4	96.6	204.5	118.9
100.5	101.9	204.5	117.4
114.8	116.5	201.7	115.6
116.5	119.3	207.0	120.7
112.9	115.4	206.6	114.9
102.0	108.5	210.6	122.0
106.0	111.5	211.1	119.6
105.3	108.8	215.0	114.6
118.8	121.8	223.9	118.4
106.1	109.6	238.2	110.9
109.3	112.2	238.9	111.6
117.2	119.6	229.6	114.6
92.5	104.1	232.2	112.1
104.2	105.3	222.1	117.4
112.5	115.0	221.6	114.8
122.4	124.1	227.3	123.4
113.3	116.8	221.0	118.1
100.0	107.5	213.6	121.9
110.7	115.6	243.4	123.3
112.8	116.2	253.8	120.1
109.8	116.3	265.3	116.7
117.3	119.0	268.2	106.3
109.1	111.9	268.5	106.8
115.9	118.6	266.9	107.5
96.0	106.9	268.4	108.5
99.8	103.2	250.8	113.6
116.8	118.6	231.2	111.6
115.7	118.7	192.0	117.1
99.4	102.8	171.4	120.0
94.3	100.6	160.0	123.0




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time1 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 & 1 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=64971&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]'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=64971&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=64971&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'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( TIP , Bouwnijverheid )0.8463071028290280
tau( TIP , IGP )0.2651965048031960.00277673539619760
tau( TIP , Energie )0.1075269505974530.225529925726653
tau( Bouwnijverheid , IGP )0.3533411099433386.8777908782458e-05
tau( Bouwnijverheid , Energie )0.1932558876115370.0296031901775597
tau( IGP , Energie )0.2066233408370350.0198975072744796

\begin{tabular}{lllllllll}
\hline
Kendall tau rank correlations for all pairs of data series \tabularnewline
pair & tau & p-value \tabularnewline
tau( TIP , Bouwnijverheid ) & 0.846307102829028 & 0 \tabularnewline
tau( TIP , IGP ) & 0.265196504803196 & 0.00277673539619760 \tabularnewline
tau( TIP , Energie
 ) & 0.107526950597453 & 0.225529925726653 \tabularnewline
tau( Bouwnijverheid , IGP ) & 0.353341109943338 & 6.8777908782458e-05 \tabularnewline
tau( Bouwnijverheid , Energie
 ) & 0.193255887611537 & 0.0296031901775597 \tabularnewline
tau( IGP , Energie
 ) & 0.206623340837035 & 0.0198975072744796 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=64971&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( TIP , Bouwnijverheid )[/C][C]0.846307102829028[/C][C]0[/C][/ROW]
[ROW][C]tau( TIP , IGP )[/C][C]0.265196504803196[/C][C]0.00277673539619760[/C][/ROW]
[ROW][C]tau( TIP , Energie
 )[/C][C]0.107526950597453[/C][C]0.225529925726653[/C][/ROW]
[ROW][C]tau( Bouwnijverheid , IGP )[/C][C]0.353341109943338[/C][C]6.8777908782458e-05[/C][/ROW]
[ROW][C]tau( Bouwnijverheid , Energie
 )[/C][C]0.193255887611537[/C][C]0.0296031901775597[/C][/ROW]
[ROW][C]tau( IGP , Energie
 )[/C][C]0.206623340837035[/C][C]0.0198975072744796[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=64971&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=64971&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( TIP , Bouwnijverheid )0.8463071028290280
tau( TIP , IGP )0.2651965048031960.00277673539619760
tau( TIP , Energie )0.1075269505974530.225529925726653
tau( Bouwnijverheid , IGP )0.3533411099433386.8777908782458e-05
tau( Bouwnijverheid , Energie )0.1932558876115370.0296031901775597
tau( IGP , Energie )0.2066233408370350.0198975072744796



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