require(lsmeans)
# Example from the function glm()
counts <- c(18,17,15,20,10,20,25,13,12)
outcome <- gl(3,1,9)
treatment <- gl(3,3)
glm.D93 <- glm(counts~outcome+treatment,family="poisson")
lsm <- lsmeans(glm.D93,~outcome)
inv.lsmeans(glm.D93,lsm)
Run the code above in your browser using DataLab