R中的仿真

Question

我需要模擬100個n = 100的數據集。 我正在嘗試類似的東西：

x1 = NULL 
for (i in 1:100){
set.seed(i)
x1[i] <- sample(1:4, 100, replace=TRUE)
}

但是我不斷收到錯誤消息：x1 [i] <-sample（1：4，100，replace = TRUE）：要替換的項目數不是替換長度的倍數

我感到難過。 任何幫助，將不勝感激。 謝謝！

Answer 1

您也可以像這樣使用“復制”功能：

x1 <- replicate(100, sample(1:4,100,replace=T))
split(x1, col(x1))

如果您有興趣，可以為R的初學者閱讀一本優秀的電子書，稱為“ R基礎和圖形”。

Answer 2

您需要做的就是將x1[i]更改為x1[[i]]

保羅