#log-likelihood function for poisson distribution like<-function(lam,n,k,m,x){ -(n-k)*lam+log(lam)*sum(x[x