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(1915,1843,1761,2858,3968,5061,4661,4269,3857,3568,3274,2987,1683,1381,1071,2772,4485,6181,5479,4782,4067,3489,2903,2330,1736,1483,1242,2334,3423,4523,3986,3462,2908,2575,2237,1904,1610,1251,941,2450,3946,5409,4741,4069,3539,3189,2960,2704,1697,1598,1456,2316,3083,4158,3469,2892,2578,2233,1947,2049) > y <- c(10144,10751,11752,13808,16203,17432,18014,16956,17982,19435,19990,20154,10327,9807,10862,13743,16458,18466,18810,17361,17411,18517,18525,17859,9499,9490,9255,10758,12375,14617,15427,14136,14308,15293,15679,16319,11196,11169,12158,14251,16237,19706,18960,18537,19103,19691,19464,17264,8957,9703,9166,9519,10535,11526,9630,7061,6021,4728,2657,1264) > x <- c(285,10144,10751,11752,13808,16203,17432,18014,16956,17982,19435,19990,20154,10327,9807,10862,13743,16458,18466,18810,17361,17411,18517,18525,17859,9499,9490,9255,10758,12375,14617,15427,14136,14308,15293,15679,16319,11196,11169,12158,14251,16237,19706,18960,18537,19103,19691,19464,17264,8957,9703,9166,9519,10535,11526,9630,7061,6021,4728,2657) > #'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.817213 > (rho23 <- cor(y, z)) [1] 0.6272954 > (rho13 <- cor(x, z)) [1] 0.4557621 > (rhoxy_z <- (rho12-(rho13*rho23))/(sqrt(1-(rho13*rho13)) * sqrt(1-(rho23*rho23)))) [1] 0.7664738 > (rhoxz_y <- (rho13-(rho12*rho23))/(sqrt(1-(rho12*rho12)) * sqrt(1-(rho23*rho23)))) [1] -0.1267086 > (rhoyz_x <- (rho23-(rho12*rho13))/(sqrt(1-(rho12*rho12)) * sqrt(1-(rho13*rho13)))) [1] 0.496768 > > #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/1z9cb1258126661.tab") > > > > proc.time() user system elapsed 0.181 0.055 0.233