我正在使用Frank Harrell的精彩rms
软件包来使lrm模型适应连续的交互。我想以图形方式显示交互。
要做到这一点,我正在使用bplot
,它会产生一个漂亮的3d回归曲面。但是,Y轴上的单位是对数。我希望显示优势比,以便读者轻松解读。
我的代码如下所示:
bplot(Predict(model,
cont.var1=seq(5,110),
cont.var2, np=50, ref.zero=TRUE),
lfun=wireframe,
adj.subtitle=FALSE,
drape=TRUE,
color.key=FALSE,
col.regions=terrain.colors(1000),
asp=1,
screen=list(z=-35, x=-60),
ylabrot=48,
xlabrot=-25,
main="Effect of cont.var2 on the relationship between cont.var1 and outcome")
我会非常感谢任何帮助!