Вы можете попробовать что-то вроде этого:
library(tidyverse)
# some sample data
Stf <- tibble(
SATISFACTION = runif(10, 1, 10),
AGE = runif(10, 15, 90),
SEVERITY = runif(10, 1, 10),
ANXIETY = runif(10, 1, 10),
SURG = sample(c("yes", "no"), size = 10, replace = T)
)
# creating separate models for SURG
Stfmodel1 <- Stf %>%
split(.$SURG) %>%
map(~ lm(SATISFACTION ~ AGE + SEVERITY + ANXIETY, data = .))