我想模拟它500次，但是这个函数给出了187个正确的值，但是在187个值之后它给出了NA值

Question

set.seed(1988)

X1<-runif(188)

n<-500

for(i in 1:n){

  logarithem[i]= ifelse (X1[i] <0.5,log (C[i] ) ,(C[i]))

  }

Answer 1

该行：

X1<-runif(188)

用188值初始化X1。如果您正在索引X1 [189]，则会出错。

试试这个：

print (X1[189])
> NA

为了纠正它，只需更改你的runif：

X1<-runif(500)