# 產生假資料
set.seed(1234)
dt <-
  data.frame(
    Num = rnorm(15),
    Group = c(rep("A", 12), rep("B", 3))
  )

Num.A <- subset(dt, Group == "A")$Num
Num.B <- subset(dt, Group == "B")$Num

res <- vector("list", 4)
# res 是一個用來記錄之後做相關檢驗的容器
# 它本身是一個 list，有 4 個「位子」

for(i in 1:4){
  Num.A.simulated <- sample(Num.A, 3)
  res[[i]] <- cor.test(Num.A.simulated, Num.B)
}

res # 所有的結果
res[[1]] # 第一次模擬的結果
sapply(res, function(x){x$estimate}) # 所有的 correlation coefficient
sapply(res, function(x){x$p.value}) # 所有的 p-value