Free Statistics

of Irreproducible Research!

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 08:04:15 -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/t1226415902f55y5uaxxlgd60f.htm/, Retrieved Sun, 19 May 2024 11:33:41 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=23565, Retrieved Sun, 19 May 2024 11:33:41 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact168
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
F       [Partial Correlation] [Q1] [2008-11-11 15:04:15] [e11d930c9e2984715c66c796cf63ef19] [Current]
Feedback Forum
2008-11-16 15:22:34 [Jan Van Riet] [reply
De correlatiewaarde ligt inderdaad het hoogst tussen y en z (vlaams en waals gewest). Ook de partiële correlatiewaarde (met de invloed van 'z' eruit gezuiverd) is hoog tussen deze 2 tijdreeksen.
2008-11-21 14:43:12 [] [reply
De 3 partiële correlaties vertonen een verschillend beeld. De correlatie tussen x en y is eerst negatief, nadat de derde variabele z erbij komt wordt de partiële correlatie positief.

Bij de correlatie tussen x en z verandert er bijna niets door het toevoegen van de derde variabele y.

Bij de correlatie tussen y en z verandert er eveneens weinig door het toevoegen van de derde variabele x. De correlatie en de partiële correlatie zijn hier ook het grootst.
2008-11-21 14:51:01 [] [reply
Door het gebruik van de partiële correlatie kan men de invloed van een derde variabele op een relatie bekijken.
2008-11-21 14:51:56 [Matthieu Blondeau] [reply
De 3 partiële correlaties vertonen een verschillend beeld. De correlatie tussen x en y is eerst negatief, nadat de derde variabele z erbij komt wordt de partiële correlatie positief.

Bij de correlatie tussen x en z verandert er bijna niets door het toevoegen van de derde variabele y.

Bij de correlatie tussen y en z verandert er eveneens weinig door het toevoegen van de derde variabele x. De correlatie en de partiële correlatie zijn hier ook het grootst.
2008-11-21 14:52:09 [Matthieu Blondeau] [reply
Door het gebruik van de partiële correlatie kan men de invloed van een derde variabele op een relatie bekijken.
2008-11-22 11:11:25 [Wim Golsteyn] [reply
Opvallend is dat door de partiele correlatie te gebruiken er niet veel verandert tussen 'x en z' en 'y en z', maar wanneer we 'z' zuiveren uit de correlatie tussen 'x en y' verandert de correlatie van negatief naar positief. De absolute waarde blijft wel ongeveer gelijk. Hier zouden conclusies uit getrokken kunnen worden, maar aangezien het voor mij niet helemaal duidelijk is welk gewest x, y en z voorstellen en ook niet wat er juist bedoeld wordt met de gewesten (% BNP? BBP? werkloosheid?) is het moeilijk hier een theorie aan vast te kleven.
2008-11-22 12:27:01 [Inge Meelberghs] [reply
Partiële correlatie kunne we omschrijven als correlatie tussen twee variabelen na correctie voor een derde variable.

Als eerst berekent het model de correlatie tussen twee variabelen(XY,XZ,UZ). Daarna wordt er een derde variable geïntroduceerd, Z. Deze dient om na te gaan of Z al dat niet een invloed heeft. Als de waarde van de simpele correlatie kort bij die van de partiële correlatie ligt, dan kunnen we zeggen dat Z geen invloedrijke variable is. Is het net andersom en liggen beide waarden relatief ver uit elkaar, dan is Z wél een invloedrijke variable en kunnen we zeggen dat de correlatie een vertekend beeld geeft.

In dit voorbeeld bedraagt de correlatie voor r(xz)-0.482775445396066 en 0.954781586139669 voor r(yz). Als we hierna de derde variabele Z toevoegen kunnen we stellen dat deze geen invloed heeft op de twee correlaties (want de waarden van de partiële correlatie zijn praktisch hetzelfde als die van de simpele correlatie) wat wil zeggen dat de correlatie geen vertekend beeld geeft.

Maar als we dan naar de correlatie van r(xy) kijken na toevoeging van de derde variabele Z, kunnen we zien dat deze wél een invloed heeft. Hier ligt de waarde van de partiële correlatie beduidend hoger dan die van de simpele correlatie ( de waarde gaat van negatief naar positief). In dit geval geeft de correlatie dus wel een vertekend beeld.
2008-11-23 16:00:27 [Nathalie Boden] [reply
Je gaat hier inderdaad 3 variabelen met elkaar gaan vergelijken. Bij de tijdreeksen x,y en z. We gaan dus de invloed van de andere variabelen uitzuiveren of wegvegen. Bijvoorbeeld het verband tussen y en z geeft een correlatie van 0.955274545743278 en krijgen we een partiële correlatie van 0.954781586139669 dus we kunnen besluiten dat de waarden hierbij dicht bij elkaar liggen en niet veel veranderen. Dus het uitzuiveren van x heeft geen invloed op op y en z. We moeten ook altijd rekening houden dat je een vertekend beeld kan hebben.
2008-11-24 16:24:16 [Bernard Femont] [reply
Zo moet je de partial correlation interpreteren:
Hierbij hebben we 3 variabelen nodig: X, Y en Z.
We gaan hierbij niet de correlatie tussen de 3 variabelen tegelijkertijd onderzoeken, maar wel 2 aan 2 correlaties en de partiële correlatie. De partiële correlatie houdt in dat we de correlatie tussen 2 variabelen gaan onderzoeken onder invloed van de 3e variabele.
Dikwijls krijg je een vertekend beeld bij een correlatie tussen 2 variabelen. Het voordeel van deze methode is dat je door het toevoegen van de 3e variabele dit effect weggezuiverd wordt.
De 3 partiële correlaties vertonen een verschillend beeld. De correlatie tussen x en y is eerst negatief, nadat de derde variabele z erbij komt wordt de partiële correlatie positief.
Bij de correlatie tussen x en z verandert er bijna niets door het toevoegen van de derde variabele y.
Bij de correlatie tussen y en z verandert er eveneens weinig door het toevoegen van de derde variabele x. De correlatie en de partiële correlatie zijn hier ook het grootst.
2008-11-24 19:21:55 [Liese Tormans] [reply
De correlatie geeft het verband weer tussen twee variabele de invloed van de derde variabele is hier nog niet uitgezuiverd. Wat soms kan leiden tot een verkeerde conclusie. De oplossing voor dit probleem is de partial correlatie, deze geeft de correlatie weer tussen twee variabelen na filtering van een derde variabele. De partial correlatie gaat dan kijken welke invloed de derde variabele heeft op de eerste en de tweede variabele.
Als de gewone correlatie dicht bij de partial correlatie ligt is de invloed van de derde variabele zeer klein
Maar ligt de gewone correlatie relatief ver van de partial correlatie dan kunnen we zeggen dat de derde variabele een grote invloed heeft op de andere variabele. Bij gevolg geeft de gewone correlatie een vertekend beeld. Deze invloed kan zowel een negatief als positief zijn.

De gewone correlatie van r(xy) = -0.368954960085345
De partial correlatie r(xy.z) (na zuivering van de variabele Z) =0.35612700684333
Z heeft een grote invloed op X en Y: de gewone correlatie geeft een vertekend beeld

De gewone correlatie van r(xz) = -0.482775445396066
De partial correlatie r(xz.y) (na zuivering van de variabele Y) = -0.474146194440677
Y heeft geen grote invloed op X en Z: de gewone correlatie geeft geen vertekend beeld

De gewone correlatie van r(yz) = 0.955274545743278
De partial correlatie r(yz.x) (na zuivering van de variabele X) = 0.954781586139669
X heeft geen grote invloed op Y en Z: de gewone correlatie geeft geen vertekend beeld
2008-11-24 19:46:23 [Jeroen Aerts] [reply
De partiële correlatie is een eenvoudige interpretatie om 3 variabele te bestuderen. Als er geen verband is tussen X en Y, wil dat niet uitsluiten dat er een verband bestaat tussen XZ en YZ. Er is dan sprake van een onrechtstreeks verband tussen X en Y.

De grootste verandering vindt hier plaats bij r(xy), waar we van een negatieve naar een positieve correlatie gaan.
2008-11-24 23:06:21 [Kristof Augustyns] [reply
Het is hier inderdaad juist wat hier gezegd wordt.
De correlatie tussen Vlaams en Waals gewest is hier duidelijk te zien.
Partiele correlatie is ook hoog tussen beide tijdreeksen.
Dit is goed om geen rekening te houden met de invloed.
De invloed wordt dus eerst uitgeveegd en dan pas het verband tussen de twee.

Post a new message
Dataseries X:
680,2
668,9
611,4
640,8
549,8
541,8
628,6
686,5
611,8
588,4
566,9
563,7
569,9
635,4
590,8
634,3
576,1
351,6
507,5
586,2
666,4
693,6
650,6
654,8
733,5
648,1
678,1
816,2
591
563,5
742,5
694,4
728,6
749
538,9
568,5
692,8
580,5
506,9
612,8
442,9
523,3
596,7
533,7
523,1
559,2
430,7
538,2
612,4
428
522,4
531,1
425,9
410,3
551
555,6
460,2
288,9
392,3
400,5
399
354,9
337,6
379,2
334,1
321,6
449,8
486,3
421,9
405,6
420
432,4
418,1
Dataseries Y:
12710,3
12120,8
12469,5
12054,6
12112,9
9617,2
12645,8
13581,3
12162,3
10969,7
11880
11887,6
12926,9
12300
12092,8
12380,8
12196,9
9455
13168
13427,9
11980,5
11884,8
11691,7
12233,8
14341,4
13130,7
12421,1
14285,8
12864,6
11160,2
14316,2
14388,7
14013,9
13419
12769,6
13315,5
15332,9
14243
13824,4
14962,9
13202,9
12199
15508,9
14199,8
15169,6
14058
13786,2
14147,9
16541,7
13587,5
15582,4
15802,8
14130,5
12923,2
15612,2
16033,7
16036,6
14037,8
15330,6
15038,3
17401,8
14992,5
16043,7
16929,6
15921,3
14417,2
15961
17851,9
16483,9
14215,5
17429,7
17839,5
17629,2
Dataseries Z:
2468,9
2469,2
2417,7
2411,1
2361
1924,1
2486,9
2674,9
2296,2
2101,5
2322
2273,8
2501,3
2435,2
2273,3
2454,7
2328,8
1897
2608,1
2712,3
2322
2282,6
2241,1
2417,2
2829
2600,6
2321
2768,5
2457,3
2142,7
2764,4
2788,9
2679,5
2536,5
2682,7
2699,6
3097,8
3015,2
2878
3010,9
2612,3
2419,3
3096,5
3013
3397,4
3423,1
3298,7
3065,7
3918,3
3154,4
3334,7
3461,6
3018,5
2832
3301,3
3342,8
3464,4
3016,6
3201,3
3135,3
3549,8
3247,2
3441,8
3535,6
3384,7
2996,6
3402,8
3900,2
3776,4
3197,5
4022,4
3845,1
3818,6




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=23565&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=23565&T=0

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







Pearson Product Moment Partial Correlation - Ungrouped Data
StatisticValue
Correlation r(xy)-0.368954960085345
Partial Correlation r(xy.z)0.35612700684333
Correlation r(xz)-0.482775445396066
Partial Correlation r(xz.y)-0.474146194440677
Correlation r(yz)0.955274545743278
Partial Correlation r(yz.x)0.954781586139669

\begin{tabular}{lllllllll}
\hline
Pearson Product Moment Partial Correlation - Ungrouped Data \tabularnewline
Statistic & Value \tabularnewline
Correlation r(xy) & -0.368954960085345 \tabularnewline
Partial Correlation r(xy.z) & 0.35612700684333 \tabularnewline
Correlation r(xz) & -0.482775445396066 \tabularnewline
Partial Correlation r(xz.y) & -0.474146194440677 \tabularnewline
Correlation r(yz) & 0.955274545743278 \tabularnewline
Partial Correlation r(yz.x) & 0.954781586139669 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=23565&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.368954960085345[/C][/ROW]
[ROW][C]Partial Correlation r(xy.z)[/C][C]0.35612700684333[/C][/ROW]
[ROW][C]Correlation r(xz)[/C][C]-0.482775445396066[/C][/ROW]
[ROW][C]Partial Correlation r(xz.y)[/C][C]-0.474146194440677[/C][/ROW]
[ROW][C]Correlation r(yz)[/C][C]0.955274545743278[/C][/ROW]
[ROW][C]Partial Correlation r(yz.x)[/C][C]0.954781586139669[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=23565&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=23565&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 Partial Correlation - Ungrouped Data
StatisticValue
Correlation r(xy)-0.368954960085345
Partial Correlation r(xy.z)0.35612700684333
Correlation r(xz)-0.482775445396066
Partial Correlation r(xz.y)-0.474146194440677
Correlation r(yz)0.955274545743278
Partial Correlation r(yz.x)0.954781586139669



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