Latest update on 2012年3月5日 (月) at 10:54:46.
2003663613 * as.bigz(2)^195000 + 1してみたらPerlよりずっと速かった。
require(survival) start <- c(2,9,3,1,1,10,6,4,3,8,7,13) end <- c(11,15,8,3,15,13,7,8,4,13,13,15) cens <- c(1,0,0,0,0,0,0,1,0,0,1,0) hasurv <- Surv(end-start,cens) print(res <- survfit(hasurv)) plot(res)とやれば,このデータにおけるγグロブリン投与後のA型肝炎発症までの中央値が9ヶ月(95%信頼区間の下限が6ヶ月,上限は無限大)とわかるし,カプラン=マイヤプロットも描ける。
▼前【714】(また締め切り(2007年1月17日) ) ▲次【716】(ともかくも(2007年1月19日) ) ●Top
△Read/Write COMMENTS