Free Statistics

of Irreproducible Research!

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, 24 Nov 2011 12:22:51 -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/Nov/24/t1322155406yq6dsg1e5thk72d.htm/, Retrieved Fri, 19 Apr 2024 03:41:27 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=147110, Retrieved Fri, 19 Apr 2024 03:41:27 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact53
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Pearson Correlation] [WS7] [2011-11-24 17:22:51] [3d1c016f2f3d502d2cccf121bc4326eb] [Current]
Feedback Forum

Post a new message
Dataseries X:
72772
45104
44525
41169
31118
28435
22162
20202
17773
17094
15153
11218
10796
9594
9309
8556
8041
7639
6884
6642
6321
6216
5865
5799
5695
5644
5446
5395
5363
5338
5160
5091
5057
5039
4880
4735
4693
4653
4586
4398
3974
3858
3826
3819
3556
3372
3193
3126
3104
2967
2848
2748
2649
2625
2572
2548
2477
2442
2392
2372
Dataseries Y:
26073
18103
15100
14738
22259
10277
6225
7663
6618
9945
7590
4293
4656
5145
2001
1779
1609
2191
1617
2554
2198
1578
3446
1380
1249
1223
834
3754
2283
3028
1100
457
1201
2192
1508
1393
952
1032
1279
1370
649
1900
666
1313
1353
1500
877
874
1133
754
695
609
696
756
670
301
630
798
436
388




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time0 seconds
R Server'Herman Ole Andreas Wold' @ wold.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 & 0 seconds \tabularnewline
R Server & 'Herman Ole Andreas Wold' @ wold.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=147110&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]0 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Herman Ole Andreas Wold' @ wold.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=147110&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=147110&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 Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time0 seconds
R Server'Herman Ole Andreas Wold' @ wold.wessa.net







Pearson Product Moment Correlation - Ungrouped Data
StatisticVariable XVariable Y
Mean10000.46666666673681.51666666667
Biased Variance166312953.68222228613982.9830556
Biased Standard Deviation12896.23796625295349.20395788528
Covariance66730383.2971751
Correlation0.951200752359902
Determination0.904782871290044
T-Test23.4762454405747
p-value (2 sided)0
p-value (1 sided)0
Degrees of Freedom58
Number of Observations60

\begin{tabular}{lllllllll}
\hline
Pearson Product Moment Correlation - Ungrouped Data \tabularnewline
Statistic & Variable X & Variable Y \tabularnewline
Mean & 10000.4666666667 & 3681.51666666667 \tabularnewline
Biased Variance & 166312953.682222 & 28613982.9830556 \tabularnewline
Biased Standard Deviation & 12896.2379662529 & 5349.20395788528 \tabularnewline
Covariance & 66730383.2971751 \tabularnewline
Correlation & 0.951200752359902 \tabularnewline
Determination & 0.904782871290044 \tabularnewline
T-Test & 23.4762454405747 \tabularnewline
p-value (2 sided) & 0 \tabularnewline
p-value (1 sided) & 0 \tabularnewline
Degrees of Freedom & 58 \tabularnewline
Number of Observations & 60 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=147110&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]10000.4666666667[/C][C]3681.51666666667[/C][/ROW]
[ROW][C]Biased Variance[/C][C]166312953.682222[/C][C]28613982.9830556[/C][/ROW]
[ROW][C]Biased Standard Deviation[/C][C]12896.2379662529[/C][C]5349.20395788528[/C][/ROW]
[ROW][C]Covariance[/C][C]66730383.2971751[/C][/ROW]
[ROW][C]Correlation[/C][C]0.951200752359902[/C][/ROW]
[ROW][C]Determination[/C][C]0.904782871290044[/C][/ROW]
[ROW][C]T-Test[/C][C]23.4762454405747[/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]58[/C][/ROW]
[ROW][C]Number of Observations[/C][C]60[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=147110&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=147110&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
StatisticVariable XVariable Y
Mean10000.46666666673681.51666666667
Biased Variance166312953.68222228613982.9830556
Biased Standard Deviation12896.23796625295349.20395788528
Covariance66730383.2971751
Correlation0.951200752359902
Determination0.904782871290044
T-Test23.4762454405747
p-value (2 sided)0
p-value (1 sided)0
Degrees of Freedom58
Number of Observations60



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