Мне известны пакеты в R (например, sandwich), которые принимают объект lm и создают устойчивые стандартные ошибки в сводке. Что мне нужно, так это манипулировать моим lm объектом так, чтобы он имел устойчивые стандартные ошибки (не просто изменяя сводку lm). В конечном счете, я хочу построить взаимодействие из моей (lm) модели, используя пакет эффектов, но вход для эффекта должен быть lm.
Все манипуляции, которые я пробовал, изменяют только вывод summary.lm(), а не саму модель lm.
Как я могу вызвать lm объект, вычислить устойчивые стандартные ошибки и сохранить переменную как lm объект?
По сути, я пытаюсь повторить функцию marginsplot в Stata после запуска модели OLS с устойчивыми стандартными ошибками.