Estou usando a biblioteca plm para executar regressões de efeito fixo e as bibliotecas lmtest sanduíche para calcular erros padrão robustos. Não tenho nenhum problema ao executar as regressões, mas em alguns casos quando vou calcular os erros padrão, recebo o seguinte erro:

library(plm)
library(sandwich)
library(lmtest)

fe_reg <- plm(y ~ x + I(x^2)+factor(date), data=reg_data, index=c("id","date"), model="within")
coeftest(fe_reg, vcov.=vcovHC(fe_reg, type="HC1"))

RRuntimeError: Error in solve.default(crossprod(demX))
  system is computationally singular: reciprocal condition number = 1.84726e-25

Não tenho problemas em calcular os coeficientes ou os erros padrão "normais" (ou seja, homoscedásticos). Além disso, não tenho problemas em calcular os erros padrão robustos quando omito o termo quadrático:

fe_reg <- plm(y ~ x +factor(date), data=reg_data, index=c("id","date"), model="within")  

Alguém tem alguma idéia do que está acontecendo? Se a matriz de design fosse singular, os coeficientes não deveriam ter sido calculados, portanto não entendo de onde vem o problema ao calcular os erros padrão.

Obrigado!