z(zは0から1の実数)を超えるまでに発生したすべての乱数を並べたベクトルを書きたい

Question

rand_q<-function(z){
u<-runif(1)
s<-0
if ( u<= z){
repeat {
t=runif(1)
if (t <= z) s<-c(u,t) ,

else break 

}}

return(s)
}

Answer 1

rand_q<-function(z){
  s<- runif(1)
repeat{
  if (max(s) <= z)
    {s<- c(s,runif(1))}
  else        break
}
  return(s)
}

Comments

[Mizugorou]

ありがとうございます。この方法は思いつきませんでした。助かります。