ryamada22 https://blog.hatena.ne.jp/ryamada22/ ryamadaの遺伝学・遺伝統計学メモ https://ryamada22.hatenablog.jp/ 生存解析 R n <- 1000 # ランダムに脱落する prob <- c(0.5,0.5) censor <- sample(0:1,n,replace=TRUE,prob=prob) # 指数乱数で死亡時刻が決まる #times <- cumsum(rexp(n)) times <- sort(rexp(n)) x <- ((n-1):0)/(n:1) v <- rep(1,length(x)) v[which(censor == 1)] <- x[which(censor ==1)] #par(mfcol = c(1,2)) #plot(times,cumprod(v),type="s",ylim … 190 <iframe src="https://hatenablog-parts.com/embed?url=https%3A%2F%2Fryamada22.hatenablog.jp%2Fentry%2F20140304%2F1393903280" title="イベント集積に関するメモ - ryamadaの遺伝学・遺伝統計学メモ" class="embed-card embed-blogcard" scrolling="no" frameborder="0" style="display: block; width: 100%; height: 190px; max-width: 500px; margin: 10px 0px;"></iframe> https://cdn-ak.f.st-hatena.com/images/fotolife/r/ryamada22/20140304/20140304122243.jpg Hatena Blog https://hatena.blog 2014-03-04 12:21:20 rich https://ryamada22.hatenablog.jp/entry/20140304/1393903280 1.0 100%