我想在R中运行固定效果回归，我为其定义了以下公式：
time.aspects＆lt; - as.formula（y~x1 + x2 + x3 + t）
time.total ＆lt; - plm（time.aspects，data = all，index = c（＆＃34; i＆＃34;，＆＃34; t＆＃34;），model =＆＃ 34;内＆＃34）
x1，x2和x3是我的自变量。我还想添加一个时间因子t来计算时间固定效应。在这方面，t代表单年1到10（包括在我的数据文件中）。
但是，如果我想以下列方式考虑强大的标准错误：
coeftest（ time.total ，vcov。= vcovSCC（ time.total ，类型=＆＃34; HC3＆＃34;）
发生以下错误：错误1 - diaghat：二元运算符的非数字参数。
有谁知道如何避免此错误消息？