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:36:58 -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/t1226417949zksof4blplo4ev2.htm/, Retrieved Sun, 19 May 2024 10:49:16 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=23606, Retrieved Sun, 19 May 2024 10:49:16 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywordsPartial correlation
Estimated Impact127
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-11 15:36:58] [0cdfeda4aa2f9e551c2e529c44a404df] [Current]
Feedback Forum
2008-11-17 07:59:20 [Laura Reussens] [reply
De correlatie tussen de 3 variabelen (xy,xz en yz) situeren zich allemaal rond de 0.6 wat op zich niet echt een sterk verband weergeeft. Bij de partial correlation wordt de controlevariable weggewerkt omdat deze een vertekenend beeld kan hebben op het verband tussen de variabelen.
Deze controlevariabele had het minste effect op het verband tussen xz, waar y werd weggewerkt.
2008-12-05 16:50:31 [a2386b643d711541400692649981f2dc] [reply
Je interpretatie is hier fout. Het doel is om na te gaan of er een verband is tussen variabele x en y maar bij partiele correlatie gaan we na of er een variabele z is, die we niet observeren, maar wel een invloed heeft op de variabele x en y. We gaan dus na of er een schijncorrelatie wordt veroorzaakt. Het verband tussen (xy) bedraagt bij jou 68% en de als we het effect van z op x en y willen wegwerken bereken we de correlatie tussen (xy,z) die dan maar 45% bedraagt. Dit wijst erop dat z een invloed had op het verband tussen x en y.

Post a new message
Dataseries X:
103,1
100,6
103,1
95,5
90,5
90,9
88,8
90,7
94,3
104,6
111,1
110,8
107,2
99
99
91
96,2
96,9
96,2
100,1
99
115,4
106,9
107,1
99,3
99,2
108,3
105,6
99,5
107,4
93,1
88,1
110,7
113,1
99,6
93,6
98,6
99,6
114,3
107,8
101,2
112,5
100,5
93,9
116,2
112
106,4
95,7
96
95,8
103
102,2
98,4
111,4
86,6
91,3
107,9
101,8
104,4
93,4
100,1
98,5
112,9
101,4
107,1
110,8
90,3
95,5
111,4
113
107,5
95,9
106,3
105,2
117,2
106,9
108,2
113
97,2
99,9
108,1
118,1
109,1
93,3
112,1
Dataseries Y:
119,5
125
145
105,3
116,9
120,1
88,9
78,4
114,6
113,3
117
99,6
99,4
101,9
115,2
108,5
113,8
121
92,2
90,2
101,5
126,6
93,9
89,8
93,4
101,5
110,4
105,9
108,4
113,9
86,1
69,4
101,2
100,5
98
106,6
90,1
96,9
125,9
112
100
123,9
79,8
83,4
113,6
112,9
104
109,9
99
106,3
128,9
111,1
102,9
130
87
87,5
117,6
103,4
110,8
112,6
102,5
112,4
135,6
105,1
127,7
137
91
90,5
122,4
123,3
124,3
120
118,1
119
142,7
123,6
129,6
151,6
110,4
99,2
130,5
136,2
129,7
128
121,6
Dataseries Z:
98,6
98
106,8
96,6
100,1
107,7
91,5
97,8
107,4
117,5
105,6
97,4
99,5
98
104,3
100,6
101,1
103,9
96,9
95,5
108,4
117
103,8
100,8
110,6
104
112,6
107,3
98,9
109,8
104,9
102,2
123,9
124,9
112,7
121,9
100,6
104,3
120,4
107,5
102,9
125,6
107,5
108,8
128,4
121,1
119,5
128,7
108,7
105,5
119,8
111,3
110,6
120,1
97,5
107,7
127,3
117,2
119,8
116,2
111
112,4
130,6
109,1
118,8
123,9
101,6
112,8
128
129,6
125,8
119,5
115,7
113,6
129,7
112
116,8
127
112,1
114,2
121,1
131,6
125
120,4
117,7




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time1 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 & 1 seconds \tabularnewline
R Server & 'Sir Ronald Aylmer Fisher' @ 193.190.124.24 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=23606&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]'Sir Ronald Aylmer Fisher' @ 193.190.124.24[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=23606&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=23606&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'Sir Ronald Aylmer Fisher' @ 193.190.124.24







Pearson Product Moment Partial Correlation - Ungrouped Data
StatisticValue
Correlation r(xy)0.580471399357484
Partial Correlation r(xy.z)0.338156687220048
Correlation r(xz)0.61998535500533
Partial Correlation r(xz.y)0.421555846971527
Correlation r(yz)0.591075502587546
Partial Correlation r(yz.x)0.361861590431523

\begin{tabular}{lllllllll}
\hline
Pearson Product Moment Partial Correlation - Ungrouped Data \tabularnewline
Statistic & Value \tabularnewline
Correlation r(xy) & 0.580471399357484 \tabularnewline
Partial Correlation r(xy.z) & 0.338156687220048 \tabularnewline
Correlation r(xz) & 0.61998535500533 \tabularnewline
Partial Correlation r(xz.y) & 0.421555846971527 \tabularnewline
Correlation r(yz) & 0.591075502587546 \tabularnewline
Partial Correlation r(yz.x) & 0.361861590431523 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=23606&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.580471399357484[/C][/ROW]
[ROW][C]Partial Correlation r(xy.z)[/C][C]0.338156687220048[/C][/ROW]
[ROW][C]Correlation r(xz)[/C][C]0.61998535500533[/C][/ROW]
[ROW][C]Partial Correlation r(xz.y)[/C][C]0.421555846971527[/C][/ROW]
[ROW][C]Correlation r(yz)[/C][C]0.591075502587546[/C][/ROW]
[ROW][C]Partial Correlation r(yz.x)[/C][C]0.361861590431523[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=23606&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=23606&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.580471399357484
Partial Correlation r(xy.z)0.338156687220048
Correlation r(xz)0.61998535500533
Partial Correlation r(xz.y)0.421555846971527
Correlation r(yz)0.591075502587546
Partial Correlation r(yz.x)0.361861590431523



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