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 computationWed, 21 Dec 2011 07:57:31 -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/21/t1324472557nhjs6id355sn8i7.htm/, Retrieved Tue, 07 May 2024 21:11:17 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=158632, Retrieved Tue, 07 May 2024 21:11:17 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact73
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Kendall tau Correlation Matrix] [PCM] [2011-12-21 12:57:31] [b625935f05df8270d3a5abfea0142dde] [Current]
Feedback Forum

Post a new message
Dataseries X:
1773	158258	90	48	20465	6200
1704	186930	58	53	33629	10265
192	7215	18	0	1423	603
2295	129098	95	51	25629	8874
3450	230632	136	76	54002	20323
6813	508313	261	128	151036	26258
1795	180745	56	62	33287	10165
1681	185559	59	83	31172	8247
1897	154581	44	55	28113	8683
2917	290658	95	67	57803	16957
1946	121844	75	50	49830	8058
2148	184039	69	77	52143	20488
1832	100324	98	46	21055	7945
3138	215855	117	79	47007	13448
1476	168265	58	56	28735	5389
1567	154647	88	54	59147	6185
1756	142018	57	81	78950	24369
1247	79030	61	6	13497	70
2779	167047	87	74	46154	17327
726	27997	24	13	53249	3878
1048	73019	59	22	10726	3149
2805	241082	100	99	83700	20517
1760	195820	72	38	40400	2570
2261	141899	53	59	33797	5162
1848	145433	86	50	36205	5299
1665	183744	32	50	30165	7233
2082	202232	161	61	58534	15657
1440	199532	93	87	44663	15329
2741	354924	118	60	92556	14881
2112	192399	44	52	40078	16318
1684	182286	44	61	34711	9556
1616	181590	45	60	31076	10462
2227	133801	105	53	74608	7192
3088	233686	123	76	58092	4362
2389	219428	53	63	42009	14349
1	0	1	0	0	0
2099	223044	63	54	36022	10881
1669	100129	51	44	23333	8022
2095	136733	48	36	53349	13073
2153	249965	64	83	92596	26641
2390	242379	71	105	49598	14426
1701	145794	59	37	44093	15604
983	96404	32	25	84205	9184
2161	195891	78	64	63369	5989
1276	117156	50	55	60132	11270
1189	157787	94	41	37403	13958
745	81293	32	23	24460	7162
2231	224049	100	67	46456	13275
2242	223789	87	54	66616	21224
2639	160344	59	68	41554	10615
658	48188	28	12	22346	2102
1917	161922	69	99	30874	12396
2489	294283	73	74	68701	18717
2026	235223	79	56	35728	9724
1911	195583	59	67	29010	9863
1716	146061	56	40	23110	8374
1852	208834	67	53	38844	8030
981	93764	24	26	27084	7509
1177	151985	66	67	35139	14146
2809	190545	95	36	57476	7768
1688	148922	60	50	33277	13823
2097	132856	80	48	31141	7230
1309	126107	60	46	61281	10170
1244	112718	37	53	25820	7573
1256	160930	35	27	23284	5753
1293	99184	40	38	35378	9791
2303	192535	70	71	74990	19365
2897	138708	65	93	29653	9422
1103	114408	38	59	64622	12310
340	31970	15	5	4157	1283
2791	225558	112	53	29245	6372
1333	137011	71	40	50008	5413
1441	113612	68	72	52338	10837
1623	108641	71	51	13310	3394
2650	162203	67	81	92901	12964
1499	100098	44	27	10956	3495
2302	174768	60	94	34241	11580
2540	158459	97	71	75043	9970
1000	80934	30	20	21152	4911
1234	84971	71	34	42249	10138
927	80545	68	54	42005	14697
2176	287191	64	49	41152	8464
957	62974	28	26	14399	4204
1551	134091	40	48	28263	10226
1014	75555	46	35	17215	3456
1771	162154	54	32	48140	8895
2613	226638	227	55	62897	22557
1203	115019	110	58	22883	6900
1337	108749	62	44	41622	8620
1524	155537	52	45	40715	7820
1829	153133	41	49	65897	12112
2229	165618	78	72	76542	13178
1233	151517	57	39	37477	7028
1365	133686	58	28	53216	6616
950	61342	40	24	40911	9570
2319	245196	117	52	57021	14612
1857	195576	70	96	73116	11219
223	19349	12	13	3895	786
2390	225371	105	38	46609	11252
1973	152796	76	41	29351	9289
700	59117	29	24	2325	593
1062	91762	24	54	31747	6562
1311	136769	54	68	32665	8208
1157	114798	61	28	19249	7488
823	85338	40	36	15292	4574
596	27676	22	2	5842	522
1545	153535	48	91	33994	12840
1130	122417	37	29	13018	1350
0	0	0	0	0	0
1082	91529	32	46	98177	10623
1135	107205	67	25	37941	5322
1366	144664	44	51	31032	7987
1452	136540	62	59	32683	10566
870	76656	60	36	34545	1900
78	3616	5	0	0	0
0	0	0	0	0	0
1127	183065	43	40	27525	10698
1582	144677	84	68	66856	14884
2034	159104	98	28	28549	6852
919	113273	38	36	38610	6873
778	43410	19	7	2781	4
1752	175774	73	70	41211	9188
957	95401	42	30	22698	5141
1875	118893	54	59	41194	4260
731	60493	40	3	32689	443
285	19764	12	10	5752	2416
1834	164062	56	46	26757	9831
1147	132696	32	34	22527	5953
1646	155367	54	54	44810	9435
256	11796	9	1	0	0
98	10674	9	0	0	0
1404	142261	57	39	100674	7642
41	6836	3	0	0	0
1824	162563	63	48	57786	6837
42	5118	3	5	0	0
528	40248	16	8	5444	775
0	0	0	0	0	0
1073	122641	47	38	28470	8191
1305	88837	38	21	61849	1661
81	7131	4	0	0	0
261	9056	14	0	2179	548
934	76611	24	15	8019	3080
1179	132697	50	50	39644	13400
1147	100681	19	17	23494	8181




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

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







Correlations for all pairs of data series (method=pearson)
TNPTTSNLBCCWNCCWNR
TNP10.890.8270.7850.6930.696
TTS0.8910.7680.7820.690.743
NL0.8270.76810.6380.6190.628
BC0.7850.7820.63810.6530.765
CWNC0.6930.690.6190.65310.722
CWNR0.6960.7430.6280.7650.7221

\begin{tabular}{lllllllll}
\hline
Correlations for all pairs of data series (method=pearson) \tabularnewline
  & TNP & TTS & NL & BC & CWNC & CWNR \tabularnewline
TNP & 1 & 0.89 & 0.827 & 0.785 & 0.693 & 0.696 \tabularnewline
TTS & 0.89 & 1 & 0.768 & 0.782 & 0.69 & 0.743 \tabularnewline
NL & 0.827 & 0.768 & 1 & 0.638 & 0.619 & 0.628 \tabularnewline
BC & 0.785 & 0.782 & 0.638 & 1 & 0.653 & 0.765 \tabularnewline
CWNC & 0.693 & 0.69 & 0.619 & 0.653 & 1 & 0.722 \tabularnewline
CWNR & 0.696 & 0.743 & 0.628 & 0.765 & 0.722 & 1 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=158632&T=1

[TABLE]
[ROW][C]Correlations for all pairs of data series (method=pearson)[/C][/ROW]
[ROW][C] [/C][C]TNP[/C][C]TTS[/C][C]NL[/C][C]BC[/C][C]CWNC[/C][C]CWNR[/C][/ROW]
[ROW][C]TNP[/C][C]1[/C][C]0.89[/C][C]0.827[/C][C]0.785[/C][C]0.693[/C][C]0.696[/C][/ROW]
[ROW][C]TTS[/C][C]0.89[/C][C]1[/C][C]0.768[/C][C]0.782[/C][C]0.69[/C][C]0.743[/C][/ROW]
[ROW][C]NL[/C][C]0.827[/C][C]0.768[/C][C]1[/C][C]0.638[/C][C]0.619[/C][C]0.628[/C][/ROW]
[ROW][C]BC[/C][C]0.785[/C][C]0.782[/C][C]0.638[/C][C]1[/C][C]0.653[/C][C]0.765[/C][/ROW]
[ROW][C]CWNC[/C][C]0.693[/C][C]0.69[/C][C]0.619[/C][C]0.653[/C][C]1[/C][C]0.722[/C][/ROW]
[ROW][C]CWNR[/C][C]0.696[/C][C]0.743[/C][C]0.628[/C][C]0.765[/C][C]0.722[/C][C]1[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=158632&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=158632&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)
TNPTTSNLBCCWNCCWNR
TNP10.890.8270.7850.6930.696
TTS0.8910.7680.7820.690.743
NL0.8270.76810.6380.6190.628
BC0.7850.7820.63810.6530.765
CWNC0.6930.690.6190.65310.722
CWNR0.6960.7430.6280.7650.7221







Correlations for all pairs of data series with p-values
pairPearson rSpearman rhoKendall tau
TNP;TTS0.89030.86570.7029
p-value(0)(0)(0)
TNP;NL0.82750.79720.6283
p-value(0)(0)(0)
TNP;BC0.78470.76820.5939
p-value(0)(0)(0)
TNP;CWNC0.69290.63370.4742
p-value(0)(0)(0)
TNP;CWNR0.69640.66950.5026
p-value(0)(0)(0)
TTS;NL0.76840.72790.5615
p-value(0)(0)(0)
TTS;BC0.78160.76870.5899
p-value(0)(0)(0)
TTS;CWNC0.68960.62530.4719
p-value(0)(0)(0)
TTS;CWNR0.74260.71530.5436
p-value(0)(0)(0)
NL;BC0.63830.66910.5033
p-value(0)(0)(0)
NL;CWNC0.61870.60510.4601
p-value(0)(0)(0)
NL;CWNR0.62790.58260.43
p-value(0)(0)(0)
BC;CWNC0.65320.63090.473
p-value(0)(0)(0)
BC;CWNR0.76540.76040.5828
p-value(0)(0)(0)
CWNC;CWNR0.72160.70840.5493
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
TNP;TTS & 0.8903 & 0.8657 & 0.7029 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
TNP;NL & 0.8275 & 0.7972 & 0.6283 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
TNP;BC & 0.7847 & 0.7682 & 0.5939 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
TNP;CWNC & 0.6929 & 0.6337 & 0.4742 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
TNP;CWNR & 0.6964 & 0.6695 & 0.5026 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
TTS;NL & 0.7684 & 0.7279 & 0.5615 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
TTS;BC & 0.7816 & 0.7687 & 0.5899 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
TTS;CWNC & 0.6896 & 0.6253 & 0.4719 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
TTS;CWNR & 0.7426 & 0.7153 & 0.5436 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
NL;BC & 0.6383 & 0.6691 & 0.5033 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
NL;CWNC & 0.6187 & 0.6051 & 0.4601 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
NL;CWNR & 0.6279 & 0.5826 & 0.43 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
BC;CWNC & 0.6532 & 0.6309 & 0.473 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
BC;CWNR & 0.7654 & 0.7604 & 0.5828 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
CWNC;CWNR & 0.7216 & 0.7084 & 0.5493 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=158632&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]TNP;TTS[/C][C]0.8903[/C][C]0.8657[/C][C]0.7029[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]TNP;NL[/C][C]0.8275[/C][C]0.7972[/C][C]0.6283[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]TNP;BC[/C][C]0.7847[/C][C]0.7682[/C][C]0.5939[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]TNP;CWNC[/C][C]0.6929[/C][C]0.6337[/C][C]0.4742[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]TNP;CWNR[/C][C]0.6964[/C][C]0.6695[/C][C]0.5026[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]TTS;NL[/C][C]0.7684[/C][C]0.7279[/C][C]0.5615[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]TTS;BC[/C][C]0.7816[/C][C]0.7687[/C][C]0.5899[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]TTS;CWNC[/C][C]0.6896[/C][C]0.6253[/C][C]0.4719[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]TTS;CWNR[/C][C]0.7426[/C][C]0.7153[/C][C]0.5436[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]NL;BC[/C][C]0.6383[/C][C]0.6691[/C][C]0.5033[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]NL;CWNC[/C][C]0.6187[/C][C]0.6051[/C][C]0.4601[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]NL;CWNR[/C][C]0.6279[/C][C]0.5826[/C][C]0.43[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]BC;CWNC[/C][C]0.6532[/C][C]0.6309[/C][C]0.473[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]BC;CWNR[/C][C]0.7654[/C][C]0.7604[/C][C]0.5828[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]CWNC;CWNR[/C][C]0.7216[/C][C]0.7084[/C][C]0.5493[/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=158632&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=158632&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
TNP;TTS0.89030.86570.7029
p-value(0)(0)(0)
TNP;NL0.82750.79720.6283
p-value(0)(0)(0)
TNP;BC0.78470.76820.5939
p-value(0)(0)(0)
TNP;CWNC0.69290.63370.4742
p-value(0)(0)(0)
TNP;CWNR0.69640.66950.5026
p-value(0)(0)(0)
TTS;NL0.76840.72790.5615
p-value(0)(0)(0)
TTS;BC0.78160.76870.5899
p-value(0)(0)(0)
TTS;CWNC0.68960.62530.4719
p-value(0)(0)(0)
TTS;CWNR0.74260.71530.5436
p-value(0)(0)(0)
NL;BC0.63830.66910.5033
p-value(0)(0)(0)
NL;CWNC0.61870.60510.4601
p-value(0)(0)(0)
NL;CWNR0.62790.58260.43
p-value(0)(0)(0)
BC;CWNC0.65320.63090.473
p-value(0)(0)(0)
BC;CWNR0.76540.76040.5828
p-value(0)(0)(0)
CWNC;CWNR0.72160.70840.5493
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')