円内に一様分布する乱数を生成する時は、俺、絶対忘れないよヤコビアンのこと - My Life as a Mock Quant

円内に一様分布する乱数を生成する時は、俺、絶対忘れないよヤコビアンのこと - My Life as a Mock Quant

頭出し ある半径Rの円内に一様分布する乱数を生成する時には注意しないといけないことがありますよとい...概要を表示 頭出し ある半径Rの円内に一様分布する乱数を生成する時には注意しないといけないことがありますよというお話。所謂「一度はやってしまうミス」系でもある。この手の話は円に限ったわけではなく、円の高次元版である球、あるいは超球(次元>3)、あるいは任意の座標変換をかませてそこにヤコビアンが出て来るときでも同じ。 使うライブラリは以下の２つ。なければinstall.packages関数でインストールしておく。 library(ggplot2) library(dplyr) また、以下のように定数を２つ定義しておく。意味はコメントにある通りだ。 #サンプル数 N <- 10^4 #半径のサイズ R <- 4 本題 さて、問題のある半径Rの円内に一様分布する乱数を生成するにはどうしたらいいのかというと、非常に単純に考えた場合、以下のように(俺は)思考する。 X方向の成分として[-R, R]の間の一様乱数

• Random
• algorithm
• R
• 数学