乱数ベクトルの任意の数の列を作成するにはどうすればよいですか？

Question

私は、多くの場合、いくつかの統計モデルをテストするために

sample_size = 10^4 
my_data <- data.frame(x1 = runif(sample_size, 0,3), x2 = runif(sample_size, 0,3), x3 = runif(sample_size, 0,3), x4 = runif(sample_size, 0,3)) 

のようなものを記述する必要があります。例えば、my_dataので

error <- rnorm(sample_size, 0, 0.1) 
y <- with(my_data, 2*x1+0.1*(x2 + x3 + x4)) + error 
my_model <- lm(y ~ ., data = my_data) 

がlmへの入力として使用され、それはデータ・フレーム（またはリスト）でなければなりません。

これを行うには、runifを4回呼び出すのが正しい方法か、より良い解決策があるのだろうかと思います。試しました

my_data <- matrix(4*runif(sample_size, 0,3), sample_size, 4, dimnames = list(NULL, paste0("x", 1:4))) 
my_data <- as.data.frame(my_data) 

しかし、私にはわかりにくいようです。

Answer 1

これを行うにはいくつかの方法があります。あなたがncol列をしたいとしましょう、ここではいくつかの良い方法があります。

ncol = 4 
sample_size = 10 

replicate(ncol, runif(sample_size, 0, 3)) 
matrix(runif(sample_size * ncol, 0, 3), ncol = ncol) 
sapply(1:ncol, function(x) runif(sample_size, 0, 3)) 

は、これらは、必要に応じて、あなたが、当然のことながら、データフレームに変換することができます行列を作成します。違いは軽微です。 replicateは本質的にsapplyの素晴らしいラッパーです。直接matrixメソッドは少し遅くであるかもしれませんが、おそらく差は数ミリ秒です。