R編程（蒙特卡羅模擬）

Question

嗨，我想問一下如何使用Monte Carlo Simulation采樣出50個虹膜數據實例（其中包含150個實例）？ 任何想法？ 非常感謝

Answer 1

我們可以使用dplyr中的 sample_n選擇50行進行替換。

# Set seed for reproducibility
set.seed(12800)

library(dplyr)
library(purrr)

iris_sub <- iris %>% sample_n(size = 50, replace = TRUE)

在這里，我將展示一種方法，可以重復這一過程1000次，使用map_dfr從purrr包。 最終結果是一個具有50000行的數據幀。 將創建一個名為“ Time ”的新列，以記錄采樣數。

iris_sample <- map_dfr(1:1000, ~iris %>%
                         sample_n(size = 50, replace = TRUE) %>%
                         mutate(Time = .x))

Answer 2

這是在基數R中執行此操作的一種方法。

您可以用替換來采樣50行

iris[sample(1:nrow(iris), size = 50, replace = TRUE), ]

要列出列表（例如，包含50行的1000個樣本進行替換），可以使用lapply 。

iris_mc_samps <- lapply(1:1000, function(x) iris[sample(1:nrow(iris), size = 50, replace = TRUE), ])