par1 <- '10' x <- x[!is.na(x)] bitmap(file='test1.png') par1 <- as.numeric(par1) myhist<-hist(x, breaks=par1, col=2, xlab=xlab, main=main) dev.off() bitmap(file='test2.png') qqnorm(x) qqline(x, col=2) dev.off() bitmap(file='test3.png') plot(density(x), xlab=xlab, main=main) dev.off() load(file='createtable') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Descriptive Statistics',3,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'mean',header=TRUE) a<-table.element(a,signif(mean(x), 6)) a<-table.element(a,hyperlink('arithmetic_mean.htm','formula','click to see the formula')) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'standard deviation',header=TRUE) a<-table.element(a,signif(sd(x), 6)) a<-table.element(a,hyperlink('unbiased1.htm','formula','click to see the formula')) a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable.tab')
|