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