ガンマ分布のハザード関数 先日調べ物をしていたときに、ガンマ分布のハザード関数が$t \to \infty$で定数$\lambda$ (rate parameter) に収束するという話を知りました。ちょっと調べものをしていたところ、ガンマ分布のハザード関数h(t)は\lim_{t→∞} h(t) =…

はじめに ちょっと色々とやっていて R の survival パッケージの tt 関数の挙動が気になるという話になりました。 survival パッケージでは時間変動共変量の解析が実施可能であり、時間変動共変量と時間変動効果（回帰係数）についての vignette には詳細な…

楽しく読ませていただきました。【新刊案内】『宇宙怪人しまりす統計よりも重要なことを学ぶ』2024/2/26発売あの宇宙怪人が帰ってきた！医療統計を学んできたしまりすに，統計よりも重要なことはあるのか。https://t.co/fnqbMEdpRn— 朝倉書店 (@AsakuraPub) …

割合の差のスコア検定 以下のツイートをしましたところ、情報提供をいただきましたので、まとめた文章を作成しました。リスク差の検定の方の式を見るとWald型の統計量で、カイ二乗検定（これはスコア型の検定に一致する）と違うものになっているので、同じに…

役に立ちそうだったので、メモを作成してみました。 library(forestplot) library(dplyr) dat1 <- read.table(textConnection(" l0 l1 l2 l3 l4 l5 hr lcl ucl Variable 1 100 50 50 0.83 (0.70, 0.99) 0.041 0.833 0.698 0.992 Variable 2 100 50 50 2.34 (…