Время обнаружения данных в 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.Спасибо