R 趋势和分段回归

我想得到分段回归的四个斜率.365天之前和365天之后，每种释放类型有两个坡度.我也知道我应该使用emmeans软件包.

这是一个虚拟数据集.

 df <- data.frame (tsr = c(0,0,9,10,19,20,20,21, 30,30,100,101,200,205,350,360, 400,401,500,501,600,605,700,710,800,801,900,902,1000,1001,1100,1105,2000,2250,2500,2501),
              release_type = c('S','H','S','H','S','S','H','S','H','S','S','H','S','H','S','S','H','S','H','S','S','H','S','H','S','S','H','S','H','S','S','H','S','H','S', 'H'),
              cond = c(250,251,250,251,300,301,351,375,250,249,216,257,264,216,250,251,250,251,300,301,351,375,250,249,216,257,264,216, 250,251,250,251,300,301,351,375),
              notch = c('A','B','C','D','A','B','C','D','A','B','C','D','E','G','E','G','A','H','J','K','L','Q','W','E','R','Y','U','I','O','P','Y','U','I','O','P', 'Z'))

 #Load libraries
 library(emmeans)
 library(lme4)

 #Set up break point manually
 bp = 365
 b1 <- function(x, bp) ifelse(x < bp, bp - x, 0)
 b2 <- function(x, bp) ifelse(x < bp, 0, x - bp)
 
 #Fit linear mixed effect model using piecewise regression
 m1 <- lmer(cond~b1(tsr, bp) + b2(tsr,bp) + b1(tsr, bp):release_type 
       + b2(tsr,bp):release_type + release_type + (1|notch), data = df)

 #Obtain slopes
 emtrends(m1, params = "bp", var = "tsr", pairwise ~ release_type)

我只得到每种版本类型的一个斜率的估计值.我做错了什么？

注意:我不能使用summary()函数来获取斜率，因为它使用上面的函数来生成这些估计值.所以这不是纯斜坡.