r-我应该如何将模拟研究中的数据存储在矩阵中

Question

Im目前正在进行蒙特卡罗模拟，其中Im测试哪种模型最适合变量选择。在这样做的时候，我需要将每个模拟的信息存储在一个矩阵中，以便以后进行分析。作为R的新手，我在执行这项任务时遇到了困难。这是相关代码：

iter <- 2000
results.matrix <- matrix(ncol = 12, nrow = iter )
for (i in iter) {
Ylink <- (1 * (age - mean(age)) + (1.5 * (gender == "1")) + 1.2 * (inkomst - mean(inkomst)))
Yprob <- 1/(1 + exp(-Ylink))
Y <- rbinom(2000, size = 1, prob = Yprob)
df <- data.frame(x1, x2, x3, x4, x5, x6, x7, age, gender, inkomst, Y)
glm0 <- glm(Y ~ x1 + x2 + x3 + x4 + x5 + x6 + x7 + age + gender + inkomst, data = df)
summary(glm0)
results.matrix[,1] <- summary(glm0)$coef[, "Pr(>|t|)"]
}

当尝试时，出现以下错误：；要替换的项目的数量不是替换长度的倍数"；

Answer 1

试试这个：

iter <- 2000
results.matrix <- matrix(ncol = 12, nrow = iter )
for (i in seq_len(iter)) {
Ylink <- (1 * (age - mean(age)) + (1.5 * (gender == "1")) + 1.2 * (inkomst - mean(inkomst)))
Yprob <- 1/(1 + exp(-Ylink))
Y <- rbinom(2000, size = 1, prob = Yprob)
df <- data.frame(x1, x2, x3, x4, x5, x6, x7, age, gender, inkomst, Y)
glm0 <- glm(Y ~ x1 + x2 + x3 + x4 + x5 + x6 + x7 + age + gender + inkomst, data = df)
results.matrix[i,] <- summary(glm0)$coef[, "Pr(>|t|)"]
}