FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

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 12:49:25 -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/t1226433014hgvx93rb8gya7p1.htm/, Retrieved Mon, 27 May 2024 05:45:26 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=23903, Retrieved Mon, 27 May 2024 05:45:26 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact192

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
F     [Univariate Data Series] [Niet werkende wer...] [2008-10-13 17:30:07] [91d2608132ba5d00ecce3524d8276757]
F RMPD    [Partial Correlation] [Partial Correlati...] [2008-11-11 19:49:25] [5e9e099b83e50415d7642e10d74756e4] [Current]
F    D      [Partial Correlation] [Partial Correlati...] [2008-11-11 19:56:04] [deb3c14ac9e4607a6d84fc9d0e0e6cc2]
F RMPD        [Hierarchical Clustering] [Hierarchical Clus...] [2008-11-11 20:56:26] [deb3c14ac9e4607a6d84fc9d0e0e6cc2]
F RM D        [Box-Cox Linearity Plot] [Box-Cox Linearity...] [2008-11-11 21:16:20] [deb3c14ac9e4607a6d84fc9d0e0e6cc2]
F RM D          [Box-Cox Normality Plot] [Box-Cox Normality...] [2008-11-11 21:26:55] [deb3c14ac9e4607a6d84fc9d0e0e6cc2]
F RM              [Maximum-likelihood Fitting - Normal Distribution] [Maximum-likelihoo...] [2008-11-11 21:28:51] [deb3c14ac9e4607a6d84fc9d0e0e6cc2]
F   PD              [Maximum-likelihood Fitting - Normal Distribution] [] [2008-11-12 13:54:54] [8d2ae74f923b31b35e9e42977c3c4399]
-   PD              [Maximum-likelihood Fitting - Normal Distribution] [Various EDA topic...] [2008-11-18 20:01:30] [82970caad4b026be9dd352fdec547fe4]
- RM D        [Box-Cox Linearity Plot] [Box-Cox Linearity...] [2008-11-11 21:18:29] [deb3c14ac9e4607a6d84fc9d0e0e6cc2]
Feedback Forum
2008-11-13 15:57:33 [Thomas Beyers] [reply] 
Reeks Omschrijving Code Partial Correlation
1 Het aantal niet werkende werkzoekenden in België. X Partial Correlation r(xy.z)
0.929732112210466
3 Het aantal vrouwelijke niet werkende werkzoekenden in het Vlaams gewest. Y Partial Correlation r(xz.y)
0.904328686919826
4 Het aantal vrouwelijke niet werkende werkzoekenden in het Waals gewest. Z Partial Correlation r(yz.x)
-0.707689052209976
2008-11-17 19:39:15 [Ruben Jacobs] [reply] 
De berekening is correct maar de interpretatie is verkeerd. Bij Partiële Correlatie ga je 3 tijdreeksen en dus 3 variabelen met elkaar vergelijken. Hierbij ga je dan de invloed van 1 variabele uitzuiveren; deze kan namelijk het beeld van de correlatie tussen andere 2 andere variabelen vertekenen.
Zo is bijvoorbeeld in jouw tabel de correlatie tussen y en z 0.8467. Maar als je de invloed van x op deze y en z eerst gaat uitzuiveren (de partial correlation). Dan kan je zien dat er zelfs een negatieve correlatie is van -0.7076 tussen y en z.

Bij de correlatie van x en y kan je hetzelfde doen, de gewone correlatie bedraagt hier 0.9607. Als je de invloed van z op x en y uitzuivert dan zie je dat er nog steeds een goede correlatie bestaat van 0.9297.

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
577992
565464
547344
554788
562325
560854
555332
543599
536662
542722
593530
610763
612613
611324
594167
595454
590865
589379
584428
573100
567456
569028
620735
628884
628232
612117
595404
597141
593408
590072
579799
574205
572775
572942
619567
625809
619916
587625
565742
557274
560576
548854
531673
525919
511038
498662
555362
564591
541657
527070
509846
514258
516922
507561
492622
490243
469357
477580
528379
533590
517945
Dataseries Y:
Download CSV
Histogram 
121148
114624
109822
112081
113534
112110
109826
107423
105540
108573
128591
139145
129700
132828
126868
128390
126830
124105
122323
119296
116822
119224
139357
144322
133676
128283
121640
122877
117284
116463
112685
113235
111692
113152
129889
131153
123770
112516
105940
104320
103582
99064
94989
92241
89752
90610
109456
110213
97694
91844
87572
89812
89050
85990
85070
83277
79586
84215
99708
100698
90861
Dataseries Z:
Download CSV
Histogram 
143119
140353
137149
138301
138709
138619
137528
134523
133593
135647
143541
144024
148607
148881
145447
144170
142148
142492
141703
139116
138461
139870
148703
150034
153670
147508
145150
143955
147608
144839
142096
138941
139093
141526
151492
155217
154530
148533
143923
141075
142261
140435
136169
136107
131063
127765
140731
145981
141834
137625
132515
131111
133603
132207
126488
125933
120927
122597
135437
138649
134576

Tables (Output of Computation)





Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time2 seconds
R Server'Herman Ole Andreas Wold' @ 193.190.124.10:1001

\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 & 'Herman Ole Andreas Wold' @ 193.190.124.10:1001 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=23903&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]'Herman Ole Andreas Wold' @ 193.190.124.10:1001[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=23903&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=23903&T=0

As an alternative you can also use a QR Code:  
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 time2 seconds
R Server'Herman Ole Andreas Wold' @ 193.190.124.10:1001






Pearson Product Moment Partial Correlation - Ungrouped Data
StatisticValue
Correlation r(xy)0.96079910027726
Partial Correlation r(xy.z)0.929732112210466
Correlation r(xz)0.9469551898897
Partial Correlation r(xz.y)0.904328686919826
Correlation r(yz)0.846780590345912
Partial Correlation r(yz.x)-0.707689052209976

\begin{tabular}{lllllllll}
\hline
Pearson Product Moment Partial Correlation - Ungrouped Data \tabularnewline
Statistic & Value \tabularnewline
Correlation r(xy) & 0.96079910027726 \tabularnewline
Partial Correlation r(xy.z) & 0.929732112210466 \tabularnewline
Correlation r(xz) & 0.9469551898897 \tabularnewline
Partial Correlation r(xz.y) & 0.904328686919826 \tabularnewline
Correlation r(yz) & 0.846780590345912 \tabularnewline
Partial Correlation r(yz.x) & -0.707689052209976 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=23903&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.96079910027726[/C][/ROW]
[ROW][C]Partial Correlation r(xy.z)[/C][C]0.929732112210466[/C][/ROW]
[ROW][C]Correlation r(xz)[/C][C]0.9469551898897[/C][/ROW]
[ROW][C]Partial Correlation r(xz.y)[/C][C]0.904328686919826[/C][/ROW]
[ROW][C]Correlation r(yz)[/C][C]0.846780590345912[/C][/ROW]
[ROW][C]Partial Correlation r(yz.x)[/C][C]-0.707689052209976[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=23903&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=23903&T=1

As an alternative you can also use a QR Code:  
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.96079910027726
Partial Correlation r(xy.z)0.929732112210466
Correlation r(xz)0.9469551898897
Partial Correlation r(xz.y)0.904328686919826
Correlation r(yz)0.846780590345912
Partial Correlation r(yz.x)-0.707689052209976


Figures (Output of Computation)

Input Parameters & R Code

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