Repository of Reproducible Computations

Free Statistics

of Irreproducible Research!

Author's title

Author

*The author of this computation has been verified*

R Software Module

rwasp_pairs.wasp

Title produced by software

Kendall tau Correlation Matrix

Date of computation

Mon, 12 Dec 2011 10:32:32 -0500

Cite this page as follows

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2011/Dec/12/t1323704643m0yyktdx2majidd.htm/, Retrieved Fri, 03 May 2024 06:55:25 +0000

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=154049, Retrieved Fri, 03 May 2024 06:55:25 +0000

QR Codes:

Paste this QR Code to cite your computation.

Original text written by user:

IsPrivate?

No (this computation is public)

User-defined keywords

Estimated Impact

123

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]
- RMPD    [Kendall tau Correlation Matrix] [] [2011-12-12 15:32:32] [5c7f48f8aac1b583ed9f086d60a1681a] [Current]

Feedback Forum

Post a new message

Dataseries X:

Download CSV

Histogram

Boxplots

1407	118540	20	18158	22622
1072	127145	38	30461	73570
192	7215	0	1423	1929
2032	112861	49	25629	36294
3033	197581	70	48758	62378
5519	377410	104	129230	167760
1321	117604	37	27376	52443
1034	120102	46	26706	57283
1388	96175	42	26505	36614
2552	253517	62	49801	93268
1735	108994	50	46580	35439
1788	156212	65	48352	72405
1292	68810	28	13899	24044
2362	149246	48	39342	55909
1150	125503	42	27465	44689
1374	125769	47	55211	49319
1503	123467	71	74098	62075
965	56088	0	13497	2341
2164	108128	50	38338	40551
633	22762	12	52505	11621
837	48554	16	10663	18741
1991	159102	73	74484	84202
1450	139108	29	28895	15334
1765	92058	38	32827	28024
1612	126311	50	36188	53306
1173	113807	33	28173	37918
1603	133994	44	54926	54819
1046	131810	59	38900	89058
2176	279065	47	88530	103354
1767	158146	40	35482	70239
1167	107352	33	26730	33045
1394	149827	50	29806	63852
1733	90912	41	41799	30905
2374	169510	73	54289	24242
1852	162204	43	36805	78907
1	0	0	0	0
1677	163310	41	33146	36005
1505	92342	44	23333	31972
1813	98357	31	47686	35853
1648	178277	71	77783	115301
1560	134927	61	36042	47689
1301	104577	27	34541	34223
784	81614	21	75620	43431
1580	119111	42	60610	52220
897	83923	40	55041	33863
959	109885	34	32087	46879
567	52851	15	16356	23228
1519	153621	46	40161	42827
1625	152572	43	55459	65765
1817	114073	38	36679	38167
658	48188	12	22346	14812
1187	90994	42	27377	32615
1863	215588	56	50273	82188
1559	176489	41	32104	51763
1340	138871	46	27016	59325
1342	104416	30	19715	48976
1505	156186	44	33629	43384
669	60368	25	27084	26692
814	95391	42	32352	53279
2189	126048	28	51845	20652
1291	96388	33	26591	38338
1378	77353	32	29677	36735
1008	83867	22	54237	42764
823	79772	30	20284	44331
789	96945	13	22741	41354
1135	88300	35	34178	47879
1875	126203	39	69551	103793
2289	110681	68	29653	52235
817	81299	32	38071	49825
340	31970	5	4157	4105
2354	185321	53	28321	58687
992	87611	33	40195	40745
984	73165	44	48158	33187
1357	82167	36	13310	14063
2048	109099	52	78474	37407
933	49164	0	6386	7190
1600	105181	49	31588	49562
1821	105922	29	61254	76324
816	60138	16	21152	21928
1121	73422	33	41272	27860
800	67727	48	34165	28078
1447	188098	33	37054	49577
750	51185	24	12368	28145
1171	84448	33	23168	36241
662	41956	16	16380	10824
1284	110827	32	41242	46892
1980	167055	41	48259	60865
879	63603	36	20790	22933
869	64995	22	34585	20787
1000	93424	26	35672	43978
1240	97229	37	52168	51305
1771	112819	57	53933	55593
945	111538	20	34474	51648
1042	102220	24	43753	30552
629	38721	18	36456	23470
1770	168764	37	51183	77530
1454	151588	70	52742	57299
222	19349	13	3895	9604
1529	125440	18	37076	34684
1303	101954	32	24079	41094
552	43803	8	2325	3439
708	47062	38	29354	25171
998	104220	41	30341	23437
872	85939	24	18992	34086
584	58660	23	15292	24649
596	27676	2	5842	2342
926	93448	52	28918	45571
576	43284	5	3738	3255
0	0	0	0	0
868	64333	43	95352	30002
736	57050	18	37478	19360
998	96933	41	26839	43320
915	70088	45	26783	35513
782	65494	29	33392	23536
78	3616	0	0	0
0	0	0	0	0
782	135104	32	25446	54438
1159	95554	41	60038	56815
1646	120307	17	28162	33838
749	84336	24	33298	32366
778	43410	7	2781	13
1335	131452	62	37121	55082
806	79015	30	22698	31334
1390	88043	49	27615	16612
680	57578	3	32689	5084
285	19764	10	5752	9927
1335	105757	42	23164	47413
840	96410	18	20304	27389
1231	105056	39	34409	30425
256	11796	1	0	0
80	7627	0	0	0
1163	117413	29	92538	33510
41	6836	0	0	0
1540	131955	45	46037	40389
42	5118	5	0	0
528	40248	8	5444	6012
0	0	0	0	0
799	77813	16	23924	22205
1086	67140	16	52230	17231
81	7131	0	0	0
61	4194	0	0	0
849	60378	15	8019	11017
970	96971	40	34542	46741
964	83484	17	21157	39869

Summary of computational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	2 seconds
R Server	'George Udny Yule' @ yule.wessa.net

\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 & 'George Udny Yule' @ yule.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=154049&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]'George Udny Yule' @ yule.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=154049&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=154049&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 Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	2 seconds
R Server	'George Udny Yule' @ yule.wessa.net

Correlations for all pairs of data series (method=pearson)
	pvw	rfc	blogs	comkar	comsec
pvw	1	0.879	0.79	0.682	0.748
rfc	0.879	1	0.793	0.687	0.863
blogs	0.79	0.793	1	0.684	0.784
comkar	0.682	0.687	0.684	1	0.722
comsec	0.748	0.863	0.784	0.722	1

\begin{tabular}{lllllllll}
\hline
Correlations for all pairs of data series (method=pearson) \tabularnewline
  & pvw & rfc & blogs & comkar & comsec \tabularnewline
pvw & 1 & 0.879 & 0.79 & 0.682 & 0.748 \tabularnewline
rfc & 0.879 & 1 & 0.793 & 0.687 & 0.863 \tabularnewline
blogs & 0.79 & 0.793 & 1 & 0.684 & 0.784 \tabularnewline
comkar & 0.682 & 0.687 & 0.684 & 1 & 0.722 \tabularnewline
comsec & 0.748 & 0.863 & 0.784 & 0.722 & 1 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=154049&T=1

[TABLE]
[ROW][C]Correlations for all pairs of data series (method=pearson)[/C][/ROW]
[ROW][C] [/C][C]pvw[/C][C]rfc[/C][C]blogs[/C][C]comkar[/C][C]comsec[/C][/ROW]
[ROW][C]pvw[/C][C]1[/C][C]0.879[/C][C]0.79[/C][C]0.682[/C][C]0.748[/C][/ROW]
[ROW][C]rfc[/C][C]0.879[/C][C]1[/C][C]0.793[/C][C]0.687[/C][C]0.863[/C][/ROW]
[ROW][C]blogs[/C][C]0.79[/C][C]0.793[/C][C]1[/C][C]0.684[/C][C]0.784[/C][/ROW]
[ROW][C]comkar[/C][C]0.682[/C][C]0.687[/C][C]0.684[/C][C]1[/C][C]0.722[/C][/ROW]
[ROW][C]comsec[/C][C]0.748[/C][C]0.863[/C][C]0.784[/C][C]0.722[/C][C]1[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=154049&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=154049&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)
	pvw	rfc	blogs	comkar	comsec
pvw	1	0.879	0.79	0.682	0.748
rfc	0.879	1	0.793	0.687	0.863
blogs	0.79	0.793	1	0.684	0.784
comkar	0.682	0.687	0.684	1	0.722
comsec	0.748	0.863	0.784	0.722	1

Correlations for all pairs of data series with p-values
pair	Pearson r	Spearman rho	Kendall tau
pvw;rfc	0.8785	0.8615	0.6937
p-value	(0)	(0)	(0)
pvw;blogs	0.7902	0.764	0.5855
p-value	(0)	(0)	(0)
pvw;comkar	0.6818	0.6503	0.4864
p-value	(0)	(0)	(0)
pvw;comsec	0.7477	0.7001	0.533
p-value	(0)	(0)	(0)
rfc;blogs	0.7931	0.7654	0.5862
p-value	(0)	(0)	(0)
rfc;comkar	0.6872	0.6409	0.4809
p-value	(0)	(0)	(0)
rfc;comsec	0.8627	0.8355	0.6709
p-value	(0)	(0)	(0)
blogs;comkar	0.6842	0.6367	0.4757
p-value	(0)	(0)	(0)
blogs;comsec	0.7844	0.7608	0.5966
p-value	(0)	(0)	(0)
comkar;comsec	0.7217	0.656	0.4975
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
pvw;rfc & 0.8785 & 0.8615 & 0.6937 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
pvw;blogs & 0.7902 & 0.764 & 0.5855 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
pvw;comkar & 0.6818 & 0.6503 & 0.4864 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
pvw;comsec & 0.7477 & 0.7001 & 0.533 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
rfc;blogs & 0.7931 & 0.7654 & 0.5862 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
rfc;comkar & 0.6872 & 0.6409 & 0.4809 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
rfc;comsec & 0.8627 & 0.8355 & 0.6709 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
blogs;comkar & 0.6842 & 0.6367 & 0.4757 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
blogs;comsec & 0.7844 & 0.7608 & 0.5966 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
comkar;comsec & 0.7217 & 0.656 & 0.4975 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=154049&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]pvw;rfc[/C][C]0.8785[/C][C]0.8615[/C][C]0.6937[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]pvw;blogs[/C][C]0.7902[/C][C]0.764[/C][C]0.5855[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]pvw;comkar[/C][C]0.6818[/C][C]0.6503[/C][C]0.4864[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]pvw;comsec[/C][C]0.7477[/C][C]0.7001[/C][C]0.533[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]rfc;blogs[/C][C]0.7931[/C][C]0.7654[/C][C]0.5862[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]rfc;comkar[/C][C]0.6872[/C][C]0.6409[/C][C]0.4809[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]rfc;comsec[/C][C]0.8627[/C][C]0.8355[/C][C]0.6709[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]blogs;comkar[/C][C]0.6842[/C][C]0.6367[/C][C]0.4757[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]blogs;comsec[/C][C]0.7844[/C][C]0.7608[/C][C]0.5966[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]comkar;comsec[/C][C]0.7217[/C][C]0.656[/C][C]0.4975[/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=154049&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=154049&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
pair	Pearson r	Spearman rho	Kendall tau
pvw;rfc	0.8785	0.8615	0.6937
p-value	(0)	(0)	(0)
pvw;blogs	0.7902	0.764	0.5855
p-value	(0)	(0)	(0)
pvw;comkar	0.6818	0.6503	0.4864
p-value	(0)	(0)	(0)
pvw;comsec	0.7477	0.7001	0.533
p-value	(0)	(0)	(0)
rfc;blogs	0.7931	0.7654	0.5862
p-value	(0)	(0)	(0)
rfc;comkar	0.6872	0.6409	0.4809
p-value	(0)	(0)	(0)
rfc;comsec	0.8627	0.8355	0.6709
p-value	(0)	(0)	(0)
blogs;comkar	0.6842	0.6367	0.4757
p-value	(0)	(0)	(0)
blogs;comsec	0.7844	0.7608	0.5966
p-value	(0)	(0)	(0)
comkar;comsec	0.7217	0.656	0.4975
p-value	(0)	(0)	(0)

Figure 1

PNG link

Postscript link

PDF link

Parameters (Session):

par1 = pearson ;

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

Free Statistics

Description of Statistical Computation

Tree of Dependent Computations

Dataset

Tables (Output of Computation)

Figures (Output of Computation)

Input Parameters & R Code