R中的时间检测数据-

Question

#simulate data
M=250
lambda=0.3
psi=0.4
Tmax=10
y<-matrix(NA, ncol=Tmax, nrow=M)
z<-rexp(n=M, rate=lambda)
for (j in 1:Tmax){
  y[,j]<-rbinom(n=M, size=1, prob=(z*psi))
}

如何模拟在R中检测数据的时间。我正在使用它并获得大量NA。 谢谢