N <- 1000 # Number of observations
Nx <- 20 # Number of independent variables
Ny <- 80 # Number of dependent variables
# Simulate outcomes that are all standard Gaussians
Y <- matrix(rnorm(N*Ny), ncol=Ny)
X <- matrix(rnorm(N*Nx), ncol=Nx)
plr(Y, X)
Run the code above in your browser using DataLab