用R语言编一个回归分析的程序,一个y变量，3个x变量，加误差项。
That is,y=a+b1*×1+b2*x2+b3*x3+b1*x4+e,a=b1=b2=b3=1，样本量N=50.考虑下面这中情况1) x1,x2,x3 are normally distributed with means0, variances 1.0, correlations .5,e~N(O,1);重复模拟1000次,每次有一个intercept和3个regression coefficients,用一千次的值算出4个回归系数的斜方差矩阵和相应的标准误差。