Я пытаюсь построить график взаимодействия с функцией interact_plot
из пакета jtools
.Моя модель - coxph
из пакета survival
.
Я получаю следующую ошибку:
Ошибка в match.arg (тип): «arg» должно бытьодин из «lp», «риск», «ожидаемый», «условия», «выживание»
Я пытался добавить type="survival"
, но продолжаю получать ту же ошибку.Любой совет, в чем может быть проблема?
Вот мой код:
library(survival)
hiitime.coxph<-coxph(formula=Surv(rep(1, 77748L), observed) ~
log(hii+1)*weekssincedisp+
strata(id) + cluster(wolfid),data = timemodtata, method = "breslow")
library(jtools)
interact_plot(hiitime.coxph, pred = "hii", modx = "weekssincedisp",
modx.values = c(0, 13, 26, 39, 52),
interval = TRUE, int.width = 0.95,
data = timemodtata,
plot.points = TRUE,
x.label = "Human impact index",
modx.labels = c("0 weeks", "13 weeks", "26 weeks", "39 weeks", "52 weeks"),
legend.main = "Weeks since dispersal"
)