R中非线性混合效应模型的White-Huber标准误差（鲁棒）

Question

在R中，lmtest和Sandwich软件包为回归对象提供了用户友好的可靠标准错误。我有一个非线性混合效应模型，我不确定该怎么做？

get_data <- function(obj, start, end) {
   get(obj)[paste(start, end, sep = "/")]
}

usagextsobjects <- with(sum_datesdf,
                            c(Map(get_data, var1, startperioddate, endperioddate),
                              Map(get_data, var2, startperioddate, endperioddate)
                             )
                       )