I found an answer for my question that others can learn from and use as well:
library(effects, quietly = TRUE) mar.ef <- Effect("x1", glm.model) mar.ef
This code gave me the same results as Stata!