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(562000,561000,555000,544000,537000,543000,594000,611000,613000,611000,594000,595000,591000,589000,584000,573000,567000,569000,621000,629000,628000,612000,595000,597000,593000,590000,580000,574000,573000,573000,620000,626000,620000,588000,566000,557000,561000,549000,532000,526000,511000,499000,555000,565000,542000,527000,510000,514000,517000,508000,493000,490000,469000,478000,528000,534000,518000,506000,502000,516000) > y <- c(52.034,58.146,79.035,63.969,34.831,42.602,23.589,21.974,29.651,39.547,27.957,10.511,63.044,50.898,63.547,55.028,40.073,56.080,28.091,29.211,33.415,35.336,31.502,1.007,75.840,69.091,88.380,61.147,59.750,50.533,33.467,28.944,32.053,42.095,27.840,-686,69.362,55.107,73.801,50.981,48.575,52.660,37.173,25.894,30.201,46.248,26.577,1.911,68.401,65.458,67.183,76.063,48.064,52.591,34.502,22.268,30.174,39.902,14.115,-2.715) > x <- c(4651.600942,4648.46719,4629.664677,4595.193403,4573.257137,4592.05965,4751.881013,4805.154801,4811.422306,4805.154801,4751.881013,4755.014766,4742.479757,4736.212252,4720.543491,4686.072217,4667.269704,4673.537208,4836.492323,4861.562341,4858.428589,4808.288553,4755.014766,4761.28227,4748.747261,4739.346005,4708.008482,4689.205969,4686.072217,4686.072217,4833.358571,4852.161084,4833.358571,4733.0785,4664.135951,4635.932181,4648.46719,4610.862164,4557.588376,4538.785863,4491.779579,4454.174553,4629.664677,4661.002199,4588.925898,4541.919615,4488.645827,4501.180836,4510.582093,4482.378323,4435.37204,4425.970783,4360.161986,4388.365756,4545.053367,4563.85588,4513.715845,4476.110818,4463.57581,4507.448341) > #'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.02640968 > (rho23 <- cor(y, z)) [1] -0.02640968 > (rho13 <- cor(x, z)) [1] 1 > (rhoxy_z <- (rho12-(rho13*rho23))/(sqrt(1-(rho13*rho13)) * sqrt(1-(rho23*rho23)))) [1] 0.01748829 > (rhoxz_y <- (rho13-(rho12*rho23))/(sqrt(1-(rho12*rho12)) * sqrt(1-(rho23*rho23)))) [1] 1 > (rhoyz_x <- (rho23-(rho12*rho13))/(sqrt(1-(rho12*rho12)) * sqrt(1-(rho13*rho13)))) [1] -0.01748829 > > #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/1485x1257355436.tab") > > > > proc.time() user system elapsed 0.195 0.026 0.244