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 computationThu, 18 Dec 2014 07:42:53 +0000
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2014/Dec/18/t14188885974n7f3u8wbgwf9zq.htm/, Retrieved Fri, 17 May 2024 13:57:00 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=270752, Retrieved Fri, 17 May 2024 13:57:00 +0000
QR Codes:

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

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Percentiles] [Intrinsic Motivat...] [2010-10-12 12:10:58] [b98453cac15ba1066b407e146608df68]
- RMPD  [Kernel Density Estimation] [] [2011-10-18 22:42:23] [b98453cac15ba1066b407e146608df68]
- RMPD    [Percentiles] [] [2011-10-18 22:46:45] [b98453cac15ba1066b407e146608df68]
- RMPD      [Notched Boxplots] [] [2011-10-18 22:58:56] [b98453cac15ba1066b407e146608df68]
- RMP         [Notched Boxplots] [] [2014-10-16 18:35:18] [bcf5edf18529a33bd1494456d2c6cb9a]
- RMPD            [Pearson Correlation] [] [2014-12-18 07:42:53] [88f8137dd67cdcd531568536dca46410] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
0
1
0
1
1
1
0
1
1
1
1
1
1
0
0
0
1
0
1
0
1
1
0
1
1
1
1
0
1
0
0
0
1
1
1
1
0
1
1
0
1
1
1
1
1
0
1
1
0
1
1
1
1
1
1
1
0
1
0
1
1
0
0
1
0
1
0
1
1
0
0
0
0
1
0
0
1
0
0
1
0
0
0
1
0
1
1
0
0
0
1
1
1
0
1
0
1
1
1
1
1
0
1
1
1
0
0
0
0
0
0
1
1
1
1
1
0
1
1
0
1
0
0
1
1
1
1
1
1
1
1
0
1
1
1
1
1
1
0
1
1
1
1
1
0
0
1
0
0
1
0
1
1
0
1
1
1
0
0
0
0
1
0
1
0
1
0
1
0
0
1
0
1
1
1
1
1
0
1
0
0
0
1
0
1
0
0
0
0
1
1
1
1
0
0
1
0
0
0
0
1
1
0
0
1
0
1
1
0
1
1
1
0
0
0
1
1
1
1
0
0
0
1
1
0
0
0
1
1
0
0
0
1
1
1
1
0
1
1
1
0
0
1
0
1
1
0
1
0
1
1
0
1
0
1
0
1
0
0
1
1
0
1
1
0
1
0
0
1
1
0
1
0
0
1
0
0
1
Dataseries Y:
Download CSV
Histogram 
91
137
92
148
131
59
128
90
83
116
42
155
96
49
104
76
66
74
99
96
108
97
116
106
80
74
114
87
140
127
74
91
98
126
98
95
133
110
70
95
86
130
96
99
68
121
131
71
102
68
89
87
49
96
100
141
102
110
100
146
147
94
52
61
98
60
118
109
68
109
115
78
118
73
162
122
65
100
82
52
115
90
121
101
104
42
96
110
108
113
57
86
88
115
85
111
102
86
114
94
64
77
105
49
95
89
78
110
117
63
131
102
63
117
57
73
77
105
31
112
139
49
56
158
128
224
105
159
167
165
159
48
119
163
153
148
188
149
63
244
150
132
161
105
162
81
97
110
104
166
88
111
145
99
162
163
109
76
109
120
91
148
108
125
119
116
149
138
148
159
164
176
202
214
188
110
48
54
50
124
121
221
150
153
154
94
156
151
157
194
159
39
100
187
105
111
162
99
186
183
138
101
177
126
101
139
114
114
162
111
75
82
159
158
67
121
32
117
165
147
165
150
154
126
149
145
109
120
172
132
169
172
114
156
167
2
113
165
165
118
173
158
155
49
220
122
151
44
141
152
103
107
175
154
110
143
131
167
137
121
168
149
94
51
168
145
140
109
66
119
164
132
126
83
142
93
117
166

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'Sir Maurice George Kendall' @ kendall.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 & 'Sir Maurice George Kendall' @ kendall.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=270752&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]'Sir Maurice George Kendall' @ kendall.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=270752&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=270752&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'Sir Maurice George Kendall' @ kendall.wessa.net






Pearson Product Moment Correlation - Ungrouped Data
StatisticVariable XVariable Y
Mean0.568345323741007116.402877697842
Biased Variance0.2453289167227371581.06790538792
Biased Standard Deviation0.49530689145492139.7626445975104
Covariance-0.774926629093837
Correlation-0.0392054033327769
Determination0.00153706365048572
T-Test-0.651830222676493
p-value (2 sided)0.515053125958899
p-value (1 sided)0.25752656297945
Degrees of Freedom276
Number of Observations278

\begin{tabular}{lllllllll}
\hline
Pearson Product Moment Correlation - Ungrouped Data \tabularnewline
Statistic & Variable X & Variable Y \tabularnewline
Mean & 0.568345323741007 & 116.402877697842 \tabularnewline
Biased Variance & 0.245328916722737 & 1581.06790538792 \tabularnewline
Biased Standard Deviation & 0.495306891454921 & 39.7626445975104 \tabularnewline
Covariance & -0.774926629093837 \tabularnewline
Correlation & -0.0392054033327769 \tabularnewline
Determination & 0.00153706365048572 \tabularnewline
T-Test & -0.651830222676493 \tabularnewline
p-value (2 sided) & 0.515053125958899 \tabularnewline
p-value (1 sided) & 0.25752656297945 \tabularnewline
Degrees of Freedom & 276 \tabularnewline
Number of Observations & 278 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=270752&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.568345323741007[/C][C]116.402877697842[/C][/ROW]
[ROW][C]Biased Variance[/C][C]0.245328916722737[/C][C]1581.06790538792[/C][/ROW]
[ROW][C]Biased Standard Deviation[/C][C]0.495306891454921[/C][C]39.7626445975104[/C][/ROW]
[ROW][C]Covariance[/C][C]-0.774926629093837[/C][/ROW]
[ROW][C]Correlation[/C][C]-0.0392054033327769[/C][/ROW]
[ROW][C]Determination[/C][C]0.00153706365048572[/C][/ROW]
[ROW][C]T-Test[/C][C]-0.651830222676493[/C][/ROW]
[ROW][C]p-value (2 sided)[/C][C]0.515053125958899[/C][/ROW]
[ROW][C]p-value (1 sided)[/C][C]0.25752656297945[/C][/ROW]
[ROW][C]Degrees of Freedom[/C][C]276[/C][/ROW]
[ROW][C]Number of Observations[/C][C]278[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=270752&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=270752&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
Mean0.568345323741007116.402877697842
Biased Variance0.2453289167227371581.06790538792
Biased Standard Deviation0.49530689145492139.7626445975104
Covariance-0.774926629093837
Correlation-0.0392054033327769
Determination0.00153706365048572
T-Test-0.651830222676493
p-value (2 sided)0.515053125958899
p-value (1 sided)0.25752656297945
Degrees of Freedom276
Number of Observations278






Normality Tests
> jarque.x
	Jarque-Bera Normality Test
data:  x
JB = 46.4005, p-value = 8.4e-11
alternative hypothesis: greater

> jarque.y
	Jarque-Bera Normality Test
data:  y
JB = 1.1364, p-value = 0.5666
alternative hypothesis: greater

> ad.x
	Anderson-Darling normality test
data:  x
A = 51.0242, p-value < 2.2e-16

> ad.y
	Anderson-Darling normality test
data:  y
A = 0.8152, p-value = 0.03473


\begin{tabular}{lllllllll}
\hline
Normality Tests \tabularnewline

> jarque.x
	Jarque-Bera Normality Test
data:  x
JB = 46.4005, p-value = 8.4e-11
alternative hypothesis: greater

 \tabularnewline

> jarque.y
	Jarque-Bera Normality Test
data:  y
JB = 1.1364, p-value = 0.5666
alternative hypothesis: greater

 \tabularnewline

> ad.x
	Anderson-Darling normality test
data:  x
A = 51.0242, p-value < 2.2e-16

 \tabularnewline

> ad.y
	Anderson-Darling normality test
data:  y
A = 0.8152, p-value = 0.03473

 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=270752&T=2

[TABLE]
[ROW][C]Normality Tests[/C][/ROW]
[ROW][C]
> jarque.x
	Jarque-Bera Normality Test
data:  x
JB = 46.4005, p-value = 8.4e-11
alternative hypothesis: greater

[/C][/ROW]
[ROW][C]
> jarque.y
	Jarque-Bera Normality Test
data:  y
JB = 1.1364, p-value = 0.5666
alternative hypothesis: greater

[/C][/ROW]
[ROW][C]
> ad.x
	Anderson-Darling normality test
data:  x
A = 51.0242, p-value < 2.2e-16

[/C][/ROW]
[ROW][C]
> ad.y
	Anderson-Darling normality test
data:  y
A = 0.8152, p-value = 0.03473

[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=270752&T=2

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

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


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

Normality Tests
> jarque.x
	Jarque-Bera Normality Test
data:  x
JB = 46.4005, p-value = 8.4e-11
alternative hypothesis: greater

> jarque.y
	Jarque-Bera Normality Test
data:  y
JB = 1.1364, p-value = 0.5666
alternative hypothesis: greater

> ad.x
	Anderson-Darling normality test
data:  x
A = 51.0242, p-value < 2.2e-16

> ad.y
	Anderson-Darling normality test
data:  y
A = 0.8152, p-value = 0.03473


Figures (Output of Computation)

Input Parameters & R Code

Parameters (Session):
Parameters (R input):
R code (references can be found in the software module):
x <- x[!is.na(y)]
y <- y[!is.na(y)]
y <- y[!is.na(x)]
x <- x[!is.na(x)]
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, sub=main)
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', na.rm = T)
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')
library(moments)
library(nortest)
jarque.x <- jarque.test(x)
jarque.y <- jarque.test(y)
if(lx>7) {
ad.x <- ad.test(x)
ad.y <- ad.test(y)
}
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Normality Tests',1,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,paste('
',RC.texteval('jarque.x'),'
',sep=''))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,paste('
',RC.texteval('jarque.y'),'
',sep=''))
a<-table.row.end(a)
if(lx>7) {
a<-table.row.start(a)
a<-table.element(a,paste('
',RC.texteval('ad.x'),'
',sep=''))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,paste('
',RC.texteval('ad.y'),'
',sep=''))
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable1.tab')
library(car)
bitmap(file='test2.png')
qq.plot(x,main='QQplot of variable x')
dev.off()
bitmap(file='test3.png')
qq.plot(y,main='QQplot of variable y')
dev.off()