我是统计学的新手,因此我很难绘制R的cox回归的交互作用。
我试图找到一个合适的程序包,并发现simPH可能有用,但是,由于我对说明/手册的理解不够,所以无法正确实现它。
以下代码显示了我的Cox回归(cox_prp)。在此示例中,我尝试检查通货膨胀(Inf_E)与另一个变量(Kammer)的相互作用是否导致政府更早崩溃。 我猜想函数“ coxsimInteract”(来自simPH包)可以工作,但老实说我不理解手册中有关单个参数的功能。
cox_prop <- coxph(Surv(Duration, termination) ~ Inf_E*Kammer, data = daten)
我也尝试了在互联网上找到的代码,但这似乎无法正常工作,实际上我在不知道参数含义的情况下将其复制了:
Sim2 <- coxsimInteract(cox_prop, b1 = "Inf_E", b2 = "Kammer", qi = "Marginal Effect", X2 = seq(2, 115, by = 2), nsim = 1000)
simGG(Sim2, ribbons = TRUE, alpha = 0.5, xlab = "Kammer\n", ylab = "Marginal effect of Inf_E\n")
哪个函数可以绘制交互关系,以及如何设置参数?