2
正しい打ち切りでおもちゃの生存(イベントまでの時間)データを作成する方法
正しく打ち切られ、比例ハザードと一定のベースラインハザードのある分布に従うおもちゃの生存(イベントまでの時間)データを作成したいと思います。 次のようにデータを作成しましたが、Cox比例ハザードモデルをシミュレーションデータにフィッティングした後、真の値に近い推定ハザード比を取得できません。 私は何を間違えましたか? Rコード: library(survival) #set parameters set.seed(1234) n = 40000 #sample size #functional relationship lambda=0.000020 #constant baseline hazard 2 per 100000 per 1 unit time b_haz <-function(t) #baseline hazard { lambda #constant hazard wrt time } x = cbind(hba1c=rnorm(n,2,.5)-2,age=rnorm(n,40,5)-40,duration=rnorm(n,10,2)-10) B = c(1.1,1.2,1.3) # hazard ratios (model coefficients) hist(x %*% …