Free Statistics

of Irreproducible Research!

Author's title

Author*Unverified author*
R Software Modulerwasp_partialcorrelation.wasp
Title produced by softwarePartial Correlation
Date of computationTue, 11 Nov 2008 04:32:10 -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/t1226403199gw2gggx0v7mvmlz.htm/, Retrieved Mon, 20 May 2024 09:45:27 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=23332, Retrieved Mon, 20 May 2024 09:45:27 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact121
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
F       [Partial Correlation] [8.3 various EDA t...] [2008-11-11 11:32:10] [0cebda6bbc99948f606f5db2560512ab] [Current]
Feedback Forum
2008-11-15 14:24:38 [Philip Van Herck] [reply
In deze berkening zijn in dataset X i.p.v. waarnemingen, periodes weergegeven. Dit geeft natuurlijk een heel vertekend beeld en is niet de bedoeling van de partieële correlatie-coëfficient. U geeft in deze berekening dat u de correlatie gaat berekenen tussen de duurzame en niet-duurzame consumptiegoederen. Dit is goed om een gewone correlatiecoëfficient te berekenen maar wanneer er een partiële correlatie moet berekend worden moet er ook een derde, een controlevariabele aanwezig zijn. Het is dan de bedoeling om aan de hand van deze derde variabele de mogelijke correlatie tussen de eerste twee te testen en dan te besluiten of deze al dan niet aan het toeval te wijten is.
2008-11-24 13:53:42 [58d427c57bd46519a715a3a7fea6a80f] [reply
r(xy,z)waarbij de 3e dimensie (z) wordt weggelaten. Het effect van z wordt dus weggewerkt en daarna wordt de correlatie tussen x en y berekend. In dit geval is deze negatief. We kunnen er niet uit afleiden of z een sterke/zwakke invloed heeft op x en y.

Post a new message
Dataseries X:
31/08/2001
30/09/2001
31/10/2001
30/11/2001
31/12/2001
31/01/2002
28/02/2002
31/03/2002
30/04/2002
31/05/2002
30/06/2002
31/07/2002
31/08/2002
30/09/2002
31/10/2002
30/11/2002
31/12/2002
31/01/2003
28/02/2003
31/03/2003
30/04/2003
31/05/2003
30/06/2003
31/07/2003
31/08/2003
30/09/2003
31/10/2003
30/11/2003
31/12/2003
31/01/2004
29/02/2004
31/03/2004
30/04/2004
31/05/2004
30/06/2004
31/07/2004
31/08/2004
30/09/2004
31/10/2004
30/11/2004
31/12/2004
31/01/2005
28/02/2005
31/03/2005
30/04/2005
31/05/2005
30/06/2005
31/07/2005
31/08/2005
30/09/2005
31/10/2005
30/11/2005
31/12/2005
31/01/2006
28/02/2006
31/03/2006
30/04/2006
31/05/2006
30/06/2006
31/07/2006
31/08/2006
30/09/2006
31/10/2006
30/11/2006
31/12/2006
31/01/2007
28/02/2007
31/03/2007
30/04/2007
31/05/2007
30/06/2007
31/07/2007
31/08/2007
30/09/2007
31/10/2007
30/11/2007
31/12/2007
31/01/2008
29/02/2008
31/03/2008
30/04/2008
31/05/2008
30/06/2008
31/07/2008
31/08/2008
Dataseries Y:
75,4
105,5
112,3
102,5
93,5
86,7
95,2
103,8
97
95,5
101
67,5
64
106,7
100,6
101,2
93,1
84,2
85,8
91,8
92,4
80,3
79,7
62,5
57,1
100,8
100,7
86,2
83,2
71,7
77,5
89,8
80,3
78,7
93,8
57,6
60,6
91
85,3
77,4
77,3
68,3
69,9
81,7
75,1
69,9
84
54,3
60
89,9
77
85,3
77,6
69,2
75,5
85,7
72,2
79,9
85,3
52,2
61,2
82,4
85,4
78,2
70,2
70,2
69,3
77,5
66,1
69
79,2
56,2
63,3
77,8
92
78,1
65,1
71,1
70,9
72
81,9
70,6
72,5
65,1
61,1
Dataseries Z:
98,4
107,4
117,7
105,7
97,5
99,9
98,2
104,5
100,8
101,5
103,9
99,6
98,4
112,7
118,4
108,1
105,4
114,6
106,9
115,9
109,8
101,8
114,2
110,8
108,4
127,5
128,6
116,6
127,4
105
108,3
125
111,6
106,5
130,3
115
116,1
134
126,5
125,8
136,4
114,9
110,9
125,5
116,8
116,8
125,5
104,2
115,1
132,8
123,3
124,8
122
117,4
117,9
137,4
114,6
124,7
129,6
109,4
120,9
134,9
136,3
133,2
127,2
122,7
120,5
137,8
119,1
124,3
134,4
121,1
122,2
127,7
137,4
132,2
129,2
124,9
124,8
128,2
134,4
118,6
132,6
123,2
112,3




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time1 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=23332&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=23332&T=0

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







Pearson Product Moment Partial Correlation - Ungrouped Data
StatisticValue
Correlation r(xy)-0.117466255250806
Partial Correlation r(xy.z)-0.125805821971128
Correlation r(xz)-0.161947159189143
Partial Correlation r(xz.y)-0.168018149124476
Correlation r(yz)-0.040608644381433
Partial Correlation r(yz.x)-0.0608509556158277

\begin{tabular}{lllllllll}
\hline
Pearson Product Moment Partial Correlation - Ungrouped Data \tabularnewline
Statistic & Value \tabularnewline
Correlation r(xy) & -0.117466255250806 \tabularnewline
Partial Correlation r(xy.z) & -0.125805821971128 \tabularnewline
Correlation r(xz) & -0.161947159189143 \tabularnewline
Partial Correlation r(xz.y) & -0.168018149124476 \tabularnewline
Correlation r(yz) & -0.040608644381433 \tabularnewline
Partial Correlation r(yz.x) & -0.0608509556158277 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=23332&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.117466255250806[/C][/ROW]
[ROW][C]Partial Correlation r(xy.z)[/C][C]-0.125805821971128[/C][/ROW]
[ROW][C]Correlation r(xz)[/C][C]-0.161947159189143[/C][/ROW]
[ROW][C]Partial Correlation r(xz.y)[/C][C]-0.168018149124476[/C][/ROW]
[ROW][C]Correlation r(yz)[/C][C]-0.040608644381433[/C][/ROW]
[ROW][C]Partial Correlation r(yz.x)[/C][C]-0.0608509556158277[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=23332&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=23332&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.117466255250806
Partial Correlation r(xy.z)-0.125805821971128
Correlation r(xz)-0.161947159189143
Partial Correlation r(xz.y)-0.168018149124476
Correlation r(yz)-0.040608644381433
Partial Correlation r(yz.x)-0.0608509556158277



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