我有50个权重指标，从x1到x50。想用循环命令，获得不同权重处理下的教育回报率系数。但下面的结果，是得到最后一个回归的单个教育回报率系数，不是所有教育回报率系数的矩阵。各位能帮忙改进程序吗？

foreach i of numlist 1(1)50{
reg wage eduy gdp experience water employ consump [pw=x`i']
matrix cc=e(b)
matrix cf=cc[1,1..4]
}
matrix edu=(cf[1..50,1])