PROSAGA码农传奇-模型可视化-如何在ggplot中可视化线性模型预测以及置信区间？

如何在ggplot中可视化线性模型预测以及置信区间？

作者: 薄情
发布时间: 2025-02-15 08:27:52 (6天前)
转自：

<div class =“post-text”itemprop =“text”>
  <P>
    如果你翻转你的轴，你可以使用
     <code>
      geom_ribbon
    </code>
    。你可以从中获得必要的数字
     <code>
      predict
    </code>
     如果你告诉它你也想要一个置信区间：
  </p>
   <pre class="lang-r prettyprint-override">
    <code>
      library(ggplot2)

my_df <- data.frame(
    res = c(2,3.1,4.5,5.1,6.5,7.1,8.5,9.11,10.1,11.8,12.3), 
    predictor1 = c(4.2, 5.3, 5.68,6.5, 7.77,8.5,9.5, 10.18,11.64,12.15,14.19), 
    predictor2 = c(3.1, 5.2, 6.3,7.1, 9.7,11.5,12.99, 14.5 ,15.5,16.41,17.6))

model1<- lm(res ~ predictor1 + predictor2, data = my_df)

#____Another dataset
observeddata <- data.frame(
    response = c(12.5,13.5,14.65,16.1,16.5,17.22,18.54,21.31,23.61,25.58,26.43), 
    predictor1 = c(14.21, 15.13, 16.25,16.5, 17.37,18.51,19.35, 22.18,23.64,25.12,26.19), 
    predictor2 = c(13.11, 15.22, 16.23,17.41, 18.72,21.5,22.99, 24.35 ,25.15,26.21,28.5))

# adds 3 columns (fit, lwr, upr), so use `cbind` instead of just adding a column
observeddata <- cbind(observeddata, 
                      predict(model1, newdata = observeddata, interval = 'confidence'))

# add extra aesthetics for geom_ribbon
ggplot(data = observeddata, aes(x = response, y = fit, ymin = lwr, ymax = upr)) + 
    geom_point() + 
    geom_ribbon(alpha = .3)    # set opacity so points are visible

</code>
  </pre>
  <P>
    如果你真的想要你的轴，请添加
     <code>
      coord_flip
    </code>
    。
  </p>
  <P>
    <img src =“https://i.imgur.com/FMXqkIl.png”alt =“”/>
  </p>
</DIV>