我不知道“某些条件”是什么意思，但如果你想通过模拟计算回归系数的不确定性，你必须考虑一般残差标准偏差的不确定性（因为 Var$(\,\hat{\beta}\,|\,X\,) = \sigma^2\,(X^TX)^{-1}$, $\sigma^2$应该$\chi^2_{n-k}$分布，但在您的情况下可能没有不确定性）以及回归系数（通常假定为多元正态分布）。

有一个关于模拟的好章节

Gelman, A. 和 Hill, J. (2006)。使用回归和多级/分层模型的数据分析。剑桥大学出版社

它还提供了一个 R 函数“sim”（在arm包中）来做你想做的事情。您可以查看他们如何在 R 中编写此函数以了解您必须做什么，或者您可以阅读他们在书中关于模拟的章节。