请求如下，希望能够得到帮助
数据如下(附件555dta文件)，共有500个公司，每个公司有若干条日度数据，每条数据后面都有一个残差（最后一列数据e），现需要每一个公司每90天的数据进行滚动回归得到残差的标准差。具体就是一个公司1-90天的90个残差算出标准差作为第91天的某一个指标，2-91天的90个残差算出的标准差作为92天的某一个指标，依次类推，共需算出500个公司的数据。请问该如何用Stata写代码呢？或者其他软件可以解决的？非常期待指导！


code Date time x1 x2 x3 x4 x5 y e
1 20aug2010 1 -.018048 -.00151 .002432 -.00529 -.002467 14.239932 19.80529
1 23aug2010 2 -.000048 .008739 -.006269 -.006542 .003935 14.009932 21.09659
1 24aug2010 3 .005953 .006162 -.002887 .000329 .002978 14.039932 22.14069
1 25aug2010 4 -.021048 .001403 -.007125 -.000121 .003066 14.089932 20.87462
1 26aug2010 5 .003953 .002926 -.004748 -.002879 .00322 14.179932 21.0652
1 27aug2010 6 .003953 .001991 .000478 .00501 -.001961 14.489932 22.17968
1 30aug2010 7 .017953 .005164 -.009715 5.000e-06 .005471 14.669931 21.88166
1 31aug2010 8 -.002048 .010774 -.009231 -.000837 .003385 14.699931 22.06408
1 01sep2010 9 -.008048 -.005091 .000284 .000896 -.00277 14.359931 20.35955
1 02sep2010 10 .013953 .005489 -.002428 -.004196 .002821 14.329931 21.40074
1 03sep2010 11 .000953 .006951 -.009174 -.005641 .003059 14.519931 20.77761
1 06sep2010 12 .010953 -.011732 .007343 .001093 .001455 14.219931 20.58143
1 07sep2010 13 .001953 .004548 -.002967 -.005159 .003449 14.399931 20.99195
1 08sep2010 14 .001953 .005514 -.01037 -.001304 .005617 14.609931 21.48178
1 09sep2010 15 -.014048 .000862 -.002873 -.000242 .000172 14.76993 20.89603
1 10sep2010 16 .003953 .006121 -.00624 -.000576 .0005 14.89993 21.538
1 13sep2010 17 .010953 .005882 -.004681 .000367 .001485 15.25993 21.93331
1 14sep2010 18 -.000048 .001958 -.00504 -.004867 .00238 14.87993 20.46824
1 15sep2010 19 -.016048 -.00152 .000663 -.000242 .0008 14.76993 20.84622
1 16sep2010 20 -.019048 -.004614 .004856 .00226 -.001649 14.87993 20.90155
1 17sep2010 21 -.001048 -.003343 -.006641 .001128 -.001371 14.76993 20.38454