R version 2.7.2 (2008-08-25)
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(2468.9,2469.2,2417.7,2411.1,2361,1924.1,2486.9,2674.9,2296.2,2101.5,2322,2273.8,2501.3,2435.2,2273.3,2454.7,2328.8,1897,2608.1,2712.3,2322,2282.6,2241.1,2417.2,2829,2600.6,2321,2768.5,2457.3,2142.7,2764.4,2788.9,2679.5,2536.5,2682.7,2699.6,3097.8,3015.2,2878,3010.9,2612.3,2419.3,3096.5,3013,3397.4,3423.1,3298.7,3065.7,3918.3,3154.4,3334.7,3461.6,3018.5,2832,3301.3,3342.8,3464.4,3016.6,3201.3,3135.3,3549.8,3247.2,3441.8,3535.6,3384.7,2996.6,3402.8,3900.2,3776.4,3197.5,4022.4,3845.1,3818.6)
> y <- c(12710.3,12120.8,12469.5,12054.6,12112.9,9617.2,12645.8,13581.3,12162.3,10969.7,11880,11887.6,12926.9,12300,12092.8,12380.8,12196.9,9455,13168,13427.9,11980.5,11884.8,11691.7,12233.8,14341.4,13130.7,12421.1,14285.8,12864.6,11160.2,14316.2,14388.7,14013.9,13419,12769.6,13315.5,15332.9,14243,13824.4,14962.9,13202.9,12199,15508.9,14199.8,15169.6,14058,13786.2,14147.9,16541.7,13587.5,15582.4,15802.8,14130.5,12923.2,15612.2,16033.7,16036.6,14037.8,15330.6,15038.3,17401.8,14992.5,16043.7,16929.6,15921.3,14417.2,15961,17851.9,16483.9,14215.5,17429.7,17839.5,17629.2)
> x <- c(680.2,668.9,611.4,640.8,549.8,541.8,628.6,686.5,611.8,588.4,566.9,563.7,569.9,635.4,590.8,634.3,576.1,351.6,507.5,586.2,666.4,693.6,650.6,654.8,733.5,648.1,678.1,816.2,591,563.5,742.5,694.4,728.6,749,538.9,568.5,692.8,580.5,506.9,612.8,442.9,523.3,596.7,533.7,523.1,559.2,430.7,538.2,612.4,428,522.4,531.1,425.9,410.3,551,555.6,460.2,288.9,392.3,400.5,399,354.9,337.6,379.2,334.1,321.6,449.8,486.3,421.9,405.6,420,432.4,418.1)
> #'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.3689550
> (rho23 <- cor(y, z))
[1] 0.9552745
> (rho13 <- cor(x, z))
[1] -0.4827754
> (rhoxy_z <- (rho12-(rho13*rho23))/(sqrt(1-(rho13*rho13)) * sqrt(1-(rho23*rho23))))
[1] 0.356127
> (rhoxz_y <- (rho13-(rho12*rho23))/(sqrt(1-(rho12*rho12)) * sqrt(1-(rho23*rho23))))
[1] -0.4741462
> (rhoyz_x <- (rho23-(rho12*rho13))/(sqrt(1-(rho12*rho12)) * sqrt(1-(rho13*rho13))))
[1] 0.9547816
> 
> #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/1vl1w1226415854.tab") 
> 
> 
> 
> proc.time()
   user  system elapsed 
  0.426   0.034   0.454