Functions to calculate predicted values and the difference between the two cases with confidence interval for lm() [linear model], glm() [general linear model], glm.nb() [negative binomial model], polr() [ordinal logistic model] and multinom() [multinomial model] using Monte Carlo simulations. Reference: Bennet A. Zelner (2009) <doi:10.1002/smj.783>.
