# Fit a linear model for each response variable using the same predictor
df <- data.frame(
y1 = rnorm(100),
y2 = rnorm(100),
x = rnorm(100)
)
# List of formulas and data
formulas <- list(y1 ~ x, y2 ~ x)
data_list <- list(df, df)
fmapn(list(formula = formulas, data = data_list), function(formula, data) {
lm(formula, data = data)
})
# Extract model summaries in parallel
models <- fmapn(list(formula = formulas, data = data_list), function(formula, data) {
summary(lm(formula, data = data))$r.squared
})
Run the code above in your browser using DataLab