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

Sun, 11 Dec 2011 07:16:49 -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/11/t1323606197ho2odoy14n2gyhw.htm/, Retrieved Sun, 28 Apr 2024 22:18:14 +0000

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=153689, Retrieved Sun, 28 Apr 2024 22:18:14 +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

114

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] [kendall tau corre...] [2011-12-11 12:16:49] [f007bbc48ca3190e286f441a6cce1887] [Current]

Feedback Forum

Post a new message

Dataseries X:

Download CSV

Histogram

Boxplots

1418	210907	79	2
869	120982	58	0
1530	176508	60	0
2172	179321	108	4
901	123185	49	0
463	52746	0	0
3201	385534	121	0
371	33170	1	0
1192	101645	20	0
1583	149061	43	0
1439	165446	69	0
1764	237213	78	0
1495	173326	86	-4
1373	133131	44	4
2187	258873	104	4
1491	180083	63	0
4041	324799	158	0
1706	230964	102	-1
2152	236785	77	0
1036	135473	82	0
1882	202925	115	0
1929	215147	101	0
2242	344297	80	1
1220	153935	50	0
1289	132943	83	0
2515	174724	123	0
2147	174415	73	3
2352	225548	81	0
1638	223632	105	-1
1222	124817	47	0
1812	221698	105	0
1677	210767	94	0
1579	170266	44	0
1731	260561	114	0
807	84853	38	0
2452	294424	107	4
829	101011	30	0
1940	215641	71	0
2662	325107	84	3
186	7176	0	0
1499	167542	59	0
865	106408	33	1
1793	96560	42	0
2527	265769	96	-2
2747	269651	106	-3
1324	149112	56	-4
2702	175824	57	0
1383	152871	59	2
1179	111665	39	0
2099	116408	34	0
4308	362301	76	2
918	78800	20	0
1831	183167	91	-4
3373	277965	115	3
1713	150629	85	0
1438	168809	76	0
496	24188	8	0
2253	329267	79	0
744	65029	21	0
1161	101097	30	0
2352	218946	76	2
2144	244052	101	2
4691	341570	94	0
1112	103597	27	0
2694	233328	92	5
1973	256462	123	0
1769	206161	75	-2
3148	311473	128	0
2474	235800	105	0
2084	177939	55	0
1954	207176	56	-2
1226	196553	41	-3
1389	174184	72	0
1496	143246	67	2
2269	187559	75	0
1833	187681	114	0
1268	119016	118	0
1943	182192	77	2
893	73566	22	0
1762	194979	66	2
1403	167488	69	0
1425	143756	105	4
1857	275541	116	4
1840	243199	88	0
1502	182999	73	0
1441	135649	99	0
1420	152299	62	2
1416	120221	53	0
2970	346485	118	0
1317	145790	30	0
1644	193339	100	2
870	80953	49	0
1654	122774	24	0
1054	130585	67	-4
937	112611	46	3
3004	286468	57	0
2008	241066	75	0
2547	148446	135	3
1885	204713	68	0
1626	182079	124	2
1468	140344	33	0
2445	220516	98	0
1964	243060	58	-1
1381	162765	68	-3
1369	182613	81	0
1659	232138	131	0
2888	265318	110	0
1290	85574	37	0
2845	310839	130	0
1982	225060	93	1
1904	232317	118	0
1391	144966	39	0
602	43287	13	0
1743	155754	74	0
1559	164709	81	0
2014	201940	109	0
2143	235454	151	0
2146	220801	51	0
874	99466	28	0
1590	92661	40	0
1590	133328	56	-3
1210	61361	27	0
2072	125930	37	0
1281	100750	83	3
1401	224549	54	0
834	82316	27	0
1105	102010	28	0
1272	101523	59	0
1944	243511	133	0
391	22938	12	0
761	41566	0	0
1605	152474	106	0
530	61857	23	0
1988	99923	44	0
1386	132487	71	3
2395	317394	116	-3
387	21054	4	0
1742	209641	62	0
620	22648	12	0
449	31414	18	0
800	46698	14	0
1684	131698	60	0
1050	91735	7	0
2699	244749	98	0
1606	184510	64	-4
1502	79863	29	0
1204	128423	32	2
1138	97839	25	-1
568	38214	16	0
1459	151101	48	0
2158	272458	100	0
1111	172494	46	3
1421	108043	45	0
2833	328107	129	0
1955	250579	130	0
2922	351067	136	0
1002	158015	59	0
1060	98866	25	0
956	85439	32	0
2186	229242	63	2
3604	351619	95	5
1035	84207	14	0
1417	120445	36	0
3261	324598	113	2
1587	131069	47	0
1424	204271	92	0
1701	165543	70	0
1249	141722	19	0
946	116048	50	0
1926	250047	41	0
3352	299775	91	0
1641	195838	111	-2
2035	173260	41	0
2312	254488	120	0
1369	104389	135	0
1577	136084	27	0
2201	199476	87	3
961	92499	25	-2
1900	224330	131	0
1254	135781	45	0
1335	74408	29	0
1597	81240	58	0
207	14688	4	0
1645	181633	47	6
2429	271856	109	-3
151	7199	7	0
474	46660	12	0
141	17547	0	0
1639	133368	37	0
872	95227	37	3
1318	152601	46	0
1018	98146	15	0
1383	79619	42	0
1314	59194	7	0
1335	139942	54	0
1403	118612	54	-2
910	72880	14	0
616	65475	16	1
1407	99643	33	0
771	71965	32	0
766	77272	21	0
473	49289	15	0
1376	135131	38	0
1232	108446	22	0
1521	89746	28	0
572	44296	10	0
1059	77648	31	0
1544	181528	32	0
1230	134019	32	0
1206	124064	43	0
1205	92630	27	0
1255	121848	37	2
613	52915	20	0
721	81872	32	0
1109	58981	0	0
740	53515	5	0
1126	60812	26	0
728	56375	10	0
689	65490	27	0
592	80949	11	0
995	76302	29	2
1613	104011	25	0
2048	98104	55	-3
705	67989	23	0
301	30989	5	-2
1803	135458	43	0
799	73504	23	0
861	63123	34	0
1186	61254	36	0
1451	74914	35	0
628	31774	0	1
1161	81437	37	0
1463	87186	28	0
742	50090	16	0
979	65745	26	0
675	56653	38	0
1241	158399	23	0
676	46455	22	0
1049	73624	30	0
620	38395	16	0
1081	91899	18	0
1688	139526	28	0
736	52164	32	0
617	51567	21	0
812	70551	23	0
1051	84856	29	0
1656	102538	50	0
705	86678	12	0
945	85709	21	0
554	34662	18	0
1597	150580	27	-4
982	99611	41	0
222	19349	13	0
1212	99373	12	0
1143	86230	21	0
435	30837	8	0
532	31706	26	0
882	89806	27	0
608	62088	13	0
459	40151	16	0
578	27634	2	0
826	76990	42	0
509	37460	5	0
717	54157	37	0
637	49862	17	0
857	84337	38	0
830	64175	37	0
652	59382	29	1
707	119308	32	0
954	76702	35	0
1461	103425	17	0
672	70344	20	0
778	43410	7	0
1141	104838	46	0
680	62215	24	0
1090	69304	40	0
616	53117	3	0
285	19764	10	0
1145	86680	37	0
733	84105	17	0
888	77945	28	0
849	89113	19	0
1182	91005	29	0
528	40248	8	0
642	64187	10	0
947	50857	15	0
819	56613	15	0
757	62792	28	0
894	72535	17	0

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	'Herman Ole Andreas Wold' @ wold.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 & 'Herman Ole Andreas Wold' @ wold.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=153689&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]'Herman Ole Andreas Wold' @ wold.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=153689&T=0

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

Correlations for all pairs of data series (method=pearson)
	page	time	blog	score
page	1	0.893	0.761	0.147
time	0.893	1	0.831	0.124
blog	0.761	0.831	1	0.122
score	0.147	0.124	0.122	1

\begin{tabular}{lllllllll}
\hline
Correlations for all pairs of data series (method=pearson) \tabularnewline
  & page & time & blog & score \tabularnewline
page & 1 & 0.893 & 0.761 & 0.147 \tabularnewline
time & 0.893 & 1 & 0.831 & 0.124 \tabularnewline
blog & 0.761 & 0.831 & 1 & 0.122 \tabularnewline
score & 0.147 & 0.124 & 0.122 & 1 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=153689&T=1

[TABLE]
[ROW][C]Correlations for all pairs of data series (method=pearson)[/C][/ROW]
[ROW][C] [/C][C]page[/C][C]time[/C][C]blog[/C][C]score[/C][/ROW]
[ROW][C]page[/C][C]1[/C][C]0.893[/C][C]0.761[/C][C]0.147[/C][/ROW]
[ROW][C]time[/C][C]0.893[/C][C]1[/C][C]0.831[/C][C]0.124[/C][/ROW]
[ROW][C]blog[/C][C]0.761[/C][C]0.831[/C][C]1[/C][C]0.122[/C][/ROW]
[ROW][C]score[/C][C]0.147[/C][C]0.124[/C][C]0.122[/C][C]1[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=153689&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=153689&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)
	page	time	blog	score
page	1	0.893	0.761	0.147
time	0.893	1	0.831	0.124
blog	0.761	0.831	1	0.122
score	0.147	0.124	0.122	1

Correlations for all pairs of data series with p-values
pair	Pearson r	Spearman rho	Kendall tau
page;time	0.8928	0.9025	0.7387
p-value	(0)	(0)	(0)
page;blog	0.7606	0.8249	0.6272
p-value	(0)	(0)	(0)
page;score	0.1472	0.0931	0.0714
p-value	(0.0123)	(0.1143)	(0.1203)
time;blog	0.8309	0.8717	0.683
p-value	(0)	(0)	(0)
time;score	0.1236	0.0948	0.0736
p-value	(0.0357)	(0.1077)	(0.1092)
blog;score	0.122	0.1148	0.0907
p-value	(0.0382)	(0.0512)	(0.0492)

\begin{tabular}{lllllllll}
\hline
Correlations for all pairs of data series with p-values \tabularnewline
pair & Pearson r & Spearman rho & Kendall tau \tabularnewline
page;time & 0.8928 & 0.9025 & 0.7387 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
page;blog & 0.7606 & 0.8249 & 0.6272 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
page;score & 0.1472 & 0.0931 & 0.0714 \tabularnewline
p-value & (0.0123) & (0.1143) & (0.1203) \tabularnewline
time;blog & 0.8309 & 0.8717 & 0.683 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
time;score & 0.1236 & 0.0948 & 0.0736 \tabularnewline
p-value & (0.0357) & (0.1077) & (0.1092) \tabularnewline
blog;score & 0.122 & 0.1148 & 0.0907 \tabularnewline
p-value & (0.0382) & (0.0512) & (0.0492) \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=153689&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]page;time[/C][C]0.8928[/C][C]0.9025[/C][C]0.7387[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]page;blog[/C][C]0.7606[/C][C]0.8249[/C][C]0.6272[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]page;score[/C][C]0.1472[/C][C]0.0931[/C][C]0.0714[/C][/ROW]
[ROW][C]p-value[/C][C](0.0123)[/C][C](0.1143)[/C][C](0.1203)[/C][/ROW]
[ROW][C]time;blog[/C][C]0.8309[/C][C]0.8717[/C][C]0.683[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]time;score[/C][C]0.1236[/C][C]0.0948[/C][C]0.0736[/C][/ROW]
[ROW][C]p-value[/C][C](0.0357)[/C][C](0.1077)[/C][C](0.1092)[/C][/ROW]
[ROW][C]blog;score[/C][C]0.122[/C][C]0.1148[/C][C]0.0907[/C][/ROW]
[ROW][C]p-value[/C][C](0.0382)[/C][C](0.0512)[/C][C](0.0492)[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=153689&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=153689&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
page;time	0.8928	0.9025	0.7387
p-value	(0)	(0)	(0)
page;blog	0.7606	0.8249	0.6272
p-value	(0)	(0)	(0)
page;score	0.1472	0.0931	0.0714
p-value	(0.0123)	(0.1143)	(0.1203)
time;blog	0.8309	0.8717	0.683
p-value	(0)	(0)	(0)
time;score	0.1236	0.0948	0.0736
p-value	(0.0357)	(0.1077)	(0.1092)
blog;score	0.122	0.1148	0.0907
p-value	(0.0382)	(0.0512)	(0.0492)

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