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 computationFri, 23 Dec 2011 05:13:10 -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/2011/Dec/23/t1324635208j2hyodqt5l4wmhj.htm/, Retrieved Mon, 29 Apr 2024 19:08:45 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=160236, Retrieved Mon, 29 Apr 2024 19:08:45 +0000
QR Codes:

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

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Multiple Regression] [] [2010-12-05 18:56:24] [b98453cac15ba1066b407e146608df68]
- R PD  [Multiple Regression] [bambix] [2011-12-09 12:38:27] [379dab8110dbf77cfcc4b7951c3a599f]
-   P     [Multiple Regression] [bambix] [2011-12-09 12:42:52] [379dab8110dbf77cfcc4b7951c3a599f]
-           [Multiple Regression] [Multiple Regression] [2011-12-09 13:03:13] [74b1e5a3104ff0b2404b2865a63336ad]
-    D        [Multiple Regression] [Multiple Linear R...] [2011-12-09 14:02:33] [74b1e5a3104ff0b2404b2865a63336ad]
-               [Multiple Regression] [Multiple linear r...] [2011-12-23 09:00:38] [74b1e5a3104ff0b2404b2865a63336ad]
- RMPD              [Pearson Correlation] [Multicolinearitei...] [2011-12-23 10:13:10] [f9bdb25068ab2a4592adc645515299ca] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
96
75
70
134
72
8
169
1
88
98
106
122
57
139
87
176
114
121
103
135
123
99
74
103
158
116
102
132
62
150
143
50
141
48
141
83
112
79
33
149
126
85
84
68
50
101
20
101
150
118
99
8
88
21
30
97
163
132
161
89
160
139
104
103
66
163
93
85
150
143
107
22
85
91
131
140
156
81
137
102
72
161
30
120
49
71
76
85
146
165
89
168
48
149
75
107
116
165
155
165
121
156
86
13
113
112
133
169
30
121
82
148
12
146
23
84
163
4
81
18
118
76
55
62
16
98
137
50
152
163
142
77
59
94
128
63
127
59
118
110
45
96
128
41
146
147
121
185
0
4
0
0
0
0
85
157
0
0
7
12
0
37
0
62
Dataseries Y:
Download CSV
Histogram 
186099
113854
99776
106194
100792
47552
250931
6853
115466
110896
169351
94853
72591
101345
113713
165354
164263
135213
111669
134163
140303
150773
111848
102509
96785
116136
158376
153990
64057
230054
184531
114198
198299
33750
189723
100826
188355
104470
58391
164808
134097
80238
133252
54518
121850
79367
56968
106314
191889
104864
160792
15049
191179
25109
45824
129711
210012
194679
197680
81180
197765
214738
96252
124527
153242
145707
113963
134904
114268
94333
102204
23824
111563
91313
89770
100125
165278
181712
80906
75881
83963
175721
68580
136323
55792
25157
100922
118845
170492
81716
115750
105590
92795
82390
135599
127667
163073
211381
189944
226168
117495
195894
80684
19630
88634
139292
128602
135848
178377
106330
178303
116938
5841
106020
24610
74151
232241
6622
127097
13155
160501
91502
24469
88229
13983
80716
157384
122975
191469
231257
258287
122531
61394
86480
195791
18284
147581
72558
147341
114651
100187
130332
134218
10901
145758
75767
134969
169216
0
7953
0
0
0
0
105406
174586
0
0
4245
21509
7670
15673
0
75882

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

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






Pearson Product Moment Correlation - Ungrouped Data
StatisticVariable XVariable Y
Mean93.4268292682927109971
Biased Variance2498.293426531833812145606.54878
Biased Standard Deviation49.982931351930861742.5753151647
Covariance2454729.66257669
Correlation0.79057115432868
Determination0.625002750056581
T-Test16.4317731267327
p-value (2 sided)0
p-value (1 sided)0
Degrees of Freedom162
Number of Observations164

\begin{tabular}{lllllllll}
\hline
Pearson Product Moment Correlation - Ungrouped Data \tabularnewline
Statistic & Variable X & Variable Y \tabularnewline
Mean & 93.4268292682927 & 109971 \tabularnewline
Biased Variance & 2498.29342653183 & 3812145606.54878 \tabularnewline
Biased Standard Deviation & 49.9829313519308 & 61742.5753151647 \tabularnewline
Covariance & 2454729.66257669 \tabularnewline
Correlation & 0.79057115432868 \tabularnewline
Determination & 0.625002750056581 \tabularnewline
T-Test & 16.4317731267327 \tabularnewline
p-value (2 sided) & 0 \tabularnewline
p-value (1 sided) & 0 \tabularnewline
Degrees of Freedom & 162 \tabularnewline
Number of Observations & 164 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=160236&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]93.4268292682927[/C][C]109971[/C][/ROW]
[ROW][C]Biased Variance[/C][C]2498.29342653183[/C][C]3812145606.54878[/C][/ROW]
[ROW][C]Biased Standard Deviation[/C][C]49.9829313519308[/C][C]61742.5753151647[/C][/ROW]
[ROW][C]Covariance[/C][C]2454729.66257669[/C][/ROW]
[ROW][C]Correlation[/C][C]0.79057115432868[/C][/ROW]
[ROW][C]Determination[/C][C]0.625002750056581[/C][/ROW]
[ROW][C]T-Test[/C][C]16.4317731267327[/C][/ROW]
[ROW][C]p-value (2 sided)[/C][C]0[/C][/ROW]
[ROW][C]p-value (1 sided)[/C][C]0[/C][/ROW]
[ROW][C]Degrees of Freedom[/C][C]162[/C][/ROW]
[ROW][C]Number of Observations[/C][C]164[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=160236&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=160236&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
Mean93.4268292682927109971
Biased Variance2498.293426531833812145606.54878
Biased Standard Deviation49.982931351930861742.5753151647
Covariance2454729.66257669
Correlation0.79057115432868
Determination0.625002750056581
T-Test16.4317731267327
p-value (2 sided)0
p-value (1 sided)0
Degrees of Freedom162
Number of Observations164


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