我正在攻读学士学位论文，我必须计算一种称为经验瘦度指标的东西。为此，我需要获取每个观察值的学生化残差。每个观察结果都是特定年份的公司。我通过将库存水平的自然对数与收入的自然对数进行回归来进行回归分析。我还按行业（SIC代码）和货币（美元或加元）对回归进行了分组。现在，我必须找到每个观测值的学生化残差，这些残差包含在一个数据框中（样本大小为62739，有494种不同的回归），这说明了每个观测值必须与其SIC的回归相匹配的事实代码和货币。
我尝试了很多具有不同想法的代码，但没有一个奏效。
这是我目前的完整代码（invt表示总库存，revt表示总收入，conm表示公司名称）：
app-launch