R version 2.7.0 (2008-04-22)
Copyright (C) 2008 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(53.3,66.9,73.7,71.9,67.9,77.5,82,40.5,76.2,75,64,77.2,63.9,71.4,71.3,50,47.9,55.2,61.8,53.6,60,72.7,63.6,62.1,58.1,68.1,65.5,59.3,65.7,68.3,77.8,49.1,96.9,64.6,67,64.6,84.3,72.6,71.3,52.8,63.9,67.9,74.1,61.7,74.1,83.5,75.6,77.9,69.2,67.8,74.8,69.7,69.9,84.3,81.7,69.3,83.8,96,74.2,70.4)
> y <- c(64.5,57.9,72.5,82.9,62.7,66.6,70.9,56.3,81.2,67.8,57.2,62.8,51.2,66.5,71.5,67.8,56.9,90.6,64.1,65.7,77.9,84.2,81.7,65.9,76.1,80.1,74.3,86.5,70.9,91.5,93.7,77.5,100.3,85.7,95,73.7,87.6,87.6,92.9,80.7,74,87.7,68.2,117.4,100.1,97.7,101.1,87.6,83.2,85.9,118,87.3,79.5,93.2,99.2,66.6,83.9,100.4,73.5,60.8)
> x <- c(744.8,672.1,666.6,760.8,756,604.4,883.9,527.9,756.2,812.9,655.6,707.6,612.6,659.2,833.4,727.8,797.2,753,762,613.7,759.2,816.4,736.8,680.1,736.5,637.2,801.9,772.3,897.3,792.1,826.8,666.8,906.6,871.4,891,739.2,833.6,715.6,871.6,751.6,1005.5,681.2,837.3,674.7,806.3,860.2,689.8,691.6,682.6,800.1,1023.7,733.5,875.3,770.2,1005.7,982.3,742.9,974.2,822.3,773.2)
> #'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.4085766
> (rho23 <- cor(y, z))
[1] 0.4352770
> (rho13 <- cor(x, z))
[1] 0.4309867
> (rhoxy_z <- (rho12-(rho13*rho23))/(sqrt(1-(rho13*rho13)) * sqrt(1-(rho23*rho23))))
[1] 0.2720097
> (rhoxz_y <- (rho13-(rho12*rho23))/(sqrt(1-(rho12*rho12)) * sqrt(1-(rho23*rho23))))
[1] 0.3080636
> (rhoyz_x <- (rho23-(rho12*rho13))/(sqrt(1-(rho12*rho12)) * sqrt(1-(rho13*rho13))))
[1] 0.3146972
> 
> #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/1wo361226269459.tab") 
> 
> 
> 
> proc.time()
   user  system elapsed 
  0.819   0.029   0.844