Repository of Reproducible Computations

Free Statistics

of Irreproducible Research!

Author's title

Author

*Unverified author*

R Software Module

rwasp_partialcorrelation.wasp

Title produced by software

Partial Correlation

Date of computation

Wed, 12 Nov 2008 10:36:31 -0700

Cite this page as follows

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2008/Nov/12/t1226511421em1pfkcpplrki9p.htm/, Retrieved Sun, 19 May 2024 12:38:34 +0000

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=24316, Retrieved Sun, 19 May 2024 12:38:34 +0000

QR Codes:

Paste this QR Code to cite your computation.

Original text written by user:

IsPrivate?

No (this computation is public)

User-defined keywords

Estimated Impact

169

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)

-     [Bivariate Kernel Density Estimation] [Bel20 en Downjones] [2008-11-12 17:23:23] [74be16979710d4c4e7c6647856088456]
F RMPD    [Partial Correlation] [kelly] [2008-11-12 17:36:31] [d41d8cd98f00b204e9800998ecf8427e] [Current]

Feedback Forum

2008-11-15 13:56:41 [Hundra Smet] [reply] 
ook hier juiste voor- en nadelen. ik zou het iets uitgebreider maken.  
 
om de werking van de partial correlation te verduidelijken zal ik een voorbeeld geven.  
bij part corr r(xz.y) noemen we de variabele y de controlevariabele. deze zal een vertekenend effect hebben op het verband tussen x en z.  
 
de correlatie tussen x en y is zeer hoog: 0,93.  
maar ook de partial correlation tussen x en y is hier groot (0,84). de beide variabelen hebben dus, ondanks de invloed van y nog steeds een positieve correlatie en het resultaat van y is niet zo vertekenend. 
2008-11-24 20:50:53 [Jonas Scheltjens] [reply] 
Bij de partiële correlatie gaat het erom dat de 3de dimensie wordt weggelaten (dit heeft de student niet besproken). Dus, er worden telkens slecht 2 reeksen vergeleken, terwijl de derde buiten beschouwing wordt gelaten. Dit had best nog vermeld moeten worden. De partiële correlatie  levert vaak verbanden op. Maar dit is echter niet steeds een oorzakelijk verband (zoals er bijvoorbeeld een positieve correlatie werd gevonden tussen de aantallen van geboortes van baby’s bij mensen en bij ooievaars, maar het is absurd om te stellen dat er hier sprake is van een oorzakelijk verband). Dit verschijnsel noemt men dan ook wel een schijncorrelatie. Waarom zo’n schijncorrelatie dan toch een positief verband kan tonen ligt aan het feit dat er tussen de gegevensreeksen onderling misschien geen verband bestaat, maar dat de relatie zo’n verloop vertoont vanwege een derde (of in veel gevallen zijn er zeer veel) variabele die niet in het onderzoek werd opgenomen. De student heeft gelijk wanneer hij zegt dat men met letters XYZ moet werken i.p.v. de namen, maar in tegenstelling tot de student vind ik dit niet echt hinderlijk. Wel is het juist dat de cijfers in een tabel van de partial correlation makkelijk vallen af te lezen. Hier is de correlatie tussen x en y , zonder z (0,84) het hoogst.

Post a new message

Dataseries X:

Download CSV

Histogram

Boxplots

Dataseries Y:

Download CSV

Histogram

Dataseries Z:

Download CSV

Histogram

Summary of computational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	2 seconds
R Server	'Sir Ronald Aylmer Fisher' @ 193.190.124.24

\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 & 2 seconds \tabularnewline
R Server & 'Sir Ronald Aylmer Fisher' @ 193.190.124.24 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=24316&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]2 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Sir Ronald Aylmer Fisher' @ 193.190.124.24[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=24316&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=24316&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 Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	2 seconds
R Server	'Sir Ronald Aylmer Fisher' @ 193.190.124.24

Pearson Product Moment Partial Correlation - Ungrouped Data
Statistic	Value
Correlation r(xy)	0.934464927146275
Partial Correlation r(xy.z)	0.835884801081939
Correlation r(xz)	0.87093574446003
Partial Correlation r(xz.y)	0.652775288916858
Correlation r(yz)	0.774757443781979
Partial Correlation r(yz.x)	-0.223482365480305

\begin{tabular}{lllllllll}
\hline
Pearson Product Moment Partial Correlation - Ungrouped Data \tabularnewline
Statistic & Value \tabularnewline
Correlation r(xy) & 0.934464927146275 \tabularnewline
Partial Correlation r(xy.z) & 0.835884801081939 \tabularnewline
Correlation r(xz) & 0.87093574446003 \tabularnewline
Partial Correlation r(xz.y) & 0.652775288916858 \tabularnewline
Correlation r(yz) & 0.774757443781979 \tabularnewline
Partial Correlation r(yz.x) & -0.223482365480305 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=24316&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.934464927146275[/C][/ROW]
[ROW][C]Partial Correlation r(xy.z)[/C][C]0.835884801081939[/C][/ROW]
[ROW][C]Correlation r(xz)[/C][C]0.87093574446003[/C][/ROW]
[ROW][C]Partial Correlation r(xz.y)[/C][C]0.652775288916858[/C][/ROW]
[ROW][C]Correlation r(yz)[/C][C]0.774757443781979[/C][/ROW]
[ROW][C]Partial Correlation r(yz.x)[/C][C]-0.223482365480305[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=24316&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=24316&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
Statistic	Value
Correlation r(xy)	0.934464927146275
Partial Correlation r(xy.z)	0.835884801081939
Correlation r(xz)	0.87093574446003
Partial Correlation r(xz.y)	0.652775288916858
Correlation r(yz)	0.774757443781979
Partial Correlation r(yz.x)	-0.223482365480305

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

Free Statistics

Description of Statistical Computation

Tree of Dependent Computations

Dataset

Tables (Output of Computation)

Figures (Output of Computation)

Input Parameters & R Code