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(91.70,74.44,100.00,74.99,76.27,105.90,74.70,71.60,78.82,88.27,76.11,82.06,78.91,84.36,99.41,93.89,109.62,106.51,79.30,94.51,99.20,90.57,85.54,109.18,109.16,96.00,116.23,93.41,102.27,106.25,87.98,90.11,95.45,100.93,78.97,88.10,90.21,80.36,97.45,80.10,78.50,97.20,81.90,87.20,72.55,105.26,80.90,72.95,83.70,91.92,83.98,87.03,78.21,91.31,74.55,72.25,83.90,76.61,78.10,91.20)
> y <- c(101.26,101.08,100.00,98.02,96.76,97.84,107.03,110.09,110.45,110.09,107.03,107.21,106.49,106.13,105.23,103.24,102.16,102.52,111.89,113.33,113.15,110.27,107.21,107.57,106.85,106.31,104.50,103.42,103.24,103.24,111.71,112.79,111.71,105.95,101.98,100.36,101.08,98.92,95.86,94.77,92.07,89.91,100.00,101.80,97.66,94.95,91.89,92.61,93.15,91.53,88.83,88.29,84.50,86.13,95.14,96.22,93.33,91.17,90.45,92.97)
> x <- c(77.21,81.52,100.00,87.41,64.23,71.41,51.22,48.64,55.06,63.65,54.83,40.83,83.69,74.18,85.44,79.01,67.92,80.67,53.28,54.31,57.99,60.25,58.05,34.12,95.26,89.41,106.63,83.90,83.21,75.88,58.09,54.10,57.26,67.79,56.72,34.51,91.31,80.10,97.17,78.26,77.23,82.16,65.17,55.56,59.97,75.42,59.51,38.84,93.18,91.74,93.80,101.36,79.60,83.25,64.49,54.04,62.03,70.43,49.78,35.72)
> #'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.3001647
> (rho23 <- cor(y, z))
[1] 0.2191293
> (rho13 <- cor(x, z))
[1] 0.3255677
> (rhoxy_z <- (rho12-(rho13*rho23))/(sqrt(1-(rho13*rho13)) * sqrt(1-(rho23*rho23))))
[1] -0.4026998
> (rhoxz_y <- (rho13-(rho12*rho23))/(sqrt(1-(rho12*rho12)) * sqrt(1-(rho23*rho23))))
[1] 0.4204802
> (rhoyz_x <- (rho23-(rho12*rho13))/(sqrt(1-(rho12*rho12)) * sqrt(1-(rho13*rho13))))
[1] 0.3513103
> 
> #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/1pzas1258230863.tab") 
> 
> 
> 
> proc.time()
   user  system elapsed 
  0.201   0.031   0.250