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(96.33,96.33,95.05,96.84,96.92,97.44,97.78,97.69,96.67,98.29,98.20,98.71,98.54,98.20,100.80,101.33,101.88,101.85,102.04,102.22,102.63,102.65,102.54,102.37,102.68,102.76,102.82,103.31,103.23,103.60,103.95,103.93,104.25,104.38,104.36,104.32,104.58,104.68,104.92,105.46,105.23,105.58,105.34,105.28,105.70,105.67,105.71,106.19,106.93,107.44,107.85,108.71,109.32,109.49,110.20,110.62,111.22,110.88,111.15,111.29,111.09,111.24,111.45,111.75,111.07,111.17,110.96,110.50,110.48,110.66) > y <- c(98.60,96.90,95.10,97.00,112.70,102.90,97.40,111.40,87.40,96.80,114.10,110.30,103.90,101.60,94.60,95.90,104.70,102.80,98.10,113.90,80.90,95.70,113.20,105.90,108.80,102.30,99.00,100.70,115.50,100.70,109.90,114.60,85.40,100.50,114.80,116.50,112.90,102.00,106.00,105.30,118.80,106.10,109.30,117.20,92.50,104.20,112.50,122.40,113.30,100.00,110.70,112.80,109.80,117.30,109.10,115.90,96.00,99.80,116.80,115.70,99.40,94.30,91.00,93.20,103.10,94.10,91.80,102.70,82.60,89.10) > x <- c(103.34,102.60,100.69,105.67,123.61,113.08,106.46,123.38,109.87,95.74,123.06,123.39,120.28,115.33,110.40,114.49,132.03,123.16,118.82,128.32,112.24,104.53,132.57,122.52,131.80,124.55,120.96,122.60,145.52,118.57,134.25,136.70,121.37,111.63,134.42,137.65,137.86,119.77,130.69,128.28,147.45,128.42,136.90,143.95,135.64,122.48,136.83,153.04,142.71,123.46,144.37,146.15,147.61,158.51,147.40,165.05,154.64,126.20,157.36,154.15,123.21,113.07,110.45,113.57,122.44,114.93,111.85,126.04,121.34,124.36) > #'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.7099101 > (rho23 <- cor(y, z)) [1] 0.03491106 > (rho13 <- cor(x, z)) [1] 0.5085509 > (rhoxy_z <- (rho12-(rho13*rho23))/(sqrt(1-(rho13*rho13)) * sqrt(1-(rho23*rho23)))) [1] 0.8043584 > (rhoxz_y <- (rho13-(rho12*rho23))/(sqrt(1-(rho12*rho12)) * sqrt(1-(rho23*rho23)))) [1] 0.687303 > (rhoyz_x <- (rho23-(rho12*rho13))/(sqrt(1-(rho12*rho12)) * sqrt(1-(rho13*rho13)))) [1] -0.5377714 > > #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/18jts1261146329.tab") > > > > proc.time() user system elapsed 0.202 0.032 0.252