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 computationFri, 07 Dec 2012 07:16:11 -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/2012/Dec/07/t1354882598hc8albu52fjjzw8.htm/, Retrieved Fri, 29 Mar 2024 05:55:18 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=197331, Retrieved Fri, 29 Mar 2024 05:55:18 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact143
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Kendall tau Correlation Matrix] [Pearson correlati...] [2012-12-07 12:16:11] [1cff1ff8324fa30b2f172cc4e03abff6] [Current]
Feedback Forum

Post a new message
Dataseries X:
1966	1	41
1966	2	39
1966	3	50
1966	4	40
1966	5	43
1966	6	38
1966	7	44
1966	8	35
1966	9	39
1966	10	35
1966	11	29
1966	12	49
1967	1	50
1967	2	59
1967	3	63
1967	4	32
1967	5	39
1967	6	47
1967	7	53
1967	8	60
1967	9	57
1967	10	52
1967	11	70
1967	12	90
1968	1	74
1968	2	62
1968	3	55
1968	4	84
1968	5	94
1968	6	70
1968	7	108
1968	8	139
1968	9	120
1968	10	97
1968	11	126
1968	12	149
1969	1	158
1969	2	124
1969	3	140
1969	4	109
1969	5	114
1969	6	77
1969	7	120
1969	8	133
1969	9	110
1969	10	92
1969	11	97
1969	12	78
1970	1	99
1970	2	107
1970	3	112
1970	4	90
1970	5	98
1970	6	125
1970	7	155
1970	8	190
1970	9	236
1970	10	189
1970	11	174
1970	12	178
1971	1	136
1971	2	161
1971	3	171
1971	4	149
1971	5	184
1971	6	155
1971	7	276
1971	8	224
1971	9	213
1971	10	279
1971	11	268
1971	12	287
1972	1	238
1972	2	213
1972	3	257
1972	4	293
1972	5	212
1972	6	246
1972	7	353
1972	8	339
1972	9	308
1972	10	247
1972	11	257
1972	12	322
1973	1	298
1973	2	273
1973	3	312
1973	4	249
1973	5	286
1973	6	279
1973	7	309
1973	8	401
1973	9	309
1973	10	328
1973	11	353
1973	12	354
1974	1	327
1974	2	324
1974	3	285
1974	4	243
1974	5	241
1974	6	287
1974	7	355
1974	8	460
1974	9	364
1974	10	487
1974	11	452
1974	12	391
1975	1	500
1975	2	451
1975	3	375
1975	4	372
1975	5	302
1975	6	316
1975	7	398
1975	8	394
1975	9	431
1975	10	431




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

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







Correlations for all pairs of data series (method=pearson)
YearMonthRobberies
Year1-0.040.928
Month-0.0410.118
Robberies0.9280.1181

\begin{tabular}{lllllllll}
\hline
Correlations for all pairs of data series (method=pearson) \tabularnewline
  & Year & Month & Robberies \tabularnewline
Year & 1 & -0.04 & 0.928 \tabularnewline
Month & -0.04 & 1 & 0.118 \tabularnewline
Robberies & 0.928 & 0.118 & 1 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=197331&T=1

[TABLE]
[ROW][C]Correlations for all pairs of data series (method=pearson)[/C][/ROW]
[ROW][C] [/C][C]Year[/C][C]Month[/C][C]Robberies[/C][/ROW]
[ROW][C]Year[/C][C]1[/C][C]-0.04[/C][C]0.928[/C][/ROW]
[ROW][C]Month[/C][C]-0.04[/C][C]1[/C][C]0.118[/C][/ROW]
[ROW][C]Robberies[/C][C]0.928[/C][C]0.118[/C][C]1[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=197331&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=197331&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)
YearMonthRobberies
Year1-0.040.928
Month-0.0410.118
Robberies0.9280.1181







Correlations for all pairs of data series with p-values
pairPearson rSpearman rhoKendall tau
Year;Month-0.04-0.0391-0.0285
p-value(0.6671)(0.6741)(0.6726)
Year;Robberies0.92760.94570.8221
p-value(0)(0)(0)
Month;Robberies0.11820.11160.0786
p-value(0.2024)(0.229)(0.2238)

\begin{tabular}{lllllllll}
\hline
Correlations for all pairs of data series with p-values \tabularnewline
pair & Pearson r & Spearman rho & Kendall tau \tabularnewline
Year;Month & -0.04 & -0.0391 & -0.0285 \tabularnewline
p-value & (0.6671) & (0.6741) & (0.6726) \tabularnewline
Year;Robberies & 0.9276 & 0.9457 & 0.8221 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
Month;Robberies & 0.1182 & 0.1116 & 0.0786 \tabularnewline
p-value & (0.2024) & (0.229) & (0.2238) \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=197331&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]Year;Month[/C][C]-0.04[/C][C]-0.0391[/C][C]-0.0285[/C][/ROW]
[ROW][C]p-value[/C][C](0.6671)[/C][C](0.6741)[/C][C](0.6726)[/C][/ROW]
[ROW][C]Year;Robberies[/C][C]0.9276[/C][C]0.9457[/C][C]0.8221[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]Month;Robberies[/C][C]0.1182[/C][C]0.1116[/C][C]0.0786[/C][/ROW]
[ROW][C]p-value[/C][C](0.2024)[/C][C](0.229)[/C][C](0.2238)[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=197331&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=197331&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
Year;Month-0.04-0.0391-0.0285
p-value(0.6671)(0.6741)(0.6726)
Year;Robberies0.92760.94570.8221
p-value(0)(0)(0)
Month;Robberies0.11820.11160.0786
p-value(0.2024)(0.229)(0.2238)



Parameters (Session):
par1 = pearson ;
Parameters (R input):
par1 = pearson ;
R code (references can be found in the software module):
par1 <- 'pearson'
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')