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, 22 Dec 2011 13:54:45 -0500
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2011/Dec/22/t13245801908lc4yrwcz1dtkdr.htm/, Retrieved Fri, 03 May 2024 06:09:03 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=159862, Retrieved Fri, 03 May 2024 06:09:03 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact111
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Univariate Data Series] [Aantal ongevallen...] [2011-12-21 22:36:19] [d67ce207bd02ca41b9162077ae11c874]
- RMPD    [Kendall tau Correlation Matrix] [Pearson Correlation] [2011-12-22 18:54:45] [e569a00cc6e8044e6afea1f18dd335a0] [Current]
Feedback Forum

Post a new message
Dataseries X:
30	112285	145	56
28	84786	101	56
38	83123	98	54
30	101193	132	89
22	38361	60	40
26	68504	38	25
25	119182	144	92
18	22807	5	18
11	17140	28	63
26	116174	84	44
25	57635	79	33
38	66198	127	84
44	71701	78	88
30	57793	60	55
40	80444	131	60
34	53855	84	66
47	97668	133	154
30	133824	150	53
31	101481	91	119
23	99645	132	41
36	114789	136	61
36	99052	124	58
30	67654	118	75
25	65553	70	33
39	97500	107	40
34	69112	119	92
31	82753	89	100
31	85323	112	112
33	72654	108	73
25	30727	52	40
33	77873	112	45
35	117478	116	60
42	74007	123	62
43	90183	125	75
30	61542	27	31
33	101494	162	77
13	27570	32	34
32	55813	64	46
36	79215	92	99
0	1423	0	17
28	55461	83	66
14	31081	41	30
17	22996	47	76
32	83122	120	146
30	70106	105	67
35	60578	79	56
20	39992	65	107
28	79892	70	58
28	49810	55	34
39	71570	39	61
34	100708	67	119
26	33032	21	42
39	82875	127	66
39	139077	152	89
33	71595	113	44
28	72260	99	66
4	5950	7	24
39	115762	141	259
18	32551	21	17
14	31701	35	64
29	80670	109	41
44	143558	133	68
21	117105	123	168
16	23789	26	43
28	120733	230	132
35	105195	166	105
28	73107	68	71
38	132068	147	112
23	149193	179	94
36	46821	61	82
32	87011	101	70
29	95260	108	57
25	55183	90	53
27	106671	114	103
36	73511	103	121
28	92945	142	62
23	78664	79	52
40	70054	88	52
23	22618	25	32
40	74011	83	62
28	83737	113	45
34	69094	118	46
33	93133	110	63
28	95536	129	75
34	225920	51	88
30	62133	93	46
33	61370	76	53
22	43836	49	37
38	106117	118	90
26	38692	38	63
35	84651	141	78
8	56622	58	25
24	15986	27	45
29	95364	91	46
20	26706	48	41
29	89691	63	144
45	67267	56	82
37	126846	144	91
33	41140	73	71
33	102860	168	63
25	51715	64	53
32	55801	97	62
29	111813	117	63
28	120293	100	32
28	138599	149	39
31	161647	187	62
52	115929	127	117
21	24266	37	34
24	162901	245	92
41	109825	87	93
33	129838	177	54
32	37510	49	144
19	43750	49	14
20	40652	73	61
31	87771	177	109
31	85872	94	38
32	89275	117	73
18	44418	60	75
23	192565	55	50
17	35232	39	61
20	40909	64	55
12	13294	26	77
17	32387	64	75
30	140867	58	72
31	120662	95	50
10	21233	25	32
13	44332	26	53
22	61056	76	42
42	101338	129	71
1	1168	11	10
9	13497	2	35
32	65567	101	65
11	25162	28	25
25	32334	36	66
36	40735	89	41
31	91413	193	86
0	855	4	16
24	97068	84	42
13	44339	23	19
8	14116	39	19
13	10288	14	45
19	65622	78	65
18	16563	14	35
33	76643	101	95
40	110681	82	49
22	29011	24	37
38	92696	36	64
24	94785	75	38
8	8773	16	34
35	83209	55	32
43	93815	131	65
43	86687	131	52
14	34553	39	62
41	105547	144	65
38	103487	139	83
45	213688	211	95
31	71220	78	29
13	23517	50	18
28	56926	39	33
31	91721	90	247
40	115168	166	139
30	111194	12	29
16	51009	57	118
37	135777	133	110
30	51513	69	67
35	74163	119	42
32	51633	119	65
27	75345	65	94
20	33416	61	64
18	83305	49	81
31	98952	101	95
31	102372	196	67
21	37238	15	63
39	103772	136	83
41	123969	89	45
13	27142	40	30
32	135400	123	70
18	21399	21	32
39	130115	163	83
14	24874	29	31
7	34988	35	67
17	45549	13	66
0	6023	5	10
30	64466	96	70
37	54990	151	103
0	1644	6	5
5	6179	13	20
1	3926	3	5
16	32755	56	36
32	34777	23	34
24	73224	57	48
17	27114	14	40
11	20760	43	43
24	37636	20	31
22	65461	72	42
12	30080	87	46
19	24094	21	33
13	69008	56	18
17	54968	59	55
15	46090	82	35
16	27507	43	59
24	10672	25	19
15	34029	38	66
17	46300	25	60
18	24760	38	36
20	18779	12	25
16	21280	29	47
16	40662	47	54
18	28987	45	53
22	22827	40	40
8	18513	30	40
17	30594	41	39
18	24006	25	14
16	27913	23	45
23	42744	14	36
22	12934	16	28
13	22574	26	44
13	41385	21	30
16	18653	27	22
16	18472	9	17
20	30976	33	31
22	63339	42	55
17	25568	68	54
18	33747	32	21
17	4154	6	14
12	19474	67	81
7	35130	33	35
17	39067	77	43
14	13310	46	46
23	65892	30	30
17	4143	0	23
14	28579	36	38
15	51776	46	54
17	21152	18	20
21	38084	48	53
18	27717	29	45
18	32928	28	39
17	11342	34	20
17	19499	33	24
16	16380	34	31
15	36874	33	35
21	48259	80	151
16	16734	32	52
14	28207	30	30
15	30143	41	31
17	41369	41	29
15	45833	51	57
15	29156	18	40
10	35944	34	44
6	36278	31	25
22	45588	39	77
21	45097	54	35
1	3895	14	11
18	28394	24	63
17	18632	24	44
4	2325	8	19
10	25139	26	13
16	27975	19	42
16	14483	11	38
9	13127	14	29
16	5839	1	20
17	24069	39	27
7	3738	5	20
15	18625	37	19
14	36341	32	37
14	24548	38	26
18	21792	47	42
12	26263	47	49
16	23686	37	30
21	49303	51	49
19	25659	45	67
16	28904	21	28
1	2781	1	19
16	29236	42	49
10	19546	26	27
19	22818	21	30
12	32689	4	22
2	5752	10	12
14	22197	43	31
17	20055	34	20
19	25272	31	20
14	82206	19	39
11	32073	34	29
4	5444	6	16
16	20154	11	27
20	36944	24	21
12	8019	16	19
15	30884	72	35
16	19540	21	14




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time2 seconds
R Server'Gwilym Jenkins' @ jenkins.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 & 'Gwilym Jenkins' @ jenkins.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=159862&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' @ jenkins.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=159862&T=0

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







Correlations for all pairs of data series (method=pearson)
compendiums_reviewedtotsizetotblogslogins
compendiums_reviewed10.7340.720.549
totsize0.73410.7980.536
totblogs0.720.79810.581
logins0.5490.5360.5811

\begin{tabular}{lllllllll}
\hline
Correlations for all pairs of data series (method=pearson) \tabularnewline
  & compendiums_reviewed & totsize & totblogs & logins \tabularnewline
compendiums_reviewed & 1 & 0.734 & 0.72 & 0.549 \tabularnewline
totsize & 0.734 & 1 & 0.798 & 0.536 \tabularnewline
totblogs & 0.72 & 0.798 & 1 & 0.581 \tabularnewline
logins & 0.549 & 0.536 & 0.581 & 1 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=159862&T=1

[TABLE]
[ROW][C]Correlations for all pairs of data series (method=pearson)[/C][/ROW]
[ROW][C] [/C][C]compendiums_reviewed[/C][C]totsize[/C][C]totblogs[/C][C]logins[/C][/ROW]
[ROW][C]compendiums_reviewed[/C][C]1[/C][C]0.734[/C][C]0.72[/C][C]0.549[/C][/ROW]
[ROW][C]totsize[/C][C]0.734[/C][C]1[/C][C]0.798[/C][C]0.536[/C][/ROW]
[ROW][C]totblogs[/C][C]0.72[/C][C]0.798[/C][C]1[/C][C]0.581[/C][/ROW]
[ROW][C]logins[/C][C]0.549[/C][C]0.536[/C][C]0.581[/C][C]1[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=159862&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=159862&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)
compendiums_reviewedtotsizetotblogslogins
compendiums_reviewed10.7340.720.549
totsize0.73410.7980.536
totblogs0.720.79810.581
logins0.5490.5360.5811







Correlations for all pairs of data series with p-values
pairPearson rSpearman rhoKendall tau
compendiums_reviewed;totsize0.73450.78990.5825
p-value(0)(0)(0)
compendiums_reviewed;totblogs0.71970.75580.5563
p-value(0)(0)(0)
compendiums_reviewed;logins0.54880.62840.4525
p-value(0)(0)(0)
totsize;totblogs0.79830.83940.6562
p-value(0)(0)(0)
totsize;logins0.53640.65880.4742
p-value(0)(0)(0)
totblogs;logins0.5810.70780.5116
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
compendiums_reviewed;totsize & 0.7345 & 0.7899 & 0.5825 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
compendiums_reviewed;totblogs & 0.7197 & 0.7558 & 0.5563 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
compendiums_reviewed;logins & 0.5488 & 0.6284 & 0.4525 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
totsize;totblogs & 0.7983 & 0.8394 & 0.6562 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
totsize;logins & 0.5364 & 0.6588 & 0.4742 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
totblogs;logins & 0.581 & 0.7078 & 0.5116 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=159862&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]compendiums_reviewed;totsize[/C][C]0.7345[/C][C]0.7899[/C][C]0.5825[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]compendiums_reviewed;totblogs[/C][C]0.7197[/C][C]0.7558[/C][C]0.5563[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]compendiums_reviewed;logins[/C][C]0.5488[/C][C]0.6284[/C][C]0.4525[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]totsize;totblogs[/C][C]0.7983[/C][C]0.8394[/C][C]0.6562[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]totsize;logins[/C][C]0.5364[/C][C]0.6588[/C][C]0.4742[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]totblogs;logins[/C][C]0.581[/C][C]0.7078[/C][C]0.5116[/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=159862&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=159862&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
compendiums_reviewed;totsize0.73450.78990.5825
p-value(0)(0)(0)
compendiums_reviewed;totblogs0.71970.75580.5563
p-value(0)(0)(0)
compendiums_reviewed;logins0.54880.62840.4525
p-value(0)(0)(0)
totsize;totblogs0.79830.83940.6562
p-value(0)(0)(0)
totsize;logins0.53640.65880.4742
p-value(0)(0)(0)
totblogs;logins0.5810.70780.5116
p-value(0)(0)(0)



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