library(MASS) library(car) chii <- x chi_df <- fitdistr(chii,'chi-squared',start=list(df=3), method='Brent',lower=0.1,upper=100) chi_k <- chi_df[[1]][1] bitmap(file='test1.png') chi_hist <- hist(chii,breaks=50,col='gray',freq=FALSE,main='Histogram scaled to density',xlab='x') curve(dchisq(x,df=chi_k),add=TRUE,col='red') dev.off() f<-chi_df f$estimate f$sd xlab <- paste('Chisq(df = ',round(f$estimate[[1]],2),sep='') xlab <- paste(xlab,')',sep='') xlab bitmap(file='test2.png') qqPlot(x, dist='chisq', df=f$estimate[[1]], ncp=0, main='QQ plot (Chi-squared 1 param.)', xlab=xlab ) 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,'df',header=TRUE) a<-table.element(a,f$estimate[1]) a<-table.element(a,f$sd[1]) a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable.tab')
|