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页