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 17:55:30 -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/t12245469738c6e7vl5urs8qij.htm/, Retrieved Sun, 19 May 2024 21:16:01 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=18308, Retrieved Sun, 19 May 2024 21:16:01 +0000
QR Codes:

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

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
F     [Pearson Correlation] [Correlatie Vlaams...] [2008-10-20 22:28:14] [0996801648c22721e57e03a7deb595f2]
F    D    [Pearson Correlation] [Q10-2] [2008-10-20 23:55:30] [9d7d6e0b01d5647e3f5dd21d8ef8600e] [Current]
Feedback Forum
2008-10-26 12:20:37 [Kevin Neelen] [reply] 
De studente heeft slechts 2 verbanden onderzocht. Misschien had dit voor alle tijdreeksen gedaan kunnen worden en zouden er eventueel verbanden ontdekt kunnen worden die op het eerste zicht niet leken te bestaan.
Dit is misschien wel gebeurd geweest door de studente, maar deze computations van die correlations staan niet vermeld in het Word-document om te bewijzen dat alle tijdreeksenn onderling wel degelijk ongezocht zijn geweest (en niet gewoon 'uit het oog verloren zijn').

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
96,1
101,7
109,9
98,1
102,4
110
71,6
74,8
110,3
109,2
104,9
100,6
85
95,9
108,9
96,1
100,1
105,8
64,4
66,9
110,4
96,1
102,5
97,6
83,6
86,5
96
91
87,2
84,5
59,1
61,6
98,8
97,9
92,8
84,2
74,6
79,8
86,7
79,8
87
91,3
58,7
62,8
87,8
90,4
80,6
73,5
71,5
70,6
78,3
76
77,4
80,9
63,5
58
88,2
81,2
84,9
76,4
71,5
76,1
82,9
78,1
82
84,7
55,7
59,5
83,2
87,6
76,2
76,4
68,2
70
76,3
70,9
72,5
80,1
57,4
62,7
82,6
88,9
80,5
72,1
69,4
69,3
77,4
79,3
78,6
76,2
64
Dataseries Y:
Download CSV
Histogram 
97,4
98
106
101,2
102,7
110
91,8
97,4
111,9
114,8
107
104,2
99,5
98,2
108,8
100,1
105
108,3
96,7
100,5
115,5
114,9
110,7
107,7
113,5
106,9
119,5
109,4
106,8
118,7
109
112,9
125,1
126,6
122,6
127,4
107,1
112
122,1
111,7
113,1
128,2
115,2
117,4
132
130,8
128
132,8
116,9
110,9
123,6
117,4
122,6
123,5
111,4
113,8
131,2
126,9
126,1
121,2
118,7
117,9
135,3
120,7
126,3
129,7
113,3
120,5
135,5
137,5
130,7
133,1
121,5
120,5
137
123,6
128,5
135
120,8
121,2
132,1
134,5
133,6
135,9
124,5
124,7
133,4
132,4
125,2
135,4
117,6

Tables (Output of Computation)





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' @ 72.249.76.132

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=18308&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 time2 seconds
R Server'George Udny Yule' @ 72.249.76.132






Pearson Product Moment Correlation - Ungrouped Data
StatisticVariable XVariable Y
Mean82.557142857143118.232967032967
Biased Variance197.936075353218131.381110976935
Biased Standard Deviation14.068975632689811.4621599612348
Covariance-33.804126984127
Correlation-0.207320333661976
Determination0.042981720749713
T-Test-1.9992945676742
p-value (2 sided)0.0486266833623537
p-value (1 sided)0.0243133416811768
Degrees of Freedom89
Number of Observations91

\begin{tabular}{lllllllll}
\hline
Pearson Product Moment Correlation - Ungrouped Data \tabularnewline
Statistic & Variable X & Variable Y \tabularnewline
Mean & 82.557142857143 & 118.232967032967 \tabularnewline
Biased Variance & 197.936075353218 & 131.381110976935 \tabularnewline
Biased Standard Deviation & 14.0689756326898 & 11.4621599612348 \tabularnewline
Covariance & -33.804126984127 \tabularnewline
Correlation & -0.207320333661976 \tabularnewline
Determination & 0.042981720749713 \tabularnewline
T-Test & -1.9992945676742 \tabularnewline
p-value (2 sided) & 0.0486266833623537 \tabularnewline
p-value (1 sided) & 0.0243133416811768 \tabularnewline
Degrees of Freedom & 89 \tabularnewline
Number of Observations & 91 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=18308&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]82.557142857143[/C][C]118.232967032967[/C][/ROW]
[ROW][C]Biased Variance[/C][C]197.936075353218[/C][C]131.381110976935[/C][/ROW]
[ROW][C]Biased Standard Deviation[/C][C]14.0689756326898[/C][C]11.4621599612348[/C][/ROW]
[ROW][C]Covariance[/C][C]-33.804126984127[/C][/ROW]
[ROW][C]Correlation[/C][C]-0.207320333661976[/C][/ROW]
[ROW][C]Determination[/C][C]0.042981720749713[/C][/ROW]
[ROW][C]T-Test[/C][C]-1.9992945676742[/C][/ROW]
[ROW][C]p-value (2 sided)[/C][C]0.0486266833623537[/C][/ROW]
[ROW][C]p-value (1 sided)[/C][C]0.0243133416811768[/C][/ROW]
[ROW][C]Degrees of Freedom[/C][C]89[/C][/ROW]
[ROW][C]Number of Observations[/C][C]91[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=18308&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=18308&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
Mean82.557142857143118.232967032967
Biased Variance197.936075353218131.381110976935
Biased Standard Deviation14.068975632689811.4621599612348
Covariance-33.804126984127
Correlation-0.207320333661976
Determination0.042981720749713
T-Test-1.9992945676742
p-value (2 sided)0.0486266833623537
p-value (1 sided)0.0243133416811768
Degrees of Freedom89
Number of Observations91


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