Я довольно новичок в сплайновой регрессии с R, и у меня возникли некоторые проблемы, надеюсь, вы поможете мне с этим. То, что я хочу сделать, это построить график сплайна HR события по дням. У меня есть 2 процедуры, и в идеале я хотел бы нанести оба этих HR на один график, чтобы противостоять им, хотя это не главная проблема.
Чтобы получить сплайн, я использую следующее код:
attach(database)
dd <-datadist(days, event)
options(datadist='dd')
model <- lrm(event~rcs(days, 4), x=TRUE, y=TRUE)
ggplot(Predict(model, days, ref.zero=TRUE, fun=exp), ylab="OR for event", xlab="Days")+ theme_classic()
detach(database)
Этот код работает для всей базы данных. Затем я устанавливаю набор данных для разделения двух обработок:
database_arm1 <- subset(database, Arm==1)
attach(database_arm1)
dd <-datadist(days, event)
options(datadist='dd')
model <- lrm(event~rcs(days, 4), x=TRUE, y=TRUE)
ggplot(Predict(model, days, ref.zero=TRUE, fun=exp), ylab="OR for event", xlab="Days")+ theme_classic()
detach(database_arm1)
И когда я запускаю команду lrm, я получаю следующую ошибку:
In lrm(event ~ rcs(days, 4), x = TRUE, :
Unable to fit model using “lrm.fit”
Если я пытаюсь используйте другую непрерывную переменную для сплайнов (es: Age), команда работает также для подгрупп.
У вас есть какие-либо подсказки, почему у меня возникает эта ошибка?
Второй вопрос, знаете ли вы, как Могу ли я построить 2 сплайна на одном графике?
Спасибо