r - 模拟以适应不同比例的逻辑

Question

我是新的 R 用户，我正在尝试进行模拟研究以拟合具有不同缺失比例的逻辑模型。我尝试了一些代码，但它没有给我任何输出，有人可以帮助我吗？

p=c(.03,0.2,0.6) #proportion of missingness
for(i in 1:length(p)){
mcar<-prodNA(dd,noNA=0.2)# function to impose missingness
mis.fit=glm(y~x1+x2+x3, family = "binomial",data=mcar)
mis.fit
mis.bias=t(coef(mis.fit))-t(coef(comp.fit))
mis.bias
mis.sum=summary(mis.fit)
mis.sum
}

标签： r

Answer 1

我们无法重现您的问题：您没有提及missForest包裹，我们也不知道是什么dd。

我的猜测是你可以从print()在你的循环中使用

例如

i <- 5
i^2

给

[1] 25

但

for (i in 1:5){
 i^2
}

没有可见的输出

尽管

for (i in 1:5){
  print(i^2)
}

给

[1] 1
[1] 4
[1] 9
[1] 16
[1] 25