fastmode function

Rapidly converging algorithm for calculating posterior modes in linear mixed models