m <- lm(y ~ x) sm <- summary(m) bitmap(file='test1.png') myscatter<-plot(x,y,main=main,xlab=xlab,ylab=ylab) abline(m) grid() dev.off() load(file='createtable') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Y = alpha + beta X',2,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Intercept (alpha)',header=TRUE) a<-table.element(a,round(m$coef[[1]],4)) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Slope (beta)',header=TRUE) a<-table.element(a,round(m$coef[[2]],4)) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'S.D. (beta)',header=TRUE) a<-table.element(a,round(sm$coef[2,2],4)) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'T-Stat (beta)',header=TRUE) a<-table.element(a,round(sm$coef[2,3],4)) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'P-Value (beta)',header=TRUE) a<-table.element(a,round(sm$coef[2,4],4)) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Correlation',header=TRUE) a<-table.element(a,round(sqrt(sm$r.sq),4)) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'R-squared',header=TRUE) a<-table.element(a,round(sm$r.sq,4)) a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable.tab')
|