R version 2.9.0 (2009-04-17)
Copyright (C) 2009 The R Foundation for Statistical Computing
ISBN 3-900051-07-0

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> z <- c(46,46,41,43,39,36,34,31,27,27,30,28,31,25,26,22,27,38,39,45,48,46,47,55,63,83,100,76,68,63,63,62,62,64,72,67,67,67,69,78,80,85,90,83,84,84,89,82,76,97,82,79,73,72,74,77,78,68,71,68,70,79,78,79,78,77,72,72,84,83,82,78,77,75,77,73,64,64,71,65,68,64,65,56,57,51,52,79,59,54,55,47,52,52,54,49,48,44,47,43,45,40,38,46,51,49,57,61,61,81,97,105,108,113,111,111,110,109,105,97,97)
> y <- c(62,64,62,64,64,69,69,65,56,58,53,62,55,60,59,58,53,57,57,53,54,53,57,57,55,49,50,49,54,58,58,52,56,52,59,53,52,53,51,50,56,52,46,48,46,48,48,49,53,48,51,48,50,55,52,53,52,55,53,53,56,54,52,55,54,59,56,56,51,53,52,51,46,49,46,55,57,53,52,53,50,54,53,50,51,52,47,51,49,53,52,45,53,51,48,48,48,48,40,43,40,39,39,36,41,39,40,39,46,40,37,37,44,41,40,36,38,43,42,45,46)
> x <- c(33,39,45,46,45,45,49,50,54,59,58,56,48,50,52,53,55,43,42,38,41,41,39,34,27,15,14,31,41,43,46,42,45,45,40,35,36,38,39,32,24,21,12,29,36,31,28,30,38,27,40,40,44,47,45,42,38,46,37,41,40,33,34,36,36,38,42,35,25,24,22,27,17,30,30,34,37,36,33,33,33,37,40,35,37,43,42,33,39,40,37,44,42,43,40,30,30,31,18,24,22,26,28,23,17,12,9,19,21,18,18,15,24,18,19,30,33,35,36,47,46)
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> (rho12 <- cor(x, y))
[1] 0.659134
> (rho23 <- cor(y, z))
[1] -0.4632186
> (rho13 <- cor(x, z))
[1] -0.4949235
> (rhoxy_z <- (rho12-(rho13*rho23))/(sqrt(1-(rho13*rho13)) * sqrt(1-(rho23*rho23))))
[1] 0.5582157
> (rhoxz_y <- (rho13-(rho12*rho23))/(sqrt(1-(rho12*rho12)) * sqrt(1-(rho23*rho23))))
[1] -0.2844811
> (rhoyz_x <- (rho23-(rho12*rho13))/(sqrt(1-(rho12*rho12)) * sqrt(1-(rho13*rho13))))
[1] -0.2096488
> 
> #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/html/rcomp/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('http://www.xycoon.com/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('http://www.xycoon.com/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('http://www.xycoon.com/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="/var/www/html/rcomp/tmp/1dbgv1257364385.tab") 
> 
> 
> 
> proc.time()
   user  system elapsed 
  0.197   0.032   0.408