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 11:35:28 -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/t12265149636i8qd9h8pwt19x9.htm/, Retrieved Sun, 19 May 2024 08:13:01 +0000

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=24366, Retrieved Sun, 19 May 2024 08:13:01 +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

Severijns Britt

Estimated Impact

191

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

-     [Bivariate Kernel Density Estimation] [Variuos EDA topic...] [2008-11-09 13:17:45] [3548296885df7a66ea8efc200c4aca50]
F RMPD    [Partial Correlation] [various EDA topic...] [2008-11-12 18:35:28] [78308c9f3efc33d1da821bcd963df161] [Current]

Feedback Forum

2008-11-23 15:31:12 [76d4ba45ca37c54e4ca3ca97939a2cd4] [reply] 
Je hebt de grafieken wel gemaakt, maar je bent vergeten er commentaar bij te geven. Je had bij deze grafiek kunnen schrijven dat het gaat om een weergave van 2 dimensies. En dat de hoogtelijnen dienen om de concentratie weer te geven.  De roze gebieden geven de kern weer, in kern is de waarschijnlijkheid het grootst dat hier gegevens liggen. Je ziet op de tekening meerdere kernen, en dus meerdere groepen van gegevens. Aangezien de kernen een ellipsvorm hebben, van links onder naar rechts boven, kunnen we zeggen dat er een positieve correlatie is tussen de 2 reeksen gegevens. Links van de rechte lijn, die het gemiddelde van alle gegevens voorstelt, zie je nog een cluster van gegevens. Deze gegevens wijken een beetje af van de andere gegevens, maar ook zij vormen een ellips, dus er is duidelijk een positief verband tussen deze twee variabelen. Wat ook logisch is, want beide variabelen gaan over het aantal werklozen. Dus het zou onlogisch zijn als er geen verband was tussen deze variabelen.
2008-11-23 15:41:19 [Aurélie Van Impe] [reply] 
Je hebt de tabel gemaakt, maar opnieuw heb je er geen uitleg bij gegeven. Je had kunnen zeggen dat de correlatie tussen variabele x en z de grootste is, omdat die van alle correlaties het dichtst bij 1 ligt. De correlatie is ook groot genoeg om van een verband te kunnen spreken. Om de invloed van variabele y te minimaliseren, werd een partiële correlatie berekend. Ook deze is voor de variabelen x en z de grootste, en het dichtst bij 1. Het is dus opnieuw duidelijk dat er tussen deze 2 variabelen een positief verband waar te nemen is.
2008-11-23 17:37:14 [7bf28d4d60530086dbc44ae6b648927e] [reply] 
Ook hier gaan we zoeken naar verbanden tussen 2 variabelen x en y. Er bestaat een variable z die we niet onderzoeken maar die wel een sterke invloed heeft op x en y. Met de partiele correlatie r(xy,z) krijgen we de correlatie tussen x en y waarbij de invloed van z is weggewerkt.

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	1 seconds
R Server	'Gwilym Jenkins' @ 72.249.127.135

\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 & 'Gwilym Jenkins' @ 72.249.127.135 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=24366&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]'Gwilym Jenkins' @ 72.249.127.135[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=24366&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=24366&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	1 seconds
R Server	'Gwilym Jenkins' @ 72.249.127.135

Pearson Product Moment Partial Correlation - Ungrouped Data
Statistic	Value
Correlation r(xy)	-0.0206031053034193
Partial Correlation r(xy.z)	-0.87861346574556
Correlation r(xz)	0.792846709181054
Partial Correlation r(xz.y)	0.95669845135911
Correlation r(yz)	0.541693279998651
Partial Correlation r(yz.x)	0.915863999673108

\begin{tabular}{lllllllll}
\hline
Pearson Product Moment Partial Correlation - Ungrouped Data \tabularnewline
Statistic & Value \tabularnewline
Correlation r(xy) & -0.0206031053034193 \tabularnewline
Partial Correlation r(xy.z) & -0.87861346574556 \tabularnewline
Correlation r(xz) & 0.792846709181054 \tabularnewline
Partial Correlation r(xz.y) & 0.95669845135911 \tabularnewline
Correlation r(yz) & 0.541693279998651 \tabularnewline
Partial Correlation r(yz.x) & 0.915863999673108 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=24366&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.0206031053034193[/C][/ROW]
[ROW][C]Partial Correlation r(xy.z)[/C][C]-0.87861346574556[/C][/ROW]
[ROW][C]Correlation r(xz)[/C][C]0.792846709181054[/C][/ROW]
[ROW][C]Partial Correlation r(xz.y)[/C][C]0.95669845135911[/C][/ROW]
[ROW][C]Correlation r(yz)[/C][C]0.541693279998651[/C][/ROW]
[ROW][C]Partial Correlation r(yz.x)[/C][C]0.915863999673108[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=24366&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=24366&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.0206031053034193
Partial Correlation r(xy.z)	-0.87861346574556
Correlation r(xz)	0.792846709181054
Partial Correlation r(xz.y)	0.95669845135911
Correlation r(yz)	0.541693279998651
Partial Correlation r(yz.x)	0.915863999673108

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