R中的行采样

Question

我使用示例数据来提出问题。

seed(1) 
X <- data.frame(matrix(rnorm(200), nrow=20))

我想每次都选择 10 个随机行而不进行替换并进行多元回归。我试过了

hi=X[sample(1:20,10),]
MR1<-lm(X10~., data=hi)
R1<-summary(MR1)$r.squared #extract the R squared

是否可以创建 25 个这样的数据集，每次采样 10 行。最后，我想存储采样数据集并进行多元回归，并从 25 个这样的模型中提取 r 平方值。

Answer 1

你可以使用lapply:

set.seed(1) 
X <- data.frame(matrix(rnorm(200), nrow=20))

n <- 25
res <- lapply(1:n, 
              function(i) {
                samples <- sample(1:20,10)
                hi=X[samples,]
                MR1<-lm(X10~., data=X)
                R1<-summary(MR1)$r.squared
                return(list(Samples=samples,Hi=hi,MR1=MR1,R1=R1))
              })