library(psych) xy <- cor(x,y,use = 'pairwise') xz <- cor(x,z,use = 'pairwise') yz <- cor(y,z,use = 'pairwise') nx <- length(na.omit(x)) ny <- length(na.omit(y)) nz <- length(na.omit(z)) nxy <- min(nx,ny) nxz <- min(nx,nz) if(par1=='paired') { r <- paired.r(xy,xz,yz,n=nxy,n2=nxz) } else { r <- paired.r(xy,xz,n=nxy,n2=nxz) } r load(file='createtable') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Testing difference between two Pearson Correlations',2,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Pearson Correlation between X and Y',header=TRUE) a<-table.element(a,xy) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Pearson Correlation between X and Z',header=TRUE) a<-table.element(a,xz) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Type of test',header=TRUE) a<-table.element(a,par1) a<-table.row.end(a) if(par1=='paired') { a<-table.row.start(a) a<-table.element(a,'t-Test for dependent correlations',header=TRUE) a<-table.element(a,r$t) a<-table.row.end(a) } else { a<-table.row.start(a) a<-table.element(a,'z-Test for independent correlations',header=TRUE) a<-table.element(a,r$z) a<-table.row.end(a) } a<-table.row.start(a) a<-table.element(a,'P-value (H0: r(xy) = r(xz))',header=TRUE) a<-table.element(a,r$p) a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable.tab')
|