我想计算1000多支股票某一年度的非系统风险，打算用该年份内股票i的日回报率与日市场回报率回归，所得残差的标准差即为非系统风险。如下式：

ri=αi+β*rm+εi

式中，rm为市场收益率，ri为股票i的收益率。

请问一下，如何批量计算呢？总不能一家公司一家公司去回归吧。请尽量详细一些，非常感谢。