pwcorr_a  var1 var2 var3 var4, star1(0.01) star5(0.05) star10(0.1)是不分年度求相关系数的代码，
而后面的回归结果需要分年度回归，回归的时候是将年度year前面加上“i.”，但是求pearson系数时不能直接加“i.”。
那么如何在求pearson系数的时候将年度变量添加进去呢？