par2 <- '' par1 <- '' library(MASS) library(car) par1 <- as.numeric(par1) if (par2 == '0') par2 = 'Sturges' else par2 <- as.numeric(par2) x <- as.ts(x) #otherwise the fitdistr function does not work properly r <- fitdistr(x,'normal') print(r) bitmap(file='test1.png') myhist<-hist(x,col=par1,breaks=par2,main=main,ylab=ylab,xlab=xlab,freq=F) curve(1/(r$estimate[2]*sqrt(2*pi))*exp(-1/2*((x-r$estimate[1])/r$estimate[2])^2),min(x),max(x),add=T) dev.off() bitmap(file='test3.png') qqPlot(x,dist='norm',main='QQ plot (Normal) with confidence intervals') grid() dev.off() load(file='createtable') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Parameter',1,TRUE) a<-table.element(a,'Estimated Value',1,TRUE) a<-table.element(a,'Standard Deviation',1,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'mean',header=TRUE) a<-table.element(a,r$estimate[1]) a<-table.element(a,r$sd[1]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'standard deviation',header=TRUE) a<-table.element(a,r$estimate[2]) a<-table.element(a,r$sd[2]) a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable.tab')
|