k <- cor.test(x,y,method='spearman') bitmap(file='test1.png') plot(x,y,main='Scatterplot',xlab=xlab,ylab=ylab) grid() dev.off() bitmap(file='test2.png') plot(rank(x),rank(y),main='Scatterplot of Ranks',xlab=xlab,ylab=ylab) grid() dev.off() load(file='createtable') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Spearman Rank Correlation',2,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'rho',header=TRUE) a<-table.element(a,k$estimate) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'2-sided p-value',header=TRUE) a<-table.element(a,k$p.value) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'S',header=TRUE) a<-table.element(a,k$statistic) a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable.tab')
|