R version 2.6.0 (2007-10-03)
Copyright (C) 2007 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(92.58,93.15,93.7,93.34,93.93,94.41,94.3,94.44,96.09,95.99,96.23,95.64,94.88,95.48,95.52,96.18,96.82,96.81,96.25,96.24,96.95,96.25,96.04,95.78,95.86,96.02,96.34,96.84,96.73,96.34,96.6,96.64,97.2,97.5,96.99,97.08,97.55,98.42,98.78,97.49,96.99,97.16,97.29,97.8,98.12,98.03,98.11,98.07,98.21,98.48,98.83,99.2,99.88,99.71,100.03,100.6,100.85,101.96,101.4,100.81,100.66,101.55,102.23,102.9,102.68,103.41,104.62,104.93,105.88,105.18,104.54,104.58,104.34,104.66,104.73,105.44,105.72,105.68,105.9,105.97,105.21,104.75,104.89,105.26,104.84,105.47,105.4,105.73,105.72,105.63,105.97,105.92,106.32)
> y <- c(89.28,89.47,89.53,90.72,90.91,91.38,91.49,90.9,90.93,90.57,91.28,90.83,91.5,91.58,92.49,94.16,95.46,95.8,95.32,95.41,95.35,95.68,95.59,94.96,96.92,96.06,96.59,96.67,97.27,96.38,96.47,96.05,96.76,96.51,96.55,95.97,97,97.46,97.9,98.42,98.54,99,98.94,99.02,100.07,98.72,98.73,98.04,99.08,99.22,99.57,100.44,100.84,100.75,100.49,99.98,99.96,99.76,100.11,99.79,100.29,101.12,102.65,102.71,103.39,102.8,102.07,102.15,101.21,101.27,101.86,101.65,101.94,102.62,102.71,103.39,104.51,104.09,104.29,104.57,105.39,105.15,106.13,105.46,106.47,106.62,106.52,108.04,107.15,107.32,107.76,107.26,107.89)
> x <- c(91.19,91.53,91.88,92.06,92.32,92.67,92.85,92.82,93.46,93.23,93.54,93.29,93.2,93.6,93.81,94.62,95.22,95.38,95.31,95.3,95.57,95.42,95.53,95.33,95.9,96.06,96.31,96.34,96.49,96.22,96.53,96.5,96.77,96.66,96.58,96.63,97.06,97.73,98.01,97.76,97.49,97.77,97.96,98.23,98.51,98.19,98.37,98.31,98.6,98.97,99.11,99.64,100.03,99.98,100.32,100.44,100.51,101,100.88,100.55,100.83,101.51,102.16,102.39,102.54,102.85,103.47,103.57,103.69,103.5,103.47,103.45,103.48,103.93,103.89,104.4,104.79,104.77,105.13,105.26,104.96,104.75,105.01,105.15,105.2,105.77,105.78,106.26,106.13,106.12,106.57,106.44,106.54)
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: Wessa P., (2007), Partial Correlation (v1.0.3) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_partialcorrelation.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.9788068
> (rho23 <- cor(y, z))
[1] 0.9311868
> (rho13 <- cor(x, z))
[1] 0.9834291
> (rhoxy_z <- (rho12-(rho13*rho23))/(sqrt(1-(rho13*rho13)) * sqrt(1-(rho23*rho23))))
[1] 0.9540256
> (rhoxz_y <- (rho13-(rho12*rho23))/(sqrt(1-(rho12*rho12)) * sqrt(1-(rho23*rho23))))
[1] 0.9641547
> (rhoyz_x <- (rho23-(rho12*rho13))/(sqrt(1-(rho12*rho12)) * sqrt(1-(rho13*rho13))))
[1] -0.8457718
> 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/1piko1194104950.tab") 
> 
> 
> 
> proc.time()
   user  system elapsed 
  0.438   0.027   0.461