R version 2.8.1 (2008-12-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(100,85.95702149,68.015992,91.25437281,104.2978511,71.11444278,94.35282359,71.91404298,74.8125937,98.4007996,115.6921539,109.5452274,110.2948526,122.0389805,91.2043978,114.8425787,118.9905047,103.1984008,124.5877061,90.85457271,108.9455272,123.1384308,107.0964518,121.0894553,110.2948526,102.1989005,90.30484758,116.4417791,118.1909045,119.6401799,104.2978511,104.4477761,123.5382309,142.0789605,122.8385807,124.4877561,125.6371814,125.5372314,103.5482259,114.0929535,127.0864568,108.8955522,122.1389305,116.5417291,106.2468766,119.6901549,167.6661669,124.3378311,132.2338831,137.6311844,90.30484758,127.986007,123.6381809,113.8430785,123.988006,76.8115942,112.6936532,137.1314343,200.049975,150.8245877,172.5137431) > y <- c(100,87.14054095,112.0054296,112.312101,109.474134,104.9746116,100.4926851,104.2154743,120.1768388,112.1028355,108.1481575,116.802197,102.1699512,95.15358705,120.6707808,111.5234277,119.9669448,113.3697401,110.0717661,111.5567342,132.424212,107.900558,122.1626615,124.3992258,110.4450505,101.5874013,122.3203962,125.2582826,125.4411543,108.9902468,118.9243879,116.7242723,134.1724901,116.8530994,124.5732995,130.9914031,123.4239103,111.4536725,124.5135991,139.2589613,129.8596099,112.3460359,131.381655,133.0004776,134.3220552,144.2379719,134.1278719,150.1891559,140.722563,114.8389975,143.1973003,140.2738676,112.1248303,102.8951536,100.5090242,103.3513901,111.4134533,104.5887587,101.7840983,114.7007441,108.7426474) > x <- c(100,99.31916948,114.7106298,114.9941367,114.6085397,116.3633855,108.0492516,112.8502449,127.2518452,116.4999655,113.8511416,126.3771815,103.8290681,108.296889,125.0969166,115.6701386,130.2545354,131.2064565,124.8665241,122.6191626,148.2368766,117.7160792,128.9232255,130.7380837,113.1192661,116.5454922,127.9782024,132.9240533,133.9732358,121.2623301,128.559702,124.692695,140.9850314,124.6823481,134.9451611,137.2642616,132.6308891,122.7122853,132.5846727,151.5410085,145.764641,134.4043595,152.9192247,143.9387459,152.9709595,162.3356557,148.1658274,167.9499207,156.9400566,140.2241843,163.13996,160.969166,135.3031662,124.4554046,119.2143202,116.21163,125.2638477,114.0987791,110.0413879,126.6317169,113.0875354) > #'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.9350637 > (rho23 <- cor(y, z)) [1] 0.1473842 > (rho13 <- cor(x, z)) [1] 0.1848110 > (rhoxy_z <- (rho12-(rho13*rho23))/(sqrt(1-(rho13*rho13)) * sqrt(1-(rho23*rho23)))) [1] 0.933937 > (rhoxz_y <- (rho13-(rho12*rho23))/(sqrt(1-(rho12*rho12)) * sqrt(1-(rho23*rho23)))) [1] 0.1340452 > (rhoyz_x <- (rho23-(rho12*rho13))/(sqrt(1-(rho12*rho12)) * sqrt(1-(rho13*rho13)))) [1] -0.07298453 > > #Note: the /var/www/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/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/rcomp/tmp/13tml1257948174.tab") > > > proc.time() user system elapsed 0.16 0.17 0.31