如何将风险比应用于R中的flexsurv模型？

我有一个flexsurv模型，我想知道是否有可能绘制一条危险比为x（其中x>=0）的曲线？例如，我想测试当我将危险比设置为不同水平时，曲线会是什么样子。我想用所有的一般分布和样条线模型来做这件事

例如，与下面的曲线相比，如何绘制HR=2的曲线

library(flexsurv)
library(survival)

data <- survival::lung
formula <- Surv(time, status) ~1
fsmodel<-flexsurvreg(formula,dist = "weibull", data = data)

plot(NULL,type='n', main = "", xlim = c(0,450), ylim = c(0,1),ylab = "Survival probability (%)",xlab = "Time (months)")
lines(fsmodel)

库（flexsurv）
图书馆（生存）
数据我不确定你想展示什么。你提到了危险比率。您需要（至少）两种条件来比较数据（和危险）。换言之，危险比将特定条件下的数据与参考数据进行比较。但是在您的模型中，
Surv（时间、状态）~1
您没有考虑不同的条件。您是否将危险比与危险函数混淆？谢谢您的回复！我只是想展示一个假设的场景，展示在不同的相对有效性水平下曲线的样子。我试图展示我正在努力实现的目标，我在我的问题中添加了一些代码。你能检查一下吗？关于这个问题有什么更新吗？我不知道你想展示什么。你提到了危险比率。您需要（至少）两种条件来比较数据（和危险）。换言之，危险比将特定条件下的数据与参考数据进行比较。但是在您的模型中，
Surv（时间、状态）~1
您没有考虑不同的条件。您是否将危险比与危险函数混淆？谢谢您的回复！我只是想展示一个假设的场景，展示在不同的相对有效性水平下曲线的样子。我试图展示我正在努力实现的目标，我在我的问题中添加了一些代码。你能检查一下吗？关于这个问题有什么更新吗？
summary<- data.frame(summary(fsmodel, t=1:450, type = "survival"))

HR <- 2.0

summary$est_HR <- summary$est^HR
lines(x=summary$time, y=summary$est_HR)