模拟以适合不同比例的物流

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
}

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