如何在R Studio中自定义2008为断点进行分段线性拟合并获取参数

针对你的需求，我们可以用两种方法实现固定2008为断点的分段线性拟合，以下是具体操作步骤：

方法1：手动构造变量拟合（无需额外包） #

这种方法通过构造分段变量，直接用基础的lm()函数拟合，直观易理解。

步骤1：整理数据 #

先把数据转换成数据框格式，方便后续处理：

year <- c(2000, 2001, 2002, 2003, 2004, 2005, 2006, 2007, 2008, 2009, 2010, 2011, 2012, 2013, 2014, 2015)
score1 <- c(35,  40,  45, 59, 62, 71, 78, 65, 78, 82, 89, 85, 78, 88, 96, 98)
score2 <- c(56,  45,  45, 59, 65, 73, 79, 69, 79, 88, 89, 88, 72, 85, 98, 92)

# 转换为宽格式数据框
df <- data.frame(year, score1, score2)

# 构造分段变量：当年份≥2008时，取值为year-2008，否则为0
df$post_2008 <- pmax(df$year - 2008, 0)

步骤2：拟合分段线性模型 #

分别对score1和score2拟合模型：

# 拟合score1的分段模型
model1 <- lm(score1 ~ year + post_2008, data = df)
# 查看结果（含截距、斜率、p值）
summary(model1)

# 拟合score2的分段模型
model2 <- lm(score2 ~ year + post_2008, data = df)
summary(model2)

结果解释 #

• 截距：输出中的(Intercept)是模型截距项，对应2008年之前线性段的截距计算基础
• 2008年前斜率：year对应的系数就是2008年之前的斜率
• 2008年后斜率：year系数 + post_2008系数，即为2008年之后的斜率
• p值：每一行的Pr(>|t|)列就是对应系数的显著性p值

方法2：用segmented包固定断点 #

如果你偏好使用segmented包，只需指定断点值即可，无需让包自动计算：

步骤1：加载包并拟合基础线性模型 #

# 首次使用时先安装包
install.packages("segmented")
library(segmented)

# 拟合基础线性模型
base_lm1 <- lm(score1 ~ year, data = df)
base_lm2 <- lm(score2 ~ year, data = df)

步骤2：指定断点拟合分段模型 #

# 对score1指定2008为断点
seg_model1 <- segmented(base_lm1, seg.Z = ~year, psi = list(year = 2008))
summary(seg_model1)

# 对score2指定2008为断点
seg_model2 <- segmented(base_lm2, seg.Z = ~year, psi = list(year = 2008))
summary(seg_model2)

结果解释 #

• Initial slope是2008年前的斜率，Break-point显示你指定的2008，Slope change是两段斜率的差值，因此2008年后的斜率=初始斜率+斜率变化值
• 各系数的p值同样在Pr(>|t|)列中查看

内容的提问来源于stack exchange，提问作者Skyer