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 computationThu, 05 Nov 2009 03:23:45 -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/Nov/05/t1257416690h4ym10itl1kaxui.htm/, Retrieved Thu, 02 May 2024 21:13:19 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=53926, Retrieved Thu, 02 May 2024 21:13:19 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact170
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Kendall tau Correlation Matrix] [3/11/2009] [2009-11-02 21:25:00] [b98453cac15ba1066b407e146608df68]
-    D    [Kendall tau Correlation Matrix] [kendall tau corre...] [2009-11-05 10:23:45] [ea241b681aafed79da4b5b99fad98471] [Current]
Feedback Forum

Post a new message
Dataseries X:
216234	258961	87130	562325
213587	258647	88620	560854
209465	256304	89563	555332
204045	250498	89056	543599
200237	247883	88542	536662
203666	249552	89504	542722
241476	262626	89428	593530
260307	264416	86040	610763
243324	273049	96240	612613
244460	272441	94423	611324
233575	267564	93028	594167
237217	265952	92285	595454
235243	263937	91685	590865
230354	264765	94260	589379
227184	263386	93858	584428
221678	258985	92437	573100
217142	257334	92980	567456
219452	257477	92099	569028
256446	271486	92803	620735
265845	274488	88551	628884
248624	281274	98334	628232
241114	272674	98329	612117
229245	269704	96455	595404
231805	268227	97109	597141
219277	276444	97687	593408
219313	272247	98512	590072
212610	268516	98673	579799
214771	263406	96028	574205
211142	263619	98014	572775
211457	265905	95580	572942
240048	281681	97838	619567
240636	287413	97760	625809
230580	289423	99913	619916
208795	281242	97588	587625
197922	273878	93942	565742
194596	269022	93656	557274
194581	272630	93365	560576
185686	270287	92881	548854
178106	260447	93120	531673
172608	262248	91063	525919
167302	252806	90930	511038
168053	238663	91946	498662
202300	258438	94624	555362
202388	266719	95484	564591
182516	263279	95862	541657
173476	258064	95530	527070
166444	248828	94574	509846
171297	248284	94677	514258
169701	253376	93845	516922
164182	251846	91533	507561
161914	239494	91214	492622
159612	239709	90922	490243
151001	228793	89563	469357
158114	229521	89945	477580
186530	249999	91850	528379
187069	254016	92505	533590
174330	251178	92437	517945
169362	242936	93876	506174
166827	241478	93561	501866
178037	243985	94119	516141
186412	246546	95264	528222
189226	247323	96089	532638
191563	247599	97160	536322
188906	248985	98644	536535
186005	241326	96266	523597
195309	242967	97938	536214
223532	263281	99757	586570
226899	268145	101550	596594
214126	263948	102449	580523




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

\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
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=53926&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]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=53926&T=0

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







Kendall tau rank correlations for all pairs of data series
pairtaup-value
tau( wl_vla , wl_wa )0.5541346973572041.65776281590979e-11
tau( wl_vla , wl_bru )0.1782516153546880.0303803313750584
tau( wl_vla , wl_bel )0.8533674339300940
tau( wl_wa , wl_bru )0.2712153764726830.000986808385882831
tau( wl_wa , wl_bel )0.6973572037510660
tau( wl_bru , wl_bel )0.2515991699982430.00224313650474284

\begin{tabular}{lllllllll}
\hline
Kendall tau rank correlations for all pairs of data series \tabularnewline
pair & tau & p-value \tabularnewline
tau( wl_vla , wl_wa ) & 0.554134697357204 & 1.65776281590979e-11 \tabularnewline
tau( wl_vla , wl_bru ) & 0.178251615354688 & 0.0303803313750584 \tabularnewline
tau( wl_vla , wl_bel ) & 0.853367433930094 & 0 \tabularnewline
tau( wl_wa , wl_bru ) & 0.271215376472683 & 0.000986808385882831 \tabularnewline
tau( wl_wa , wl_bel ) & 0.697357203751066 & 0 \tabularnewline
tau( wl_bru , wl_bel ) & 0.251599169998243 & 0.00224313650474284 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=53926&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( wl_vla , wl_wa )[/C][C]0.554134697357204[/C][C]1.65776281590979e-11[/C][/ROW]
[ROW][C]tau( wl_vla , wl_bru )[/C][C]0.178251615354688[/C][C]0.0303803313750584[/C][/ROW]
[ROW][C]tau( wl_vla , wl_bel )[/C][C]0.853367433930094[/C][C]0[/C][/ROW]
[ROW][C]tau( wl_wa , wl_bru )[/C][C]0.271215376472683[/C][C]0.000986808385882831[/C][/ROW]
[ROW][C]tau( wl_wa , wl_bel )[/C][C]0.697357203751066[/C][C]0[/C][/ROW]
[ROW][C]tau( wl_bru , wl_bel )[/C][C]0.251599169998243[/C][C]0.00224313650474284[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=53926&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=53926&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( wl_vla , wl_wa )0.5541346973572041.65776281590979e-11
tau( wl_vla , wl_bru )0.1782516153546880.0303803313750584
tau( wl_vla , wl_bel )0.8533674339300940
tau( wl_wa , wl_bru )0.2712153764726830.000986808385882831
tau( wl_wa , wl_bel )0.6973572037510660
tau( wl_bru , wl_bel )0.2515991699982430.00224313650474284



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