library(MASS) (f<-fitdistr(x, 'logistic')) xlab <- paste('Logistic(location=',round(f$estimate[[1]],2)) xlab <- paste(xlab,', scale=') xlab <- paste(xlab,round(f$estimate[[2]],2)) xlab <- paste(xlab,')') bitmap(file='test2.png') qqplot(qlogis(ppoints(x), location=f$estimate[[1]], scale=f$estimate[[2]]), x, main='QQ plot (Logistic)', 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,'location',header=TRUE) a<-table.element(a,f$estimate[1]) a<-table.element(a,f$sd[1]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'scale',header=TRUE) a<-table.element(a,f$estimate[2]) a<-table.element(a,f$sd[2]) a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable.tab')
|