R语言多元正态分布采样中的浮点陷阱
2025-05-22
本文讲述了作者帮助同事解决R语言代码中多元正态分布采样结果不可重复的问题。问题根源并非R语言或MASS包本身的错误，而是浮点运算的固有特性导致的。即使使用`set.seed()`设定随机数种子，不同机器上运行相同的代码，由于浮点计算的舍入误差，会导致`MASS::mvrnorm()`函数产生不同的结果。作者深入分析了`MASS::mvrnorm()`和`mvtnorm::rmvnorm()`函数的底层实现，指出`MASS::mvrnorm()`使用特征值分解，对微小的输入扰动非常敏感，可能导致特征向量符号反转，从而影响结果的可重复性。而`mvtnorm::rmvnorm()`使用Cholesky分解，更加鲁棒，解决了这个问题。作者建议使用`mvtnorm::rmvnorm()`并设置`method = "chol"`来提高代码的可重复性。
阅读更多
开发