统计绘图
之前请教大佬们,说是用R语言的RCS分析可以做下面的图(第二张),但是我查阅相关的资料,只能做出另一张图 (第一张),第二张的纵坐标是血管内治疗相对药物治疗的OR值,但是我只能分别做出血管内治疗和药物治疗的OR值,我的算法是这样的,哪位老师帮忙看一下,怎么改可以出现上面的图:
dd <- datadist(ASPECTS_omit)
options(datadist='dd')
model1<- lrm(mRS03 ~ rcs(ASPECTS,3)+group, x=T, y=T, data=ASPECTS_omit)
an <- anova(model1)
ggplot(Predict(model1), anova=an, pval=TRUE)
dd$limits$ASPECTS[2] <-7
model1=update(model1)
OR <- Predict(model1, ASPECTS=seq(1,10), group, ref.zero = T, fun=exp)
ggplot(OR, ylab='Odds Ratio', add layer=geom_hline (yintercept=1, linetype="dashed")) +theme_bw()+scale_colour_discrete(name="group",breaks=c("0", "1"),labels=c("SMT", "EVT"))