如何为每个行业年度运行线性回归模型，不包括R中的公司i观察值？

Question

这是我的数据集在R ......中的dput输出...

data1 Firm name= c（“ A”，“ A”，“ B”，“ B”，“ C”，“ C”，“ D”，“ D”，“ D”，“ E”，“ E”，“ F”，“ F”，“ G”，“ G”，“ H”，“ H”），行业= c（“ AUTO”，“ AUTO”，“ AUTO”，“ AUTO”，“ AUTO”，“ AUTO”，“ AUTO”，“ AUTO”，“ AUTO”，“ Pharma”，“ Pharma”，“ Pharma”，“ Pharma”，“ Pharma”，“ Pharma”，“ Pharma”，“ Pharma”），X = c（1、2、5、6、7、9、10、11、12，13，15，16，17，17，18，19，20，21），Y = c（30，31，34，35，36，38、39、40、41、42、44、45、46、47、48、49、50），Z = c（23，29、47、53、59、71、77、83、89、95、107、113、119、125、131，137，143）），row.names = c（NA，-17L），class = c（“ tbl_df”，“ tbl”，“ data.frame”），na.action =结构（c（1 = 1L），class =“省略”））37、38、39、40、41、42、43、44、45、46、47、48、49、50），Z = c（23，29、35、41、47、53、59、65、71、77、83、89、95、101、107、113，119，125，131，137，143）），row.names = c（NA，-21L），class = c（“ tbl_df”，“ tbl”，“ data.frame”），na.action = structure（c（1 = 1L），class =“ omit”））

在这里，我尝试对每个行业年度进行回归[[Y〜X + Z，但不包括公司i的观察值。对于每个公司，我想估计线性回归模型使用所有行业同行公司的观察值，但不包括公司的观察值。例如，对于公司A，我想通过使用跨行业的所有同行公司（B，C和D）的所有观察结果来回归Y〜X + Z，但不包括公司A的观察结果。同样，我想通过使用企业A，C和D（企业B与B属于同一行业的一部分）的所有观察值（企业B观察值除外）运行企业B的模型。同样，对于公司的C＆D也是如此。我想对每个行业的每个公司进行此练习。请帮助。

这是我的数据集在R ...... data1中的dput输出，>

Answer 1

显示的代码不好看。但是从您写的内容来看，我建议您使用嵌套循环，例如：

for(y in year){
    for(comp in FirmName){
      # transform data : select only companys in this industry, but exclude comp
       lm(..)
     }
 }