library(rms)
patient_data <- data.frame(
PatientID = 1:1000,
Death = sample(c(0, 1), size = 1000, replace = TRUE),
Age = round(runif(1000, min = 20, max = 90), 1)
)
lrm(Death~Age, data = patient_data) #works fine
glm(Death ~ Age), family = "binomial", data = patient_data) #works fine
glm(Death ~ offset(2*Age), family = "binomial", data = patient_data) #works fine
lrm(Death ~ offset(2*Age),data = patient_data) #ERROR
rms
作者似乎仍然是真的:
offset
。我们可以通过创建一个玩具协变量并将其添加到您的
lrm
公式:来检查一下:
patient_data$x = runif(1000) lrm(Death ~ x + offset(2*Age), data = patient_data)
,现在似乎有效。