Questions tagged «piecewise-linear»

当我还需要估计其他随机效应时，有人可以告诉我如何让R估计分段线性模型中的断裂点（作为固定或随机参数）吗？ 我在下面提供了一个玩具示例，该示例适合曲棍球杆/折断的杆回归，其随机点的斜率变化和y轴截距的随机变化为4的断裂点。我想估算断裂点而不是指定断裂点。它可以是随机效果（最好）或固定效果。 library(lme4) str(sleepstudy) #Basis functions bp = 4 b1 <- function(x, bp) ifelse(x < bp, bp - x, 0) b2 <- function(x, bp) ifelse(x < bp, 0, x - bp) #Mixed effects model with break point = 4 (mod <- lmer(Reaction ~ b1(Days, bp) + b2(Days, bp) + (b1(Days, …

14 r  mixed-model  lme4-nlme  change-point  piecewise-linear 

情况 我有一个具有一个因变量和一个独立变量的数据集。我想用出现在已知/固定断点拟合连续分段线性回归。已知breakpoins没有不确定性，所以我不想估计它们。然后，我拟合以下形式的回归（OLS）： 这是ÿÿyXXxķķk（一个1个，一个2，… ，一个ķ）（一个1个，一个2，…，一个ķ）(a_{1}, a_{2}, \ldots, a_{k})ÿ一世=β0+β1个X一世+β2最高（X一世-一个1个，0 ）+β3最高（X一世-一个2，0 ）+ … +βk + 1最高（X一世-一个ķ，0）+ϵ一世ÿ一世=β0+β1个X一世+β2最高⁡（X一世-一个1个，0）+β3最高⁡（X一世-一个2，0）+…+βķ+1个最高⁡（X一世-一个ķ，0）+ϵ一世 y_{i} = \beta_{0} + \beta_{1}x_{i} + \beta_{2}\operatorname{max}(x_{i}-a_{1},0) + \beta_{3}\operatorname{max}(x_{i}-a_{2},0) +\ldots+ \beta_{k+1}\operatorname{max}(x_{i}-a_{k},0) +\epsilon_{i} R set.seed(123) x <- c(1:10, 13:22) y <- numeric(20) y[1:10] <- 20:11 + rnorm(10, 0, 1.5) y[11:20] <- seq(11, 15, len=10) + rnorm(10, 0, 2) …

9 r  regression  standard-error  piecewise-linear