par2 <- as.numeric(par2) par3 <- as.numeric(par3) par4 <- as.numeric(par4) if (par5 == 'TRUE') par5 <- TRUE else par5 <- FALSE xc <- x if (par1 == 'TRUE') { xc <- cut(x, quantile(x, probs=seq(par2, par3, par4), na.rm=par5)) } myt <- as.matrix(table(xc)) myt cmyt <- cumsum(myt) smyt <- sum(myt) breaks <- round((par3-par2)/par4) bitmap(file='table.png') hist(x,breaks=breaks) dev.off() bitmap(file='table1.png') hist(cumsum(x),breaks=breaks) dev.off() load(file='createtable') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Frequency Table 1',5,TRUE) a<-table.row.end(a) a<-table.row.start(a) if (par1 == 'TRUE') { a<-table.element(a,'Intervals',header=TRUE) } else { a<-table.element(a,'Level',header=TRUE) } a<-table.element(a,'Non-cumul.',2,TRUE) a<-table.element(a,'Cumulative',2,TRUE) a<-table.row.end(a) a<-table.row.start(a) if (par1 == 'TRUE') { a<-table.element(a,'Lower,Upper',header=TRUE) } else { a<-table.element(a,'',header=TRUE) } a<-table.element(a,'Abs.',1,TRUE) a<-table.element(a,'Rel.',1,TRUE) a<-table.element(a,'Abs.',1,TRUE) a<-table.element(a,'Rel.',1,TRUE) a<-table.row.end(a) for (i in 1:length(myt)) { a<-table.row.start(a) a<-table.element(a,rownames(myt)[i],header=TRUE) a<-table.element(a,myt[i]) a<-table.element(a,round(myt[i]/smyt,4)) a<-table.element(a,cmyt[i]) a<-table.element(a,round(cmyt[i]/smyt,4)) a<-table.row.end(a) } a<-table.end(a) table.save(a,file='mytable.tab')
|