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, 13 Dec 2011 12:21:22 -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/13/t1323796954hpg9oforo0qsezc.htm/, Retrieved Thu, 31 Oct 2024 23:13:15 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=154560, Retrieved Thu, 31 Oct 2024 23:13:15 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact107
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Kendall tau Correlation Matrix] [workshop 10: (1)] [2011-12-13 17:21:22] [d7127d50f40450f0f3837a0965e389eb] [Current]
Feedback Forum

Post a new message
Dataseries X:
2050	2650	13	7	1	0	1639
2080	2600	*	4	1	0	1088
2150	2664	6	5	1	0	1193
2150	2921	3	6	1	0	1635
1999	2580	4	4	1	0	1732
1900	2580	4	4	1	0	1534
1800	2774	2	4	1	0	1765
1560	1920	1	5	1	0	1161
1450	2150	*	4	1	0	*
1449	1710	1	3	1	0	1010
1375	1837	4	5	1	0	1191
1270	1880	8	6	1	0	930
1250	2150	15	3	1	0	984
1235	1894	14	5	1	0	1112
1170	1928	18	8	1	0	600
1180	1830	*	3	1	0	733
1155	1767	16	4	1	0	794
1110	1630	15	3	1	1	867
1139	1680	17	4	1	1	750
995	1725	*	3	1	0	923
995	1500	15	4	1	0	743
975	1430	*	3	1	0	752
975	1360	*	4	1	0	696
900	1400	16	2	1	1	731
960	1573	17	6	1	0	768
860	1385	*	2	1	0	653
1695	2931	28	3	1	1	1142
1553	2200	28	4	1	0	1035
1250	2277	*	4	1	0	*
1300	2000	*	3	1	0	1076
1020	1478	53	3	1	1	626
1020	1713	30	4	1	1	600
922	1326	*	4	1	0	668
925	1050	*	2	1	1	553
899	1464	*	2	1	0	566
850	1190	41	1	1	0	600
876	1156	*	1	1	0	*
890	1746	*	2	1	0	591
870	1280	*	1	1	0	599
700	1215	*	3	1	0	477
720	1121	46	4	1	0	398
720	1050	*	1	1	0	*
749	1733	43	6	1	0	656
731	1299	*	6	1	0	585
725	1140	*	3	1	1	490
670	1181	*	4	1	0	440
2150	2848	4	6	1	0	1487
1599	2440	*	5	1	0	1265
1350	2253	23	4	1	0	939
1299	2743	25	5	1	1	1232
1250	2180	17	4	1	1	1141
1239	1706	14	4	1	0	810
1200	1948	*	4	1	0	899
1125	1710	16	4	1	0	800
1100	1657	*	4	1	0	865
1080	2200	26	4	1	0	1076
1050	1680	13	4	1	0	875
1049	1900	34	3	1	0	690
955	1565	*	3	1	0	648
934	1543	20	3	1	0	820
875	1173	6	4	1	0	456
889	1549	*	4	1	0	723
855	1900	*	3	1	0	780
835	1560	*	5	1	1	638
810	1365	*	2	1	0	673
805	1258	7	4	1	1	821
799	1314	*	2	1	0	671
750	1338	*	3	1	1	649
759	997	4	4	1	0	461
755	1275	*	5	1	0	*
750	1030	*	1	1	0	486
730	1027	*	3	1	0	427
729	1007	19	6	1	0	513
710	1083	22	4	1	0	504
773	1320	*	5	1	0	*
690	1348	15	2	1	0	*
670	1350	*	2	1	0	622
619	837	*	2	1	0	342
1295	3750	*	4	0	1	1200
975	1500	7	3	0	1	700
939	1428	40	2	0	0	701
820	1375	*	1	0	0	585
780	1080	*	3	0	0	600
770	900	*	3	0	0	391
700	1505	*	2	0	1	591
620	1480	*	4	0	0	*
540	1142	*	0	0	0	223
1070	1464	*	2	0	0	376
2100	2116	25	3	0	0	1209
725	1280	*	3	0	0	447
660	1159	*	0	0	0	225
600	1198	*	4	0	0	*
580	1051	15	2	0	0	426
1844	2250	40	6	0	0	915
1580	2563	*	2	0	0	1189
699	1400	45	1	0	1	481
1330	1850	5	5	0	1	*
1160	1720	5	4	0	0	867
1109	1740	4	3	0	0	816
1129	1700	6	4	0	0	725
1050	1620	6	4	0	0	800
1045	1630	6	4	0	0	750
1050	1920	8	4	0	0	944
1020	1606	5	4	0	0	811
1000	1535	7	5	0	1	668
1030	1540	6	2	0	1	826
975	1739	13	3	0	0	880
950	1715	*	3	0	0	900
940	1305	5	3	0	0	647
920	1415	7	4	0	0	866
945	1580	9	3	0	0	810
874	1236	3	4	0	0	707
872	1229	6	3	0	0	721
870	1273	4	4	0	0	638
869	1165	7	4	0	0	694
766	1200	7	4	0	1	634
739	970	4	4	0	1	541




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time0 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 & 0 seconds \tabularnewline
R Server & 'Herman Ole Andreas Wold' @ wold.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=154560&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]0 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=154560&T=0

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



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