我正在尝试创建1000个数据集,每个数据集具有4个变量,X1是正常(0,1),X2是正常(0,1),X3是expit(0.5X1 + 0.8X2),A是bernouli(X3 )。所有变量都有2500个观测值。
我想将1000个数据集中的每一个都放入glm模型,但无济于事。
那是我所做的
sim_list = replicate(n =3,expr = {data.frame(X1 = X1 <- rnorm(2500,mean = 0,sd = 1),X2 = X2 <- rnorm(2500,X3 = X3 <- rnorm(2500,X4 = X4 <- expit(0.5*X1 + 0.8*X3),A = A <- rbern(2500,X4)},simplify = F)
for (i in 1:3){
glm1[i] <- glm(A~X1,family = binomial,data = sim_list[[i]])
score[i] <- glm1[i]$fitted.values
}
在此方面,我真的很感谢帮助