使用使用group_by函数的回归模型的系数

Question

我运行以下代码来生成回归模型

library (dplyr)
fitted_models <- df %>% 
  group_by(sic, fyear) %>% 
  do (model = lm (TACCdTA ~  Inverse_TA + DeL_RevRec + PPEdTA , data = .)) 

然后要获取每个sic和fyear的系数，我运行以下代码

library(broom)
fitted_models %>% tidy(model)

我获得了每个 fyear 的每个 sic 的系数。

现在我的问题是-在每个 fyear 的每个 sic 下，有很多观测值（例如1000个观测值）-如何计算每个观测值的拟合值使用上面的模型生成的系数，在每个原文和 fyear 下。

另一个小问题-对于我运行模型的第一个代码，如何确保仅针对每个 sic 和< strong> fyear 组合（sic-fyear）至少有10个观测值。