diff --git a/src/GLM.jl b/src/GLM.jl
index 807019fd..91e47688 100644
--- a/src/GLM.jl
+++ b/src/GLM.jl
@@ -12,12 +12,12 @@ module GLM
     import Statistics: cor
     import StatsBase: coef, coeftable, confint, deviance, nulldeviance, dof, dof_residual,
                       loglikelihood, nullloglikelihood, nobs, stderror, vcov, residuals, predict,
-                      fit, model_response, r2, r², adjr2, adjr², PValue
+                      fit, model_response, response, modelmatrix, r2, r², adjr2, adjr², PValue
     import StatsFuns: xlogy
     import SpecialFunctions: erfc, erfcinv, digamma, trigamma
     export coef, coeftable, confint, deviance, nulldeviance, dof, dof_residual,
            loglikelihood, nullloglikelihood, nobs, stderror, vcov, residuals, predict,
-           fit, fit!, model_response, r2, r², adjr2, adjr²
+           fit, fit!, model_response, response, modelmatrix, r2, r², adjr2, adjr²
 
     export                              # types
         Bernoulli,
diff --git a/src/linpred.jl b/src/linpred.jl
index d8d0b79b..6f15e3d8 100644
--- a/src/linpred.jl
+++ b/src/linpred.jl
@@ -228,7 +228,7 @@ end
 
 modelframe(obj::LinPredModel) = obj.fr
 modelmatrix(obj::LinPredModel) = obj.pp.X
-model_response(obj::LinPredModel) = obj.rr.y
+response(obj::LinPredModel) = obj.rr.y
 
 fitted(m::LinPredModel) = m.rr.mu
 predict(mm::LinPredModel) = fitted(mm)