Без предоставления каких-либо данных или воспроизводимого кода, нелегко точно знать, что вы хотите.Помогает ли это?
# Create some mock data:
n <- 100
x <- rnorm(n)
y <- rnorm(n, mean = x)
grp <- sample(LETTERS[1:3], replace = TRUE, size = n)
# Do linear regressions for each group
lin_regs <- lapply(unique(grp), function(g) lm(y ~ x, subset = g==grp))
# Plot lines
plot(x, y, type = "n") # Initialize plot
for (i in seq_along(lin_regs)) {
abline(lin_regs[[i]], col = i) # Plot each linear reg. with abline()
}