「mvrnorm」の記事一覧

• 再現性が壊れるとき：疑似乱数と浮動小数点の罠

  2025-05-19 07:17

  科学・技術

  Rでset.seed()を用いても、マルチバリアント正規分布の乱数生成が異なるマシン間で再現されない事象についての詳細な分析記事です。原因は、浮動小数点演算における微小な差異とそれに起因する固有ベクトルの符号反転にありました。特にMASS::mvrnorm()での固有分解によって大きな違いが生じ、わずかな行列の違いが大きな結果の差異を生むことが実証されています。この記事は、再現性を求める科学者にとっての教訓を豊富に含んでいます。