Let's say I have a dataset library(tidyverse) set.seed(123) df = tibble( x = rnorm(1000), y = -0.1 * x + rnorm(1000, sd = 0.1), truncate = (y > (0.25x -0.4)) & (y < 0.25x+0.4))

...