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 15:54:10 -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/t13245872670mkgkx4jexebvy5.htm/, Retrieved Fri, 03 May 2024 07:08:57 +0000

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

-     [Multiple Regression] [] [2010-12-05 18:56:24] [b98453cac15ba1066b407e146608df68]
- RMPD    [Kendall tau Correlation Matrix] [misser] [2011-12-22 20:54:10] [5c55e7d277583a4a66c326a86fdb470e] [Current]
- R         [Kendall tau Correlation Matrix] [error] [2011-12-22 21:04:19] [9a01a2fa0c801b8d05d8c9909bdf1c0f] 

Feedback Forum

Post a new message

Dataseries X:

Download CSV

Histogram

Boxplots

1845	162687	95	595	115	0	48	21
1797	201906	63	545	76	1	58	20
192	7215	18	72	1	0	0	0
2444	146367	97	679	155	0	67	27
3567	257045	139	1201	125	0	83	31
6953	528532	268	1975	278	1	137	36
1873	191582	60	602	93	1	65	26
1740	195674	60	496	59	0	86	30
2079	177020	45	670	87	0	62	30
3120	330255	100	1047	130	1	72	27
1946	121844	75	634	158	2	50	24
2370	203938	72	743	120	0	88	30
1956	116374	107	689	87	0	62	22
3198	220751	120	1086	264	4	79	28
1496	173259	65	420	51	4	56	18
1574	156326	89	474	85	3	54	22
1808	145178	59	442	100	0	81	37
1309	89171	61	373	72	5	13	15
2820	172624	88	899	147	0	74	34
799	43391	28	253	49	0	18	18
1162	87927	62	399	40	0	31	15
2818	241285	103	850	99	0	99	30
1780	200429	75	648	127	1	38	25
2316	146946	57	717	164	1	59	34
1994	159763	89	619	41	1	54	21
1806	207078	34	657	160	0	63	21
2153	212394	167	691	92	0	66	25
1458	201536	96	366	59	0	90	31
3000	394662	121	994	89	0	72	31
2236	217892	46	929	90	0	61	20
1685	182286	45	490	76	0	61	28
1667	188748	48	574	116	2	63	22
2257	137978	107	738	92	4	53	17
3402	259627	132	1039	361	0	120	25
2571	236489	55	844	85	1	73	25
1	0	1	0	0	0	0	0
2143	230761	65	1000	63	0	54	31
1878	132807	54	629	138	3	54	14
2197	158599	51	533	270	9	46	35
2186	253254	68	811	64	0	83	34
2533	269329	72	837	96	2	106	22
1823	161273	61	682	62	0	44	34
1095	107181	33	400	35	2	27	23
2303	213097	81	831	66	1	73	24
1365	139667	51	419	56	2	71	26
1245	171101	99	334	41	2	44	23
756	81407	33	216	49	1	23	35
2419	247596	106	786	121	0	78	24
2327	239807	90	752	113	1	60	31
2787	172743	60	964	190	8	73	30
658	48188	28	205	37	0	12	22
2013	169355	71	506	52	0	104	23
2616	325322	77	830	89	0	86	27
2072	241518	80	694	73	0	57	30
1912	195583	60	691	49	1	67	33
1775	159913	57	547	77	8	44	12
1943	223936	71	547	58	0	53	26
1047	101694	26	329	75	1	26	26
1190	157258	68	427	32	0	67	23
2891	202536	100	972	59	10	36	38
1868	181076	66	564	71	6	56	32
2255	150518	84	836	91	0	52	21
1392	141491	64	376	87	11	54	22
1355	130108	40	471	48	3	61	26
1318	166049	37	430	63	0	27	28
1526	124197	43	483	41	0	58	33
2335	195043	71	504	86	8	76	36
2898	138708	66	887	152	2	93	25
1118	116552	40	271	49	0	59	25
340	31970	15	101	40	0	5	21
2977	258158	115	1097	135	3	57	19
1452	151194	79	470	83	1	42	12
1550	135926	68	528	62	2	88	30
1685	119629	73	475	91	1	53	21
2728	171518	71	698	95	0	81	39
1574	108949	45	425	82	2	35	32
2413	183471	60	709	112	1	102	28
2563	159966	98	824	70	0	71	29
1080	93786	35	336	78	0	28	21
1235	84971	72	395	105	0	34	31
980	88882	76	234	49	0	54	26
2246	304603	65	830	60	0	49	29
1076	75101	30	334	49	1	30	23
1638	145043	41	524	132	0	57	25
1208	95827	48	393	49	0	54	22
1866	173924	59	574	71	0	38	26
2727	241957	238	672	102	0	63	33
1208	115367	115	284	74	0	58	24
1427	118689	65	452	49	7	46	24
1610	164078	54	653	74	0	46	21
1865	158931	42	684	59	5	51	28
2413	184139	83	706	91	1	87	28
1238	152856	58	417	68	0	39	25
1468	146159	61	551	81	0	28	15
974	62535	43	394	33	0	26	13
2319	245196	117	730	166	0	52	36
1890	199841	71	571	97	0	96	27
223	19349	12	67	15	0	13	1
2527	247280	109	877	105	3	43	24
2105	164457	86	865	61	0	42	31
778	72128	30	306	11	0	30	4
1194	104253	26	382	45	0	59	21
1424	151090	57	435	89	0	73	27
1386	147990	68	348	72	1	40	26
839	87448	42	227	27	1	36	12
596	27676	22	194	59	0	2	16
1684	170326	52	413	127	0	103	29
1168	132148	38	273	48	1	30	26
0	0	0	0	0	0	0	0
1106	95778	34	343	58	0	46	25
1149	109001	68	376	57	0	25	21
1485	158833	46	495	60	0	59	24
1529	150013	66	448	77	1	60	21
962	89887	63	313	71	0	36	21
78	3616	5	14	5	0	0	0
0	0	0	0	0	0	0	0
1184	199005	45	410	70	0	45	23
1672	160930	93	606	76	0	79	33
2142	177948	102	593	124	2	30	32
1016	136061	40	312	56	0	43	23
778	43410	19	292	63	0	7	1
1857	184277	75	547	92	1	80	29
1084	109873	45	315	58	0	32	20
2297	151030	59	660	64	8	84	33
731	60493	40	174	29	3	3	12
285	19764	12	75	19	1	10	2
1872	177559	56	572	64	3	47	21
1181	140281	35	389	79	0	35	28
1725	164249	54	562	104	0	54	35
256	11796	9	79	22	0	1	2
98	10674	9	33	7	0	0	0
1435	151322	59	487	37	0	46	18
41	6836	3	11	5	0	0	1
1931	174712	68	664	48	6	51	21
42	5118	3	6	1	0	5	0
528	40248	16	183	34	1	8	4
0	0	0	0	0	0	0	0
1122	127628	51	342	53	0	38	29
1305	88837	38	269	44	0	21	26
81	7131	4	27	0	1	0	0
262	9056	15	99	18	0	0	4
1104	88589	29	306	52	1	18	19
1290	144470	53	327	56	0	53	22
1248	111408	20	459	50	1	17	22

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
R Engine error message	Error in dimnames(data) <- dimnames : length of 'dimnames' [1] not equal to array extent Calls: array Execution halted

\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
R Engine error message & Error in dimnames(data) <- dimnames : 
  length of 'dimnames' [1] not equal to array extent
Calls: array
Execution halted
 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=159972&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]
[ROW][C]R Engine error message[/C][C]Error in dimnames(data) <- dimnames : 
  length of 'dimnames' [1] not equal to array extent
Calls: array
Execution halted
[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=159972&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=159972&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
R Engine error message	Error in dimnames(data) <- dimnames : length of 'dimnames' [1] not equal to array extent Calls: array Execution halted

Parameters (Session):

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