次の式を定義する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 : 二項演算子の引数が数値ではありません。
このエラーメッセージを回避する方法を知っている人はいますか?