#下面是用R 写的，用来做300 次模拟。如果你数实际数据，则需要将变量读进x,y中。（参加R 的函数read.table）可以参考下。
n = 100                                               ---- 样本量
sim = 300                                            ---- 重复次数
beta = matrix(0, ncol=2, nrow=sim)    ---- 矩阵的每一列分别用来装估计的系数
for(i in 1:p)
{
   x = rnorm(n)
  intercept = 0.5
   beta0 = 2
   y = intercept + x*beta0 + rnorm(n)
   beta[i,] = lm (y~x)$coefficients
}
  beta