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_correlation.wasp

Title produced by software

Pearson Correlation

Date of computation

Mon, 28 Dec 2009 07:57:00 -0700

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/2009/Dec/28/t12620123414ag4fbu96o2dqmt.htm/, Retrieved Sun, 05 May 2024 13:05:38 +0000

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=70987, Retrieved Sun, 05 May 2024 13:05:38 +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

Paper

Estimated Impact

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)

-     [Pearson Correlation] [SHW_WS5_Correlation] [2009-10-29 18:55:36] [8b1aef4e7013bd33fbc2a5833375c5f5]
- RMPD    [Pearson Correlation] [e[t]_vergelijken] [2009-12-28 14:57:00] [5b5bced41faf164488f2c271c918b21f] [Current]

Feedback Forum

Post a new message

Dataseries X:

Download CSV

Histogram

Boxplots

-50,74706514
-48,93785253
-47,15852374
-45,17543128
-46,94276659
-47,29735115
-42,25443623
-45,636011
-40,5706558
-34,94664303
-32,38276206
-32,7877051
-26,14063339
-22,09709979
-21,64764034
-16,98759233
-17,90735002
-15,42376615
-11,53425643
-13,32206639
-9,014273186
-4,253762759
-0,425986315
2,843751224
5,304276145
3,352359662
4,868775794
2,749804351
0,226872673
0,617183142
-0,16183631
0,436478416
2,548934314
6,596089819
10,3692672
9,720903334
12,32020087
12,14404607
12,94839943
13,76345412
16,31332119
16,81965467
15,68801627
17,19533256
18,78454744
24,92140923
27,36016736
27,69815372
30,48854798
30,43069115
30,29093572
28,97275052
27,2951215
25,22718201
26,17030799
25,74716525
29,830356
38,0397072
56,54216423
42,59868849

Dataseries Y:

Download CSV

Histogram

-7,384479833
-348,812296
656,3387963
168,4268193
163,9341261
-219,1238071
131,7351533
68,87952107
-499,4813227
88,460018
-3,676899151
-250,2346886
174,1014888
-231,4249098
544,4571966
124,3024887
245,2234366
98,23681742
372,6587033
36,0872455
-337,1878898
-72,72145156
-121,8431801
-137,6524521
19,966092
-447,5428034
769,4438158
243,4453361
-164,2602175
454,3943467
-69,44942508
405,8522526
-358,8682021
-83,29905889
50,36917358
-94,12503992
-307,9726466
-171,4972374
581,7501227
244,4978459
74,07234434
218,1732414
41,97296758
230,9245802
-379,8087485
-230,0475763
-14,66583292
-134,2125562
-277,3250263
-389,3945325
682,2209027
-261,6970045
148,0023313
158,047103
-258,1716876
153,2005099
-667,3263955
-162,5924918
-58,96570282
-588,4092144

Summary of computational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	3 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 & 3 seconds \tabularnewline
R Server & 'Gwilym Jenkins' @ 72.249.127.135 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=70987&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]3 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=70987&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=70987&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	3 seconds
R Server	'Gwilym Jenkins' @ 72.249.127.135

Pearson Product Moment Correlation - Ungrouped Data
Statistic	Variable X	Variable Y
Mean	-0.275474166933333	-3.06666614324769e-09
Biased Variance	728.100939911876	95422.97835522
Biased Standard Deviation	26.9833456026468	308.906099575939
Covariance	-1356.47724206378
Correlation	-0.160026163488448
Determination	0.0256083730008315
T-Test	-1.23463395999701
p-value (2 sided)	0.221945447664239
p-value (1 sided)	0.110972723832120
Degrees of Freedom	58
Number of Observations	60

\begin{tabular}{lllllllll}
\hline
Pearson Product Moment Correlation - Ungrouped Data \tabularnewline
Statistic & Variable X & Variable Y \tabularnewline
Mean & -0.275474166933333 & -3.06666614324769e-09 \tabularnewline
Biased Variance & 728.100939911876 & 95422.97835522 \tabularnewline
Biased Standard Deviation & 26.9833456026468 & 308.906099575939 \tabularnewline
Covariance & -1356.47724206378 \tabularnewline
Correlation & -0.160026163488448 \tabularnewline
Determination & 0.0256083730008315 \tabularnewline
T-Test & -1.23463395999701 \tabularnewline
p-value (2 sided) & 0.221945447664239 \tabularnewline
p-value (1 sided) & 0.110972723832120 \tabularnewline
Degrees of Freedom & 58 \tabularnewline
Number of Observations & 60 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=70987&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]-0.275474166933333[/C][C]-3.06666614324769e-09[/C][/ROW]
[ROW][C]Biased Variance[/C][C]728.100939911876[/C][C]95422.97835522[/C][/ROW]
[ROW][C]Biased Standard Deviation[/C][C]26.9833456026468[/C][C]308.906099575939[/C][/ROW]
[ROW][C]Covariance[/C][C]-1356.47724206378[/C][/ROW]
[ROW][C]Correlation[/C][C]-0.160026163488448[/C][/ROW]
[ROW][C]Determination[/C][C]0.0256083730008315[/C][/ROW]
[ROW][C]T-Test[/C][C]-1.23463395999701[/C][/ROW]
[ROW][C]p-value (2 sided)[/C][C]0.221945447664239[/C][/ROW]
[ROW][C]p-value (1 sided)[/C][C]0.110972723832120[/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=70987&T=1

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

As an alternative you can also use a QR Code:

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

Pearson Product Moment Correlation - Ungrouped Data
Statistic	Variable X	Variable Y
Mean	-0.275474166933333	-3.06666614324769e-09
Biased Variance	728.100939911876	95422.97835522
Biased Standard Deviation	26.9833456026468	308.906099575939
Covariance	-1356.47724206378
Correlation	-0.160026163488448
Determination	0.0256083730008315
T-Test	-1.23463395999701
p-value (2 sided)	0.221945447664239
p-value (1 sided)	0.110972723832120
Degrees of Freedom	58
Number of Observations	60

Figure 1

PNG link

Postscript link

PDF link

Parameters (Session):

par1 = 0 ; par2 = 36 ;

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

Free Statistics

Description of Statistical Computation

Tree of Dependent Computations

Dataset

Tables (Output of Computation)

Figures (Output of Computation)

Input Parameters & R Code