ggeffects::ggpredict是否有fullrange参数控制拟合线仅覆盖数据范围
实现方法
核心逻辑是限制ggpredict()生成预测值的范围仅匹配每个分组下自变量的实际数据区间,有两种常用实现方案:
方案1:用内置limit_range参数(推荐,ggeffects 1.3.0+版本支持)
ggpredict已经内置了对应功能的参数,直接开启即可完全复刻geom_smooth(method = "lm", fullrange = F)的效果:
library(tidyverse) library(ggeffects) # 准备数据 data(mtcars) mtcars <- mutate(mtcars, cyl = as.factor(cyl)) # 拟合模型 m <- lm(mpg ~ wt * cyl, data = mtcars) # 生成预测时开启limit_range = TRUE data_m <- ggeffects::ggpredict(m, terms = c("wt", "cyl"), limit_range = TRUE) # 绘图,和原有逻辑一致 ggplot() + geom_point(data = mtcars, aes(x = wt, y = mpg, color = cyl)) + geom_ribbon(data = data_m, aes(x = x, y = predicted, group = group, ymin = conf.low, ymax = conf.high, fill = group), alpha = .3) + geom_line(data = data_m, aes(x = x, y = predicted, color = group))
方案2:自定义预测数据集(兼容所有版本)
如果你的ggeffects版本较低没有limit_range参数,可以手动为每个分组生成对应自变量范围的预测数据集,再传入ggpredict即可:
# 按cyl分组生成每个分组wt实际区间内的均匀取值点 new_df <- mtcars %>% group_by(cyl) %>% summarise(wt = seq(min(wt), max(wt), length.out = 100), .groups = "drop") # 基于自定义数据集生成预测 data_m <- ggeffects::ggpredict(m, newdata = new_df, terms = c("wt", "cyl"))
后续绘图逻辑和之前完全相同,输出效果和方案1一致。
内容的提问来源于stack exchange,提问作者Markus Bauer
相关产品推荐
相关产品推荐

