1

我想使用 runif 在 R 中创建一系列随机二维点。有什么好的方法可以做到这一点?

这是我当前的实现(不起作用)

getpoints <- function(n)
  {
    t <- 1:n;
    for (i in 1:n)
      {
        t[i] <- runif(2,-1,1)
      }
    return (t)
  }
4

4 回答 4

9
n <- 200
o <- matrix(runif(2*n, min = -1, max = 1), ncol = 2)
于 2013-01-13T10:57:45.473 回答
3
 n <- 5
 cbind(runif(n,-1,1),runif(n,-1,1))
            [,1]       [,2]
[1,] -0.68434317 -0.7772889
[2,] -0.91200792 -0.6408075
[3,] -0.01888610  0.7350491
[4,]  0.01782097  0.9674700
[5,] -0.56707264 -0.9991566
于 2013-01-13T10:57:37.910 回答
2

只需添加您在这里的众多答案。这个,像@Arun's 很容易在必要时处理超过 2 列。:

replicate(2, runif(100))

它可能几乎是基准测试时间。

于 2013-01-13T11:27:19.403 回答
-1

为了概述这个线程,首先我将展示与所使用的几种方法相关的速度(这样您就可以得出自己的结论来判断哪种方法最好),然后我将向您展示如何使用函数;即您可能希望您的功能的方式。

n <- 1000000

> system.time(cbind(runif(n, -1, 1), runif(n, -1, 1)))
   user  system elapsed 
   0.28    0.03    0.31 
> system.time(matrix(runif(2*n, min = -1, max = 1), ncol = 2))
   user  system elapsed 
   0.32    0.00    0.33 
> system.time(replicate(2, runif(n, -1, 1)))
   user  system elapsed 
   0.26    0.02    0.35 

这是函数和时间基准:

getPoints <- function (N) {
  pointMatrix <- matrix(1, nrow = 1000000, ncol = 2) 
  for (i in 1:N) {
    pointMatrix[i, ] <- runif(2, -1, 1) # it needs to run for all columns
  }
}

> system.time(getPoints(n))
   user  system elapsed 
  17.59    0.04   19.94 
> 
于 2013-01-13T13:44:47.547 回答