如何在lrm模型(rms包)的bplot中显示优势比而不是Log Odds

时间:2018-04-01 03:57:16

标签: r interaction rms

我正在使用Frank Harrell的精彩rms软件包来使lrm模型适应连续的交互。我想以图形方式显示交互。

要做到这一点,我正在使用bplot,它会产生一个漂亮的3d回归曲面。但是,Y轴上的单位是对数。我希望显示优势比,以便读者轻松解读。

3d bplot of continuous-by-continuous interaction

我的代码如下所示:

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")

我会非常感谢任何帮助!

0 个答案:

没有答案