FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_correlation.wasp
Title produced by softwarePearson Correlation
Date of computationMon, 20 Oct 2008 16:39:37 -0600
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2008/Oct/21/t1224542428e412qy8a8ot7mpu.htm/, Retrieved Sun, 19 May 2024 17:12:38 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=18274, Retrieved Sun, 19 May 2024 17:12:38 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact160

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Pearson Correlation] [Investigating ass...] [2007-10-22 22:08:56] [8cd6641b921d30ebe00b648d1481bba0]
F    D    [Pearson Correlation] [Associatie bouwve...] [2008-10-20 22:39:37] [de3f0516a1536f7c4a656924d8bc8d07] [Current]
Feedback Forum
2008-10-25 13:27:17 [Davy De Nef] [reply] 
Het onderling verband is hier al wat lager. Er is wel nog steeds een positieve correlatie van 16%.
Als we een rechte trekken, ligt deze wat lager dan bij de voorgaande grafieken. Je ziet hier wel een paar duidelijke outliers.
2008-10-25 13:49:49 [Astrid Sniekers] [reply] 
De student heeft geen vergelijking gemaakt tussen tijdreeks 2 en 3, 2 en 4 en 3 en 4!!!
De uitleg van de student is beperkt tot het minimum. Beknopt. De uitleg had iets uitgebreider gekund.
2008-10-26 23:21:15 [Kristof Augustyns] [reply] 
Hetzelfde als vorige twee, maar deze keer is de correlatie maar 0,16 dat dus helemaal niet veel is.
Er is geen groot verband tussen bouwvergunningen van België en van het Brussels Hoofdstedelijke Gewest.
De student heeft niet veel uitleg gegeven.

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
2266
1878
2267
2069
1746
2299
2360
2214
2825
2355
2333
3016
2155
2172
2150
2533
2058
2160
2259
2498
2695
2799
2945
2930
2318
2540
2570
2669
2450
2842
3439
2677
2979
2257
2842
2546
2455
2293
2379
2478
2054
2272
2351
2271
2542
2304
2194
2722
2395
2146
1894
2548
2087
2063
2481
2476
2212
2834
2148
2598
Dataseries Y:
Download CSV
Histogram 
41
15
34
46
24
19
50
28
40
38
29
73
36
35
13
17
21
21
33
45
57
26
22
19
21
25
35
21
14
21
24
20
22
27
32
24
26
28
35
23
28
20
26
29
11
15
22
25
21
13
24
22
18
20
21
25
25
20
26
35

Tables (Output of Computation)





Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time1 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135

\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 & 'Gwilym Jenkins' @ 72.249.127.135 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=18274&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]'Gwilym Jenkins' @ 72.249.127.135[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=18274&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=18274&T=0

As an alternative you can also use a QR Code:  
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 time1 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135






Pearson Product Moment Correlation - Ungrouped Data
StatisticVariable XVariable Y
Mean2421.827.1
Biased Variance99145.9266666667120.123333333333
Biased Standard Deviation314.87446175685110.9600790751405
Covariance567.596610169492
Correlation0.161729568815679
Determination0.0261564534293053
T-Test1.24812714949686
p-value (2 sided)0.216999730021595
p-value (1 sided)0.108499865010797
Degrees of Freedom58
Number of Observations60

\begin{tabular}{lllllllll}
\hline
Pearson Product Moment Correlation - Ungrouped Data \tabularnewline
Statistic & Variable X & Variable Y \tabularnewline
Mean & 2421.8 & 27.1 \tabularnewline
Biased Variance & 99145.9266666667 & 120.123333333333 \tabularnewline
Biased Standard Deviation & 314.874461756851 & 10.9600790751405 \tabularnewline
Covariance & 567.596610169492 \tabularnewline
Correlation & 0.161729568815679 \tabularnewline
Determination & 0.0261564534293053 \tabularnewline
T-Test & 1.24812714949686 \tabularnewline
p-value (2 sided) & 0.216999730021595 \tabularnewline
p-value (1 sided) & 0.108499865010797 \tabularnewline
Degrees of Freedom & 58 \tabularnewline
Number of Observations & 60 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=18274&T=1

[TABLE]
[ROW][C]Pearson Product Moment Correlation - Ungrouped Data[/C][/ROW]
[ROW][C]Statistic[/C][C]Variable X[/C][C]Variable Y[/C][/ROW]
[ROW][C]Mean[/C][C]2421.8[/C][C]27.1[/C][/ROW]
[ROW][C]Biased Variance[/C][C]99145.9266666667[/C][C]120.123333333333[/C][/ROW]
[ROW][C]Biased Standard Deviation[/C][C]314.874461756851[/C][C]10.9600790751405[/C][/ROW]
[ROW][C]Covariance[/C][C]567.596610169492[/C][/ROW]
[ROW][C]Correlation[/C][C]0.161729568815679[/C][/ROW]
[ROW][C]Determination[/C][C]0.0261564534293053[/C][/ROW]
[ROW][C]T-Test[/C][C]1.24812714949686[/C][/ROW]
[ROW][C]p-value (2 sided)[/C][C]0.216999730021595[/C][/ROW]
[ROW][C]p-value (1 sided)[/C][C]0.108499865010797[/C][/ROW]
[ROW][C]Degrees of Freedom[/C][C]58[/C][/ROW]
[ROW][C]Number of Observations[/C][C]60[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=18274&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=18274&T=1

As an alternative you can also use a QR Code:  
QR Code


The GUIDs for individual cells are displayed in the table below:

Pearson Product Moment Correlation - Ungrouped Data
StatisticVariable XVariable Y
Mean2421.827.1
Biased Variance99145.9266666667120.123333333333
Biased Standard Deviation314.87446175685110.9600790751405
Covariance567.596610169492
Correlation0.161729568815679
Determination0.0261564534293053
T-Test1.24812714949686
p-value (2 sided)0.216999730021595
p-value (1 sided)0.108499865010797
Degrees of Freedom58
Number of Observations60


Figures (Output of Computation)

Input Parameters & R Code

Parameters (Session):
Parameters (R input):
R code (references can be found in the software module):
bitmap(file='test1.png')
histx <- hist(x, plot=FALSE)
histy <- hist(y, plot=FALSE)
maxcounts <- max(c(histx$counts, histx$counts))
xrange <- c(min(x),max(x))
yrange <- c(min(y),max(y))
nf <- layout(matrix(c(2,0,1,3),2,2,byrow=TRUE), c(3,1), c(1,3), TRUE)
par(mar=c(4,4,1,1))
plot(x, y, xlim=xrange, ylim=yrange, xlab=xlab, ylab=ylab)
par(mar=c(0,4,1,1))
barplot(histx$counts, axes=FALSE, ylim=c(0, maxcounts), space=0)
par(mar=c(4,0,1,1))
barplot(histy$counts, axes=FALSE, xlim=c(0, maxcounts), space=0, horiz=TRUE)
dev.off()
lx = length(x)
makebiased = (lx-1)/lx
varx = var(x)*makebiased
vary = var(y)*makebiased
corxy <- cor.test(x,y,method='pearson')
cxy <- as.matrix(corxy$estimate)[1,1]
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Pearson Product Moment Correlation - Ungrouped Data',3,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Statistic',1,TRUE)
a<-table.element(a,'Variable X',1,TRUE)
a<-table.element(a,'Variable Y',1,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,hyperlink('arithmetic_mean.htm','Mean',''),header=TRUE)
a<-table.element(a,mean(x))
a<-table.element(a,mean(y))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,hyperlink('biased.htm','Biased Variance',''),header=TRUE)
a<-table.element(a,varx)
a<-table.element(a,vary)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,hyperlink('biased1.htm','Biased Standard Deviation',''),header=TRUE)
a<-table.element(a,sqrt(varx))
a<-table.element(a,sqrt(vary))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,hyperlink('covariance.htm','Covariance',''),header=TRUE)
a<-table.element(a,cov(x,y),2)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,hyperlink('pearson_correlation.htm','Correlation',''),header=TRUE)
a<-table.element(a,cxy,2)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,hyperlink('coeff_of_determination.htm','Determination',''),header=TRUE)
a<-table.element(a,cxy*cxy,2)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,hyperlink('ttest_statistic.htm','T-Test',''),header=TRUE)
a<-table.element(a,as.matrix(corxy$statistic)[1,1],2)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'p-value (2 sided)',header=TRUE)
a<-table.element(a,(p2 <- as.matrix(corxy$p.value)[1,1]),2)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'p-value (1 sided)',header=TRUE)
a<-table.element(a,p2/2,2)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Degrees of Freedom',header=TRUE)
a<-table.element(a,lx-2,2)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Number of Observations',header=TRUE)
a<-table.element(a,lx,2)
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable.tab')