ガリガリ計算して確率の対数をとって、もっともらしいpとqの値をそれぞれ0.01～0.05まで0.0050刻みで選べる関数をつくった。 p<-runif(1,0.01,0.05) q<-runif(1,0.01,0.05) Surv.time<-function(p,q,l){ cnt<-1 repeat{ if (runif(1)< p) {x<- c(cnt,1);brea…

生存曲線のシミュレーション 4 昨日のコード、ところどころ修正したミスを書き換えてなかった部分があったので、ちゃんとしたものを記録として残しておく。 Surv.time<-function(p,q,l){ cnt<-1 repeat{ if (cnt==l){x<- c(cnt,0);break} else if (runif(1)<…