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 computationTue, 20 Dec 2011 19:28:26 -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/20/t1324427338apw4h5zwfbmgo34.htm/, Retrieved Mon, 06 May 2024 08:05:30 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=158340, Retrieved Mon, 06 May 2024 08:05:30 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact118
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Kendall tau Correlation Matrix] [paper statistiek,...] [2011-12-19 15:59:25] [4b648d52023f19d55c572f0eddd72b1f]
-    D  [Kendall tau Correlation Matrix] [PAPER - DEEL 3 - ...] [2011-12-20 22:38:04] [da10aa57c5e54f8a2ad733cadd93c4c3]
-    D      [Kendall tau Correlation Matrix] [PAPER - DEEL 3 - ...] [2011-12-21 00:28:26] [e524eb56e6915a531809c7eb50783bc6] [Current]
-             [Kendall tau Correlation Matrix] [PAPER - DEEL 3 - ...] [2011-12-21 00:30:23] [da10aa57c5e54f8a2ad733cadd93c4c3]
-   P           [Kendall tau Correlation Matrix] [PAPER - DEEL 3 - ...] [2011-12-21 08:03:17] [da10aa57c5e54f8a2ad733cadd93c4c3]
Feedback Forum

Post a new message
Dataseries X:
2	210907	30	94	112285	146283	144	145
4	179321	30	103	101193	96933	135	132
0	149061	26	93	116174	95757	84	84
0	237213	38	123	66198	143983	130	127
-4	173326	44	148	71701	75851	82	78
4	133131	30	90	57793	59238	60	60
4	258873	40	124	80444	93163	131	131
0	324799	47	168	97668	151511	140	133
-1	230964	30	115	133824	136368	151	150
0	236785	31	71	101481	112642	91	91
1	344297	30	108	67654	127766	119	118
0	174724	34	120	69112	85646	123	119
3	174415	31	114	82753	98579	90	89
-1	223632	33	120	72654	131741	113	108
4	294424	33	124	101494	171975	175	162
3	325107	36	126	79215	159676	96	92
1	106408	14	37	31081	58391	41	41
0	96560	17	38	22996	31580	47	47
-2	265769	32	120	83122	136815	126	120
-3	269651	30	93	70106	120642	105	105
-4	149112	35	95	60578	69107	80	79
2	152871	28	90	79892	108016	73	70
2	362301	34	110	100708	79336	68	67
-4	183167	39	138	82875	93176	127	127
3	277965	39	133	139077	161632	154	152
2	218946	29	96	80670	102996	112	109
2	244052	44	164	143558	160604	137	133
0	341570	21	78	117105	158051	135	123
5	233328	28	102	120733	162647	230	230
-2	206161	28	99	73107	60622	71	68
0	311473	38	129	132068	179566	147	147
-2	207176	32	114	87011	96144	105	101
-3	196553	29	99	95260	129847	107	108
2	143246	27	104	106671	71180	116	114
2	182192	40	138	70054	86767	89	88
2	194979	40	151	74011	93487	84	83
0	167488	28	72	83737	82981	113	113
4	143756	34	120	69094	73815	120	118
4	275541	33	115	93133	94552	110	110
2	152299	33	98	61370	67808	78	76
2	193339	35	71	84651	106175	145	141
-4	130585	29	107	95364	76669	91	91
3	112611	20	73	26706	57283	48	48
3	148446	37	129	126846	72413	150	144
2	182079	33	118	102860	96971	181	168
-1	243060	29	104	111813	120336	121	117
-3	162765	28	107	120293	93913	99	100
0	85574	21	36	24266	32036	40	37
1	225060	41	139	109825	102255	87	87
-3	133328	20	56	40909	63506	66	64
3	100750	30	93	140867	68370	58	58
0	101523	22	87	61056	50517	77	76
0	243511	42	110	101338	103950	130	129
0	152474	32	83	65567	84396	101	101
3	132487	36	98	40735	55515	120	89
-3	317394	31	82	91413	209056	195	193
0	244749	33	115	76643	142775	106	101
-4	184510	40	140	110681	68847	83	82
2	128423	38	120	92696	20112	37	36
-1	97839	24	66	94785	61023	77	75
3	172494	43	139	86687	112494	144	131
2	229242	31	119	91721	78876	95	90
5	351619	40	141	115168	170745	169	166
2	324598	37	133	135777	122037	134	133
-2	195838	31	98	102372	112283	197	196
0	254488	39	117	103772	120691	140	136
3	199476	32	105	135400	122422	125	123
-2	92499	18	55	21399	25899	21	21
0	224330	39	132	130115	139296	167	163
6	181633	30	73	64466	89455	96	96
-3	271856	37	86	54990	147866	151	151
3	95227	32	48	34777	14336	23	23
0	98146	17	48	27114	30059	21	14
-2	118612	12	43	30080	41907	90	87
1	65475	13	46	69008	35885	60	56
0	108446	17	65	46300	55764	26	25
2	121848	17	52	30594	35619	41	41
2	76302	20	68	30976	40557	35	33
-3	98104	17	47	25568	44197	68	68
-2	30989	17	41	4154	4103	6	6
1	31774	17	47	4143	4694	0	0
-4	150580	22	71	45588	62991	41	39
0	54157	15	30	18625	24261	38	37
1	59382	12	24	26263	21425	47	47
0	84105	17	63	20055	27184	34	34




\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 & 'AstonUniversity' @ aston.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 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=158340&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]'AstonUniversity' @ aston.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] [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=158340&T=0

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



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