生存曲線のシミュレーション 2 について、 ri<-function(d,t){ if (length(d[,1][d[,1]>=t])==0)return(1) else return(length(d[,1][d[,1]>=t]) }としておかないと、length(d[,1][d[,1]>=t])==0のときにS.hatがNaNとなってしまうことに気が付いた。（「1」…

生存曲線のシミュレーション 時刻を0,1,2…とし、時間1の間の死亡率を、脱落率をとするモデルを考えると、 Surv.time<-function(p,q,l){ cnt<-1 repeat{ if (cnt==l){x<- c(cnt,0);break} else if (runif(1)< p) {x<- c(cnt,1);break} else if (runif(1)< q) …