R:glmmLasso中的REML估计到底是怎么回事。。。

R的glmmLasso中的control函数里面的method如是提到：
two methods for the computation of the random-effects variance-covariance pa-
rameter estimates can be chosen, an EM-type estimate and an REML-typees-
timate. The REML-type estimate uses  bobyqa function for optimization.
Defaultis EM.
其中，EM我倒是能明白，但是广义线性混合模型，怎么做REML-type估计？
见附件第5页

文章里面好像说的特别自然。。。。 但是不是只有普通的高斯分布才有REML一说吗？其它分布怎么做？

想要收回论坛币。。。