The following R code may be used for constructing two-sided likelihood based intervals for the predicted probabilities of a logistic regression model. These likelihood based intervals are also referred to as likelihood ratio bounds, or profile likelihood intervals.