Replace a full with a block-diagonal variance covariance matrix Check if a full random effects variance covariance matrix is specified for a single variable. In that case, it is identical to a block-diagonal matrix. Change the rd_vcov specification to blockdiag for clarity (because then the variable name is used in the name of b, D, invD, ...)