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_partialcorrelation.wasp
Title produced by softwarePartial Correlation
Date of computationTue, 11 Nov 2008 10:30:13 -0700
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/Nov/11/t1226424671mulwuauisoseoxi.htm/, Retrieved Tue, 28 May 2024 04:05:20 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=23758, Retrieved Tue, 28 May 2024 04:05:20 +0000
QR Codes:

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

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
F       [Partial Correlation] [Various EDA topic...] [2008-11-11 17:30:13] [e4cb5a8878d0401c2e8d19a1768b515b] [Current]
Feedback Forum
2008-11-22 17:22:21 [Kenny Simons] [reply] 
Hier was het echter niet de bedoeling om enkel te zien tussen welke variabelen de correlatie het grootst was.

Met deze techniek ga je op zoek naar de partiële correlatie tussen 3 variabelen. Je moet uitzoeken of de derde variabele (Z), de relatie van de eerste 2 variabelen (X, Y) al dan niet beïnvloedt.

Als de relatie tussen X,Y relatief groot is, en de partiële relatie tussen X,Y en Z is veel kleiner, dan kan je veronderstellen dat Z de relatie tussen de Xen Y beïvnloedt. Z kan dus voor een gedeelte de relatie uitleggen tussen X & Y. Het zal de relatie uitleggen maar we zullen niet te weten komen wat de relatie veroorzaakt. Als de partiële correlatie dicht bij de gewone correlatie ligt, kan je stellen dat de Z-variabele geen vertekend beeld geeft tussen de relatie X en Y.

Correlation r(xy) 0.9270237189657
Partial Correlation r(xy.z) 0.933373766144216

Hier is de correlatie bijna gelijk aan de gewone correlatie, dus kunnen we stellen dat de variabele Z de relatie tussen X en Y niet beïnvloedt.

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
101,5
101,3
99,3
100,6
101,2
99,8
100,6
101,1
101,2
101,5
102,2
102,5
101,4
103,8
105,2
105,3
104,4
104,9
106,9
107,6
106,7
106,1
106,3
105,8
104,4
103,8
102,4
103,3
103,5
104,5
103,5
103,9
103,1
102,2
104,7
105,9
106,6
106,6
107,5
107,2
109
108,4
107
108
110,8
110,9
109,7
111
111,5
111
111,8
111,4
110,8
111,9
112,9
111,8
111
112,3
112,4
111,1
Dataseries Y:
Download CSV
Histogram 
110,4
112,9
109,4
111,9
108,9
113,8
114,5
113,2
111
114,6
113,1
113,2
115,1
117,6
117,8
115,7
115,7
118,3
117,9
117,3
119,4
117,1
119
120
118,9
116
115,6
119,7
119,7
120,8
120
120,2
121,7
116,9
122,4
122,6
123,7
120,9
124,2
122,6
125,7
123,1
122,2
126,2
124,4
127,8
124,2
126,7
126,1
128,2
130,4
130,2
129,2
129,7
131
129,2
131,1
132,9
135,2
132,3
Dataseries Z:
Download CSV
Histogram 
92,1
88,5
84,6
87
83,6
84,8
84,1
84,1
80,5
82,6
85,6
83,3
86,1
84,7
85,7
84,9
84,2
85,2
86,1
86
84,5
87,2
83,5
81,9
78,5
81,1
79,2
80,9
81,8
79,4
83,4
81,1
79,8
79,7
84
83,7
83,5
83,6
86
86,8
86,9
89
87,8
86,8
88,8
85,9
86,7
87,9
88,5
88,7
88,1
85,7
86,1
85,7
84,3
86,4
85,4
86,9
85,6
86,4

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

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






Pearson Product Moment Partial Correlation - Ungrouped Data
StatisticValue
Correlation r(xy)0.9270237189657
Partial Correlation r(xy.z)0.933373766144216
Correlation r(xz)0.4301472411727
Partial Correlation r(xz.y)0.503477773338069
Correlation r(yz)0.267778247802182
Partial Correlation r(yz.x)-0.386895555160939

\begin{tabular}{lllllllll}
\hline
Pearson Product Moment Partial Correlation - Ungrouped Data \tabularnewline
Statistic & Value \tabularnewline
Correlation r(xy) & 0.9270237189657 \tabularnewline
Partial Correlation r(xy.z) & 0.933373766144216 \tabularnewline
Correlation r(xz) & 0.4301472411727 \tabularnewline
Partial Correlation r(xz.y) & 0.503477773338069 \tabularnewline
Correlation r(yz) & 0.267778247802182 \tabularnewline
Partial Correlation r(yz.x) & -0.386895555160939 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=23758&T=1

[TABLE]
[ROW][C]Pearson Product Moment Partial Correlation - Ungrouped Data[/C][/ROW]
[ROW][C]Statistic[/C][C]Value[/C][/ROW]
[ROW][C]Correlation r(xy)[/C][C]0.9270237189657[/C][/ROW]
[ROW][C]Partial Correlation r(xy.z)[/C][C]0.933373766144216[/C][/ROW]
[ROW][C]Correlation r(xz)[/C][C]0.4301472411727[/C][/ROW]
[ROW][C]Partial Correlation r(xz.y)[/C][C]0.503477773338069[/C][/ROW]
[ROW][C]Correlation r(yz)[/C][C]0.267778247802182[/C][/ROW]
[ROW][C]Partial Correlation r(yz.x)[/C][C]-0.386895555160939[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=23758&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=23758&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 Partial Correlation - Ungrouped Data
StatisticValue
Correlation r(xy)0.9270237189657
Partial Correlation r(xy.z)0.933373766144216
Correlation r(xz)0.4301472411727
Partial Correlation r(xz.y)0.503477773338069
Correlation r(yz)0.267778247802182
Partial Correlation r(yz.x)-0.386895555160939


Figures (Output of Computation)

Input Parameters & R Code

Parameters (Session):
Parameters (R input):
R code (references can be found in the software module):
(rho12 <- cor(x, y))
(rho23 <- cor(y, z))
(rho13 <- cor(x, z))
(rhoxy_z <- (rho12-(rho13*rho23))/(sqrt(1-(rho13*rho13)) * sqrt(1-(rho23*rho23))))
(rhoxz_y <- (rho13-(rho12*rho23))/(sqrt(1-(rho12*rho12)) * sqrt(1-(rho23*rho23))))
(rhoyz_x <- (rho23-(rho12*rho13))/(sqrt(1-(rho12*rho12)) * sqrt(1-(rho13*rho13))))
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Pearson Product Moment Partial Correlation - Ungrouped Data',2,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Statistic',1,TRUE)
a<-table.element(a,'Value',1,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Correlation r(xy)',header=TRUE)
a<-table.element(a,rho12)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,hyperlink('partial_correlation1.htm','Partial Correlation r(xy.z)',''),header=TRUE)
a<-table.element(a,rhoxy_z)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Correlation r(xz)',header=TRUE)
a<-table.element(a,rho13)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,hyperlink('partial_correlation1.htm','Partial Correlation r(xz.y)',''),header=TRUE)
a<-table.element(a,rhoxz_y)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Correlation r(yz)',header=TRUE)
a<-table.element(a,rho23)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,hyperlink('partial_correlation1.htm','Partial Correlation r(yz.x)',''),header=TRUE)
a<-table.element(a,rhoyz_x)
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable.tab')