Adjust model for unequal variances
.updateUnequalVariances(input, fit, num_iter)
input
data.table
fit
lm
num_iter
number of iterations
merMod