我想在 R 中运行固定效应回归，为此我定义了以下公式：
time.aspects <- as.formula(y ~ x1 + x2 + x3 + t)
time.total <- plm(time.aspects, data=all, index=c("i","t"), model = "内”）
x1、x2 和 x3 是我的自变量。我还想添加一个时间因子 t 来说明时间固定效应。在这方面，t 代表单年 1 到 10（包含在我的数据文件中）。
但是，如果我想通过以下方式考虑稳健的标准错误：
coeftest( time.total , vcov.= vcovSCC( time.total , type = "HC3"))
发生以下错误： 1 - diaghat 中的错误：二元运算符的非数字参数。
有谁知道如何避免这个错误信息？