我想申请的甩采样方法来模拟随机矢量Y=(Y_1, Y_2)从一个单元盘状的均匀分布的D = { (X_1 , X_2) \in R^2: \sqrt{x^2_1 + x^2_2} ? 1}，使得X = (X_1 , X_ 2)处于正方形的均匀分布的随机矢量S = [?1, 1]^2和接头密度f(y_1,y_2) = \frac{1}{\pi} 1_{D(y_1,y_2)}.
在拒绝方法中，我们一般接受一个样本如果f(x) \leq C * g(x)。我正在使用以下代码：
x=runif(100,-1,1)
y=runif(100,-1,1)
d=data.frame(x=x,y=y)
disc_sample=d[(d$x^2+d$y^2)<1,]
plot(disc_sample)
我有两个问题：
{使用上面的代码，从逻辑上讲，的大小d应该大于的大小，disc_sample但是当我调用它们时，我看到它们中的每一个都有 100 个元素。这怎么可能。为什么尺寸相同。}这部分已解决，感谢下面的评论。
现在的问题
另外，我怎样才能重新编写我的代码，以便为我提供按照条件获取 100 个样本所需的样本总数。即给我拒绝的样品数量，直到我得到 100 个所需的样品？
感谢r2evans的回答，但我希望写一些更简单的东西，一个 while 循环将所有可能的样本存储在矩阵或数据帧而不是列表中，然后从该数据帧调用，只是样本遵循条件。我在没有使用列表和 sapply 函数的情况下修改了答案中的代码，但它没有给出所需的结果，它只产生一行。
i=0
samps <- data.frame()
goods <- data.frame()
nr <- 0L
sampsize …