FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

Author's title

Author*Unverified author*
R Software Modulerwasp_pairs.wasp
Title produced by softwareKendall tau Correlation Matrix
Date of computationThu, 29 Nov 2007 12:42:25 -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/2007/Nov/29/t1196364724lunse4nwg2umbh7.htm/, Retrieved Fri, 03 May 2024 04:42:22 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=7583, Retrieved Fri, 03 May 2024 04:42:22 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact162

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Kendall tau Correlation Matrix] [kendall tau ] [2007-11-29 19:42:25] [0ad9b3c11abcaba622af0629507f53fa] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
8,10	4,2	0	-0,05
7,30	0,8	0	-0,02
-2,30	1,7	0	0
-1,40	4,6	0	-0,01
-4,40	0,6	0	0,01
13,90	0,2	0	-0,01
11,70	14,2	0	0
5,40	2,6	0	0,02
-6,30	6,5	0	0,42
2,30	10,3	0	0,03
3,90	-0,4	0	-0,1
4,00	5,8	0	-0,07
5,30	9,2	0	-0,04
-10,10	-4,4	0	-0,01
-10,90	10,2	0	-0,01
4,20	12,6	0	-0,02
2,30	5,2	0,09	0,01
9,40	1,1	0,01	-0,01
-1,80	-0,3	0	0,02
-14,60	4,3	0	0,02
15,80	6,1	0	0,32
0,80	3,4	0,01	-0,06
4,80	-6,3	0,01	-0,13
0,40	5,5	0	-0,14
-20,10	1,9	0,01	-0,1
1,60	0,7	0	-0,03
13,00	6,3	0	-0,02
5,70	0,4	0,01	0
-6,90	6,2	0,01	0
-14,00	4,5	0	0
20,10	9,6	0,01	0
12,60	-0,2	0	0,01
-16,50	10,8	0	0,14
6,90	0,1	0	0,07
-13,90	6,0	0	-0,02
-22,60	5,7	0,01	-0,04
0,80	0,0	0	-0,02
9,00	14,1	0	0,01
3,00	5,8	0	0
-33,10	-2,3	0	1,01
13,90	4,6	0,01	-0,02
5,70	7,6	0	0,03
-6,30	5,7	0,01	0,02
-4,80	9,4	0,01	-0,02
10,90	1,9	0	0,01
-19,00	4,9	0	0,01

Tables (Output of Computation)





Summary of compuational 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

\begin{tabular}{lllllllll}
\hline
Summary of compuational 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
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=7583&T=0

[TABLE]
[ROW][C]Summary of compuational 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]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=7583&T=0

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

As an alternative you can also use a QR Code:  
QR Code


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

Summary of compuational 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


Figures (Output of Computation)

Input Parameters & R Code

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