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

Thu, 22 Dec 2011 19:00:47 -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/22/t132459903394laqyejtsokkvi.htm/, Retrieved Fri, 03 May 2024 04:16:10 +0000

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=160107, Retrieved Fri, 03 May 2024 04:16:10 +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

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] [Pearson correlatie] [2011-12-23 00:00:47] [de50302416ae5d0bdedd77e4c0468c33] [Current]

Feedback Forum

Post a new message

Dataseries X:

Download CSV

Histogram

Boxplots

162687	95	595	115
223169	66	575	79
7215	18	72	1
164587	99	737	158
283430	141	1255	127
546996	275	2021	278
192501	61	606	95
213538	64	533	64
182282	46	687	92
336547	102	1074	130
122275	77	637	158
203938	72	743	120
119300	110	701	87
220796	122	1087	264
174005	67	422	51
156326	89	474	85
164063	60	483	100
90025	63	375	72
179987	90	929	147
47066	29	262	49
109572	64	437	40
241285	103	850	99
208339	77	652	127
164166	59	754	166
159763	89	619	41
207078	34	657	160
217028	169	695	92
201536	96	366	59
408960	124	1015	89
250260	48	1029	104
216527	46	576	81
212949	51	656	116
164248	110	812	105
276436	135	1089	380
238654	59	852	88
0	1	0	0
233971	66	1009	63
149649	55	658	138
161703	52	547	270
254893	70	826	64
269492	73	838	96
169526	62	704	62
107893	35	404	35
229714	83	848	66
139667	51	419	56
175983	102	349	46
81407	33	216	49
251259	110	796	121
239807	90	752	113
172743	60	964	190
48188	28	205	37
169355	71	506	52
335398	78	841	89
244729	81	699	73
208286	62	746	61
159913	57	547	77
232137	72	561	63
101694	26	329	75
157258	68	427	32
211586	101	993	59
181076	66	564	71
158024	86	858	92
141491	64	376	87
130108	40	471	48
166420	39	432	63
125440	44	492	41
195043	72	504	86
138708	66	887	152
116552	40	271	49
31970	15	101	40
277661	116	1135	135
167825	82	506	86
135926	69	528	62
124349	76	491	96
171518	71	698	95
108980	46	426	83
183471	61	709	112
167426	101	847	77
112510	49	367	78
92421	77	413	114
116704	83	271	55
304603	65	830	60
75101	30	334	49
145043	41	524	132
95827	48	393	49
173931	60	574	71
250424	252	695	102
115367	116	284	74
125839	66	462	49
164078	54	653	74
158931	42	684	59
190382	85	714	91
155226	59	420	68
146159	61	551	81
62641	44	396	33
258585	121	741	166
199841	71	571	97
19349	12	67	15
247280	109	877	105
164984	87	870	61
72128	30	306	11
104253	26	382	45
151090	57	435	89
147990	68	348	72
87448	42	227	27
27676	22	194	59
170326	52	413	127
132148	38	273	48
0	0	0	0
133868	36	390	58
109001	68	376	57
158833	46	495	60
150013	66	448	77
89887	63	313	71
3616	5	14	5
0	0	0	0
216479	48	445	78
177323	102	637	76
177948	102	593	124
140106	41	326	67
43410	19	292	63
206059	76	573	92
109873	45	315	58
157084	61	683	65
60493	40	174	29
19764	12	75	19
177559	57	572	64
154169	36	414	79
164249	54	562	104
11796	9	79	22
10674	9	33	7
151322	59	487	37
6836	3	11	5
174712	68	664	48
5118	3	6	1
40248	16	183	34
0	0	0	0
127628	51	342	53
88837	38	269	44
7131	4	27	0
9056	15	99	18
97191	31	322	52
157478	59	367	60
121330	22	502	50

Summary of computational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	1 seconds
R Server	'Gertrude Mary Cox' @ cox.wessa.net
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 & 'Gertrude Mary Cox' @ cox.wessa.net \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=160107&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]'Gertrude Mary Cox' @ cox.wessa.net[/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=160107&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=160107&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	1 seconds
R Server	'Gertrude Mary Cox' @ cox.wessa.net
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'.

Correlations for all pairs of data series (method=pearson)
	A	B	C	D
A	1	0.76	0.886	0.613
B	0.76	1	0.755	0.593
C	0.886	0.755	1	0.705
D	0.613	0.593	0.705	1

\begin{tabular}{lllllllll}
\hline
Correlations for all pairs of data series (method=pearson) \tabularnewline
  & A & B & C & D
 \tabularnewline
A & 1 & 0.76 & 0.886 & 0.613 \tabularnewline
B & 0.76 & 1 & 0.755 & 0.593 \tabularnewline
C & 0.886 & 0.755 & 1 & 0.705 \tabularnewline
D
 & 0.613 & 0.593 & 0.705 & 1 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=160107&T=1

[TABLE]
[ROW][C]Correlations for all pairs of data series (method=pearson)[/C][/ROW]
[ROW][C] [/C][C]A[/C][C]B[/C][C]C[/C][C]D
[/C][/ROW]
[ROW][C]A[/C][C]1[/C][C]0.76[/C][C]0.886[/C][C]0.613[/C][/ROW]
[ROW][C]B[/C][C]0.76[/C][C]1[/C][C]0.755[/C][C]0.593[/C][/ROW]
[ROW][C]C[/C][C]0.886[/C][C]0.755[/C][C]1[/C][C]0.705[/C][/ROW]
[ROW][C]D
[/C][C]0.613[/C][C]0.593[/C][C]0.705[/C][C]1[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=160107&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=160107&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)
	A	B	C	D
A	1	0.76	0.886	0.613
B	0.76	1	0.755	0.593
C	0.886	0.755	1	0.705
D	0.613	0.593	0.705	1

Correlations for all pairs of data series with p-values
pair	Pearson r	Spearman rho	Kendall tau
A;B	0.7598	0.7242	0.5596
p-value	(0)	(0)	(0)
A;C	0.8865	0.8563	0.6891
p-value	(0)	(0)	(0)
A;D	0.6126	0.6738	0.5019
p-value	(0)	(0)	(0)
B;C	0.7546	0.7297	0.5686
p-value	(0)	(0)	(0)
B;D	0.5928	0.6306	0.4687
p-value	(0)	(0)	(0)
C;D	0.7047	0.7366	0.5573
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
A;B & 0.7598 & 0.7242 & 0.5596 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
A;C & 0.8865 & 0.8563 & 0.6891 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
A;D
 & 0.6126 & 0.6738 & 0.5019 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
B;C & 0.7546 & 0.7297 & 0.5686 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
B;D
 & 0.5928 & 0.6306 & 0.4687 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
C;D
 & 0.7047 & 0.7366 & 0.5573 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=160107&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]A;B[/C][C]0.7598[/C][C]0.7242[/C][C]0.5596[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]A;C[/C][C]0.8865[/C][C]0.8563[/C][C]0.6891[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]A;D
[/C][C]0.6126[/C][C]0.6738[/C][C]0.5019[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]B;C[/C][C]0.7546[/C][C]0.7297[/C][C]0.5686[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]B;D
[/C][C]0.5928[/C][C]0.6306[/C][C]0.4687[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]C;D
[/C][C]0.7047[/C][C]0.7366[/C][C]0.5573[/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=160107&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=160107&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
A;B	0.7598	0.7242	0.5596
p-value	(0)	(0)	(0)
A;C	0.8865	0.8563	0.6891
p-value	(0)	(0)	(0)
A;D	0.6126	0.6738	0.5019
p-value	(0)	(0)	(0)
B;C	0.7546	0.7297	0.5686
p-value	(0)	(0)	(0)
B;D	0.5928	0.6306	0.4687
p-value	(0)	(0)	(0)
C;D	0.7047	0.7366	0.5573
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