如何在 R 中的特定范围内用扁平线拟合数据

使用dpseg计算分段断点.每行dp$segments定义一个具有start和end的段 列是该段的开始和结束点的索引，并且x1和x2列是相应的时间.最后画出每一段数据的平均值.

library(dpseg)

plot(Value ~ Time, Data)

dp <- with(Data, dpseg(x = Time, y = Value, jumps = TRUE))
segs <- dp$segments

xx <- with(segs, c(rbind(x1, x2)))
yy <- mapply(\(st,en) mean(Data$Value[st:en]), segs$start, segs$end) |> rep(each=2)
lines(xx, yy, col = "red", lwd = 3)

如果不需要垂直线段，也可以将上面xx<-...行以上的所有内容替换为以下内容.

nr <- nrow(segs)
for(i in 1:nr) with(segs[i, ],
  segments(x0 = x1, 
           y0 = mean(Data$Value[start:end]),
           x1 = x2, col = "red", lwd = 3))

年长的

但这并不奏效.

使用Strucchange中的breakpoints计算断点，然后绘制数据和拟合值.

library(strucchange)
bp <- with(Data, breakpoints(Value ~ 1))
plot(Value ~ Time, Data)
lines(fitted(bp) ~ Time, Data, col = "red", lwd = 2)

注意事项

由于在问题中使用随机数来生成数据，因此我们首先使用了set.Seed以使其可重现.

set.seed(123)
Data <- data.frame(Time = c(1:100, 101:500, 501:1000), 
                   Value = c(runif(100, 0, 500), 
                             runif(400, 500, 1000), 
                             runif(500, 0, 100)),
                   Desired_fit = c(rep_len(250,100), 
                                   rep_len(750,400),
                                   rep_len(50,500)))