r - 円と正方形に乱数を描く方法は?

Question

Polar メソッドのアルゴリズムを 2 つの視点 (鳥瞰図とカエルの視点) から視覚化したいと考えています。

メソッドの最初のステップを描画するには、正方形と円の一様分布の乱数が必要です。

すでに円を正方形にプロットできましたがz1、この構成内に乱数 ( ) をプロットするにはどうすればよいですか?

require (plotrix)
require (grid)
z1 = runif (100)

plot (c(-1,1), c(-1,1), type="n", asp=1)
rect(-1,-1,1,1)
draw.circle (0,0,1)

そして、視点を変える方法は？

Answer 1

ポイントを「拒否」せずに行うことができます。次の R 関数は3*n乱数を必要とnし、radius の円内にランダムに選択された点を生成しrます。

randp <- function(n = 1, r = 1) {
    if (n < 1 || r < 0) return(c())
    x <- rnorm(n)
    y <- rnorm(n)
    r <- r * sqrt(runif(n)/(x^2 + y^2))
    if (n == 1) U <- c(x, y)
    else        U <- cbind(r*x, r*y)
    return(U)
}