Estoy usando la biblioteca plm para ejecutar regresiones de efectos fijos y las bibliotecas sandwich, lmtest para calcular errores estándar robustos. No tengo problemas para ejecutar las regresiones, pero en algunos casos cuando voy a calcular los errores estándar me sale el siguiente error:

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

No tengo ningún problema para calcular los coeficientes o los errores estándar "normales" (es decir, homoscedastic). Además, no tengo problemas para calcular los errores estándar robustos cuando omito el término cuadrático:

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

Alguien tiene alguna idea de lo que está pasando? Si la matriz de diseño fuera singular, entonces los coeficientes no deberían haberse calculado, por lo que no entiendo de dónde proviene el problema al calcular los errores estándar.

¡Gracias!