I have fitted nested linear mixed-effects models (two levels) using the nlme package in R and I am now interested in bootstrap procedures to obtain robust estimators and confidence intervals.

A further alternative to obtain robust estimators mentioned in the literature referred to robust standard errors (e.g. sandwich estimator). I have found the "clubSandwich" package and was not sure whether it is appropriate for my case.

Any advice is greatly appreciated!