我知道R中的软件包（例如sandwich）采用lm对象并在摘要中创建可靠的标准错误。我需要操纵lm对象，使其具有可靠的标准错误（而不仅仅是更改lm的摘要）。最终，我想使用效果包从我的（lm）模型中绘制一个交互，但是效果的输入必须是lm。
我尝试过的所有操作都只会更改summary.lm()输出，而不会更改lm模型本身。
如何调用lm对象，计算健壮的标准误差，并将变量保持为lm对象？
基本上，在运行具有可靠标准错误的OLS模型后，我试图在Stata中复制marginsplot函数。