1poisson.sim<-function(data, mle=NULL) 2{ 3 if(is.null(mle)) 4 data$Observed<-rpois(length(data$Observed), lambda=data$Expected) 5 else 6 data$Observed<-rpois(n=mle$n, lambda=mle$lambda) 7 return(data) 8} 9