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(100,79.53631927,77.14613839,94.99712891,98.11943727,85.46511628,82.48636233,92.01119724,100.97617,99.94257824,120.1909274,109.804766,84.76169968,92.07579673,109.4028137,109.4315245,109.1085271,135.0775194,98.70801034,95.93023256,118.712317,100.796727,124.2750502,95.5211025,116.0134941,107.931381,105.9933965,112.3743899,113.5371806,101.0407695,100.0358886,119.3224232,138.4582257,97.10020098,96.16709733,89.80763709,107.0700546,110.5871375,99.31811657,131.1728395,114.219064,92.53517083,92.736147,98.21274763,91.06373816,104.3927649,89.89376974,93.69796153,91.50875682,78.22279644,114.9511915,101.8446741,86.76428366,82.23514212,81.04364054,77.06000574,87.80505312,81.50301464,78.32328452,80.12489233) > y <- c(100,96.21064363,96.31280765,107.1793443,114.9066592,92.56060184,114.9995356,107.1236185,117.7765394,107.3650971,106.2970187,114.5072908,98.0031578,103.0649206,100.2879168,104.6066685,111.1544534,104.9874617,109.9284852,111.5352466,132.4974459,100.3436426,123.0983561,114.2379493,104.569518,109.0833101,106.9843039,133.6769759,124.8537197,122.5132349,116.8013374,116.0118882,129.7575926,125.1973623,143.7912139,127.9465032,130.2962757,108.4424631,129.3675118,143.6797622,131.8844618,117.6186496,118.9560695,104.8202842,134.624315,140.401226,143.8005015,153.4317823,153.2924677,127.3149438,153.5525216,136.9276493,131.7730101,144.3391845,107.4208229,113.6249652,124.2221603,102.0618557,96.36853348,111.6838488) > x <- c(100,97.82226485,94.04971502,91.12460521,93.13202153,93.88342812,92.55349954,94.43494835,96.25017563,100.4355715,101.5036685,99.39789728,99.68990733,101.6895041,103.6652759,103.0532766,100.9500712,102.345366,101.6472299,99.56809393,95.67727392,96.58494865,96.32604937,95.37109101,96.00056203,96.88367859,94.85280372,92.46943974,93.99180173,93.45262168,92.26698759,90.39653498,90.43001228,91.04995327,89.07845784,89.69314509,87.92459054,85.8789319,83.20612366,83.85722053,83.01393462,82.84508195,78.68864276,77.56959675,78.53689529,78.55717715,77.4761291,81.58931659,85.02428326,91.71290159,95.96293061,90.84689022,92.28788036,95.56511274,93.62452884,92.63071726,89.50914211,87.17171779,86.72624975,85.63212844) > #'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.4882211 > (rho23 <- cor(y, z)) [1] 0.1757351 > (rho13 <- cor(x, z)) [1] 0.1832597 > (rhoxy_z <- (rho12-(rho13*rho23))/(sqrt(1-(rho13*rho13)) * sqrt(1-(rho23*rho23)))) [1] -0.5377606 > (rhoxz_y <- (rho13-(rho12*rho23))/(sqrt(1-(rho12*rho12)) * sqrt(1-(rho23*rho23)))) [1] 0.3131711 > (rhoyz_x <- (rho23-(rho12*rho13))/(sqrt(1-(rho12*rho12)) * sqrt(1-(rho13*rho13)))) [1] 0.3091199 > > #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/1klec1258721532.tab") > > > > proc.time() user system elapsed 0.205 0.030 0.231