Free Statistics

of Irreproducible Research!

Author's title

Author*The author of this computation has been verified*
R Software ModulePatrick.Wessarwasp_pairs.wasp
Title produced by softwareKendall tau Correlation Matrix
Date of computationSat, 18 Dec 2010 21:19:58 +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/18/t1292707110p0r6kpio7lpf73k.htm/, Retrieved Tue, 30 Apr 2024 01:51:57 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=112201, Retrieved Tue, 30 Apr 2024 01:51:57 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact148
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Kendall tau Correlation Matrix] [] [2010-12-05 17:44:33] [b98453cac15ba1066b407e146608df68]
-   PD  [Kendall tau Correlation Matrix] [pearson correlati...] [2010-12-10 10:32:01] [74deae64b71f9d77c839af86f7c687b5]
-   PD      [Kendall tau Correlation Matrix] [pearson correlati...] [2010-12-18 21:19:58] [e665313c9926a9f4bdf6ad1ee5aefad6] [Current]
Feedback Forum

Post a new message
Dataseries X:
101,82	107,34	93,63	99,85	101,76
101,68	107,34	93,63	99,91	102,37
101,68	107,34	93,63	99,87	102,38
102,45	107,34	96,13	99,86	102,86
102,45	107,34	96,13	100,10	102,87
102,45	107,34	96,13	100,10	102,92
102,45	107,34	96,13	100,12	102,95
102,45	107,34	96,13	99,95	103,02
102,45	112,60	96,13	99,94	104,08
102,52	112,60	96,13	100,18	104,16
102,52	112,60	96,13	100,31	104,24
102,85	112,60	96,13	100,65	104,33
102,85	112,61	96,13	100,65	104,73
102,85	112,61	96,13	100,69	104,86
103,25	112,61	96,13	101,26	105,03
103,25	112,61	98,73	101,26	105,62
103,25	112,61	98,73	101,38	105,63
103,25	112,61	98,73	101,38	105,63
104,45	112,61	98,73	101,38	105,94
104,45	112,61	98,73	101,44	106,61
104,45	118,65	98,73	101,40	107,69
104,80	118,65	98,73	101,40	107,78
104,80	118,65	98,73	100,58	107,93
105,29	118,65	98,73	100,58	108,48
105,29	114,29	98,73	100,58	108,14
105,29	114,29	98,73	100,59	108,48
105,29	114,29	98,73	100,81	108,48
106,04	114,29	101,67	100,75	108,89
105,94	114,29	101,67	100,75	108,93
105,94	114,29	101,67	100,96	109,21
105,94	114,29	101,67	101,31	109,47
106,28	114,29	101,67	101,64	109,80
106,48	123,33	101,67	101,46	111,73
107,19	123,33	101,67	101,73	111,85
108,14	123,33	101,67	101,73	112,12
108,22	123,33	101,67	101,64	112,15
108,22	123,33	101,67	101,77	112,17
108,61	123,33	101,67	101,74	112,67
108,61	123,33	101,67	101,89	112,80
108,61	123,33	107,94	101,89	113,44
108,61	123,33	107,94	101,93	113,53
109,06	123,33	107,94	101,93	114,53
109,06	123,33	107,94	102,32	114,51
112,93	123,33	107,94	102,41	115,05
115,84	129,03	107,94	103,58	116,67
118,57	128,76	107,94	104,12	117,07
118,57	128,76	107,94	104,10	116,92
118,86	128,76	107,94	104,15	117,00
118,98	128,76	107,94	104,15	117,02
119,27	128,76	107,94	104,16	117,35
119,39	128,76	107,94	102,94	117,36
119,49	128,76	110,30	103,07	117,82
119,59	128,76	110,30	103,04	117,88
120,12	128,76	110,30	103,06	118,24
120,14	128,76	110,30	103,05	118,50
120,14	128,76	110,30	102,95	118,80
120,14	132,63	110,30	102,95	119,76
120,14	132,63	110,30	103,05	120,09
120,62	132,63	110,30	102,65	120,16




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

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







Correlations for all pairs of data series (method=pearson)
biosschouwburgeedagsacttractiehuurDVDvrijetijdsbesteding
bios10.920.9290.9170.949
schouwburg0.9210.9240.9030.976
eedagsacttractie0.9290.92410.8920.972
huurDVD0.9170.9030.89210.909
vrijetijdsbesteding0.9490.9760.9720.9091

\begin{tabular}{lllllllll}
\hline
Correlations for all pairs of data series (method=pearson) \tabularnewline
  & bios & schouwburg & eedagsacttractie & huurDVD & vrijetijdsbesteding \tabularnewline
bios & 1 & 0.92 & 0.929 & 0.917 & 0.949 \tabularnewline
schouwburg & 0.92 & 1 & 0.924 & 0.903 & 0.976 \tabularnewline
eedagsacttractie & 0.929 & 0.924 & 1 & 0.892 & 0.972 \tabularnewline
huurDVD & 0.917 & 0.903 & 0.892 & 1 & 0.909 \tabularnewline
vrijetijdsbesteding & 0.949 & 0.976 & 0.972 & 0.909 & 1 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=112201&T=1

[TABLE]
[ROW][C]Correlations for all pairs of data series (method=pearson)[/C][/ROW]
[ROW][C] [/C][C]bios[/C][C]schouwburg[/C][C]eedagsacttractie[/C][C]huurDVD[/C][C]vrijetijdsbesteding[/C][/ROW]
[ROW][C]bios[/C][C]1[/C][C]0.92[/C][C]0.929[/C][C]0.917[/C][C]0.949[/C][/ROW]
[ROW][C]schouwburg[/C][C]0.92[/C][C]1[/C][C]0.924[/C][C]0.903[/C][C]0.976[/C][/ROW]
[ROW][C]eedagsacttractie[/C][C]0.929[/C][C]0.924[/C][C]1[/C][C]0.892[/C][C]0.972[/C][/ROW]
[ROW][C]huurDVD[/C][C]0.917[/C][C]0.903[/C][C]0.892[/C][C]1[/C][C]0.909[/C][/ROW]
[ROW][C]vrijetijdsbesteding[/C][C]0.949[/C][C]0.976[/C][C]0.972[/C][C]0.909[/C][C]1[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=112201&T=1

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

As an alternative you can also use a QR Code:  

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

Correlations for all pairs of data series (method=pearson)
biosschouwburgeedagsacttractiehuurDVDvrijetijdsbesteding
bios10.920.9290.9170.949
schouwburg0.9210.9240.9030.976
eedagsacttractie0.9290.92410.8920.972
huurDVD0.9170.9030.89210.909
vrijetijdsbesteding0.9490.9760.9720.9091







Correlations for all pairs of data series with p-values
pairPearson rSpearman rhoKendall tau
bios;schouwburg0.91990.97210.8875
p-value(0)(0)(0)
bios;eedagsacttractie0.92940.98070.9199
p-value(0)(0)(0)
bios;huurDVD0.91660.92580.7827
p-value(0)(0)(0)
bios;vrijetijdsbesteding0.94920.9980.9769
p-value(0)(0)(0)
schouwburg;eedagsacttractie0.92370.93990.8644
p-value(0)(0)(0)
schouwburg;huurDVD0.90270.92480.7811
p-value(0)(0)(0)
schouwburg;vrijetijdsbesteding0.97570.9740.8832
p-value(0)(0)(0)
eedagsacttractie;huurDVD0.89190.91420.7841
p-value(0)(0)(0)
eedagsacttractie;vrijetijdsbesteding0.97220.98190.9108
p-value(0)(0)(0)
huurDVD;vrijetijdsbesteding0.90940.92690.7835
p-value(0)(0)(0)

\begin{tabular}{lllllllll}
\hline
Correlations for all pairs of data series with p-values \tabularnewline
pair & Pearson r & Spearman rho & Kendall tau \tabularnewline
bios;schouwburg & 0.9199 & 0.9721 & 0.8875 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
bios;eedagsacttractie & 0.9294 & 0.9807 & 0.9199 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
bios;huurDVD & 0.9166 & 0.9258 & 0.7827 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
bios;vrijetijdsbesteding & 0.9492 & 0.998 & 0.9769 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
schouwburg;eedagsacttractie & 0.9237 & 0.9399 & 0.8644 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
schouwburg;huurDVD & 0.9027 & 0.9248 & 0.7811 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
schouwburg;vrijetijdsbesteding & 0.9757 & 0.974 & 0.8832 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
eedagsacttractie;huurDVD & 0.8919 & 0.9142 & 0.7841 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
eedagsacttractie;vrijetijdsbesteding & 0.9722 & 0.9819 & 0.9108 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
huurDVD;vrijetijdsbesteding & 0.9094 & 0.9269 & 0.7835 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=112201&T=2

[TABLE]
[ROW][C]Correlations for all pairs of data series with p-values[/C][/ROW]
[ROW][C]pair[/C][C]Pearson r[/C][C]Spearman rho[/C][C]Kendall tau[/C][/ROW]
[ROW][C]bios;schouwburg[/C][C]0.9199[/C][C]0.9721[/C][C]0.8875[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]bios;eedagsacttractie[/C][C]0.9294[/C][C]0.9807[/C][C]0.9199[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]bios;huurDVD[/C][C]0.9166[/C][C]0.9258[/C][C]0.7827[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]bios;vrijetijdsbesteding[/C][C]0.9492[/C][C]0.998[/C][C]0.9769[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]schouwburg;eedagsacttractie[/C][C]0.9237[/C][C]0.9399[/C][C]0.8644[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]schouwburg;huurDVD[/C][C]0.9027[/C][C]0.9248[/C][C]0.7811[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]schouwburg;vrijetijdsbesteding[/C][C]0.9757[/C][C]0.974[/C][C]0.8832[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]eedagsacttractie;huurDVD[/C][C]0.8919[/C][C]0.9142[/C][C]0.7841[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]eedagsacttractie;vrijetijdsbesteding[/C][C]0.9722[/C][C]0.9819[/C][C]0.9108[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]huurDVD;vrijetijdsbesteding[/C][C]0.9094[/C][C]0.9269[/C][C]0.7835[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=112201&T=2

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

As an alternative you can also use a QR Code:  

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

Correlations for all pairs of data series with p-values
pairPearson rSpearman rhoKendall tau
bios;schouwburg0.91990.97210.8875
p-value(0)(0)(0)
bios;eedagsacttractie0.92940.98070.9199
p-value(0)(0)(0)
bios;huurDVD0.91660.92580.7827
p-value(0)(0)(0)
bios;vrijetijdsbesteding0.94920.9980.9769
p-value(0)(0)(0)
schouwburg;eedagsacttractie0.92370.93990.8644
p-value(0)(0)(0)
schouwburg;huurDVD0.90270.92480.7811
p-value(0)(0)(0)
schouwburg;vrijetijdsbesteding0.97570.9740.8832
p-value(0)(0)(0)
eedagsacttractie;huurDVD0.89190.91420.7841
p-value(0)(0)(0)
eedagsacttractie;vrijetijdsbesteding0.97220.98190.9108
p-value(0)(0)(0)
huurDVD;vrijetijdsbesteding0.90940.92690.7835
p-value(0)(0)(0)



Parameters (Session):
par1 = 1 ; par2 = Include Monthly Dummies ; par3 = Linear Trend ;
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=par1)
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')
n <- length(y[,1])
n
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,paste('Correlations for all pairs of data series (method=',par1,')',sep=''),n+1,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,' ',header=TRUE)
for (i in 1:n) {
a<-table.element(a,dimnames(t(x))[[2]][i],header=TRUE)
}
a<-table.row.end(a)
for (i in 1:n) {
a<-table.row.start(a)
a<-table.element(a,dimnames(t(x))[[2]][i],header=TRUE)
for (j in 1:n) {
r <- cor.test(y[i,],y[j,],method=par1)
a<-table.element(a,round(r$estimate,3))
}
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable.tab')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Correlations for all pairs of data series with p-values',4,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'pair',1,TRUE)
a<-table.element(a,'Pearson r',1,TRUE)
a<-table.element(a,'Spearman rho',1,TRUE)
a<-table.element(a,'Kendall tau',1,TRUE)
a<-table.row.end(a)
cor.test(y[1,],y[2,],method=par1)
for (i in 1:(n-1))
{
for (j in (i+1):n)
{
a<-table.row.start(a)
dum <- paste(dimnames(t(x))[[2]][i],';',dimnames(t(x))[[2]][j],sep='')
a<-table.element(a,dum,header=TRUE)
rp <- cor.test(y[i,],y[j,],method='pearson')
a<-table.element(a,round(rp$estimate,4))
rs <- cor.test(y[i,],y[j,],method='spearman')
a<-table.element(a,round(rs$estimate,4))
rk <- cor.test(y[i,],y[j,],method='kendall')
a<-table.element(a,round(rk$estimate,4))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'p-value',header=T)
a<-table.element(a,paste('(',round(rp$p.value,4),')',sep=''))
a<-table.element(a,paste('(',round(rs$p.value,4),')',sep=''))
a<-table.element(a,paste('(',round(rk$p.value,4),')',sep=''))
a<-table.row.end(a)
}
}
a<-table.end(a)
table.save(a,file='mytable1.tab')