FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

Author's title

Author*Unverified author*
R Software Modulerwasp_partialcorrelation.wasp
Title produced by softwarePartial Correlation
Date of computationSun, 09 Nov 2008 06:21:03 -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/09/t1226237349sxkxy93616cz8v5.htm/, Retrieved Sun, 19 May 2024 10:11:18 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=22745, Retrieved Sun, 19 May 2024 10:11:18 +0000
QR Codes:

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

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] [Partial Correlation] [2008-11-09 13:21:03] [3b916296c2d2371d528ff188880e3d2b] [Current]
Feedback Forum
2008-11-22 13:39:35 [An Knapen] [reply] 
De partiële correlatie is inderdaad een methode om de invloed van een derde variabele(controle variabele) te onderzoeken op het verband tussen 2 variabelen. Die derde variabele zou er immers voor kunnen zorgen dat er een vertekend beeld is van het verband tussen de 2 variabelen.
(x,y.z) = effect van de variabele z op het verband tussen x en y.
Dit effect is= 0.78 wat meteen ook de grootste waarde is van de 3.
Het kleinste effect wordt veroorzaakt door de invloed van x op het verband tussen y en z. Dit is ook de enige waarde die negatief is.
2008-11-23 13:38:45 [Nathalie Daneels] [reply] 
Evaluatie opdracht 3 - Blok 8 (Q1)

De student zou wel bij de conclusie moeten vermelden welke variabele precies de derde variabele (z) is.
De student zou de evaluatie nog kunnen aanvullen met het volgende:
In deze tabel wordt telkens de correlatie berekend tussen 2 variabelen en vervolgens de invloed van een derde variabele op de correlatie van de 2 (andere) variabelen. Bij de bivariate density zijn we op zoek gegaan naar een verband tussen de variabelen x en y. We kwamen daar tot de conclusie dat dit verband afgerond 0,82 bedraagt. (Dit zien we ook in de tabel van de partiële correlatie). Nu zou het eventueel mogelijk kunnen zijn dat deze correlatie een schijncorrelatie is: Dit betekent dat er een verband lijkt te zijn tussen twee variabelen, maar dat dit verband grotendeels wordt bepaald door de invloed van een derde variabele. Bijvoorbeeld: Er bleek een positief verband te bestaan tussen het aantal ooievaarsnesten en het aantal pasgeborenen in de Westerse wereld. Hier kunnen we dus spreken van een schijncorrelatie. Het zou zeer vreemd zijn dat er een verband bestaat tussen deze 2 variabelen. Na verder onderzoek bleek er inderdaad een derde (vooraf onbekende) variabele te bestaan die een invloed had op de correlatie tussen het aantal ooievaarsnesten en het aantal pasgeborenen. Deze derde variabele is 'de vervuiling door industrialisatie', die ervoor zorgde dat de andere 2 variabelen een verband leken te hebben.
Om na te gaan of deze correlatie inderdaad een schijncorrelatie is, gaan we de invloed van variabele z op de correlatie tussen x en y berekenen. Uit de grafiek kunnen we afleiden dat deze correlatie (xy, z) 0,77 afgerond is. We kunnen dus inderdaad stellen dat de variabele z een invloed heeft op de correlatie tussen de variabelen x en y. We kunnen geen uitspraak doen over het feit of deze invloed groot of klein is. We kunnen niet precies zeggen wanneer een invloed groot of klein is.
Er treedt perfecte correlatie op als het -1 of 1 is en er is geen correlatie als het 0 is. Als correlatie van (xy,z) bv. 0,03 is, wel dan kunnen we eigenlijk wel stellen dat de derde variabele geen invloed heeft op de correlatie van x en y en die correlatie dus geen schijncorrelatie is.

De student heeft wel goed uitgelegd waarom we precies de partiële correlatie gebruiken.

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
97,3
101
113,2
101
105,7
113,9
86,4
96,5
103,3
114,9
105,8
94,2
98,4
99,4
108,8
112,6
104,4
112,2
81,1
97,1
112,6
113,8
107,8
103,2
103,3
101,2
107,7
110,4
101,9
115,9
89,9
88,6
117,2
123,9
100
103,6
94,1
98,7
119,5
112,7
104,4
124,7
89,1
97
121,6
118,8
114
111,5
97,2
102,5
113,4
109,8
104,9
126,1
80
96,8
117,2
112,3
117,3
111,1
102,2
104,3
122,9
107,6
121,3
131,5
89
104,4
128,9
135,9
133,3
121,3
120,5
120,4
137,9
126,1
133,2
151,1
105
119
140,4
156,6
137,1
122,7
Dataseries Y:
Download CSV
Histogram 
93,5
94,7
112,9
99,2
105,6
113
83,1
81,1
96,9
104,3
97,7
102,6
89,9
96
112,7
107,1
106,2
121
101,2
83,2
105,1
113,3
99,1
100,3
93,5
98,8
106,2
98,3
102,1
117,1
101,5
80,5
105,9
109,5
97,2
114,5
93,5
100,9
121,1
116,5
109,3
118,1
108,3
105,4
116,2
111,2
105,8
122,7
99,5
107,9
124,6
115
110,3
132,7
99,7
96,5
118,7
112,9
130,5
137,9
115
116,8
140,9
120,7
134,2
147,3
112,4
107,1
128,4
137,7
135
151
137,4
132,4
161,3
139,8
146
166,5
143,3
121
152,6
154,4
154,6
158
Dataseries Z:
Download CSV
Histogram 
105,5
106,4
117,9
89,7
88,5
106,4
61,4
92,3
95,5
92,5
89,6
84,3
76,3
80,7
96,3
81
82,9
90,3
74,8
70,1
86,7
86,4
89,9
88,1
78,8
81,1
85,4
82,6
80,3
81,2
68
67,4
91,3
94,9
82,8
88,6
73,1
76,7
93,2
84,9
83,8
93,5
91,9
69,6
87
90,2
82,7
91,4
74,6
76,1
87,1
78,4
81,3
99,3
71
73,2
95,6
84
90,8
93,6
80,9
84,4
97,3
83,5
88,8
100,7
69,4
74,6
96,6
96,6
93,1
91,8
85,7
79,1
91,3
84,2
85,8
94,6
77,1
76,5
89,7
103,6
100
96,6

Tables (Output of Computation)





Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time27 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 & 27 seconds \tabularnewline
R Server & 'George Udny Yule' @ 72.249.76.132 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=22745&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]27 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=22745&T=0

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






Pearson Product Moment Partial Correlation - Ungrouped Data
StatisticValue
Correlation r(xy)0.821575855124337
Partial Correlation r(xy.z)0.778798066662208
Correlation r(xz)0.578217403445815
Partial Correlation r(xz.y)0.44059095810883
Correlation r(yz)0.427389934930984
Partial Correlation r(yz.x)-0.102464046850072

\begin{tabular}{lllllllll}
\hline
Pearson Product Moment Partial Correlation - Ungrouped Data \tabularnewline
Statistic & Value \tabularnewline
Correlation r(xy) & 0.821575855124337 \tabularnewline
Partial Correlation r(xy.z) & 0.778798066662208 \tabularnewline
Correlation r(xz) & 0.578217403445815 \tabularnewline
Partial Correlation r(xz.y) & 0.44059095810883 \tabularnewline
Correlation r(yz) & 0.427389934930984 \tabularnewline
Partial Correlation r(yz.x) & -0.102464046850072 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=22745&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.821575855124337[/C][/ROW]
[ROW][C]Partial Correlation r(xy.z)[/C][C]0.778798066662208[/C][/ROW]
[ROW][C]Correlation r(xz)[/C][C]0.578217403445815[/C][/ROW]
[ROW][C]Partial Correlation r(xz.y)[/C][C]0.44059095810883[/C][/ROW]
[ROW][C]Correlation r(yz)[/C][C]0.427389934930984[/C][/ROW]
[ROW][C]Partial Correlation r(yz.x)[/C][C]-0.102464046850072[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=22745&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=22745&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.821575855124337
Partial Correlation r(xy.z)0.778798066662208
Correlation r(xz)0.578217403445815
Partial Correlation r(xz.y)0.44059095810883
Correlation r(yz)0.427389934930984
Partial Correlation r(yz.x)-0.102464046850072


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