ধরা যাক আমার একটি "কিডনি ক্যাথেটার" ডেটা সেট রয়েছে। আমি একটি কক্স মডেল ব্যবহার করে বেঁচে থাকার কার্ভকে মডেল করার চেষ্টা করছি। যদি আমি একটি কক্স মডেল বিবেচনা করি:
survival
প্যাকেজ আর ফাংশনটি ব্যবহার করে basehaz()
, আমি এটি সহজেই এটি করতে পারি:
library(survival)
data(kidney)
fit <- coxph(Surv(time, status) ~ age , kidney)
basehaz(fit)
তবে যদি আমি প্যারামিটারের একটি নির্ধারিত অনুমানের জন্য বেসলাইন বিপদের এক ধাপে ধাপে লিখতে চাই তবে আমি b
কীভাবে এগিয়ে যেতে পারি? আমি চেষ্টা করেছিলাম:
bhaz <- function(beta, time, status, x) {
data <- data.frame(time,status,x)
data <- data[order(data$time), ]
dt <- data$time
k <- length(dt)
risk <- exp(data.matrix(data[,-c(1:2)]) %*% beta)
h <- rep(0,k)
for(i in 1:k) {
h[i] <- data$status[data$time==dt[i]] / sum(risk[data$time>=dt[i]])
}
return(data.frame(h, dt))
}
h0 <- bhaz(fit$coef, kidney$time, kidney$status, kidney$age)
কিন্তু এটি হিসাবে একই ফলাফল দেয় না basehaz(fit)
। সমস্যাটা কি?