r - 分段 R 的固定断点 (psi) 值

Question

嗨，我正在使用 R 中的分段。在这样做之后，我获得了 10 个断点，所以现在我想获取这个断点的一个子集并再次进行分段，因为对于我的 porpouse 来说，10 个断点太多了。此外，我想从我获得的十个中获取一些特定的断点。

# Segmented does not allow referenced to the dataframe class
lx <- life_table$lx
x <- life_table$x
#First we create a linear model
LinearMod <- lm(lx~x)
seg_5 <- segmented(LinearMod, seg.Z = ~x)
vec_breakp = as.numeric(seg_5$psi[,2])
vec_psi <- vec_breakp[4:6]
seg_mod <- segmented(LinearMod, seg.Z = ~x, psi = vec_breakp)

The output is:
    ***Regression Model with Segmented Relationship(s)***

Call: 
segmented.lm(obj = LinearMod, seg.Z = ~x, psi = vec_psi)

Estimated Break-Point(s):

           Est. St.Err

psi1.x  21.448  0.572
psi2.x 134.938  3.214
psi3.x 525.438 10.565

断点并不完全是我选择的。有可能修复它们吗？

标签： r