Ошибка вычисления ошибок Robust Standard в регрессионной модели Panel (plm, R)

Я использую библиотеку plm для запуска регрессий с фиксированными эффектами и сэндвич, библиотеки lmtest для вычисления устойчивых стандартных ошибок. У меня нет проблем с запуском регрессий, но в некоторых случаях, когда я иду, чтобы вычислить стандартные ошибки, я получаю следующую ошибку:

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

У меня нет проблем с вычислением коэффициентов или «нормальных» стандартных ошибок (т. Е. Гомоскедастических). Кроме того, у меня нет проблем с вычислением устойчивых стандартных ошибок, когда я опускаю квадратный термин:

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

У кого-нибудь есть идеи, что происходит? Если бы матрица проекта была единственной, то коэффициенты не должны были быть вычислены, поэтому я не понимаю, откуда возникает проблема при вычислении стандартных ошибок.

Спасибо!