library(MASS) library(gamlss.dist) par2 <- as.numeric(par2) #length par8 <- as.numeric(par8) #instances xarr <- array(NA,dim=c(par2,par8)) if (par1 == 'Poisson (lambda)') { if (par3=='') par3 <- '0.3' par3 <- as.numeric(par3) for(i in 1:par8) { xarr[,i] <- rpois(par2, lambda=par3) } x <- xarr[,1] } if (par1 == 'Normal (mu; sigma)') { if (par3=='') par3 <- '0' if (par4=='') par4 <- '1' par3 <- as.numeric(par3) par4 <- as.numeric(par4) for(i in 1:par8) { xarr[,i] <- rnorm(par2, mean=par3, sd=par4) } x <- xarr[,1] } if (par1 == 'Lognormal (mu; sigma)') { if (par3=='') par3 <- '0' if (par4=='') par4 <- '1' par3 <- as.numeric(par3) par4 <- as.numeric(par4) for(i in 1:par8) { xarr[,i] <- rlnorm(par2, mean=par3, sd=par4) } x <- xarr[,1] } if (par1 == 'Negative Binomial (mu; theta)') { if (par3=='') par3 <- '1' if (par4=='') par4 <- '4' par3 <- as.numeric(par3) par4 <- as.numeric(par4) for(i in 1:par8) { xarr[,i] <- rnegbin(par2, mu=par3, theta=par4) } x <- xarr[,1] } if (par1 == 'Zero-Inflated Poisson (mu; sigma)') { if (par3=='') par3 <- '5' if (par4=='') par4 <- '0.1' par3 <- as.numeric(par3) par4 <- as.numeric(par4) for(i in 1:par8) { xarr[,i] <- rZIP(par2, mu=par3, sigma=par4) } x <- xarr[,1] } if (par1 == 'Sichel - Generalized Inverse Gaussian (mu; sigma; nu)') { if (par3=='') par3 <- '0.5' if (par4=='') par4 <- '0.02' if (par5=='') par5 <- '-0.5' par3 <- as.numeric(par3) par4 <- as.numeric(par4) par5 <- as.numeric(par5) for(i in 1:par8) { xarr[,i] <- rSI(par2, mu=par3, sigma=par4, nu=par5) } x <- xarr[,1] } if (par1 == 'Sichel2 - Generalized Inverse Gaussian (mu; sigma; nu; max.value)') { if (par3=='') par3 <- '0.5' if (par4=='') par4 <- '0.02' if (par5=='') par5 <- '-0.5' if (par6=='') par6 <- '10000' par3 <- as.numeric(par3) par4 <- as.numeric(par4) par5 <- as.numeric(par5) par6 <- as.numeric(par6) for(i in 1:par8) { xarr[,i] <- rSICHEL(par2, mu=par3, sigma=par4, nu=par5, max.value=par6) } x <- xarr[,1] } x write.csv(xarr,file=paste('tmp/',par9,'.csv',sep='')) bitmap(file='test0.png') plot(x,type='b') dev.off() bitmap(file='test1.png') myhist<-hist(x,col=2) dev.off() load(file='createtable') a<-table.start() a<-table.row.start(a) a<-table.element(a,'What is next?',1,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,hyperlink(paste('https://automated.biganalytics.eu/Patrick.Wessa/rwasp_demand_forecasting_simulate.wasp#output',sep=''),'Simulate Time Series','',target='')) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,hyperlink(paste('https://automated.biganalytics.eu/Patrick.Wessa/rwasp_demand_forecasting_croston.wasp',sep=''),'Generate Forecasts','',target='')) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,hyperlink(paste('https://automated.biganalytics.eu/Patrick.Wessa/rwasp_demand_forecasting_analysis.wasp',sep=''),'Forecast Analysis','',target='')) a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable.tab') -SERVER-wessa.org
|