Fitting a linear model with multiple LHS
If I don’t get you wrong, you are working with a dataset like this: set.seed(0) dat <- data.frame(y1 = rnorm(30), y2 = rnorm(30), y3 = rnorm(30), x1 = rnorm(30), x2 = rnorm(30), x3 = rnorm(30)) x1, x2 and x3 are covariates, and y1, y2, y3 are three independent response. You are trying to fit three … Read more