I am trying to run poisson regression with multiple DVs and a single IV:
DVs<-names(counts)[1:n] summary(fits <- sapply(DVs, function(x) {glm(substitute