(rho12 <- cor(x, y)) (rho23 <- cor(y, z)) (rho13 <- cor(x, z)) (rhoxy_z <- (rho12-(rho13*rho23))/(sqrt(1-(rho13*rho13)) * sqrt(1-(rho23*rho23)))) (rhoxz_y <- (rho13-(rho12*rho23))/(sqrt(1-(rho12*rho12)) * sqrt(1-(rho23*rho23)))) (rhoyz_x <- (rho23-(rho12*rho13))/(sqrt(1-(rho12*rho12)) * sqrt(1-(rho13*rho13)))) load(file='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('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('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('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='mytable.tab')
|